adding MATRIX ttbar runcards

README.md

@@ -1,4 +1,16 @@
-# Runcards for NNPDF
+These runcards were used to compute top pair production processes for NNPDF4.1 with MATRIX at NNLO except for the 
+following changes:
+* The top mass was varied up (mt_+ = 175.0 GeV) and down (mt_- = 170.0 GeV) with respect to the central value mt_0 = 172.5 GeV. 
+* The renormalisation and factorisation scales were dynamical but evaluated at the central value mt_0. That is, 
+  the scales were set to `HT(mt=mt_0, pT)/4`, with `HT = (ET_t + ET_tx)` and `ET_t = sqrt(mt_0 ** 2 + pT_t ** 2)`, 
+  `ET_tx = sqrt(mt_0 ** 2 + pT_tx ** 2)`. This choice is reflected in the `paramter.dat` card,
+```
+: dynamic_scale     =  11
+```
+where `dynamic_scale = 11` corresponds to the definition written in `specify.scales.cxx`
+
+To reproduce the NNPDF4.1 predictions, these runcards should be modified by setting in `model.dat`
+```
+ 6 172.5 # M_t
+```
 
-This repository stores all pinecards needed to generate PineAPPL grids for the
-processes included in NNPDF.

TTBAR_13TEV_DIF/13tev_dddistribution.dat

@@ -0,0 +1,41 @@
+################ LHC TTBAR 13 TEV ##############
+
+# keys: gm- grids to be merged, bm- bins to be merged
+
+########### 2d distributions
+
+# atlas_l+j (merge the bins according to exp. data)
+
+dddistributionname = ATLAS_TTBAR_13TEV_LJ_DIF_PTT-YT_bm
+distributionname_1 = atl_lj_2d_yt_ptt-yt
+distributionname_2 = atl_lj_2d_ptt_ptt-yt
+
+dddistributionname = ATLAS_TTBAR_13TEV_LJ_DIF_MTTBAR-PTT_bm
+distributionname_1 = atl_lj_2d_mttb_mttb-ptt
+distributionname_2 = atl_lj_2d_ptt_mttb-ptt
+
+# atlas hadr. (merge the grids)
+
+dddistributionname      =       ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR_gm1
+distributionname_1      =       atl_hadr_2d_m1_ttx
+distributionname_2      =       atl_hadr_2d_ym1_ttx
+
+dddistributionname      =       ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR_gm2
+distributionname_1      =       atl_hadr_2d_m2_ttx
+distributionname_2      =       atl_hadr_2d_ym2_ttx
+
+dddistributionname      =       ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR_gm3
+distributionname_1      =       atl_hadr_2d_m3_ttx
+distributionname_2      =       atl_hadr_2d_ym3_ttx
+
+# cms l+j
+
+dddistributionname      =       CMS_TTBAR_13TEV_LJ_DIF_MTTBAR-YTTBAR
+distributionname_1      =       cms_lj_2d_mttbar
+distributionname_2      =       cms_lj_2d_yttbar
+
+# cms 2l
+
+dddistributionname      =       CMS_TTBAR_13TEV_2L_138FB-1_DIF_MTTBAR-YTTBAR
+distributionname_1      =       cms_2l_2d_mttbar
+distributionname_2      =       cms_2l_2d_yttbar

TTBAR_13TEV_DIF/13tev_distribution.dat

@@ -0,0 +1,219 @@
+################ LHC TTBAR 13 TEV ##############
+
+# keys: catr- convert absolute to regular
+
+########### 1d distributions
+
+# atlas_l+j
+
+distributionname = ATLAS_TTBAR_13TEV_LJ_DIF_MTTBAR
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 325:400:480:580:700:860:1020:1250:1500:2000
+
+distributionname = ATLAS_TTBAR_13TEV_LJ_DIF_PTT
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:50:100:160:225:300:360:475:1000
+
+distributionname = ATLAS_TTBAR_13TEV_LJ_DIF_YT
+distributiontype = absy
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:0.4:0.8:1.2:1.7:2.5
+
+distributionname = ATLAS_TTBAR_13TEV_LJ_DIF_YTTBAR
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0:0.25:0.5:0.8:1.1:1.4:1.8:2.5
+
+# atlas hadr.
+
+distributionname = ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 325.0:458.75:565.75:646.0:753.0:860.0:967.0:1100.75:1261.25:3000.0
+
+distributionname = ATLAS_TTBAR_13TEV_HADR_DIF_YTTBAR
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.12:0.24:0.36:0.49:0.62:0.76:0.91:1.06:1.21:1.39:1.59:2.4
+
+
+# cms 2l
+
+distributionname = CMS_TTBAR_13TEV_2L_DIF_MTTBAR
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 300:380:470:620:820:1100:1500:2500
+
+distributionname = CMS_TTBAR_13TEV_2L_138FB-1_DIF_PTT
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:55:100:165:240:330:440:600
+
+
+# (convert absolute to regular)
+distributionname = CMS_TTBAR_13TEV_2L_DIF_YTTBAR_catr
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0:0.4:0.8:1.2:1.6:2.6
+
+distributionname = CMS_TTBAR_13TEV_2L_DIF_YT
+distributiontype = y
+particle 1       = top 1
+binningtype      = irregular
+edges            = -2.6:-1.8:-1.35:-0.9:-0.45:0:0.45:0.9:1.35:1.8:2.6
+
+distributionname = CMS_TTBAR_13TEV_2L_DIF_PTT
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:65:125:200:290:400:550
+
+# cms l+j
+
+distributionname = CMS_TTBAR_13TEV_LJ_DIF_MTTBAR
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 250:400:480:560:640:720:800:900:1000:1150:1300:1500:1700:2000:2300:3500
+
+distributionname = CMS_TTBAR_13TEV_LJ_DIF_YTTBAR
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0:0.2:0.4:0.6:0.8:1:1.2:1.4:1.6:1.8:2.4
+
+distributionname = CMS_TTBAR_13TEV_LJ_DIF_YT
+distributiontype = absy
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:0.2:0.4:0.6:0.8:1:1.2:1.4:1.6:1.8:2:2.5
+
+distributionname = CMS_TTBAR_13TEV_LJ_DIF_PTT
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:40:80:120:160:200:250:300:350:400:450:500:600:700:850:1000:1500
+
+########### 1d distributions of the 2d distributions
+
+# atlas l+j
+
+distributionname = atl_lj_2d_yt_ptt-yt
+distributiontype = absy
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:0.75:1.5:2.5
+
+distributionname = atl_lj_2d_ptt_ptt-yt
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:85:175:285:300:390:1000
+
+distributionname = atl_lj_2d_mttb_mttb-ptt
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 325:500:700:1000:2000
+
+distributionname = atl_lj_2d_ptt_mttb-ptt
+distributiontype = pT
+particle 1       = top 1
+binningtype      = irregular
+edges            = 0:80:90:170:180:270:280:370:1000
+
+# atlas hadr.
+
+distributionname = atl_hadr_2d_m1_ttx
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:700.0
+
+distributionname = atl_hadr_2d_m2_ttx
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 700.0:970.0
+
+distributionname = atl_hadr_2d_m3_ttx
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 970.0:3000.0
+
+distributionname = atl_hadr_2d_ym1_ttx
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.46:0.91:1.55:2.5
+
+distributionname = atl_hadr_2d_ym2_ttx
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.33:0.62:0.95:2.5
+
+distributionname = atl_hadr_2d_ym3_ttx
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.4:0.79:2.5
+
+# cms l+j
+
+distributionname = cms_lj_2d_mttbar
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 250:420:520:620:800:1000:3500
+
+distributionname = cms_lj_2d_yttbar
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.3:0.6:0.9:1.2:1.5:2.5
+
+# cms 2l
+
+distributionname = cms_2l_2d_mttbar
+distributiontype = m
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 300:400:500:650:1500
+
+distributionname = cms_2l_2d_yttbar
+distributiontype = absy
+particle 1       = top 1
+particle 1       = atop 1
+binningtype      = irregular
+edges            = 0.0:0.35:0.75:1.15:2.50

TTBAR_13TEV_DIF/13tev_parameter.dat

@@ -0,0 +1,112 @@
+##########################
+# MATRIX input parameter #
+##########################
+
+#----------------------\
+# General run settings |
+#----------------------/
+process_class    =  pp-tt~+X   #  process id
+E                =  6500.      #  energy per beam
+coll_choice      =  1          #  (1) PP collider; (2) PPbar collider
+switch_PineAPPL  =  1          #  switch to turn on (1) or off(0) the generation of PineAppl grids
+use_TSV          =  1          #  switch to use new TSV scale variation instead of CV (1) or in addition to CV (2)
+                               #  (TSV variation that has more flexibility and features, and will become the new
+                               #   standard in future releases; the final result output differs slightly from CV,
+                               #   but contains at the moment the identical information, unless obove switch
+                               #   power_corrections or extrapolate_binwise is turned on, then TSV variation will
+                               #   be used regardless of setting use_TSV and those features included in the results)
+extrapolate_binwise = 1        #  switch for bin-wise r_cut extrapolation of distributions
+                               #  (note: increases written output for distributions by factor of 8)
+enhance_tails    =  1          #  switch to improve statistics in tail of distributions (a factor of two slower)
+
+
+#----------------\
+# Scale settings |
+#----------------/
+scale_ren         =  172.5     #  renormalization (muR) scale
+scale_fact        =  172.5     #  factorization (muF) scale
+dynamic_scale     =  11        #  dynamic ren./fac. scale
+                               #  0: fixed scale above
+                               #  1: invariant mass (Q) of system (of the colourless final states)
+                               #  2: transverse mass (mT^2=Q^2+pT^2) of system (of the colourless final states)
+factor_central_scale = 1       #  relative factor for central scale (important for dynamic scales)
+scale_variation   =  1         #  switch for muR/muF uncertainties (0) off; (1) 7-point (default); (2) 9-point variation
+variation_factor  =  2         #  symmetric variation factor; usually a factor of 2 up and down (default)
+
+
+#------------------------------\
+# Order-dependent run settings |
+#------------------------------/
+# LO-run
+run_LO             =  0        #  switch for LO cross section (1) on; (0) off
+LHAPDF_LO          =  NNPDF30_lo_as_0118  #  LO LHAPDF set
+PDFsubset_LO       =  0        #  member of LO PDF set
+precision_LO       =  1.e-2    #  precision of LO cross section
+
+# NLO-run
+run_NLO_QCD        =  0        #  switch for NLO QCD cross section (1) on; (0) off
+LHAPDF_NLO         =  NNPDF30_nlo_as_0118  #  NLO LHAPDF set
+PDFsubset_NLO      =  0        #  member of NLO PDF set
+precision_NLO_QCD  =  1.e-2    #  precision of NLO QCD cross section
+NLO_subtraction_method = 1     #  switch to use (2) qT subtraction (1) Catani-Seymour at NLO
+
+# NNLO-run
+run_NNLO_QCD       =  1        #  switch for NNLO QCD cross section (1) on; (0) off 
+LHAPDF_NNLO        =  NNPDF40_nnlo_as_01180  #  NNLO LHAPDF set
+PDFsubset_NNLO     =  0        #  member of NNLO PDF set
+precision_NNLO_QCD =  2.e-4    #  precision of NNLO QCD cross section
+switch_qT_accuracy =  0        #  switch to improve qT-subtraction accuracy (slower numerical convergence)
+                               #  0: lowest value of r_cut = 0.0015 varied up to 0.01 (default)
+                               #  1: lowest value of r_cut = 0.0005 varied up to 0.01
+                               #  2: lowest value of r_cut = 0.0001 varied up to 0.01 (only if extrapolate_binwise=1;
+                               #     output of fixed-r_cut result remains 0.0005, while 0.0001 used for extrapolation)
+                               #     for Drell-Yann it is recommended to turn on power_corrections rather than use switch_qT_accuracy
+power_corrections  =  0        #  switch to include leading power corrections in qT-subtraction through recoil
+                               #  (not recommended for processes involving photons and heavy quarks)
+power_corrections_pT0 = 25.    #  characteristic transverse momentum pT0 used to optimise the generation
+                               #  of the phase space for the integration of the power corrections. It should be set
+                               #  to the minimum requirement on the transverse momentum of the 2-body final state
+                               #  (for Drell-Yan for instance this should be the minimum transverse momentum of the leptons)
+
+
+#----------------------------\
+# Settings for fiducial cuts |
+#----------------------------/
+# Jet algorithm
+jet_algorithm = 3              #  (1) Cambridge-Aachen (2) kT (3) anti-kT
+jet_R_definition = 0           #  (0) pseudo-rapidity (1) rapidity
+jet_R = 0.4                    #  DeltaR
+
+# Jet cuts
+define_pT jet = 25.            #  requirement on jet transverse momentum (lower cut)
+define_eta jet = 4.5           #  requirement on jet pseudo-rapidity (upper cut)
+define_y jet = 1.e99           #  requirement on jet rapidity (upper cut)
+n_observed_min jet = 0         #  minimal number of observed jets (with cuts above)
+n_observed_max jet = 99        #  maximal number of observed jets (with cuts above)
+
+####################
+# User-defined cuts
+# (only used if defined in 'MATRIX/prc/$process/user/specify.cuts.cxx')
+#
+# no user-defined cuts implemented at the moment for this process
+####################
+# Fiducial cuts
+# (defined via general interface)
+#
+
+#-----------------\
+# MATRIX behavior |
+#-----------------/
+max_time_per_job = 8           #  very rough time(hours) one main run job shall take (default: 24h)
+                               #  unreliable when < 1h, use as tuning parameter for degree of parallelization
+                               #  note: becomes ineffective when job number > max_nr_parallel_jobs
+                               #        which is set in MATRIX_configuration file
+switch_distribution = 1        #  switch to turn on (1) and off (0) distributions
+save_previous_result = 1       #  switch to save previous result of this run (in result/"run"/saved_result_$i)
+save_previous_log = 0          #  switch to save previous log of this run (in log/"run"/saved_result_$i)
+#include_pre_in_results = 0     #  switch to (0) only include main run in results; (1) also all extrapolation (pre) runs;
+                               #  crucial to set to 0 if re-running main with different inputs (apart from precision)
+                               #  note: if missing (default) pre runs used if important for precision
+                               #  (separately for each contribution)
+reduce_workload = 0            #  switch to keep full job output (0), reduce (1) or minimize (2) workload on slow clusters
+random_seed = 0                #  specify integer value (grid-/pre-run reproducible)