diff --git a/n3fit/runcards/examples/runcard_OtG-nlo.yml b/n3fit/runcards/examples/runcard_OtG-nlo.yml
new file mode 100644
index 000000000..2d46d6dfd
--- /dev/null
+++ b/n3fit/runcards/examples/runcard_OtG-nlo.yml
@@ -0,0 +1,202 @@
+#
+# Configuration file for n3fit
+#
+############################################################
+description: "SIMUnet OtG-nlo benchmark with ttbar data."
+
+############################################################
+# frac: training fraction
+# ewk: apply ewk k-factors
+# sys: systematics treatment (see systypes)
+dataset_inputs:
+- {dataset: NMCPD_dw_ite, frac: 0.75}
+- {dataset: NMC, frac: 0.75}
+- {dataset: SLACP_dwsh, frac: 0.75}
+- {dataset: SLACD_dw_ite, frac: 0.75}
+- {dataset: BCDMSP_dwsh, frac: 0.75}
+- {dataset: BCDMSD_dw_ite, frac: 0.75}
+- {dataset: CHORUSNUPb_dw_ite, frac: 0.75}
+- {dataset: CHORUSNBPb_dw_ite, frac: 0.75}
+- {dataset: NTVNUDMNFe_dw_ite, frac: 0.75, cfac: [MAS]}
+- {dataset: NTVNBDMNFe_dw_ite, frac: 0.75, cfac: [MAS]}
+- {dataset: HERACOMBNCEM, frac: 0.75}
+- {dataset: HERACOMBNCEP460, frac: 0.75}
+- {dataset: HERACOMBNCEP575, frac: 0.75}
+- {dataset: HERACOMBNCEP820, frac: 0.75}
+- {dataset: HERACOMBNCEP920, frac: 0.75}
+- {dataset: HERACOMBCCEM, frac: 0.75}
+- {dataset: HERACOMBCCEP, frac: 0.75}
+- {dataset: HERACOMB_SIGMARED_C, frac: 0.75} 
+- {dataset: HERACOMB_SIGMARED_B, frac: 0.75}   
+- {dataset: DYE886R_dw_ite, frac: 0.75, cfac: [QCD]}
+- {dataset: DYE886P, frac: 0.75, cfac: [QCD]}
+- {dataset: DYE605_dw_ite, frac: 0.75, cfac: [QCD]}
+- {dataset: DYE906R_dw_ite, frac: 0.75, cfac: [ACC,QCD]}
+- {dataset: CDFZRAP_NEW, frac: 0.75, cfac: [QCD]}
+- {dataset: D0ZRAP_40, frac: 0.75, cfac: [QCD]}
+- {dataset: D0WMASY, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLASWZRAP36PB, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLASZHIGHMASS49FB, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLASLOMASSDY11EXT, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLASWZRAP11CC, frac: 0.75, cfac: [QCD]}               
+- {dataset: ATLASWZRAP11CF, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLASDY2D8TEV, frac: 0.75, cfac: [QCDEWK]}
+- {dataset: ATLAS_DY_2D_8TEV_LOWMASS, frac: 0.75, cfac: [QCD]} 
+- {dataset: ATLAS_WZ_TOT_13TEV, frac: 0.75, cfac: [NRM, QCD]}       
+- {dataset: ATLAS_WP_JET_8TEV_PT, frac: 0.75, cfac: [QCD]}         
+- {dataset: ATLAS_WM_JET_8TEV_PT, frac: 0.75, cfac: [QCD]}         
+- {dataset: ATLASZPT8TEVMDIST, frac: 0.75, cfac: [QCD], sys: 10}
+- {dataset: ATLASZPT8TEVYDIST, frac: 0.75, cfac: [QCD], sys: 10}
+- {dataset: ATLAS_1JET_8TEV_R06_DEC, frac: 0.75, cfac: [QCD]}
+- {dataset: ATLAS_2JET_7TEV_R06, frac: 0.75, cfac: [QCD]}     
+- {dataset: ATLASPHT15_SF, frac: 0.75, cfac: [QCD, EWK]}         
+- {dataset: CMSWEASY840PB, frac: 0.75, cfac: [QCD]}
+- {dataset: CMSWMASY47FB, frac: 0.75, cfac: [QCD]}
+- {dataset: CMSDY2D11, frac: 0.75, cfac: [QCD]}
+- {dataset: CMSWMU8TEV, frac: 0.75, cfac: [QCD]}
+- {dataset: CMSZDIFF12, frac: 0.75, cfac: [QCD, NRM], sys: 10}
+- {dataset: CMS_2JET_7TEV, frac: 0.75, cfac: [QCD]}       
+- {dataset: CMS_1JET_8TEV, frac: 0.75, cfac: [QCD]}
+- {dataset: LHCBZ940PB, frac: 0.75, cfac: [QCD]}
+- {dataset: LHCBZEE2FB_40, frac: 0.75, cfac: [QCD]}
+- {dataset: LHCBWZMU7TEV, frac: 0.75, cfac: [NRM, QCD]}
+- {dataset: LHCBWZMU8TEV, frac: 0.75, cfac: [NRM, QCD]}
+- {dataset: LHCB_Z_13TEV_DIMUON, frac: 0.75, cfac: [QCD]}
+- {dataset: LHCB_Z_13TEV_DIELECTRON, frac: 0.75, cfac: [QCD]} # HERE STARTS THE TOP SECTOR 
+#- {dataset: ATLASTTBARTOT7TEV, cfac: [QCD]} # No C-facs
+- {dataset: ATLASTTBARTOT8TEV, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TOPDIFF_DILEPT_8TEV_TTMNORM, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_8TEV_LJETS_TOTAL, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTB_DIFF_8TEV_LJ_TRAPNORM, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTB_DIFF_8TEV_LJ_TTRAPNORM, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_DILEPTON_TOTAL, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_HADRONIC_TOTAL, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_HADRONIC_2D_TTM_ABSYTTNORM, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_LJETS_TOTAL, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_TTMNORM, cfac: [QCD], fit_cfac: True} 
+#- {dataset: CMSTTBARTOT5TEV, frac: 0.75, cfac: [QCD]} No C-facs 
+#- {dataset: CMSTTBARTOT7TEV, frac: 0.75, cfac: [QCD]} # No C-facs
+- {dataset: CMSTTBARTOT8TEV, frac: 0.75, cfac: [QCD], fit_cfac: True} 
+- {dataset: CMS_TTBAR_2D_DIFF_MTT_TTRAP_NORM, frac: 0.75, cfac: [QCD], fit_cfac: True} 
+- {dataset: CMSTOPDIFF8TEVTTRAPNORM, frac: 0.75, cfac: [QCD], fit_cfac: True} 
+- {dataset: CMSTTBARTOT13TEV, frac: 0.75, cfac: [QCD], fit_cfac: True}
+- {dataset: CMS_TTB_DIFF_13TEV_2016_2L_TTMNORM, frac: 0.75, fit_cfac: True}
+- {dataset: CMS_TTBAR_13TEV_LJETS_TOTAL, frac: 0.75, cfac: [QCD], fit_cfac: True} 
+- {dataset: CMS_TTBAR_13TEV_LJETS_2D_TTM_ABSYTTNORM, frac: 0.75, cfac: [QCD], fit_cfac: True} 
+- {dataset: ATLAS_CMS_TTBAR_8TEV_ASY, frac: 0.75, fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_13TEV_ASY, frac: 0.75, fit_cfac: True} 
+- {dataset: ATLAS_TTBAR_8TEV_ASY, frac: 0.75, fit_cfac: True} 
+- {dataset: CMS_TTBAR_13TEV_ASY, frac: 0.75, fit_cfac: True} 
+- {dataset: CMS_TTBAR_8TEV_ASY, frac: 0.75, fit_cfac: True}        # FROM HERE DOWNWARDS THEY DON'T ENTER THE BENCHMARK
+#- {dataset: CMS_TTBARZ_8TEV_TOTAL, frac: 0.75}
+#- {dataset: CMS_TTBARZ_13TEV_TOTAL, frac: 0.75, cfac: [QCD]}
+#- {dataset: ATLAS_TTBARZ_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_TTBARZ_13TEV_TOTAL, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_TTBARW_8TEV_TOTAL, frac: 0.75}
+#- {dataset: CMS_TTBARW_13TEV_TOTAL, frac: 0.75, cfac: [QCD]}
+#- {dataset: ATLAS_TTBARW_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_TTBARW_13TEV_TOTAL, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_TCH_8TEV_T, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_TCH_8TEV_TB, frac: 0.75, cfac: [QCD]}
+#- {dataset: ATLAS_SINGLETOP_TCH_DIFF_8TEV_T_RAP, frac: 0.75, cfac: [QCD]}
+#- {dataset: ATLAS_SINGLETOP_TCH_DIFF_8TEV_TBAR_RAP, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_SCH_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_SINGLETOP_SCH_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_SINGLETOP_TCH_13TEV_T, frac: 0.75, cfac: [QCD]}
+#- {dataset: ATLAS_SINGLETOP_TCH_13TEV_TB, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_TCH_13TEV_T, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_TCH_13TEV_TB, frac: 0.75, cfac: [QCD]}
+#- {dataset: CMS_SINGLETOP_TCH_13TEV_YTNORM, frac: 0.75}
+#- {dataset: ATLAS_SINGLETOPW_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_SINGLETOPW_8TEV_SLEP_TOTAL, frac: 0.75}
+#- {dataset: CMS_SINGLETOPW_8TEV_TOTAL, frac: 0.75}
+#- {dataset: ATLAS_SINGLETOPW_13TEV_TOTAL, frac: 0.75}
+#- {dataset: CMS_SINGLETOPW_13TEV_TOTAL, frac: 0.75}
+#- {dataset: CMS_SINGLETOPW_13TEV_SLEP_TOTAL, frac: 0.75}
+
+############################################################
+fit_cfactors: ["OtG-nlo"]
+
+############################################################
+datacuts:
+  t0pdfset: 210219-02-rs-nnpdf40-1000 # PDF set to generate t0 covmat
+  q2min: 3.49                        # Q2 minimum
+  w2min: 12.5                        # W2 minimum
+
+############################################################
+theory:
+  theoryid: 200     # database id
+
+norm_threshold: 4
+
+############################################################
+cfactorscale: 1.0
+trvlseed: 475038818
+nnseed: 2394641471
+mcseed: 1831662593
+save: false
+genrep: true      # true = generate MC replicas, false = use real data
+
+############################################################
+
+parameters: # This defines the parameter dictionary that is passed to the Model Trainer
+  nodes_per_layer: [25, 20, 8]
+  activation_per_layer: [tanh, tanh, linear]
+  initializer: glorot_normal
+  optimizer:
+    clipnorm: 6.073e-6
+    learning_rate: 2.621e-3
+    optimizer_name: Nadam
+  epochs: 17000
+  positivity:
+    initial: 184.8
+    multiplier:
+  integrability:
+    initial: 10
+    multiplier:
+  stopping_patience: 0.1
+  layer_type: dense
+  dropout: 0.0
+  threshold_chi2: 3.5
+
+fitting:
+  savepseudodata: True
+  fitbasis: EVOL  # EVOL (7), EVOLQED (8), etc.
+  basis:
+  - {fl: sng, trainable: false, smallx: [1.093, 1.121], largex: [1.486, 3.287]}
+  - {fl: g, trainable: false, smallx: [0.8329, 1.071], largex: [3.084, 6.767]}
+  - {fl: v, trainable: false, smallx: [0.5202, 0.7431], largex: [1.556, 3.639]}
+  - {fl: v3, trainable: false, smallx: [0.1205, 0.4839], largex: [1.736, 3.622]}
+  - {fl: v8, trainable: false, smallx: [0.5864, 0.7987], largex: [1.559, 3.569]}
+  - {fl: t3, trainable: false, smallx: [-0.5019, 1.126], largex: [1.754, 3.479]}
+  - {fl: t8, trainable: false, smallx: [0.6305, 0.8806], largex: [1.544, 3.481]}
+  - {fl: t15, trainable: false, smallx: [1.087, 1.139], largex: [1.48, 3.365]}
+
+############################################################
+positivity:
+  posdatasets:
+  - {dataset: POSF2U, maxlambda: 1e6}        # Positivity Lagrange Multiplier
+  - {dataset: POSF2DW, maxlambda: 1e6}
+  - {dataset: POSF2S, maxlambda: 1e6}
+  - {dataset: POSFLL, maxlambda: 1e6}
+  - {dataset: POSDYU, maxlambda: 1e10}
+  - {dataset: POSDYD, maxlambda: 1e10}
+  - {dataset: POSDYS, maxlambda: 1e10}
+  - {dataset: POSF2C, maxlambda: 1e6}
+  - {dataset: POSXUQ, maxlambda: 1e6}        # Positivity of MSbar PDFs
+  - {dataset: POSXUB, maxlambda: 1e6}
+  - {dataset: POSXDQ, maxlambda: 1e6}
+  - {dataset: POSXDB, maxlambda: 1e6}
+  - {dataset: POSXSQ, maxlambda: 1e6}
+  - {dataset: POSXSB, maxlambda: 1e6}
+  - {dataset: POSXGL, maxlambda: 1e6}
+
+############################################################
+integrability:
+  integdatasets:
+  - {dataset: INTEGXT8, maxlambda: 1e2}
+  - {dataset: INTEGXT3, maxlambda: 1e2}
+
+############################################################
+debug: false
+maxcores: 4