diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_customizecards.dat new file mode 100644 index 00000000000..d43837c52b1 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_customizecards.dat @@ -0,0 +1,8 @@ +set param_card mass 9900012 @@n1mass@@ +set param_card mass 9900014 @@n2mass@@ +set param_card mass 9900016 999999 +set param_card mass 34 @@wrmass@@ +set decay 9900012 @@n1width@@ +set decay 9900014 @@n2width@@ +set decay 9900016 10 +set decay 34 AUTO diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_extramodels.dat new file mode 100644 index 00000000000..9ec50e402fd --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_extramodels.dat @@ -0,0 +1 @@ +EffLRSM_NLO.tgz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_proc_card.dat new file mode 100644 index 00000000000..82d5abf1398 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_proc_card.dat @@ -0,0 +1,8 @@ +import model EffLRSM_NLO +define p = u d c s b u~ d~ c~ s~ b~ g +define j = p +@@definell@@ +@@definenn@@ +generate p p > nn ll, (nn > ll t b~, (t > w+ b, w+ > j j)) +add process p p > nn ll, (nn > ll t~ b, (t~ > w- b~, w- > j j)) +output @@output@@ diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_run_card.dat new file mode 100644 index 00000000000..2966ad13e0f --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/base/base_run_card.dat @@ -0,0 +1,224 @@ +#********************************************************************* +# 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 + @@beame@@ = ebeam1 ! beam 1 total energy in GeV + @@beame@@ = 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) +#********************************************************************* + 0 = ickkw ! 0 no matching, 1 MLM + 1.0 = alpsfact ! scale factor for QCD emission vx + False = chcluster ! cluster only according to channel diag + 5 = 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) + 0 = 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) * +#********************************************************************* + 5 = 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 + + False = cut_decays ! Cut decay products + False = pdfwgt ! for ickkw=1, perform pdf reweighting diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_customizecards.dat new file mode 100644 index 00000000000..cf934cf9887 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_customizecards.dat @@ -0,0 +1,8 @@ +set param_card mass 9900012 999999 +set param_card mass 9900014 700 +set param_card mass 9900016 999999 +set param_card mass 34 1000 +set decay 9900012 10 +set decay 9900014 AUTO +set decay 9900016 10 +set decay 34 AUTO diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_extramodels.dat new file mode 100644 index 00000000000..9ec50e402fd --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_extramodels.dat @@ -0,0 +1 @@ +EffLRSM_NLO.tgz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_proc_card.dat new file mode 100644 index 00000000000..ddbcd2d330c --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_proc_card.dat @@ -0,0 +1,8 @@ +import model EffLRSM_NLO +define p = u d c s b u~ d~ c~ s~ b~ g +define j = p +define ll = mu+ mu- +define nn = n2 +generate p p > nn ll, (nn > ll t b~, (t > w+ b, w+ > j j)) +add process p p > nn ll, (nn > ll t~ b, (t~ > w- b~, w- > j j)) +output WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_run_card.dat new file mode 100644 index 00000000000..2c3b31b2c73 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/cards/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV/WRtoNMutoMuMuTB-HadTop_MWR-1000_MN-700_13TeV_run_card.dat @@ -0,0 +1,224 @@ +#********************************************************************* +# 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 + 6500 = ebeam1 ! beam 1 total energy in GeV + 6500 = 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) +#********************************************************************* + 0 = ickkw ! 0 no matching, 1 MLM + 1.0 = alpsfact ! scale factor for QCD emission vx + False = chcluster ! cluster only according to channel diag + 5 = 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) + 0 = 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) * +#********************************************************************* + 5 = 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 + + False = cut_decays ! Cut decay products + False = pdfwgt ! for ickkw=1, perform pdf reweighting diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/write.py b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/write.py new file mode 100644 index 00000000000..2df6e27d823 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WRtoNLtoLLTB-HadTop/write.py @@ -0,0 +1,61 @@ +import os + +flavors = { + "Mu" : ["define ll = mu+ mu-", "define nn = n2"], + "El" : ["define ll = e+ e-", "define nn = n1"] +} +run = { + "13TeV" : "6500", + "13p6TeV" : "6800" +} +masses = {} +count = 0 +for i_wr in range(2, 3):# 13): + m_wr = i_wr*500 + m_n = 300 + masses[str(m_wr)] = [] + while (m_n < m_wr): + if m_n >= m_wr*0.6: + masses[str(m_wr)].append(str(m_n)) + count = count + 1 + m_n = m_n + 200 + if m_n-10 > m_wr: + masses[str(m_wr)].append(str(m_wr-10)) + else: + masses[str(m_wr)].append(str(m_n-10)) + +def CopyAtoB(dataset_name, card_name): + os.system(f"cp base/base_{card_name} cards/{dataset_name}/{dataset_name}_{card_name}") + +def SedAtoB(dataset_name, card_name, A, B): + os.system(f"sed -i 's|@@{A}@@|{B}|g' cards/{dataset_name}/{dataset_name}_{card_name}") + +for flavor, defines in flavors.items(): + for com, ebeam in run.items(): + for m_wr, m_ns in masses.items(): + for m_n in m_ns: + dataset_name = f"WRtoN{flavor}to{flavor}{flavor}TB-HadTop_MWR-{m_wr}_MN-{m_n}_{com}" + os.system(f"mkdir -p cards/{dataset_name}") + CopyAtoB(dataset_name, "proc_card.dat") + CopyAtoB(dataset_name, "run_card.dat") + CopyAtoB(dataset_name, "customizecards.dat") + CopyAtoB(dataset_name, "extramodels.dat") + SedAtoB(dataset_name, "proc_card.dat", "output", dataset_name) + SedAtoB(dataset_name, "proc_card.dat", "definell", defines[0]) + SedAtoB(dataset_name, "proc_card.dat", "definenn", defines[1]) + SedAtoB(dataset_name, "run_card.dat", "beame", ebeam) + dummymass = 999999 + dummywidth = 10 + auto = "AUTO" + if flavor == "El": + SedAtoB(dataset_name, "customizecards.dat", "n1mass", m_n) + SedAtoB(dataset_name, "customizecards.dat", "n1width", auto) + SedAtoB(dataset_name, "customizecards.dat", "n2mass", dummymass) + SedAtoB(dataset_name, "customizecards.dat", "n2width", dummywidth) + else: + SedAtoB(dataset_name, "customizecards.dat", "n2mass", m_n) + SedAtoB(dataset_name, "customizecards.dat", "n2width", auto) + SedAtoB(dataset_name, "customizecards.dat", "n1mass", dummymass) + SedAtoB(dataset_name, "customizecards.dat", "n1width", dummywidth) + SedAtoB(dataset_name, "customizecards.dat", "wrmass", m_wr) + print (f"./gridpack_generation.sh {dataset_name} cards/{dataset_name}")