From 31e5e49682634977142315d1aea3d58e82ac415f Mon Sep 17 00:00:00 2001 From: jacoterh <54140851+jacoterh@users.noreply.github.com> Date: Fri, 7 Feb 2025 11:44:48 +0000 Subject: [PATCH 1/2] adding MATRIX ttbar runcards --- ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat | 17 ++ ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat | 163 ++++++++++++ ATLAS_TTBAR_13TEV_HADR_DIF/model.dat | 61 +++++ ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat | 111 +++++++++ ATLAS_TTBAR_13TEV_HADR_DIF/readme.md | 16 ++ ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx | 106 ++++++++ .../dddistribution.dat | 25 ++ .../distribution.dat | 177 +++++++++++++ .../model.dat | 61 +++++ .../parameter.dat | 111 +++++++++ .../readme.md | 16 ++ .../specify.scales.cxx | 106 ++++++++ CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat | 14 ++ CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat | 135 ++++++++++ CMS_TTBAR_5TEV_TOT_X-SEC/model.dat | 61 +++++ CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat | 111 +++++++++ CMS_TTBAR_5TEV_TOT_X-SEC/readme.md | 16 ++ CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx | 106 ++++++++ README.md | 18 +- TTBAR_13TEV_DIF/dddistribution.dat | 18 ++ TTBAR_13TEV_DIF/distribution.dat | 233 +++++++++++++++++ TTBAR_13TEV_DIF/model.dat | 61 +++++ TTBAR_13TEV_DIF/parameter.dat | 111 +++++++++ TTBAR_13TEV_DIF/readme.md | 16 ++ TTBAR_13TEV_DIF/specify.scales.cxx | 106 ++++++++ TTBAR_7TEV_TOT_X-SEC/dddistribution.dat | 14 ++ TTBAR_7TEV_TOT_X-SEC/distribution.dat | 135 ++++++++++ TTBAR_7TEV_TOT_X-SEC/model.dat | 61 +++++ TTBAR_7TEV_TOT_X-SEC/parameter.dat | 111 +++++++++ TTBAR_7TEV_TOT_X-SEC/readme.md | 16 ++ TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx | 106 ++++++++ TTBAR_8TEV_DIF/dddistribution.dat | 28 +++ TTBAR_8TEV_DIF/distribution.dat | 235 ++++++++++++++++++ TTBAR_8TEV_DIF/model.dat | 61 +++++ TTBAR_8TEV_DIF/parameter.dat | 111 +++++++++ TTBAR_8TEV_DIF/readme.md | 16 ++ TTBAR_8TEV_DIF/specify.scales.cxx | 106 ++++++++ 37 files changed, 2973 insertions(+), 3 deletions(-) create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/model.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/readme.md create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/model.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md create mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/model.dat create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/readme.md create mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx create mode 100644 TTBAR_13TEV_DIF/dddistribution.dat create mode 100644 TTBAR_13TEV_DIF/distribution.dat create mode 100644 TTBAR_13TEV_DIF/model.dat create mode 100644 TTBAR_13TEV_DIF/parameter.dat create mode 100644 TTBAR_13TEV_DIF/readme.md create mode 100644 TTBAR_13TEV_DIF/specify.scales.cxx create mode 100644 TTBAR_7TEV_TOT_X-SEC/dddistribution.dat create mode 100644 TTBAR_7TEV_TOT_X-SEC/distribution.dat create mode 100644 TTBAR_7TEV_TOT_X-SEC/model.dat create mode 100644 TTBAR_7TEV_TOT_X-SEC/parameter.dat create mode 100644 TTBAR_7TEV_TOT_X-SEC/readme.md create mode 100644 TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx create mode 100644 TTBAR_8TEV_DIF/dddistribution.dat create mode 100644 TTBAR_8TEV_DIF/distribution.dat create mode 100644 TTBAR_8TEV_DIF/model.dat create mode 100644 TTBAR_8TEV_DIF/parameter.dat create mode 100644 TTBAR_8TEV_DIF/readme.md create mode 100644 TTBAR_8TEV_DIF/specify.scales.cxx diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat new file mode 100644 index 00000000..0e2c32a8 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat @@ -0,0 +1,17 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + +dddistributionname = dd_mttbar-yttbar +distributionname_1 = dd_m_ttx +distributionname_2 = dd_y_ttx + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat new file mode 100644 index 00000000..38bcc1f2 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat @@ -0,0 +1,163 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +distributionname = m_ttx +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 = y_ttx +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 + +distributionname = dd_m_ttx +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0.0:700.0:970.0:3000.0 + +distributionname = dd_y_ttx +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0.0:0.46:0.91:1.55:2.5 + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat b/ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat new file mode 100644 index 00000000..ae3eda34 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat @@ -0,0 +1,111 @@ +########################## +# 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 | +#----------------/ +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 = 2 # 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 = 1.e-3 # 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 = 12 # 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 = 1 # 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) + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md b/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx b/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat new file mode 100644 index 00000000..49069504 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat @@ -0,0 +1,25 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + +dddistributionname = dd1 +distributionname_1 = m1_ttx +distributionname_2 = ym1_ttx + +dddistributionname = dd2 +distributionname_1 = m2_ttx +distributionname_2 = ym2_ttx + +dddistributionname = dd3 +distributionname_1 = m3_ttx +distributionname_2 = ym3_ttx + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat new file mode 100644 index 00000000..526422dd --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat @@ -0,0 +1,177 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +distributionname = m1_ttx +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0.0:700.0 + +distributionname = m2_ttx +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 700.0:970.0 + +distributionname = m3_ttx +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 970.0:3000.0 + +distributionname = 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 = 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 = ym3_ttx +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0.0:0.4:0.79:2.5 + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat new file mode 100644 index 00000000..4ef9159a --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat @@ -0,0 +1,111 @@ +########################## +# 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 | +#----------------/ +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 = 2 # 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-1 # 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 = 1.e-3 # 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 = 12 # 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 = 1 # 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) + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat new file mode 100644 index 00000000..efacdc7d --- /dev/null +++ b/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat @@ -0,0 +1,14 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + + diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat new file mode 100644 index 00000000..dd181867 --- /dev/null +++ b/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat @@ -0,0 +1,135 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat new file mode 100644 index 00000000..13ce4f2a --- /dev/null +++ b/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat @@ -0,0 +1,111 @@ +########################## +# MATRIX input parameter # +########################## + +#----------------------\ +# General run settings | +#----------------------/ +process_class = pp-tt~+X # process id +E = 2510. # 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 = 0 # switch to improve statistics in tail of distributions (a factor of two slower) + + +#----------------\ +# Scale settings | +#----------------/ +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 = 2 # 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 = 1.e-3 # 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 = 12 # 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) + diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md b/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx b/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ diff --git a/README.md b/README.md index abb36d85..b71278d5 100644 --- a/README.md +++ b/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. diff --git a/TTBAR_13TEV_DIF/dddistribution.dat b/TTBAR_13TEV_DIF/dddistribution.dat new file mode 100644 index 00000000..6fbb39f1 --- /dev/null +++ b/TTBAR_13TEV_DIF/dddistribution.dat @@ -0,0 +1,18 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + + +dddistributionname = cms_lj_mttbar-yttbar +distributionname_1 = cms_lj_2d_mttbar +distributionname_2 = cms_lj_2d_yttbar + diff --git a/TTBAR_13TEV_DIF/distribution.dat b/TTBAR_13TEV_DIF/distribution.dat new file mode 100644 index 00000000..26315dfa --- /dev/null +++ b/TTBAR_13TEV_DIF/distribution.dat @@ -0,0 +1,233 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +# atlas-lj + +distributionname = atlas_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_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:50:100:160:225:300:360:475:1000 + +distributionname = atlas_yt +distributiontype = absy +particle 1 = top 1 +binningtype = irregular +edges = 0:0.4:0.8:1.2:1.7:2.5 + +distributionname = atlas_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 + +# cms-2l + +distributionname = cms_2l_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 300:380:470:620:820:1100:1500:2500 + +distributionname = cms_2l_yttbar +distributiontype = y +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = -2.6:-1.6:-1.2:-0.8:-0.4:0:0.4:0.8:1.2:1.6:2.6 + +distributionname = cms_2l_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_2l_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:65:125:200:290:400:550 + +# cms-lj + +distributionname = cms_lj_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_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_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_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 + +distributionname = cms_lj_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_lj_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 + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/TTBAR_13TEV_DIF/parameter.dat b/TTBAR_13TEV_DIF/parameter.dat new file mode 100644 index 00000000..457fff1f --- /dev/null +++ b/TTBAR_13TEV_DIF/parameter.dat @@ -0,0 +1,111 @@ +########################## +# 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 | +#----------------/ +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 = 2 # 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 = 1.e-3 # 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 = 12 # 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) + diff --git a/TTBAR_13TEV_DIF/readme.md b/TTBAR_13TEV_DIF/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/TTBAR_13TEV_DIF/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/TTBAR_13TEV_DIF/specify.scales.cxx b/TTBAR_13TEV_DIF/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/TTBAR_13TEV_DIF/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ diff --git a/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat b/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat new file mode 100644 index 00000000..efacdc7d --- /dev/null +++ b/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat @@ -0,0 +1,14 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + + diff --git a/TTBAR_7TEV_TOT_X-SEC/distribution.dat b/TTBAR_7TEV_TOT_X-SEC/distribution.dat new file mode 100644 index 00000000..dd181867 --- /dev/null +++ b/TTBAR_7TEV_TOT_X-SEC/distribution.dat @@ -0,0 +1,135 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/TTBAR_7TEV_TOT_X-SEC/parameter.dat b/TTBAR_7TEV_TOT_X-SEC/parameter.dat new file mode 100644 index 00000000..2331d4d5 --- /dev/null +++ b/TTBAR_7TEV_TOT_X-SEC/parameter.dat @@ -0,0 +1,111 @@ +########################## +# MATRIX input parameter # +########################## + +#----------------------\ +# General run settings | +#----------------------/ +process_class = pp-tt~+X # process id +E = 3500. # 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 = 0 # switch to improve statistics in tail of distributions (a factor of two slower) + + +#----------------\ +# Scale settings | +#----------------/ +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 = 2 # 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 = 1.e-3 # 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 = 12 # 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) + diff --git a/TTBAR_7TEV_TOT_X-SEC/readme.md b/TTBAR_7TEV_TOT_X-SEC/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/TTBAR_7TEV_TOT_X-SEC/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx b/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ diff --git a/TTBAR_8TEV_DIF/dddistribution.dat b/TTBAR_8TEV_DIF/dddistribution.dat new file mode 100644 index 00000000..cd7fb948 --- /dev/null +++ b/TTBAR_8TEV_DIF/dddistribution.dat @@ -0,0 +1,28 @@ +####################################################### +# MATRIX double differential distributions definition # +####################################################### +# +# In this file you can customize the distributions created during the run (examples below) +# +# They are simply defined as a combination of two single differential distributions that +# must be defined in the file 'distribution.dat' +# +#----------------------\ +# define distributions | +#----------------------/ + + +#pT distribution of the ttx system in rapidity bins. Binning defined in 'distribution.dat' + +dddistributionname = pTt-yt +distributionname_1 = cms_2l_yt +distributionname_2 = cms_2l_pTt + +dddistributionname = mttbar-yt +distributionname_1 = cms_2l_mttbar +distributionname_2 = cms_2l_yt + +dddistributionname = mttbar-yttbar +distributionname_1 = cms_2l_mttbar +distributionname_2 = cms_2l_yttbar + diff --git a/TTBAR_8TEV_DIF/distribution.dat b/TTBAR_8TEV_DIF/distribution.dat new file mode 100644 index 00000000..a04a35c9 --- /dev/null +++ b/TTBAR_8TEV_DIF/distribution.dat @@ -0,0 +1,235 @@ +################################## +# MATRIX distribution definition # +################################## +# +# In this file you can customize the distributions created during the run (examples below) +# please read the INSTRUCTIONS at the END OF THIS FILE... +# +#------\ +# Info | +#------/ +# Total rates and jet multiplicities (possibly within cuts) will automatically be included +# Add/remove arbitrary distribution-blocks, but always add/remove a full block. +# +#----------------------\ +# define distributions | +#----------------------/ + +# atlas_2l + +distributionname = atlas_2l_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 250:450:570:700:850:1000:2700 + +distributionname = atlas_2l_yttbar +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.4:0.8:1.2:2:2.8 + +# atlas_lj + +distributionname = atlas_lj_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 345:400:470:550:650:800:1100:1600 + +distributionname = atlas_lj_yttbar +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.3:0.6:0.9:1.3:2.5 + +distributionname = atlas_lj_yt +distributiontype = absy +particle 1 = top 1 +binningtype = irregular +edges = 0:0.4:0.8:1.2:1.6:2.5 + +distributionname = atlas_lj_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:60:100:150:200:260:320:400:500 + +# cms_2l + +distributionname = cms_2l_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 340:400:500:650:1500 + +distributionname = cms_2l_yttbar +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.35:0.75:1.15:2.5 + +distributionname = cms_2l_yt +distributiontype = absy +particle 1 = top 1 +binningtype = irregular +edges = 0:0.35:0.85:1.45:2.5 + +distributionname = cms_2l_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:80:150:250:600 + +# cms_lj + +distributionname = cms_lj_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 345:400:470:550:650:800:1100:1600 + +distributionname = cms_lj_yttbar +distributiontype = y +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = -2.5:-1.3:-0.9:-0.6:-0.3:0:0.3:0.6:0.9:1.3:2.5 + +distributionname = cms_lj_yt +distributiontype = y +particle 1 = top 1 +binningtype = irregular +edges = -2.5:-1.6:-1.2:-0.8:-0.4:0:0.4:0.8:1.2:1.6:2.5 + +distributionname = cms_lj_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:60:100:150:200:260:320:400:500 + +#--------\ +# Syntax | +#--------/ +# "distributionname" +# starts new distribution-block; user-defined, unique label for identification at end of run +# +# +# "distributiontype" +# specifies observable to be binned (some require more than one particle) +# +# pre-defined types: (custom definitions not yet supported) +# pT -- transverse momentum of "particle 1" +# (scalar sum of pTs if more than one "particle" given) +# m -- invariant mass of "particle 1" +# dm -- invariant-mass difference between "particle 1" and "particle 2" +# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" +# mmin -- minimal invariant mass of "particle 1" and "particle 2" +# mmax -- maximal invariant mass of "particle 1" and "particle 2" +# y -- rapidy of "particle 1" +# absy -- absolute rapidy of "particle 1" +# dy -- rapidy difference between "particle 1" and "particle 2" +# absdy -- absolute rapidy difference between "particle 1" and "particle 2" +# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" +# absdabsy -- absolute difference between absolute rapidities of "particle 1" +# and "particle 2" +# eta -- pseudo-rapidy of "particle 1" +# abseta -- absolute pseudo-rapidy of "particle 1" +# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" +# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" +# dabseta -- difference between absolute pseudo-rapidities of"particle 1" +# and "particle 2" +# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" +# and "particle 2" +# phi -- transverse-angle (phi) of "particle 1", or phi difference between +# "particle 1" and "particle 2" if two particles are given +# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" +# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" +# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" +# (scalar sum of ETs if more than one "particle" given) +# mT -- transverse mass: ET of "particle 1", if only one particle given +# experimental definition with missing energy, if more than one particle +# given; all neutrinos (missing energie) go into "particle 1" all other +# particles in "particle 2,3,4,..." +# +# special distribution types: +# muR -- Distribution of scale muR set in computation (no "particle" definition) +# muF -- Distribution of scale muF set in computation (no "particle" definition) +# pTveto -- Cumulative cross section sigma(pT0) + +#------------------\ +# Yukawa couplings | +#------------------/ +#Block YUKAWA +# 5 4.750000e+00 # M_YB +# 6 1.730000e+02 # M_YT +# 15 1.777000e+00 # M_YTAU + +#---------------\ +# decays widths | +#---------------/ +DECAY 6 0.000000 # WT +DECAY 23 0.000000 # WZ +DECAY 24 0.000000 # WW +DECAY 25 4.070000e-03 # WH + +#-----------\ +# EW inputs | +#-----------/ +Block EWINPUTS + 1 1 # ew_scheme - determines scheme used for EW inputs + # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) + # 1: G_mu scheme (G_F above used to determine inputs) + # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) + 3 1 # use_cms - switch for the complex mass scheme + # 0: off + # 1: on + # 2: on, but alpha_e is determined through real parameters + diff --git a/TTBAR_8TEV_DIF/parameter.dat b/TTBAR_8TEV_DIF/parameter.dat new file mode 100644 index 00000000..eba664f9 --- /dev/null +++ b/TTBAR_8TEV_DIF/parameter.dat @@ -0,0 +1,111 @@ +########################## +# MATRIX input parameter # +########################## + +#----------------------\ +# General run settings | +#----------------------/ +process_class = pp-tt~+X # process id +E = 4000. # 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 | +#----------------/ +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 = 2 # 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 = NNPDF40_nnlo_as_01180 # 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 = NNPDF40_nnlo_as_01180 # 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 = 1.e-3 # 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 = 12 # 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) + diff --git a/TTBAR_8TEV_DIF/readme.md b/TTBAR_8TEV_DIF/readme.md new file mode 100644 index 00000000..b71278d5 --- /dev/null +++ b/TTBAR_8TEV_DIF/readme.md @@ -0,0 +1,16 @@ +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 +``` + diff --git a/TTBAR_8TEV_DIF/specify.scales.cxx b/TTBAR_8TEV_DIF/specify.scales.cxx new file mode 100644 index 00000000..47cf79de --- /dev/null +++ b/TTBAR_8TEV_DIF/specify.scales.cxx @@ -0,0 +1,106 @@ +{ + if (sd == 1){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 2){ + // transverse mass of the ttx system + double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); + temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); + } + else if (sd == 3){ + // geometric average of transverse masses of and tx - ET + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = sqrt(ET_t * ET_tx); + } + else if (sd == 4){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 5){ + // sum of transverse masses of t and tx and of jets - HTp + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ + temp_mu_central += PARTICLE("all")[i_l].ET; + } + // Alternatively: + // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ + } + else if (sd == 6){ + // sum of transverse masses of t and tx - HTint + double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); + double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); + temp_mu_central = xET_t + xET_tx; + } + else if (sd == 7){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 8){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + + else if (sd == 9){ + // transverse mass of the leading top-jet - mT_t1 + temp_mu_central = PARTICLE("tjet")[0].ET; + } + else if (sd == 10){ + // transverse mass of the subleading top-jet - mT_t2 + temp_mu_central = PARTICLE("tjet")[1].ET; + } + else if (sd == 11){ + // sum of transverse energy of t and tx at central mtop- HT / 4 + double pT2_t = PARTICLE("top")[0].pT2; + double pT2_tx = PARTICLE("atop")[0].pT2; + double ET_t = sqrt(pow(172.5, 2) + pT2_t); + double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + else if (sd == 70){ + // sum of transverse energies of t and tx varying mtop - HT / 4 + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = 0.25 * (ET_t + ET_tx); + } + + + /* + // these were only required because of some former bug in the implementation: + else if (sd == 11){ + // invariant mass m of the ttx system + temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); + } + else if (sd == 14){ + // sum of transverse masses of t and tx - HT + double ET_t = PARTICLE("top")[0].ET; + double ET_tx = PARTICLE("atop")[0].ET; + temp_mu_central = ET_t + ET_tx; + } + else if (sd == 17){ + // transverse mass of t - mT_t + temp_mu_central = PARTICLE("top")[0].ET; + } + else if (sd == 18){ + // transverse mass of t - mT_tx + temp_mu_central = PARTICLE("atop")[0].ET; + } + */ + else{ + assert(false && "Scale setting not defined. Reset dynamical scale to different value."); + } +} + + /* +if (sd == 1){ + fourvector Q; + for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} + temp_mu_central = Q.m(); +} + */ From cab0cd30b55d3d14d3139705e1984ef2513712dc Mon Sep 17 00:00:00 2001 From: jacoterh <54140851+jacoterh@users.noreply.github.com> Date: Thu, 24 Jul 2025 13:41:03 +0100 Subject: [PATCH 2/2] Updating MATRIX runcards to include the CMS TTBAR 2L legacy measurement that got added after --- ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat | 17 -- ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat | 163 ------------ ATLAS_TTBAR_13TEV_HADR_DIF/model.dat | 61 ----- ATLAS_TTBAR_13TEV_HADR_DIF/readme.md | 16 -- ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx | 106 -------- .../dddistribution.dat | 25 -- .../distribution.dat | 177 ------------- .../model.dat | 61 ----- .../parameter.dat | 111 --------- .../readme.md | 16 -- .../specify.scales.cxx | 106 -------- CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat | 14 -- CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat | 135 ---------- CMS_TTBAR_5TEV_TOT_X-SEC/model.dat | 61 ----- CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat | 111 --------- CMS_TTBAR_5TEV_TOT_X-SEC/readme.md | 16 -- CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx | 106 -------- TTBAR_13TEV_DIF/13tev_dddistribution.dat | 41 +++ TTBAR_13TEV_DIF/13tev_distribution.dat | 219 ++++++++++++++++ .../{parameter.dat => 13tev_parameter.dat} | 15 +- TTBAR_13TEV_DIF/dddistribution.dat | 18 -- TTBAR_13TEV_DIF/distribution.dat | 233 ----------------- TTBAR_13TEV_DIF/model.dat | 17 +- TTBAR_13TEV_DIF/readme.md | 16 -- TTBAR_7TEV_TOT_X-SEC/dddistribution.dat | 14 -- TTBAR_7TEV_TOT_X-SEC/distribution.dat | 135 ---------- TTBAR_7TEV_TOT_X-SEC/model.dat | 61 ----- TTBAR_7TEV_TOT_X-SEC/parameter.dat | 111 --------- TTBAR_7TEV_TOT_X-SEC/readme.md | 16 -- TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx | 106 -------- TTBAR_8TEV_DIF/8tev_dddistribution.dat | 17 ++ TTBAR_8TEV_DIF/8tev_distribution.dat | 108 ++++++++ .../8tev_parameter.dat | 17 +- TTBAR_8TEV_DIF/dddistribution.dat | 28 --- TTBAR_8TEV_DIF/distribution.dat | 235 ------------------ TTBAR_8TEV_DIF/model.dat | 17 +- TTBAR_8TEV_DIF/parameter.dat | 111 --------- TTBAR_8TEV_DIF/readme.md | 5 +- 38 files changed, 420 insertions(+), 2422 deletions(-) delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/model.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/readme.md delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/model.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md delete mode 100644 ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/model.dat delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/readme.md delete mode 100644 CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx create mode 100644 TTBAR_13TEV_DIF/13tev_dddistribution.dat create mode 100644 TTBAR_13TEV_DIF/13tev_distribution.dat rename TTBAR_13TEV_DIF/{parameter.dat => 13tev_parameter.dat} (92%) delete mode 100644 TTBAR_13TEV_DIF/dddistribution.dat delete mode 100644 TTBAR_13TEV_DIF/distribution.dat delete mode 100644 TTBAR_7TEV_TOT_X-SEC/dddistribution.dat delete mode 100644 TTBAR_7TEV_TOT_X-SEC/distribution.dat delete mode 100644 TTBAR_7TEV_TOT_X-SEC/model.dat delete mode 100644 TTBAR_7TEV_TOT_X-SEC/parameter.dat delete mode 100644 TTBAR_7TEV_TOT_X-SEC/readme.md delete mode 100644 TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx create mode 100644 TTBAR_8TEV_DIF/8tev_dddistribution.dat create mode 100644 TTBAR_8TEV_DIF/8tev_distribution.dat rename ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat => TTBAR_8TEV_DIF/8tev_parameter.dat (91%) delete mode 100644 TTBAR_8TEV_DIF/dddistribution.dat delete mode 100644 TTBAR_8TEV_DIF/distribution.dat delete mode 100644 TTBAR_8TEV_DIF/parameter.dat diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat deleted file mode 100644 index 0e2c32a8..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF/dddistribution.dat +++ /dev/null @@ -1,17 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - -dddistributionname = dd_mttbar-yttbar -distributionname_1 = dd_m_ttx -distributionname_2 = dd_y_ttx - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat deleted file mode 100644 index 38bcc1f2..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF/distribution.dat +++ /dev/null @@ -1,163 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -distributionname = m_ttx -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 = y_ttx -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 - -distributionname = dd_m_ttx -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0.0:700.0:970.0:3000.0 - -distributionname = dd_y_ttx -distributiontype = absy -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0.0:0.46:0.91:1.55:2.5 - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) - -#------------------\ -# Yukawa couplings | -#------------------/ -#Block YUKAWA -# 5 4.750000e+00 # M_YB -# 6 1.730000e+02 # M_YT -# 15 1.777000e+00 # M_YTAU - -#---------------\ -# decays widths | -#---------------/ -DECAY 6 0.000000 # WT -DECAY 23 0.000000 # WZ -DECAY 24 0.000000 # WW -DECAY 25 4.070000e-03 # WH - -#-----------\ -# EW inputs | -#-----------/ -Block EWINPUTS - 1 1 # ew_scheme - determines scheme used for EW inputs - # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) - # 1: G_mu scheme (G_F above used to determine inputs) - # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) - 3 1 # use_cms - switch for the complex mass scheme - # 0: off - # 1: on - # 2: on, but alpha_e is determined through real parameters - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md b/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md deleted file mode 100644 index b71278d5..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF/readme.md +++ /dev/null @@ -1,16 +0,0 @@ -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 -``` - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx b/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx deleted file mode 100644 index 47cf79de..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF/specify.scales.cxx +++ /dev/null @@ -1,106 +0,0 @@ -{ - if (sd == 1){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 2){ - // transverse mass of the ttx system - double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); - temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); - } - else if (sd == 3){ - // geometric average of transverse masses of and tx - ET - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = sqrt(ET_t * ET_tx); - } - else if (sd == 4){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 5){ - // sum of transverse masses of t and tx and of jets - HTp - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ - temp_mu_central += PARTICLE("all")[i_l].ET; - } - // Alternatively: - // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ - } - else if (sd == 6){ - // sum of transverse masses of t and tx - HTint - double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); - double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); - temp_mu_central = xET_t + xET_tx; - } - else if (sd == 7){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 8){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - - else if (sd == 9){ - // transverse mass of the leading top-jet - mT_t1 - temp_mu_central = PARTICLE("tjet")[0].ET; - } - else if (sd == 10){ - // transverse mass of the subleading top-jet - mT_t2 - temp_mu_central = PARTICLE("tjet")[1].ET; - } - else if (sd == 11){ - // sum of transverse energy of t and tx at central mtop- HT / 4 - double pT2_t = PARTICLE("top")[0].pT2; - double pT2_tx = PARTICLE("atop")[0].pT2; - double ET_t = sqrt(pow(172.5, 2) + pT2_t); - double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - else if (sd == 70){ - // sum of transverse energies of t and tx varying mtop - HT / 4 - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - - - /* - // these were only required because of some former bug in the implementation: - else if (sd == 11){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 14){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 17){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 18){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - */ - else{ - assert(false && "Scale setting not defined. Reset dynamical scale to different value."); - } -} - - /* -if (sd == 1){ - fourvector Q; - for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} - temp_mu_central = Q.m(); -} - */ diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat deleted file mode 100644 index 49069504..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/dddistribution.dat +++ /dev/null @@ -1,25 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - -dddistributionname = dd1 -distributionname_1 = m1_ttx -distributionname_2 = ym1_ttx - -dddistributionname = dd2 -distributionname_1 = m2_ttx -distributionname_2 = ym2_ttx - -dddistributionname = dd3 -distributionname_1 = m3_ttx -distributionname_2 = ym3_ttx - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat deleted file mode 100644 index 526422dd..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/distribution.dat +++ /dev/null @@ -1,177 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -distributionname = m1_ttx -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0.0:700.0 - -distributionname = m2_ttx -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 700.0:970.0 - -distributionname = m3_ttx -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 970.0:3000.0 - -distributionname = 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 = 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 = ym3_ttx -distributiontype = absy -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0.0:0.4:0.79:2.5 - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) - -#------------------\ -# Yukawa couplings | -#------------------/ -#Block YUKAWA -# 5 4.750000e+00 # M_YB -# 6 1.730000e+02 # M_YT -# 15 1.777000e+00 # M_YTAU - -#---------------\ -# decays widths | -#---------------/ -DECAY 6 0.000000 # WT -DECAY 23 0.000000 # WZ -DECAY 24 0.000000 # WW -DECAY 25 4.070000e-03 # WH - -#-----------\ -# EW inputs | -#-----------/ -Block EWINPUTS - 1 1 # ew_scheme - determines scheme used for EW inputs - # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) - # 1: G_mu scheme (G_F above used to determine inputs) - # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) - 3 1 # use_cms - switch for the complex mass scheme - # 0: off - # 1: on - # 2: on, but alpha_e is determined through real parameters - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat deleted file mode 100644 index 4ef9159a..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/parameter.dat +++ /dev/null @@ -1,111 +0,0 @@ -########################## -# 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 | -#----------------/ -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 = 2 # 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-1 # 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 = 1.e-3 # 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 = 12 # 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 = 1 # 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) - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md deleted file mode 100644 index b71278d5..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/readme.md +++ /dev/null @@ -1,16 +0,0 @@ -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 -``` - diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx b/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx deleted file mode 100644 index 47cf79de..00000000 --- a/ATLAS_TTBAR_13TEV_HADR_DIF_MTTBAR-YTTBAR/specify.scales.cxx +++ /dev/null @@ -1,106 +0,0 @@ -{ - if (sd == 1){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 2){ - // transverse mass of the ttx system - double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); - temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); - } - else if (sd == 3){ - // geometric average of transverse masses of and tx - ET - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = sqrt(ET_t * ET_tx); - } - else if (sd == 4){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 5){ - // sum of transverse masses of t and tx and of jets - HTp - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ - temp_mu_central += PARTICLE("all")[i_l].ET; - } - // Alternatively: - // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ - } - else if (sd == 6){ - // sum of transverse masses of t and tx - HTint - double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); - double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); - temp_mu_central = xET_t + xET_tx; - } - else if (sd == 7){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 8){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - - else if (sd == 9){ - // transverse mass of the leading top-jet - mT_t1 - temp_mu_central = PARTICLE("tjet")[0].ET; - } - else if (sd == 10){ - // transverse mass of the subleading top-jet - mT_t2 - temp_mu_central = PARTICLE("tjet")[1].ET; - } - else if (sd == 11){ - // sum of transverse energy of t and tx at central mtop- HT / 4 - double pT2_t = PARTICLE("top")[0].pT2; - double pT2_tx = PARTICLE("atop")[0].pT2; - double ET_t = sqrt(pow(172.5, 2) + pT2_t); - double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - else if (sd == 70){ - // sum of transverse energies of t and tx varying mtop - HT / 4 - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - - - /* - // these were only required because of some former bug in the implementation: - else if (sd == 11){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 14){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 17){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 18){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - */ - else{ - assert(false && "Scale setting not defined. Reset dynamical scale to different value."); - } -} - - /* -if (sd == 1){ - fourvector Q; - for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} - temp_mu_central = Q.m(); -} - */ diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat deleted file mode 100644 index efacdc7d..00000000 --- a/CMS_TTBAR_5TEV_TOT_X-SEC/dddistribution.dat +++ /dev/null @@ -1,14 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - - diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat deleted file mode 100644 index dd181867..00000000 --- a/CMS_TTBAR_5TEV_TOT_X-SEC/distribution.dat +++ /dev/null @@ -1,135 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) - -#------------------\ -# Yukawa couplings | -#------------------/ -#Block YUKAWA -# 5 4.750000e+00 # M_YB -# 6 1.730000e+02 # M_YT -# 15 1.777000e+00 # M_YTAU - -#---------------\ -# decays widths | -#---------------/ -DECAY 6 0.000000 # WT -DECAY 23 0.000000 # WZ -DECAY 24 0.000000 # WW -DECAY 25 4.070000e-03 # WH - -#-----------\ -# EW inputs | -#-----------/ -Block EWINPUTS - 1 1 # ew_scheme - determines scheme used for EW inputs - # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) - # 1: G_mu scheme (G_F above used to determine inputs) - # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) - 3 1 # use_cms - switch for the complex mass scheme - # 0: off - # 1: on - # 2: on, but alpha_e is determined through real parameters - diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat b/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat deleted file mode 100644 index 13ce4f2a..00000000 --- a/CMS_TTBAR_5TEV_TOT_X-SEC/parameter.dat +++ /dev/null @@ -1,111 +0,0 @@ -########################## -# MATRIX input parameter # -########################## - -#----------------------\ -# General run settings | -#----------------------/ -process_class = pp-tt~+X # process id -E = 2510. # 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 = 0 # switch to improve statistics in tail of distributions (a factor of two slower) - - -#----------------\ -# Scale settings | -#----------------/ -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 = 2 # 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 = 1.e-3 # 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 = 12 # 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) - diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md b/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md deleted file mode 100644 index b71278d5..00000000 --- a/CMS_TTBAR_5TEV_TOT_X-SEC/readme.md +++ /dev/null @@ -1,16 +0,0 @@ -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 -``` - diff --git a/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx b/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx deleted file mode 100644 index 47cf79de..00000000 --- a/CMS_TTBAR_5TEV_TOT_X-SEC/specify.scales.cxx +++ /dev/null @@ -1,106 +0,0 @@ -{ - if (sd == 1){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 2){ - // transverse mass of the ttx system - double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); - temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); - } - else if (sd == 3){ - // geometric average of transverse masses of and tx - ET - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = sqrt(ET_t * ET_tx); - } - else if (sd == 4){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 5){ - // sum of transverse masses of t and tx and of jets - HTp - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ - temp_mu_central += PARTICLE("all")[i_l].ET; - } - // Alternatively: - // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ - } - else if (sd == 6){ - // sum of transverse masses of t and tx - HTint - double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); - double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); - temp_mu_central = xET_t + xET_tx; - } - else if (sd == 7){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 8){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - - else if (sd == 9){ - // transverse mass of the leading top-jet - mT_t1 - temp_mu_central = PARTICLE("tjet")[0].ET; - } - else if (sd == 10){ - // transverse mass of the subleading top-jet - mT_t2 - temp_mu_central = PARTICLE("tjet")[1].ET; - } - else if (sd == 11){ - // sum of transverse energy of t and tx at central mtop- HT / 4 - double pT2_t = PARTICLE("top")[0].pT2; - double pT2_tx = PARTICLE("atop")[0].pT2; - double ET_t = sqrt(pow(172.5, 2) + pT2_t); - double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - else if (sd == 70){ - // sum of transverse energies of t and tx varying mtop - HT / 4 - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - - - /* - // these were only required because of some former bug in the implementation: - else if (sd == 11){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 14){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 17){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 18){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - */ - else{ - assert(false && "Scale setting not defined. Reset dynamical scale to different value."); - } -} - - /* -if (sd == 1){ - fourvector Q; - for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} - temp_mu_central = Q.m(); -} - */ diff --git a/TTBAR_13TEV_DIF/13tev_dddistribution.dat b/TTBAR_13TEV_DIF/13tev_dddistribution.dat new file mode 100644 index 00000000..bef37f25 --- /dev/null +++ b/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 diff --git a/TTBAR_13TEV_DIF/13tev_distribution.dat b/TTBAR_13TEV_DIF/13tev_distribution.dat new file mode 100644 index 00000000..69154c74 --- /dev/null +++ b/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 \ No newline at end of file diff --git a/TTBAR_13TEV_DIF/parameter.dat b/TTBAR_13TEV_DIF/13tev_parameter.dat similarity index 92% rename from TTBAR_13TEV_DIF/parameter.dat rename to TTBAR_13TEV_DIF/13tev_parameter.dat index 457fff1f..51e01852 100644 --- a/TTBAR_13TEV_DIF/parameter.dat +++ b/TTBAR_13TEV_DIF/13tev_parameter.dat @@ -23,12 +23,14 @@ enhance_tails = 1 # switch to improve statistics in tail of distri #----------------\ # Scale settings | #----------------/ -dynamic_scale = 11 # dynamic ren./fac. scale +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 = 2 # switch for muR/muF uncertainties (0) off; (1) 7-point (default); (2) 9-point variation +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) @@ -49,10 +51,10 @@ 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 +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 = 1.e-3 # precision of NNLO QCD cross section +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 @@ -95,17 +97,16 @@ n_observed_max jet = 99 # maximal number of observed jets (with cuts abo #-----------------\ # MATRIX behavior | #-----------------/ -max_time_per_job = 12 # very rough time(hours) one main run job shall take (default: 24h) +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; +#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) - diff --git a/TTBAR_13TEV_DIF/dddistribution.dat b/TTBAR_13TEV_DIF/dddistribution.dat deleted file mode 100644 index 6fbb39f1..00000000 --- a/TTBAR_13TEV_DIF/dddistribution.dat +++ /dev/null @@ -1,18 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - - -dddistributionname = cms_lj_mttbar-yttbar -distributionname_1 = cms_lj_2d_mttbar -distributionname_2 = cms_lj_2d_yttbar - diff --git a/TTBAR_13TEV_DIF/distribution.dat b/TTBAR_13TEV_DIF/distribution.dat deleted file mode 100644 index 26315dfa..00000000 --- a/TTBAR_13TEV_DIF/distribution.dat +++ /dev/null @@ -1,233 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -# atlas-lj - -distributionname = atlas_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_pTt -distributiontype = pT -particle 1 = top 1 -binningtype = irregular -edges = 0:50:100:160:225:300:360:475:1000 - -distributionname = atlas_yt -distributiontype = absy -particle 1 = top 1 -binningtype = irregular -edges = 0:0.4:0.8:1.2:1.7:2.5 - -distributionname = atlas_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 - -# cms-2l - -distributionname = cms_2l_mttbar -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 300:380:470:620:820:1100:1500:2500 - -distributionname = cms_2l_yttbar -distributiontype = y -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = -2.6:-1.6:-1.2:-0.8:-0.4:0:0.4:0.8:1.2:1.6:2.6 - -distributionname = cms_2l_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_2l_pTt -distributiontype = pT -particle 1 = top 1 -binningtype = irregular -edges = 0:65:125:200:290:400:550 - -# cms-lj - -distributionname = cms_lj_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_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_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_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 - -distributionname = cms_lj_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_lj_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 - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) + 1 1.321003e+02 # 1/alpha_e(MZ) + 2 1.166378e-05 # G_F +111 1.370359e+02 # 1/alpha_e(mu->0) #------------------\ # Yukawa couplings | @@ -58,4 +58,3 @@ Block EWINPUTS # 0: off # 1: on # 2: on, but alpha_e is determined through real parameters - diff --git a/TTBAR_13TEV_DIF/readme.md b/TTBAR_13TEV_DIF/readme.md index b71278d5..e69de29b 100644 --- a/TTBAR_13TEV_DIF/readme.md +++ b/TTBAR_13TEV_DIF/readme.md @@ -1,16 +0,0 @@ -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 -``` - diff --git a/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat b/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat deleted file mode 100644 index efacdc7d..00000000 --- a/TTBAR_7TEV_TOT_X-SEC/dddistribution.dat +++ /dev/null @@ -1,14 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - - diff --git a/TTBAR_7TEV_TOT_X-SEC/distribution.dat b/TTBAR_7TEV_TOT_X-SEC/distribution.dat deleted file mode 100644 index dd181867..00000000 --- a/TTBAR_7TEV_TOT_X-SEC/distribution.dat +++ /dev/null @@ -1,135 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) - -#------------------\ -# Yukawa couplings | -#------------------/ -#Block YUKAWA -# 5 4.750000e+00 # M_YB -# 6 1.730000e+02 # M_YT -# 15 1.777000e+00 # M_YTAU - -#---------------\ -# decays widths | -#---------------/ -DECAY 6 0.000000 # WT -DECAY 23 0.000000 # WZ -DECAY 24 0.000000 # WW -DECAY 25 4.070000e-03 # WH - -#-----------\ -# EW inputs | -#-----------/ -Block EWINPUTS - 1 1 # ew_scheme - determines scheme used for EW inputs - # 0: alpha_e_0 scheme (alpha_e(mu->0) above used to determine inputs) - # 1: G_mu scheme (G_F above used to determine inputs) - # 2: alpha_e_MZ scheme (alpha_e(MZ) above used to determine inputs) - 3 1 # use_cms - switch for the complex mass scheme - # 0: off - # 1: on - # 2: on, but alpha_e is determined through real parameters - diff --git a/TTBAR_7TEV_TOT_X-SEC/parameter.dat b/TTBAR_7TEV_TOT_X-SEC/parameter.dat deleted file mode 100644 index 2331d4d5..00000000 --- a/TTBAR_7TEV_TOT_X-SEC/parameter.dat +++ /dev/null @@ -1,111 +0,0 @@ -########################## -# MATRIX input parameter # -########################## - -#----------------------\ -# General run settings | -#----------------------/ -process_class = pp-tt~+X # process id -E = 3500. # 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 = 0 # switch to improve statistics in tail of distributions (a factor of two slower) - - -#----------------\ -# Scale settings | -#----------------/ -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 = 2 # 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 = 1.e-3 # 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 = 12 # 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) - diff --git a/TTBAR_7TEV_TOT_X-SEC/readme.md b/TTBAR_7TEV_TOT_X-SEC/readme.md deleted file mode 100644 index b71278d5..00000000 --- a/TTBAR_7TEV_TOT_X-SEC/readme.md +++ /dev/null @@ -1,16 +0,0 @@ -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 -``` - diff --git a/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx b/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx deleted file mode 100644 index 47cf79de..00000000 --- a/TTBAR_7TEV_TOT_X-SEC/specify.scales.cxx +++ /dev/null @@ -1,106 +0,0 @@ -{ - if (sd == 1){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 2){ - // transverse mass of the ttx system - double m = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - double pT = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).pT(); - temp_mu_central = sqrt(pow(m, 2) + pow(pT, 2)); - } - else if (sd == 3){ - // geometric average of transverse masses of and tx - ET - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = sqrt(ET_t * ET_tx); - } - else if (sd == 4){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 5){ - // sum of transverse masses of t and tx and of jets - HTp - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - for (int i_l = 5; i_l < PARTICLE("all").size(); i_l++){ - temp_mu_central += PARTICLE("all")[i_l].ET; - } - // Alternatively: - // for (int i_l = 3; i_l < PARTICLE("all").size(); i_l++){ - } - else if (sd == 6){ - // sum of transverse masses of t and tx - HTint - double xET_t = sqrt(msi->M2_t / 4 + PARTICLE("top")[0].pT2); - double xET_tx = sqrt(msi->M2_t / 4 + PARTICLE("atop")[0].pT2); - temp_mu_central = xET_t + xET_tx; - } - else if (sd == 7){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 8){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - - else if (sd == 9){ - // transverse mass of the leading top-jet - mT_t1 - temp_mu_central = PARTICLE("tjet")[0].ET; - } - else if (sd == 10){ - // transverse mass of the subleading top-jet - mT_t2 - temp_mu_central = PARTICLE("tjet")[1].ET; - } - else if (sd == 11){ - // sum of transverse energy of t and tx at central mtop- HT / 4 - double pT2_t = PARTICLE("top")[0].pT2; - double pT2_tx = PARTICLE("atop")[0].pT2; - double ET_t = sqrt(pow(172.5, 2) + pT2_t); - double ET_tx = sqrt(pow(172.5, 2) + pT2_tx); - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - else if (sd == 70){ - // sum of transverse energies of t and tx varying mtop - HT / 4 - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = 0.25 * (ET_t + ET_tx); - } - - - /* - // these were only required because of some former bug in the implementation: - else if (sd == 11){ - // invariant mass m of the ttx system - temp_mu_central = (PARTICLE("top")[0].momentum + PARTICLE("atop")[0].momentum).m(); - } - else if (sd == 14){ - // sum of transverse masses of t and tx - HT - double ET_t = PARTICLE("top")[0].ET; - double ET_tx = PARTICLE("atop")[0].ET; - temp_mu_central = ET_t + ET_tx; - } - else if (sd == 17){ - // transverse mass of t - mT_t - temp_mu_central = PARTICLE("top")[0].ET; - } - else if (sd == 18){ - // transverse mass of t - mT_tx - temp_mu_central = PARTICLE("atop")[0].ET; - } - */ - else{ - assert(false && "Scale setting not defined. Reset dynamical scale to different value."); - } -} - - /* -if (sd == 1){ - fourvector Q; - for (int i_p = 3; i_p < 5; i_p++){Q = Q + esi->p_parton[i_a][i_p];} - temp_mu_central = Q.m(); -} - */ diff --git a/TTBAR_8TEV_DIF/8tev_dddistribution.dat b/TTBAR_8TEV_DIF/8tev_dddistribution.dat new file mode 100644 index 00000000..a65f10ae --- /dev/null +++ b/TTBAR_8TEV_DIF/8tev_dddistribution.dat @@ -0,0 +1,17 @@ +################ LHC TTBAR 8 TEV ############## + +########### 2d distributions + +# cms_2l + +dddistributionname = CMS_TTBAR_8TEV_2L_DIF_PTT-YT +distributionname_1 = cms_2l_yt +distributionname_2 = cms_2l_pTt + +dddistributionname = CMS_TTBAR_8TEV_2L_DIF_MTTBAR-YT +distributionname_1 = cms_2l_mttbar +distributionname_2 = cms_2l_yt + +dddistributionname = CMS_TTBAR_8TEV_2L_DIF_MTTBAR-YTTBAR +distributionname_1 = cms_2l_mttbar +distributionname_2 = cms_2l_yttbar diff --git a/TTBAR_8TEV_DIF/8tev_distribution.dat b/TTBAR_8TEV_DIF/8tev_distribution.dat new file mode 100644 index 00000000..c8c1a8e4 --- /dev/null +++ b/TTBAR_8TEV_DIF/8tev_distribution.dat @@ -0,0 +1,108 @@ +################ LHC TTBAR 8 TEV ############## + +# keys: catr- convert absolute to regular + +########### 1d distributions + +# atlas_2l + +distributionname = ATLAS_TTBAR_8TEV_2L_DIF_MTTBAR +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 250:450:570:700:850:1000:2700 + +distributionname = ATLAS_TTBAR_8TEV_2L_DIF_YTTBAR +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.4:0.8:1.2:2:2.8 + +# atlas_lj + +distributionname = ATLAS_TTBAR_8TEV_LJ_DIF_MTTBAR +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 345:400:470:550:650:800:1100:1600 + +distributionname = ATLAS_TTBAR_8TEV_LJ_DIF_YTTBAR +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.3:0.6:0.9:1.3:2.5 + +distributionname = ATLAS_TTBAR_8TEV_LJ_DIF_YT +distributiontype = absy +particle 1 = top 1 +binningtype = irregular +edges = 0:0.4:0.8:1.2:1.6:2.5 + +distributionname = ATLAS_TTBAR_8TEV_LJ_DIF_PTT +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:60:100:150:200:260:320:400:500 + +# cms_lj + +distributionname = CMS_TTBAR_8TEV_LJ_DIF_MTTBAR +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 345:400:470:550:650:800:1100:1600 + +# (convert absolute to regular) +distributionname = CMS_TTBAR_8TEV_LJ_DIF_YTTBAR_catr +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.3:0.6:0.9:1.3:2.5 + +distributionname = CMS_TTBAR_8TEV_LJ_DIF_YT +distributiontype = y +particle 1 = top 1 +binningtype = irregular +edges = -2.5:-1.6:-1.2:-0.8:-0.4:0:0.4:0.8:1.2:1.6:2.5 + +distributionname = CMS_TTBAR_8TEV_LJ_DIF_PTT +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:60:100:150:200:260:320:400:500 + +########### 1d distributions of the 2d distributions + +# cms_2l + +distributionname = cms_2l_mttbar +distributiontype = m +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 340:400:500:650:1500 + +distributionname = cms_2l_yttbar +distributiontype = absy +particle 1 = top 1 +particle 1 = atop 1 +binningtype = irregular +edges = 0:0.35:0.75:1.15:2.5 + +distributionname = cms_2l_yt +distributiontype = absy +particle 1 = top 1 +binningtype = irregular +edges = 0:0.35:0.85:1.45:2.5 + +distributionname = cms_2l_pTt +distributiontype = pT +particle 1 = top 1 +binningtype = irregular +edges = 0:80:150:250:600 diff --git a/ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat b/TTBAR_8TEV_DIF/8tev_parameter.dat similarity index 91% rename from ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat rename to TTBAR_8TEV_DIF/8tev_parameter.dat index ae3eda34..26f8150b 100644 --- a/ATLAS_TTBAR_13TEV_HADR_DIF/parameter.dat +++ b/TTBAR_8TEV_DIF/8tev_parameter.dat @@ -6,7 +6,7 @@ # General run settings | #----------------------/ process_class = pp-tt~+X # process id -E = 6500. # energy per beam +E = 4000. # 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) @@ -23,12 +23,14 @@ enhance_tails = 1 # switch to improve statistics in tail of distri #----------------\ # Scale settings | #----------------/ -dynamic_scale = 11 # dynamic ren./fac. scale +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 = 2 # switch for muR/muF uncertainties (0) off; (1) 7-point (default); (2) 9-point variation +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) @@ -52,7 +54,7 @@ NLO_subtraction_method = 1 # switch to use (2) qT subtraction (1) Catani-Se 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 = 1.e-3 # precision of NNLO QCD cross section +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 @@ -95,17 +97,16 @@ n_observed_max jet = 99 # maximal number of observed jets (with cuts abo #-----------------\ # MATRIX behavior | #-----------------/ -max_time_per_job = 12 # very rough time(hours) one main run job shall take (default: 24h) +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 = 1 # 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; +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) - diff --git a/TTBAR_8TEV_DIF/dddistribution.dat b/TTBAR_8TEV_DIF/dddistribution.dat deleted file mode 100644 index cd7fb948..00000000 --- a/TTBAR_8TEV_DIF/dddistribution.dat +++ /dev/null @@ -1,28 +0,0 @@ -####################################################### -# MATRIX double differential distributions definition # -####################################################### -# -# In this file you can customize the distributions created during the run (examples below) -# -# They are simply defined as a combination of two single differential distributions that -# must be defined in the file 'distribution.dat' -# -#----------------------\ -# define distributions | -#----------------------/ - - -#pT distribution of the ttx system in rapidity bins. Binning defined in 'distribution.dat' - -dddistributionname = pTt-yt -distributionname_1 = cms_2l_yt -distributionname_2 = cms_2l_pTt - -dddistributionname = mttbar-yt -distributionname_1 = cms_2l_mttbar -distributionname_2 = cms_2l_yt - -dddistributionname = mttbar-yttbar -distributionname_1 = cms_2l_mttbar -distributionname_2 = cms_2l_yttbar - diff --git a/TTBAR_8TEV_DIF/distribution.dat b/TTBAR_8TEV_DIF/distribution.dat deleted file mode 100644 index a04a35c9..00000000 --- a/TTBAR_8TEV_DIF/distribution.dat +++ /dev/null @@ -1,235 +0,0 @@ -################################## -# MATRIX distribution definition # -################################## -# -# In this file you can customize the distributions created during the run (examples below) -# please read the INSTRUCTIONS at the END OF THIS FILE... -# -#------\ -# Info | -#------/ -# Total rates and jet multiplicities (possibly within cuts) will automatically be included -# Add/remove arbitrary distribution-blocks, but always add/remove a full block. -# -#----------------------\ -# define distributions | -#----------------------/ - -# atlas_2l - -distributionname = atlas_2l_mttbar -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 250:450:570:700:850:1000:2700 - -distributionname = atlas_2l_yttbar -distributiontype = absy -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0:0.4:0.8:1.2:2:2.8 - -# atlas_lj - -distributionname = atlas_lj_mttbar -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 345:400:470:550:650:800:1100:1600 - -distributionname = atlas_lj_yttbar -distributiontype = absy -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0:0.3:0.6:0.9:1.3:2.5 - -distributionname = atlas_lj_yt -distributiontype = absy -particle 1 = top 1 -binningtype = irregular -edges = 0:0.4:0.8:1.2:1.6:2.5 - -distributionname = atlas_lj_pTt -distributiontype = pT -particle 1 = top 1 -binningtype = irregular -edges = 0:60:100:150:200:260:320:400:500 - -# cms_2l - -distributionname = cms_2l_mttbar -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 340:400:500:650:1500 - -distributionname = cms_2l_yttbar -distributiontype = absy -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 0:0.35:0.75:1.15:2.5 - -distributionname = cms_2l_yt -distributiontype = absy -particle 1 = top 1 -binningtype = irregular -edges = 0:0.35:0.85:1.45:2.5 - -distributionname = cms_2l_pTt -distributiontype = pT -particle 1 = top 1 -binningtype = irregular -edges = 0:80:150:250:600 - -# cms_lj - -distributionname = cms_lj_mttbar -distributiontype = m -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = 345:400:470:550:650:800:1100:1600 - -distributionname = cms_lj_yttbar -distributiontype = y -particle 1 = top 1 -particle 1 = atop 1 -binningtype = irregular -edges = -2.5:-1.3:-0.9:-0.6:-0.3:0:0.3:0.6:0.9:1.3:2.5 - -distributionname = cms_lj_yt -distributiontype = y -particle 1 = top 1 -binningtype = irregular -edges = -2.5:-1.6:-1.2:-0.8:-0.4:0:0.4:0.8:1.2:1.6:2.5 - -distributionname = cms_lj_pTt -distributiontype = pT -particle 1 = top 1 -binningtype = irregular -edges = 0:60:100:150:200:260:320:400:500 - -#--------\ -# Syntax | -#--------/ -# "distributionname" -# starts new distribution-block; user-defined, unique label for identification at end of run -# -# -# "distributiontype" -# specifies observable to be binned (some require more than one particle) -# -# pre-defined types: (custom definitions not yet supported) -# pT -- transverse momentum of "particle 1" -# (scalar sum of pTs if more than one "particle" given) -# m -- invariant mass of "particle 1" -# dm -- invariant-mass difference between "particle 1" and "particle 2" -# absdm -- absolute invariant-mass difference between "particle 1" and "particle 2" -# mmin -- minimal invariant mass of "particle 1" and "particle 2" -# mmax -- maximal invariant mass of "particle 1" and "particle 2" -# y -- rapidy of "particle 1" -# absy -- absolute rapidy of "particle 1" -# dy -- rapidy difference between "particle 1" and "particle 2" -# absdy -- absolute rapidy difference between "particle 1" and "particle 2" -# dabsy -- difference between absolute rapidities of"particle 1" and "particle 2" -# absdabsy -- absolute difference between absolute rapidities of "particle 1" -# and "particle 2" -# eta -- pseudo-rapidy of "particle 1" -# abseta -- absolute pseudo-rapidy of "particle 1" -# deta -- pseudo-rapidy difference between "particle 1" and "particle 2" -# absdeta -- absolute pseud-rapidy difference between "particle 1" and "particle 2" -# dabseta -- difference between absolute pseudo-rapidities of"particle 1" -# and "particle 2" -# absdabseta -- absolute difference between absolute pseudo-rapidities of "particle 1" -# and "particle 2" -# phi -- transverse-angle (phi) of "particle 1", or phi difference between -# "particle 1" and "particle 2" if two particles are given -# dR -- distance (in y-phi-plane) between "particle 1" and "particle 2" -# dReta -- distance (in eta-phi-plane) between "particle 1" and "particle 2" -# ET -- transverse mass sqrt(m^2+pT^2) of "particle 1" -# (scalar sum of ETs if more than one "particle" given) -# mT -- transverse mass: ET of "particle 1", if only one particle given -# experimental definition with missing energy, if more than one particle -# given; all neutrinos (missing energie) go into "particle 1" all other -# particles in "particle 2,3,4,..." -# -# special distribution types: -# muR -- Distribution of scale muR set in computation (no "particle" definition) -# muF -- Distribution of scale muF set in computation (no "particle" definition) -# pTveto -- Cumulative cross section sigma(pT0) + 1 1.321003e+02 # 1/alpha_e(MZ) + 2 1.166378e-05 # G_F +111 1.370359e+02 # 1/alpha_e(mu->0) #------------------\ # Yukawa couplings | @@ -58,4 +58,3 @@ Block EWINPUTS # 0: off # 1: on # 2: on, but alpha_e is determined through real parameters - diff --git a/TTBAR_8TEV_DIF/parameter.dat b/TTBAR_8TEV_DIF/parameter.dat deleted file mode 100644 index eba664f9..00000000 --- a/TTBAR_8TEV_DIF/parameter.dat +++ /dev/null @@ -1,111 +0,0 @@ -########################## -# MATRIX input parameter # -########################## - -#----------------------\ -# General run settings | -#----------------------/ -process_class = pp-tt~+X # process id -E = 4000. # 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 | -#----------------/ -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 = 2 # 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 = NNPDF40_nnlo_as_01180 # 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 = NNPDF40_nnlo_as_01180 # 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 = 1.e-3 # 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 = 12 # 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) - diff --git a/TTBAR_8TEV_DIF/readme.md b/TTBAR_8TEV_DIF/readme.md index b71278d5..d0868424 100644 --- a/TTBAR_8TEV_DIF/readme.md +++ b/TTBAR_8TEV_DIF/readme.md @@ -1,5 +1,4 @@ -These runcards were used to compute top pair production processes for NNPDF4.1 with MATRIX at NNLO except for the -following changes: +These runcards were used to compute top pair production processes for NNPDF4.1 with MATRIX at NNLO. * 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)`, @@ -9,7 +8,7 @@ following changes: ``` 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` +To reproduce the NNPDF4.1 predictions, the only thing to do in `model.dat`is to make sure the top mass is set to ``` 6 172.5 # M_t ```