add flat syst for EIcC and include NLO & NNLO variants

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/data_nlo.yaml

@@ -0,0 +1,21 @@
+data_central:
+- 0.028210496511775596
+- 0.03096850302069775
+- 0.028405214684745425
+- 0.0031665906441841023
+- 0.03813398704512427
+- 0.04244941465988825
+- 0.03395072969615667
+- -0.005556107683203207
+- 0.01825684505621425
+- 0.0063719846735700286
+- 0.006106962503821987
+- -0.006758105867902573
+- 0.02560844071973423
+- -0.008759632019291676
+- 0.04107561191577555
+- 0.035788966870756375
+- 0.06851722723586647
+- -0.03574387083780241
+- 0.07755847113097429
+- 0.004509355592291467

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/data_nnlo.yaml

@@ -0,0 +1,21 @@
+data_central:
+- 0.022150955423378773
+- 0.11064663829869174
+- 0.19055515624028524
+- 0.23223430658861016
+- 0.30739887788814296
+- 0.35414350508835424
+- 0.4308815798230844
+- 0.447739711257262
+- 0.465814133492104
+- 0.5389832976785235
+- 0.5720400569186944
+- 0.6070984096401357
+- 0.6403726085326178
+- 0.652022714098933
+- 0.7049357934993605
+- 0.6962202977759955
+- 0.7345431060084674
+- 0.7329213926296617
+- 0.6371952442378601
+- 0.1503768622983805

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/filter.py

@@ -15,6 +15,18 @@
 if __name__ == "__main__":
     input_txt = Path("./rawdata/EICC_3_20_A1c_100fb-1.txt")
     df = read_txt_data(input_txt)
-    cv_preds = read_central_values(Path("./rawdata/EICC_NC_15GEV_EP.yaml"))
-    fluctuated_cv = fluctuate_data(cv_preds, df["abs"].values)
-    write_data(df, asym=False, abserr=fluctuated_cv, add_fluctuate=True)
+
+    for pto in ["NLO", "NNLO"]:
+        cv_preds = read_central_values(Path(f"./rawdata/EICC_NC_15GEV_EP_{pto}.yaml"))
+        fluctuated_cv = fluctuate_data(cv_preds, df["abs"].values)
+        # Flat systematic error of 3.2% due to beam polarization
+        # from Yuxiang suggestion
+        sys_error = 0.050
+        write_data(
+            df,
+            asym=False,
+            abserr=fluctuated_cv,
+            add_fluctuate=True,
+            sys_error=sys_error,
+            suffix=pto.lower(),
+        )

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/metadata.yaml

@@ -17,9 +17,9 @@ nnpdf_metadata:
   experiment: "EIC"
 
 implemented_observables:
-  - observable_name: "G1F1CHARMRATIO"
+  - observable_name: "G1F1CHARMRATIO-NLO"
     observable:
-      description: "EICC projection measurements on polarized ep scattering at 15GeV"
+      description: "EICC projection measurements on polarized ep scattering at 15GeV using NLO theory"
       label: "$g_1^c / F_1^c$"
       units: ""
     process_type: "DIS_POL"
@@ -43,9 +43,43 @@ implemented_observables:
         Q2: { description: "virtuality", label: "$Q^2$", units: "$GeV^2$" }
       file: kinematics.yaml
 
-    data_central: data.yaml
+    data_central: data_nlo.yaml
     data_uncertainties:
-      - uncertainties.yaml
+      - uncertainties_nlo.yaml
+
+    theory:
+      FK_tables:
+        - - EIcC_NC_15GEV_EP_G1
+      operation: "null"
+  - observable_name: "G1F1CHARMRATIO-NNLO"
+    observable:
+      description: "EICC projection measurements on polarized ep scattering at 15GeV using NNLO theory"
+      label: "$g_1^c / F_1^c$"
+      units: ""
+    process_type: "DIS_POL"
+    ndata: 20
+    tables: []
+    npoints: [] # number of datapoints in each table
+
+    # Plotting information
+    plotting:
+      kinematics_override: dis_sqrt_scale
+      dataset_label: "EICC 15GeV (g1c/F1c)"
+      plot_x: x
+      line_by: [Q2]
+      figure_by: []
+      y_label: "$g_{1,p}^c(x, Q^2) / F_{1,p}^c(x, Q^2)$"
+    kinematic_coverage: [x, Q2]
+
+    kinematics:
+      variables:
+        x: { description: "momentum fraction", label: "$x$", units: "" }
+        Q2: { description: "virtuality", label: "$Q^2$", units: "$GeV^2$" }
+      file: kinematics.yaml
+
+    data_central: data_nnlo.yaml
+    data_uncertainties:
+      - uncertainties_nnlo.yaml
 
     theory:
       FK_tables:

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/rawdata/EICC_NC_15GEV_EP_NLO.yaml

@@ -0,0 +1,21 @@
+predictions_central:
+- 0.010555049103444576
+- 0.007086527909554543
+- 0.006547566243491389
+- 0.006567545892765427
+- 0.00670800108690537
+- 0.006995637438592223
+- 0.007368994860621027
+- 0.007764406824890421
+- 0.008182083942964291
+- 0.008567462446159062
+- 0.00898084519247325
+- 0.009348265543498852
+- 0.009656663912778592
+- 0.009828543306697114
+- 0.009769402693172113
+- 0.009361906063754795
+- 0.008283342645634911
+- 0.006555636418967348
+- 0.002896542313283272
+- -0.013171498061207662

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/rawdata/EICC_NC_15GEV_EP_NNLO.yaml

@@ -0,0 +1,21 @@
+predictions_central:
+- 0.03897899400867568
+- 0.12181248092748868
+- 0.18619698924227923
+- 0.24465145008218833
+- 0.29854301454322285
+- 0.34851807502101745
+- 0.3964844353133432
+- 0.44183763775436
+- 0.4847243439040604
+- 0.5197528272036527
+- 0.5571224773329402
+- 0.592577237802076
+- 0.6235427919188653
+- 0.6580416532417295
+- 0.6846068194659319
+- 0.7062840853754252
+- 0.7185271667037855
+- 0.6930629357128841
+- 0.6108911277189878
+- 0.15523776948470808

nnpdf_data/nnpdf_data/new_commondata/EICC_NC_15GEV_EP/uncertainties_nlo.yaml

@@ -0,0 +1,50 @@
+definitions:
+  stat:
+    description: statistical uncertainty
+    treatment: ADD
+    type: UNCORR
+  sys:
+    description: systematic uncertainty
+    treatment: MULT
+    type: CORR
+bins:
+- stat: 0.0131741
+  sys: 0.0014105248255887799
+- stat: 0.0137443
+  sys: 0.0015484251510348877
+- stat: 0.0145816
+  sys: 0.0014202607342372713
+- stat: 0.0148745
+  sys: 0.00015832953220920511
+- stat: 0.0154729
+  sys: 0.0019066993522562133
+- stat: 0.0163357
+  sys: 0.0021224707329944126
+- stat: 0.0166505
+  sys: 0.0016975364848078338
+- stat: 0.0173117
+  sys: -0.00027780538416016033
+- stat: 0.0177981
+  sys: 0.0009128422528107125
+- stat: 0.0184076
+  sys: 0.00031859923367850144
+- stat: 0.0196155
+  sys: 0.0003053481251910994
+- stat: 0.0201606
+  sys: -0.00033790529339512867
+- stat: 0.0206674
+  sys: 0.0012804220359867117
+- stat: 0.0224443
+  sys: -0.0004379816009645838
+- stat: 0.0235975
+  sys: 0.0020537805957887775
+- stat: 0.0249343
+  sys: 0.0017894483435378188
+- stat: 0.028093
+  sys: 0.0034258613617933233
+- stat: 0.0294843
+  sys: -0.0017871935418901205
+- stat: 0.03264
+  sys: 0.003877923556548715
+- stat: 0.0412479
+  sys: 0.00022546777961457337