Fix bug in DIS evolution

NNPDF · Aug 6, 2024 · e8dba2b · e8dba2b
1 parent 2416f20
commit e8dba2b
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/pineappl/src/evolution.rs b/pineappl/src/evolution.rs
@@ -500,7 +500,7 @@ pub(crate) fn evolve_slice_with_one(
         }));
     }
 
-    let pid = grid.channels()[0].entry()[0].0[index];
+    let pid = grid.channels()[0].entry()[0].0[1 - index];
 
     Ok((
         Array1::from_iter(sub_fk_tables)

diff --git a/pineappl/src/grid.rs b/pineappl/src/grid.rs
@@ -1582,6 +1582,9 @@ impl Grid {
         // TODO: convert this unwrap into error
         let grid = lhs.unwrap();
 
+        // check that the convolutions are unchanged
+        assert_eq!(self.convolutions(), grid.convolutions());
+
         // UNWRAP: merging evolved slices should be a proper FkTable again
         Ok(FkTable::try_from(grid).unwrap_or_else(|_| unreachable!()))
     }