From e6923607ee9c70c6e4867d2b6bef6b20e289e492 Mon Sep 17 00:00:00 2001 From: Andre Govinda Stahl Leiton Date: Tue, 17 Dec 2024 14:48:41 +0100 Subject: [PATCH] Add Madgraph for HINPbPbSpring23 and HINppWinter24 central campaigns --- .../HINPbPbSpring23/HINPbPbSpring23.json | 16 ++ .../MadGraph5_aMCatNLO/ModelParams/massW.dat | 22 ++ .../MadGraph5_aMCatNLO/ModelParams/massZ.dat | 22 ++ .../Templates/LO_run_card.dat | 227 ++++++++++++++++++ .../Templates/NLO_run_card.dat | 174 ++++++++++++++ .../Templates/NLO_run_card_madstr.dat | 188 +++++++++++++++ Campaigns/HINppWinter24/HINppWinter24.json | 16 ++ .../MadGraph5_aMCatNLO/ModelParams/massW.dat | 22 ++ .../MadGraph5_aMCatNLO/ModelParams/massZ.dat | 22 ++ .../Templates/LO_run_card.dat | 227 ++++++++++++++++++ .../Templates/NLO_run_card.dat | 174 ++++++++++++++ .../Templates/NLO_run_card_madstr.dat | 188 +++++++++++++++ .../TT-2Jets_amcatnloFXFX-pythia8.json | 26 ++ ...Jets_amcatnloFXFX-pythia8_madspin_card.dat | 13 + ...T-2Jets_amcatnloFXFX-pythia8_proc_card.dat | 12 + 15 files changed, 1349 insertions(+) create mode 100644 Campaigns/HINPbPbSpring23/HINPbPbSpring23.json create mode 100644 Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massW.dat create mode 100644 Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massZ.dat create mode 100644 Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/LO_run_card.dat create mode 100644 Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat create mode 100644 Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat create mode 100644 Campaigns/HINppWinter24/HINppWinter24.json create mode 100644 Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massW.dat create mode 100644 Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massZ.dat create mode 100644 Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/LO_run_card.dat create mode 100644 Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat create mode 100644 Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat create mode 100644 Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8.json create mode 100644 Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_madspin_card.dat create mode 100644 Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_proc_card.dat diff --git a/Campaigns/HINPbPbSpring23/HINPbPbSpring23.json b/Campaigns/HINPbPbSpring23/HINPbPbSpring23.json new file mode 100644 index 00000000..4c81b133 --- /dev/null +++ b/Campaigns/HINPbPbSpring23/HINPbPbSpring23.json @@ -0,0 +1,16 @@ +{ + "chain": "chain_HINPbPbSpring23wmLHEGSHIMix_flowHINPbPbSpring23Digi_flowHINPbPbSpring23Reco_flowHINPbPbSpring23MiniAOD", + "model_params_vars": {}, + "template_vars": { + + }, + "fragment_vars": { + "concurrentLHEforMadGraph5_aMCatNLO": "False", + "concurrentLHEforPowheg": "True", + "concurrentPS": "Concurrent" + }, + "tune": "CP5", + "beam": 2681, + "gridpack_directory": "/eos/cms/store/group/phys_generator/cvmfs/gridpacks/PdmV/HINPbPbSpring23" +} + diff --git a/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massW.dat b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massW.dat new file mode 100644 index 00000000..41a8eb41 --- /dev/null +++ b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massW.dat @@ -0,0 +1,22 @@ +# These inputs are based on fixing the PDG vlaues of mZ, mW, and G_fermi +# The MC programs use a constant width scheme, whereas PDG values are sqrt(s) dependent. +# We correct for this following Sect. 3.1 in https://arxiv.org/pdf/1606.02330.pdf + +set param_card SMINPUTS 1 132.30648214177197 +set param_card SMINPUTS 2 1.1663787e-05 + +# Z boson +set param_card mass 23 91.153509740726733 +set param_card DECAY 23 2.4932018986110700 + +# W boson +set param_card mass 24 80.351812293789408 +set param_card DECAY 24 2.0904310808144846 + +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCTopWGSummaryPlots#Top_Quark_Mass +# CMS combined result (Sep 2015) + +# top quark +set param_card mass 6 172.5 +set param_card DECAY 6 1.311 +set param_card yukawa 6 172.5 diff --git a/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massZ.dat b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massZ.dat new file mode 100644 index 00000000..3dfd7695 --- /dev/null +++ b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/ModelParams/massZ.dat @@ -0,0 +1,22 @@ +# These inputs are based on fixing the PDG vlaues of mZ, mW, and G_fermi +# The MC programs use a constant width scheme, whereas PDG values are sqrt(s) dependent. +# We correct for this following Sect. 3.1 in https://arxiv.org/pdf/1606.02330.pdf + +set param_card SMINPUTS 1 128.82531590804655 +set param_card SMINPUTS 2 1.1663787e-05 + +# Z boson +set param_card mass 23 91.153509740726733 +set param_card DECAY 23 2.4932018986110700 + +# W boson +set param_card mass 24 79.906853549493746 +set param_card DECAY 24 2.085 + +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCTopWGSummaryPlots#Top_Quark_Mass +# CMS combined result (Sep 2015) + +# top quark +set param_card mass 6 172.5 +set param_card DECAY 6 1.311 +set param_card yukawa 6 172.5 diff --git a/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/LO_run_card.dat b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/LO_run_card.dat new file mode 100644 index 00000000..b64eee3e --- /dev/null +++ b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/LO_run_card.dat @@ -0,0 +1,227 @@ +#********************************************************************* +# MadGraph5_aMC@NLO * +# * +# run_card.dat MadEvent * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a '# ' are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# To display more options, you can type the command: * +# update full_run_card * +#********************************************************************* +# +#********************************************************************* +# Tag name for the run (one word) * +#********************************************************************* + tag_1 = run_tag ! name of the run +#********************************************************************* +# Number of events and rnd seed * +# Warning: Do not generate more than 1M events in a single run * +#********************************************************************* + 3000 = nevents ! Number of unweighted events requested + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron, 4=photon from muon * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + $ebeam1 = ebeam1 ! beam 1 total energy in GeV + $ebeam2 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_nPDF_SETS = lhaid ! if pdlabel=lhapdf, this is the lhapdf number + $DEFAULT_nPDF_MEMBERS = reweight_PDF +# To see heavy ion options: type "update ion_pdf" +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.188 = scale ! fixed ren scale + 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 + 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 + -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices + 1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Type and output format +#********************************************************************* + True = gridpack !True = setting up the grid pack + -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written) + average = event_norm ! average/sum. Normalization of the weight in the LHEF +#********************************************************************* +# Matching parameter (MLM only) +#********************************************************************* + $ickkw = ickkw ! 0 no matching, 1 MLM + 1.0 = alpsfact ! scale factor for QCD emission vx + False = chcluster ! cluster only according to channel diag + $maxjetflavor = asrwgtflavor ! highest quark flavor for a_s reweight + True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 + ! (turn off for VBF and single top processes) + $xqcut = xqcut ! minimum kt jet measure between partons + +#*********************************************************************** +# Turn on either the ktdurham or ptlund cut to activate * +# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] * +#*********************************************************************** + -1.0 = ktdurham + 0.4 = dparameter + -1.0 = ptlund + 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above + +#********************************************************************* +# +#********************************************************************* +# Phase-Space Optimization strategy (basic options) +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. + ! 0: sum over helicity, 1: importance sampling + 1 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx) + ! 1 is old strategy (using amp square) + ! 2 is new strategy (using only the denominator) +# To see advanced option for Phase-Space optimization: type "update psoptim" +#********************************************************************* +# Generation bias, check the wiki page below for more information: * +# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' * +#********************************************************************* + None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-] + {} = bias_parameters ! Specifies the parameters of the module. +# +#******************************* +# Parton level cuts definition * +#******************************* +# +# +#********************************************************************* +# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay +#********************************************************************* + 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 0.0 = pta ! minimum pt for the photons + 0.0 = ptj ! minimum pt for the jets + 0.0 = ptl ! minimum pt for the charged leptons + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptlmax ! maximum pt for the charged leptons + {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle + {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) +# +# For display option for energy cut in the partonic center of mass frame type 'update ecut' +# +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + -1.0 = etaa ! max rap for the photons + -1.0 = etaj ! max rap for the jets + -1.0 = etal ! max rap for the charged leptons + 0.0 = etalmin ! main rap for the charged leptons + {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle + {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5}) +#********************************************************************* +# Minimum and maximum DeltaR distance * +#********************************************************************* + 0.0 = draj ! min distance between gamma and jet + 0.0 = dral ! min distance between gamma and lepton + 0.0 = drjj ! min distance between jets + 0.0 = drll ! min distance between leptons + 0.0 = drjl ! min distance between jet and lepton + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = dralmax ! maxdistance between gamma and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drllmax ! max distance between leptons + -1.0 = drjlmax ! max distance between jet and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 0.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair + -1.0 = mmjjmax ! max invariant mass of a jet pair + -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair + {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250}) + {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only + ! to pairs of particle/antiparticle and not to pairs of the same pdg codes. + #********************************************************************* + # Minimum and maximum invariant mass for all letpons * + #********************************************************************* + 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) + -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) + #********************************************************************* + # Minimum and maximum pt for 4-momenta sum of leptons / neutrino * + # for pair of lepton includes only same flavor, opposite charge + #********************************************************************* + 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) + -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) +#********************************************************************* +# Inclusive cuts * +#********************************************************************* + 0.0 = xptj ! minimum pt for at least one jet + 0.0 = xptl ! minimum pt for at least one charged lepton + #********************************************************************* + # Control the pt's of the jets sorted by pt * + #********************************************************************* + 0.0 = ptj1min ! minimum pt for the leading jet in pt + 0.0 = ptj2min ! minimum pt for the second jet in pt + 0.0 = ptj3min ! minimum pt for the third jet in pt + 0.0 = ptj4min ! minimum pt for the fourth jet in pt + -1.0 = ptj1max ! maximum pt for the leading jet in pt + -1.0 = ptj2max ! maximum pt for the second jet in pt + -1.0 = ptj3max ! maximum pt for the third jet in pt + -1.0 = ptj4max ! maximum pt for the fourth jet in pt + 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts + #********************************************************************* + # Control the pt's of leptons sorted by pt * + #********************************************************************* + 0.0 = ptl1min ! minimum pt for the leading lepton in pt + 0.0 = ptl2min ! minimum pt for the second lepton in pt + -1.0 = ptl1max ! maximum pt for the leading lepton in pt + -1.0 = ptl2max ! maximum pt for the second lepton in pt + #********************************************************************* + # Control the Ht(k)=Sum of k leading jets * + #********************************************************************* + 0.0 = htjmin ! minimum jet HT=Sum(jet pt) + -1.0 = htjmax ! maximum jet HT=Sum(jet pt) + 0.0 = ihtmin !inclusive Ht for all partons (including b) + -1.0 = ihtmax !inclusive Ht for all partons (including b) + 0.0 = ht2min ! minimum Ht for the two leading jets + 0.0 = ht3min ! minimum Ht for the three leading jets + 0.0 = ht4min ! minimum Ht for the four leading jets + -1.0 = ht2max ! maximum Ht for the two leading jets + -1.0 = ht3max ! maximum Ht for the three leading jets + -1.0 = ht4max ! maximum Ht for the four leading jets + #********************************************************************* + # WBF cuts * + #********************************************************************* + 0.0 = xetamin ! minimum rapidity for two jets in the WBF case + 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* + $maxjetflavor = maxjetflavor ! Maximum jet pdg code +#********************************************************************* +# +#********************************************************************* +# Store info for systematics studies * +# WARNING: Do not use for interference type of computation * +#********************************************************************* + True = use_syst ! Enable systematics studies + +#************************************** +# Manually added by Sihyun Jeon +#************************************** + False = cut_decays ! Cut decay products + False = pdfwgt ! for ickkw=1, perform pdf reweighting diff --git a/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat new file mode 100644 index 00000000..b1adb7be --- /dev/null +++ b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat @@ -0,0 +1,174 @@ +#*********************************************************************** +# MadGraph5_aMC@NLO * +# * +# run_card.dat aMC@NLO * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a hash (#) are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# Some of the values of variables can be list. These can either be * +# comma or space separated. * +# * +# To display additional parameter, you can use the command: * +# update to_full * +#*********************************************************************** +# +#******************* +# Running parameters +#******************* +# +#*********************************************************************** +# Tag name for the run (one word) * +#*********************************************************************** + tag_1 = run_tag ! name of the run +#*********************************************************************** +# Number of LHE events (and their normalization) and the required * +# (relative) accuracy on the Xsec. * +# These values are ignored for fixed order runs * +#*********************************************************************** + 5000 = nevents ! Number of unweighted events requested + 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents) + -1 = nevt_job! Max number of events per job in event generation. + ! (-1= no split). +#*********************************************************************** +# Normalize the weights of LHE events such that they sum or average to * +# the total cross section * +#*********************************************************************** + average = event_norm ! valid settings: average, sum, bias +#*********************************************************************** +# Number of points per itegration channel (ignored for aMC@NLO runs) * +#*********************************************************************** + 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the + ! number of points and iter. below) +# These numbers are ignored except if req_acc_FO is equal to -1 + 5000 = npoints_FO_grid ! number of points to setup grids + 4 = niters_FO_grid ! number of iter. to setup grids + 10000 = npoints_FO ! number of points to compute Xsec + 6 = niters_FO ! number of iter. to compute Xsec +#*********************************************************************** +# Random number seed * +#*********************************************************************** + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#*********************************************************************** +# Collider type and energy * +#*********************************************************************** + 1 = lpp1 ! beam 1 type (0 = no PDF) + 1 = lpp2 ! beam 2 type (0 = no PDF) + $ebeam1 = ebeam1 ! beam 1 energy in GeV + $ebeam2 = ebeam2 ! beam 2 energy in GeV +#*********************************************************************** +# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * +#*********************************************************************** + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_nPDF_SETS = lhaid +#*********************************************************************** +# Include the NLO Monte Carlo subtr. terms for the following parton * +# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * +# WARNING: PYTHIA6PT works only for processes without FSR!!!! * +#*********************************************************************** + $parton_shower = parton_shower + 1.0 = shower_scale_factor ! multiply default shower starting + ! scale by this factor +#*********************************************************************** +# Renormalization and factorization scales * +# (Default functional form for the non-fixed scales is the sum of * +# the transverse masses divided by two of all final state particles * +# and partons. This can be changed in SubProcesses/set_scales.f or via * +# dynamical_scale_choice option) * +#*********************************************************************** + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.118 = muR_ref_fixed ! fixed ren reference scale + 91.118 = muF_ref_fixed ! fixed fact reference scale + -1 = dynamical_scale_choice ! Choose one (or more) of the predefined + ! dynamical choices. Can be a list; scale choices beyond the + ! first are included via reweighting + 1.0 = muR_over_ref ! ratio of current muR over reference muR + 1.0 = muF_over_ref ! ratio of current muF over reference muF +#*********************************************************************** +# Reweight variables for scale dependence and PDF uncertainty * +#*********************************************************************** + 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting + 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting + True = reweight_scale ! Reweight to get scale variation using the + ! rw_rscale and rw_fscale factors. Should be a list of + ! booleans of equal length to dynamical_scale_choice to + ! specify for which choice to include scale dependence. +#*********************************************************************** +# Store reweight information in the LHE file for off-line model- * +# parameter reweighting at NLO+PS accuracy * +#*********************************************************************** + True = store_rwgt_info ! Store info for reweighting in LHE file +#*********************************************************************** +# ickkw parameter: * +# 0: No merging * +# 3: FxFx Merging - WARNING! Applies merging only at the hard-event * +# level. After showering an MLM-type merging should be applied as * +# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * +# 4: UNLOPS merging (with pythia8 only). No interface from within * +# MG5_aMC available, but available in Pythia8. * +# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * +#*********************************************************************** + $ickkw = ickkw +#*********************************************************************** +# +#*********************************************************************** +# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are * +# written in the LHE event file * +#*********************************************************************** + 15.0 = bwcutoff +#*********************************************************************** +# Cuts on the jets. Jet clustering is performed by FastJet. * +# - When matching to a parton shower, these generation cuts should be * +# considerably softer than the analysis cuts. * +# - More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) + 1.0 = jetradius ! The radius parameter for the jet algorithm + 0.0 = ptj ! Min jet transverse momentum + -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) +#*********************************************************************** +# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptl ! Min lepton transverse momentum + -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) + 0.0 = drll ! Min distance between opposite sign lepton pairs + 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs + 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs + 0.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, * +# all the other parameters are ignored. * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + -1.0 = etagamma ! Max photon abs(pseudo-rap) + 0.15 = R0gamma ! Radius of isolation code + 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 + 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 + True = isoEM ! isolate photons from EM energy (photons and leptons) +#*********************************************************************** +# Cuts associated to MASSIVE particles identified by their PDG codes. * +# All cuts are applied to both particles and anti-particles, so use * +# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or * +# {6:100, 25:200} for multiple particles * +#*********************************************************************** + {} = pt_min_pdg ! Min pT for a massive particle + {} = pt_max_pdg ! Max pT for a massive particle + {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles +#*********************************************************************** +# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* +#*********************************************************************** + 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) +#*********************************************************************** +## Maximal PDG code for quark to be considered a jet when applying cuts.* +## At least all massless quarks of the model should be included here. * +##*********************************************************************** + $maxjetflavor = maxjetflavor +# #*********************************************************************** diff --git a/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat new file mode 100644 index 00000000..8188c726 --- /dev/null +++ b/Campaigns/HINPbPbSpring23/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat @@ -0,0 +1,188 @@ +#*********************************************************************** +# MadGraph5_aMC@NLO * +# * +# run_card.dat aMC@NLO * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a hash (#) are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# Some of the values of variables can be list. These can either be * +# comma or space separated. * +# * +# To display additional parameter, you can use the command: * +# update to_full * +#*********************************************************************** +# +#******************* +# Running parameters +#******************* +# +#*********************************************************************** +# Tag name for the run (one word) * +#*********************************************************************** + tag_1 = run_tag ! name of the run +#*********************************************************************** +# Number of LHE events (and their normalization) and the required * +# (relative) accuracy on the Xsec. * +# These values are ignored for fixed order runs * +#*********************************************************************** + 5000 = nevents ! Number of unweighted events requested + 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents) + -1 = nevt_job! Max number of events per job in event generation. + ! (-1= no split). +#*********************************************************************** +# Normalize the weights of LHE events such that they sum or average to * +# the total cross section * +#*********************************************************************** + average = event_norm ! valid settings: average, sum, bias +#*********************************************************************** +# Number of points per itegration channel (ignored for aMC@NLO runs) * +#*********************************************************************** + 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the + ! number of points and iter. below) +# These numbers are ignored except if req_acc_FO is equal to -1 + 5000 = npoints_FO_grid ! number of points to setup grids + 4 = niters_FO_grid ! number of iter. to setup grids + 10000 = npoints_FO ! number of points to compute Xsec + 6 = niters_FO ! number of iter. to compute Xsec +#*********************************************************************** +# Random number seed * +#*********************************************************************** + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#*********************************************************************** +# Collider type and energy * +#*********************************************************************** + 1 = lpp1 ! beam 1 type (0 = no PDF) + 1 = lpp2 ! beam 2 type (0 = no PDF) + $ebeam1 = ebeam1 ! beam 1 energy in GeV + $ebeam2 = ebeam2 ! beam 2 energy in GeV +#*********************************************************************** +# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * +#*********************************************************************** + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_nPDF_SETS = lhaid +#*********************************************************************** +# Include the NLO Monte Carlo subtr. terms for the following parton * +# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * +# WARNING: PYTHIA6PT works only for processes without FSR!!!! * +#*********************************************************************** + $parton_shower = parton_shower + 1.0 = shower_scale_factor ! multiply default shower starting + ! scale by this factor +#*********************************************************************** +# Renormalization and factorization scales * +# (Default functional form for the non-fixed scales is the sum of * +# the transverse masses divided by two of all final state particles * +# and partons. This can be changed in SubProcesses/set_scales.f or via * +# dynamical_scale_choice option) * +#*********************************************************************** + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.118 = muR_ref_fixed ! fixed ren reference scale + 91.118 = muF_ref_fixed ! fixed fact reference scale + -1 = dynamical_scale_choice ! Choose one (or more) of the predefined + ! dynamical choices. Can be a list; scale choices beyond the + ! first are included via reweighting + 1.0 = muR_over_ref ! ratio of current muR over reference muR + 1.0 = muF_over_ref ! ratio of current muF over reference muF +#*********************************************************************** +# Reweight variables for scale dependence and PDF uncertainty * +#*********************************************************************** + 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting + 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting + True = reweight_scale ! Reweight to get scale variation using the + ! rw_rscale and rw_fscale factors. Should be a list of + ! booleans of equal length to dynamical_scale_choice to + ! specify for which choice to include scale dependence. +#*********************************************************************** +# Store reweight information in the LHE file for off-line model- * +# parameter reweighting at NLO+PS accuracy * +#*********************************************************************** + True = store_rwgt_info ! Store info for reweighting in LHE file +#*********************************************************************** +# Parameters relevant for the MasSTR plugin: * +# iSTR controls the strategy for the resonance treatment * +# istr = 1 -> DR without interferece * +# istr = 2 -> DR with interferece * +# istr = 3 -> DS with reshuffling on initial state, standard BW * +# istr = 4 -> DS with reshuffling on initial state, running BW * +# istr = 5 -> DS with reshuffling on all FS particles, standard BW * +# istr = 6 -> DS with reshuffling on all FS particles, running BW * +#*********************************************************************** + $istr = istr ! strategy to be used to remove resonances + ! appearing in real emissions + True = str_include_pdf ! compensate for PDFs when doing reshuffling + True = str_include_flux ! compensate for flux when doing reshuffling +#*********************************************************************** +# ickkw parameter: * +# 0: No merging * +# 3: FxFx Merging - WARNING! Applies merging only at the hard-event * +# level. After showering an MLM-type merging should be applied as * +# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * +# 4: UNLOPS merging (with pythia8 only). No interface from within * +# MG5_aMC available, but available in Pythia8. * +# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * +#*********************************************************************** + $ickkw = ickkw +#*********************************************************************** +# +#*********************************************************************** +# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are * +# written in the LHE event file * +#*********************************************************************** + 15.0 = bwcutoff +#*********************************************************************** +# Cuts on the jets. Jet clustering is performed by FastJet. * +# - When matching to a parton shower, these generation cuts should be * +# considerably softer than the analysis cuts. * +# - More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) + 1.0 = jetradius ! The radius parameter for the jet algorithm + 0.0 = ptj ! Min jet transverse momentum + -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) +#*********************************************************************** +# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptl ! Min lepton transverse momentum + -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) + 0.0 = drll ! Min distance between opposite sign lepton pairs + 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs + 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs + 0.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, * +# all the other parameters are ignored. * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + -1.0 = etagamma ! Max photon abs(pseudo-rap) + 0.15 = R0gamma ! Radius of isolation code + 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 + 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 + True = isoEM ! isolate photons from EM energy (photons and leptons) +#*********************************************************************** +# Cuts associated to MASSIVE particles identified by their PDG codes. * +# All cuts are applied to both particles and anti-particles, so use * +# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or * +# {6:100, 25:200} for multiple particles * +#*********************************************************************** + {} = pt_min_pdg ! Min pT for a massive particle + {} = pt_max_pdg ! Max pT for a massive particle + {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles +#*********************************************************************** +# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* +#*********************************************************************** + 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) +#*********************************************************************** +## Maximal PDG code for quark to be considered a jet when applying cuts.* +## At least all massless quarks of the model should be included here. * +##*********************************************************************** + $maxjetflavor = maxjetflavor +# #*********************************************************************** diff --git a/Campaigns/HINppWinter24/HINppWinter24.json b/Campaigns/HINppWinter24/HINppWinter24.json new file mode 100644 index 00000000..10f99ba5 --- /dev/null +++ b/Campaigns/HINppWinter24/HINppWinter24.json @@ -0,0 +1,16 @@ +{ + "chain": "chain_HINppWinter24wmLHEGS_flowHINppWinter24DRPremix_flowHINppWinter24MiniAOD_flowHINppWinter24NanoAOD", + "model_params_vars": {}, + "template_vars": { + + }, + "fragment_vars": { + "concurrentLHEforMadGraph5_aMCatNLO": "False", + "concurrentLHEforPowheg": "True", + "concurrentPS": "Concurrent" + }, + "tune": "CP5", + "beam": 2680, + "gridpack_directory": "/eos/cms/store/group/phys_generator/cvmfs/gridpacks/PdmV/HINppWinter24" +} + diff --git a/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massW.dat b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massW.dat new file mode 100644 index 00000000..41a8eb41 --- /dev/null +++ b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massW.dat @@ -0,0 +1,22 @@ +# These inputs are based on fixing the PDG vlaues of mZ, mW, and G_fermi +# The MC programs use a constant width scheme, whereas PDG values are sqrt(s) dependent. +# We correct for this following Sect. 3.1 in https://arxiv.org/pdf/1606.02330.pdf + +set param_card SMINPUTS 1 132.30648214177197 +set param_card SMINPUTS 2 1.1663787e-05 + +# Z boson +set param_card mass 23 91.153509740726733 +set param_card DECAY 23 2.4932018986110700 + +# W boson +set param_card mass 24 80.351812293789408 +set param_card DECAY 24 2.0904310808144846 + +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCTopWGSummaryPlots#Top_Quark_Mass +# CMS combined result (Sep 2015) + +# top quark +set param_card mass 6 172.5 +set param_card DECAY 6 1.311 +set param_card yukawa 6 172.5 diff --git a/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massZ.dat b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massZ.dat new file mode 100644 index 00000000..3dfd7695 --- /dev/null +++ b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/ModelParams/massZ.dat @@ -0,0 +1,22 @@ +# These inputs are based on fixing the PDG vlaues of mZ, mW, and G_fermi +# The MC programs use a constant width scheme, whereas PDG values are sqrt(s) dependent. +# We correct for this following Sect. 3.1 in https://arxiv.org/pdf/1606.02330.pdf + +set param_card SMINPUTS 1 128.82531590804655 +set param_card SMINPUTS 2 1.1663787e-05 + +# Z boson +set param_card mass 23 91.153509740726733 +set param_card DECAY 23 2.4932018986110700 + +# W boson +set param_card mass 24 79.906853549493746 +set param_card DECAY 24 2.085 + +# https://twiki.cern.ch/twiki/bin/view/LHCPhysics/LHCTopWGSummaryPlots#Top_Quark_Mass +# CMS combined result (Sep 2015) + +# top quark +set param_card mass 6 172.5 +set param_card DECAY 6 1.311 +set param_card yukawa 6 172.5 diff --git a/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/LO_run_card.dat b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/LO_run_card.dat new file mode 100644 index 00000000..4814e9e1 --- /dev/null +++ b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/LO_run_card.dat @@ -0,0 +1,227 @@ +#********************************************************************* +# MadGraph5_aMC@NLO * +# * +# run_card.dat MadEvent * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a '# ' are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# To display more options, you can type the command: * +# update full_run_card * +#********************************************************************* +# +#********************************************************************* +# Tag name for the run (one word) * +#********************************************************************* + tag_1 = run_tag ! name of the run +#********************************************************************* +# Number of events and rnd seed * +# Warning: Do not generate more than 1M events in a single run * +#********************************************************************* + 3000 = nevents ! Number of unweighted events requested + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron, 4=photon from muon * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + $ebeam1 = ebeam1 ! beam 1 total energy in GeV + $ebeam2 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid ! if pdlabel=lhapdf, this is the lhapdf number + $DEFAULT_PDF_MEMBERS = reweight_PDF +# To see heavy ion options: type "update ion_pdf" +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.188 = scale ! fixed ren scale + 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 + 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 + -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices + 1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Type and output format +#********************************************************************* + True = gridpack !True = setting up the grid pack + -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written) + average = event_norm ! average/sum. Normalization of the weight in the LHEF +#********************************************************************* +# Matching parameter (MLM only) +#********************************************************************* + $ickkw = ickkw ! 0 no matching, 1 MLM + 1.0 = alpsfact ! scale factor for QCD emission vx + False = chcluster ! cluster only according to channel diag + $maxjetflavor = asrwgtflavor ! highest quark flavor for a_s reweight + True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 + ! (turn off for VBF and single top processes) + $xqcut = xqcut ! minimum kt jet measure between partons + +#*********************************************************************** +# Turn on either the ktdurham or ptlund cut to activate * +# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] * +#*********************************************************************** + -1.0 = ktdurham + 0.4 = dparameter + -1.0 = ptlund + 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above + +#********************************************************************* +# +#********************************************************************* +# Phase-Space Optimization strategy (basic options) +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. + ! 0: sum over helicity, 1: importance sampling + 1 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx) + ! 1 is old strategy (using amp square) + ! 2 is new strategy (using only the denominator) +# To see advanced option for Phase-Space optimization: type "update psoptim" +#********************************************************************* +# Generation bias, check the wiki page below for more information: * +# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' * +#********************************************************************* + None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-] + {} = bias_parameters ! Specifies the parameters of the module. +# +#******************************* +# Parton level cuts definition * +#******************************* +# +# +#********************************************************************* +# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay +#********************************************************************* + 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 0.0 = pta ! minimum pt for the photons + 0.0 = ptj ! minimum pt for the jets + 0.0 = ptl ! minimum pt for the charged leptons + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptlmax ! maximum pt for the charged leptons + {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle + {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) +# +# For display option for energy cut in the partonic center of mass frame type 'update ecut' +# +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + -1.0 = etaa ! max rap for the photons + -1.0 = etaj ! max rap for the jets + -1.0 = etal ! max rap for the charged leptons + 0.0 = etalmin ! main rap for the charged leptons + {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle + {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5}) +#********************************************************************* +# Minimum and maximum DeltaR distance * +#********************************************************************* + 0.0 = draj ! min distance between gamma and jet + 0.0 = dral ! min distance between gamma and lepton + 0.0 = drjj ! min distance between jets + 0.0 = drll ! min distance between leptons + 0.0 = drjl ! min distance between jet and lepton + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = dralmax ! maxdistance between gamma and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drllmax ! max distance between leptons + -1.0 = drjlmax ! max distance between jet and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 0.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair + -1.0 = mmjjmax ! max invariant mass of a jet pair + -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair + {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250}) + {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only + ! to pairs of particle/antiparticle and not to pairs of the same pdg codes. + #********************************************************************* + # Minimum and maximum invariant mass for all letpons * + #********************************************************************* + 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) + -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) + #********************************************************************* + # Minimum and maximum pt for 4-momenta sum of leptons / neutrino * + # for pair of lepton includes only same flavor, opposite charge + #********************************************************************* + 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) + -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) +#********************************************************************* +# Inclusive cuts * +#********************************************************************* + 0.0 = xptj ! minimum pt for at least one jet + 0.0 = xptl ! minimum pt for at least one charged lepton + #********************************************************************* + # Control the pt's of the jets sorted by pt * + #********************************************************************* + 0.0 = ptj1min ! minimum pt for the leading jet in pt + 0.0 = ptj2min ! minimum pt for the second jet in pt + 0.0 = ptj3min ! minimum pt for the third jet in pt + 0.0 = ptj4min ! minimum pt for the fourth jet in pt + -1.0 = ptj1max ! maximum pt for the leading jet in pt + -1.0 = ptj2max ! maximum pt for the second jet in pt + -1.0 = ptj3max ! maximum pt for the third jet in pt + -1.0 = ptj4max ! maximum pt for the fourth jet in pt + 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts + #********************************************************************* + # Control the pt's of leptons sorted by pt * + #********************************************************************* + 0.0 = ptl1min ! minimum pt for the leading lepton in pt + 0.0 = ptl2min ! minimum pt for the second lepton in pt + -1.0 = ptl1max ! maximum pt for the leading lepton in pt + -1.0 = ptl2max ! maximum pt for the second lepton in pt + #********************************************************************* + # Control the Ht(k)=Sum of k leading jets * + #********************************************************************* + 0.0 = htjmin ! minimum jet HT=Sum(jet pt) + -1.0 = htjmax ! maximum jet HT=Sum(jet pt) + 0.0 = ihtmin !inclusive Ht for all partons (including b) + -1.0 = ihtmax !inclusive Ht for all partons (including b) + 0.0 = ht2min ! minimum Ht for the two leading jets + 0.0 = ht3min ! minimum Ht for the three leading jets + 0.0 = ht4min ! minimum Ht for the four leading jets + -1.0 = ht2max ! maximum Ht for the two leading jets + -1.0 = ht3max ! maximum Ht for the three leading jets + -1.0 = ht4max ! maximum Ht for the four leading jets + #********************************************************************* + # WBF cuts * + #********************************************************************* + 0.0 = xetamin ! minimum rapidity for two jets in the WBF case + 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* + $maxjetflavor = maxjetflavor ! Maximum jet pdg code +#********************************************************************* +# +#********************************************************************* +# Store info for systematics studies * +# WARNING: Do not use for interference type of computation * +#********************************************************************* + True = use_syst ! Enable systematics studies + +#************************************** +# Manually added by Sihyun Jeon +#************************************** + False = cut_decays ! Cut decay products + False = pdfwgt ! for ickkw=1, perform pdf reweighting diff --git a/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat new file mode 100644 index 00000000..a0336423 --- /dev/null +++ b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card.dat @@ -0,0 +1,174 @@ +#*********************************************************************** +# MadGraph5_aMC@NLO * +# * +# run_card.dat aMC@NLO * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a hash (#) are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# Some of the values of variables can be list. These can either be * +# comma or space separated. * +# * +# To display additional parameter, you can use the command: * +# update to_full * +#*********************************************************************** +# +#******************* +# Running parameters +#******************* +# +#*********************************************************************** +# Tag name for the run (one word) * +#*********************************************************************** + tag_1 = run_tag ! name of the run +#*********************************************************************** +# Number of LHE events (and their normalization) and the required * +# (relative) accuracy on the Xsec. * +# These values are ignored for fixed order runs * +#*********************************************************************** + 5000 = nevents ! Number of unweighted events requested + 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents) + -1 = nevt_job! Max number of events per job in event generation. + ! (-1= no split). +#*********************************************************************** +# Normalize the weights of LHE events such that they sum or average to * +# the total cross section * +#*********************************************************************** + average = event_norm ! valid settings: average, sum, bias +#*********************************************************************** +# Number of points per itegration channel (ignored for aMC@NLO runs) * +#*********************************************************************** + 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the + ! number of points and iter. below) +# These numbers are ignored except if req_acc_FO is equal to -1 + 5000 = npoints_FO_grid ! number of points to setup grids + 4 = niters_FO_grid ! number of iter. to setup grids + 10000 = npoints_FO ! number of points to compute Xsec + 6 = niters_FO ! number of iter. to compute Xsec +#*********************************************************************** +# Random number seed * +#*********************************************************************** + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#*********************************************************************** +# Collider type and energy * +#*********************************************************************** + 1 = lpp1 ! beam 1 type (0 = no PDF) + 1 = lpp2 ! beam 2 type (0 = no PDF) + $ebeam1 = ebeam1 ! beam 1 energy in GeV + $ebeam2 = ebeam2 ! beam 2 energy in GeV +#*********************************************************************** +# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * +#*********************************************************************** + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid +#*********************************************************************** +# Include the NLO Monte Carlo subtr. terms for the following parton * +# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * +# WARNING: PYTHIA6PT works only for processes without FSR!!!! * +#*********************************************************************** + $parton_shower = parton_shower + 1.0 = shower_scale_factor ! multiply default shower starting + ! scale by this factor +#*********************************************************************** +# Renormalization and factorization scales * +# (Default functional form for the non-fixed scales is the sum of * +# the transverse masses divided by two of all final state particles * +# and partons. This can be changed in SubProcesses/set_scales.f or via * +# dynamical_scale_choice option) * +#*********************************************************************** + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.118 = muR_ref_fixed ! fixed ren reference scale + 91.118 = muF_ref_fixed ! fixed fact reference scale + -1 = dynamical_scale_choice ! Choose one (or more) of the predefined + ! dynamical choices. Can be a list; scale choices beyond the + ! first are included via reweighting + 1.0 = muR_over_ref ! ratio of current muR over reference muR + 1.0 = muF_over_ref ! ratio of current muF over reference muF +#*********************************************************************** +# Reweight variables for scale dependence and PDF uncertainty * +#*********************************************************************** + 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting + 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting + True = reweight_scale ! Reweight to get scale variation using the + ! rw_rscale and rw_fscale factors. Should be a list of + ! booleans of equal length to dynamical_scale_choice to + ! specify for which choice to include scale dependence. +#*********************************************************************** +# Store reweight information in the LHE file for off-line model- * +# parameter reweighting at NLO+PS accuracy * +#*********************************************************************** + True = store_rwgt_info ! Store info for reweighting in LHE file +#*********************************************************************** +# ickkw parameter: * +# 0: No merging * +# 3: FxFx Merging - WARNING! Applies merging only at the hard-event * +# level. After showering an MLM-type merging should be applied as * +# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * +# 4: UNLOPS merging (with pythia8 only). No interface from within * +# MG5_aMC available, but available in Pythia8. * +# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * +#*********************************************************************** + $ickkw = ickkw +#*********************************************************************** +# +#*********************************************************************** +# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are * +# written in the LHE event file * +#*********************************************************************** + 15.0 = bwcutoff +#*********************************************************************** +# Cuts on the jets. Jet clustering is performed by FastJet. * +# - When matching to a parton shower, these generation cuts should be * +# considerably softer than the analysis cuts. * +# - More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) + 1.0 = jetradius ! The radius parameter for the jet algorithm + 0.0 = ptj ! Min jet transverse momentum + -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) +#*********************************************************************** +# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptl ! Min lepton transverse momentum + -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) + 0.0 = drll ! Min distance between opposite sign lepton pairs + 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs + 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs + 0.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, * +# all the other parameters are ignored. * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + -1.0 = etagamma ! Max photon abs(pseudo-rap) + 0.15 = R0gamma ! Radius of isolation code + 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 + 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 + True = isoEM ! isolate photons from EM energy (photons and leptons) +#*********************************************************************** +# Cuts associated to MASSIVE particles identified by their PDG codes. * +# All cuts are applied to both particles and anti-particles, so use * +# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or * +# {6:100, 25:200} for multiple particles * +#*********************************************************************** + {} = pt_min_pdg ! Min pT for a massive particle + {} = pt_max_pdg ! Max pT for a massive particle + {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles +#*********************************************************************** +# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* +#*********************************************************************** + 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) +#*********************************************************************** +## Maximal PDG code for quark to be considered a jet when applying cuts.* +## At least all massless quarks of the model should be included here. * +##*********************************************************************** + $maxjetflavor = maxjetflavor +# #*********************************************************************** diff --git a/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat new file mode 100644 index 00000000..a4693780 --- /dev/null +++ b/Campaigns/HINppWinter24/MadGraph5_aMCatNLO/Templates/NLO_run_card_madstr.dat @@ -0,0 +1,188 @@ +#*********************************************************************** +# MadGraph5_aMC@NLO * +# * +# run_card.dat aMC@NLO * +# * +# This file is used to set the parameters of the run. * +# * +# Some notation/conventions: * +# * +# Lines starting with a hash (#) are info or comments * +# * +# mind the format: value = variable ! comment * +# * +# Some of the values of variables can be list. These can either be * +# comma or space separated. * +# * +# To display additional parameter, you can use the command: * +# update to_full * +#*********************************************************************** +# +#******************* +# Running parameters +#******************* +# +#*********************************************************************** +# Tag name for the run (one word) * +#*********************************************************************** + tag_1 = run_tag ! name of the run +#*********************************************************************** +# Number of LHE events (and their normalization) and the required * +# (relative) accuracy on the Xsec. * +# These values are ignored for fixed order runs * +#*********************************************************************** + 5000 = nevents ! Number of unweighted events requested + 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents) + -1 = nevt_job! Max number of events per job in event generation. + ! (-1= no split). +#*********************************************************************** +# Normalize the weights of LHE events such that they sum or average to * +# the total cross section * +#*********************************************************************** + average = event_norm ! valid settings: average, sum, bias +#*********************************************************************** +# Number of points per itegration channel (ignored for aMC@NLO runs) * +#*********************************************************************** + 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the + ! number of points and iter. below) +# These numbers are ignored except if req_acc_FO is equal to -1 + 5000 = npoints_FO_grid ! number of points to setup grids + 4 = niters_FO_grid ! number of iter. to setup grids + 10000 = npoints_FO ! number of points to compute Xsec + 6 = niters_FO ! number of iter. to compute Xsec +#*********************************************************************** +# Random number seed * +#*********************************************************************** + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#*********************************************************************** +# Collider type and energy * +#*********************************************************************** + 1 = lpp1 ! beam 1 type (0 = no PDF) + 1 = lpp2 ! beam 2 type (0 = no PDF) + $ebeam1 = ebeam1 ! beam 1 energy in GeV + $ebeam2 = ebeam2 ! beam 2 energy in GeV +#*********************************************************************** +# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. * +#*********************************************************************** + 'lhapdf' = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid +#*********************************************************************** +# Include the NLO Monte Carlo subtr. terms for the following parton * +# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) * +# WARNING: PYTHIA6PT works only for processes without FSR!!!! * +#*********************************************************************** + $parton_shower = parton_shower + 1.0 = shower_scale_factor ! multiply default shower starting + ! scale by this factor +#*********************************************************************** +# Renormalization and factorization scales * +# (Default functional form for the non-fixed scales is the sum of * +# the transverse masses divided by two of all final state particles * +# and partons. This can be changed in SubProcesses/set_scales.f or via * +# dynamical_scale_choice option) * +#*********************************************************************** + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.118 = muR_ref_fixed ! fixed ren reference scale + 91.118 = muF_ref_fixed ! fixed fact reference scale + -1 = dynamical_scale_choice ! Choose one (or more) of the predefined + ! dynamical choices. Can be a list; scale choices beyond the + ! first are included via reweighting + 1.0 = muR_over_ref ! ratio of current muR over reference muR + 1.0 = muF_over_ref ! ratio of current muF over reference muF +#*********************************************************************** +# Reweight variables for scale dependence and PDF uncertainty * +#*********************************************************************** + 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting + 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting + True = reweight_scale ! Reweight to get scale variation using the + ! rw_rscale and rw_fscale factors. Should be a list of + ! booleans of equal length to dynamical_scale_choice to + ! specify for which choice to include scale dependence. +#*********************************************************************** +# Store reweight information in the LHE file for off-line model- * +# parameter reweighting at NLO+PS accuracy * +#*********************************************************************** + True = store_rwgt_info ! Store info for reweighting in LHE file +#*********************************************************************** +# Parameters relevant for the MasSTR plugin: * +# iSTR controls the strategy for the resonance treatment * +# istr = 1 -> DR without interferece * +# istr = 2 -> DR with interferece * +# istr = 3 -> DS with reshuffling on initial state, standard BW * +# istr = 4 -> DS with reshuffling on initial state, running BW * +# istr = 5 -> DS with reshuffling on all FS particles, standard BW * +# istr = 6 -> DS with reshuffling on all FS particles, running BW * +#*********************************************************************** + $istr = istr ! strategy to be used to remove resonances + ! appearing in real emissions + True = str_include_pdf ! compensate for PDFs when doing reshuffling + True = str_include_flux ! compensate for flux when doing reshuffling +#*********************************************************************** +# ickkw parameter: * +# 0: No merging * +# 3: FxFx Merging - WARNING! Applies merging only at the hard-event * +# level. After showering an MLM-type merging should be applied as * +# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. * +# 4: UNLOPS merging (with pythia8 only). No interface from within * +# MG5_aMC available, but available in Pythia8. * +# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. * +#*********************************************************************** + $ickkw = ickkw +#*********************************************************************** +# +#*********************************************************************** +# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are * +# written in the LHE event file * +#*********************************************************************** + 15.0 = bwcutoff +#*********************************************************************** +# Cuts on the jets. Jet clustering is performed by FastJet. * +# - When matching to a parton shower, these generation cuts should be * +# considerably softer than the analysis cuts. * +# - More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT) + 1.0 = jetradius ! The radius parameter for the jet algorithm + 0.0 = ptj ! Min jet transverse momentum + -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut) +#*********************************************************************** +# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptl ! Min lepton transverse momentum + -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut) + 0.0 = drll ! Min distance between opposite sign lepton pairs + 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs + 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs + 0.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs +#*********************************************************************** +# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, * +# all the other parameters are ignored. * +# More specific cuts can be specified in SubProcesses/cuts.f * +#*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + -1.0 = etagamma ! Max photon abs(pseudo-rap) + 0.15 = R0gamma ! Radius of isolation code + 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442 + 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442 + True = isoEM ! isolate photons from EM energy (photons and leptons) +#*********************************************************************** +# Cuts associated to MASSIVE particles identified by their PDG codes. * +# All cuts are applied to both particles and anti-particles, so use * +# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or * +# {6:100, 25:200} for multiple particles * +#*********************************************************************** + {} = pt_min_pdg ! Min pT for a massive particle + {} = pt_max_pdg ! Max pT for a massive particle + {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles +#*********************************************************************** +# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* +#*********************************************************************** + 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) +#*********************************************************************** +## Maximal PDG code for quark to be considered a jet when applying cuts.* +## At least all massless quarks of the model should be included here. * +##*********************************************************************** + $maxjetflavor = maxjetflavor +# #*********************************************************************** diff --git a/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8.json b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8.json new file mode 100644 index 00000000..b6463386 --- /dev/null +++ b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8.json @@ -0,0 +1,26 @@ +{ + "template": "NLO_run_card.dat", + "template_user": [], + "template_vars":{ + "parton_shower": "pythia8", + "ickkw": "3", + "maxjetflavor": "5" + }, + "model_params": "massW.dat", + "model_params_user": [ + "set run_card ptj 20" + ], + "model_params_vars": {}, + "fragment": ["Generator/ExternalLHEProducer_MadGraph5_aMCatNLO.dat", "PartonShower/amcatnloFXFX-pythia8.dat"], + "fragment_user": [], + "fragment_vars": { + "processParameters": [ + "'JetMatching:qCut = 40.'", + "'JetMatching:qCutME = 20.'", + "'JetMatching:nQmatch = 5'", + "'JetMatching:nJetMax = 2'", + "'TimeShower:mMaxGamma = 4.0'" + ] + } +} + diff --git a/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_madspin_card.dat b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_madspin_card.dat new file mode 100644 index 00000000..a0642dc2 --- /dev/null +++ b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_madspin_card.dat @@ -0,0 +1,13 @@ +#directory for gridpack mode +set ms_dir ./madspingrid + +#initialization parameters +set Nevents_for_max_weight 250 # number of events for the estimate of the max. weight +set max_weight_ps_point 400 # number of PS to estimate the maximum for each event + +set max_running_process 1 + +decay t > w+ b, w+ > all all +decay t~ > w- b~, w- > all all + +launch diff --git a/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_proc_card.dat b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_proc_card.dat new file mode 100644 index 00000000..0a5945eb --- /dev/null +++ b/Cards/MadGraph5_aMCatNLO/TT/TT-2Jets_amcatnloFXFX-pythia8/TT-2Jets_amcatnloFXFX-pythia8_proc_card.dat @@ -0,0 +1,12 @@ +import model loop_sm-ckm_no_b_mass + +define ell+ = e+ mu+ ta+ +define ell- = e- mu- ta- +define vl = ve vm vt +define vl~ = ve~ vm~ vt~ + +generate p p > t t~ [QCD] @0 +add process p p > t t~ j [QCD] @1 +add process p p > t t~ j j [QCD] @2 + +output TT-2Jets_amcatnloFXFX-pythia8 -nojpeg