Add Reference struct

pineappl/src/grid.rs

@@ -25,6 +25,7 @@ use std::io::{self, BufRead, BufReader, BufWriter, Read, Write};
 use std::ops::Range;
 use std::{iter, mem};
 use thiserror::Error;
+use super::reference::Reference;
 
 const BIN_AXIS: Axis = Axis(1);
 
@@ -133,6 +134,7 @@ pub struct Grid {
     scales: Scales,
     metadata: BTreeMap<String, String>,
     more_members: MoreMembers,
+    reference: Reference,
 }
 
 impl Grid {
@@ -193,6 +195,7 @@ impl Grid {
             scales,
             metadata: default_metadata(),
             more_members: MoreMembers::V4(Mmv4),
+            reference: Reference::default(),
         }
     }
 
@@ -1384,6 +1387,8 @@ impl Grid {
                 },
                 metadata: self.metadata.clone(),
                 more_members: self.more_members.clone(),
+                // TODO: transform the reference result to match the FKTable structure
+                reference: self.reference.clone(),
             };
 
             if let Some(lhs) = &mut lhs {

pineappl/src/lib.rs

@@ -47,4 +47,5 @@ pub mod interp_subgrid;
 pub mod interpolation;
 pub mod packed_array;
 pub mod pids;
+pub mod reference;
 pub mod subgrid;

pineappl/src/reference.rs

@@ -0,0 +1,32 @@
+//! TODO
+
+use serde::{Deserialize, Serialize};
+
+#[derive(Clone, Deserialize, Serialize)]
+enum RefRes{
+    ByBin(f64),
+    ByBinOrder(Vec<f64>),
+    ByBinChannel(Vec<f64>),
+    ByBinOrderChannel(Vec<Vec<f64>>),
+}
+
+#[derive(Clone, Deserialize, Serialize)]
+enum RefRelUnc{
+    ByBin(f64),
+    ByBinOrder(Vec<f64>, CombOp),
+    ByBinChannel(Vec<f64>, CombOp),
+    ByBinOrderChannel(Vec<Vec<f64>>, CombOp),
+}
+
+#[derive(Clone, Deserialize, Serialize)]
+enum CombOp{
+    Sum,
+    Quadrature,
+}
+
+/// TODO
+#[derive(Clone, Default, Deserialize, Serialize)]
+pub struct Reference{
+    ref_res_unc: Vec<(RefRes, RefRelUnc)>,
+    ref_conv_fun: Vec<String>,
+}