diff --git a/bin/MadGraph5_aMCatNLO/Utilities/source_condor.sh b/bin/MadGraph5_aMCatNLO/Utilities/source_condor.sh index ae18a8105eb1..ecdbe7801db8 100644 --- a/bin/MadGraph5_aMCatNLO/Utilities/source_condor.sh +++ b/bin/MadGraph5_aMCatNLO/Utilities/source_condor.sh @@ -1,7 +1,13 @@ #!/bin/bash # HTCondor python bindings are lost after cmsenv/scram -PYTHON_BINDINGS="$(python3 -c 'import htcondor; import os; print(os.path.dirname(htcondor.__path__[0]))' 2>/dev/null)" -if [ -z "$PYTHON_BINDINGS" ]; then - pip3 install htcondor --user #FIXME need better way to interface HTCondor Python API for Python3.9 +#PYTHON_BINDINGS="$(python3 -c 'import htcondor; import os; print(os.path.dirname(htcondor.__path__[0]))' 2>/dev/null)" +#if [ -z "$PYTHON_BINDINGS" ]; then +if [ -d $CMSSW_BASE/venv ]; +then + echo "venv exists" +else + scram-venv + cmsenv + pip3 install htcondor --upgrade #FIXME need better way to interface HTCondor Python API for Python3.9 fi diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_proc_card.dat new file mode 100644 index 000000000000..8b54df96e466 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > w+{0} z{0} j j QED=4 QCD=0, w+ > j j, z > l- l+ @1 +add process p p > w-{0} z{0} j j QED=4 QCD=0, w- > j j, z > l- l+ @2 +output VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP/VBS_WZ_PolarizLL_WLtoHAD_ZLtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_proc_card.dat new file mode 100644 index 000000000000..6d020ea2aeec --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > w+{0} z{0} j j QED=4 QCD=0, w+ > l+ vl, z > j j @1 +add process p p > w-{0} z{0} j j QED=4 QCD=0, w- > l- vl~ , z > j j @2 +output VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizLL_semilep/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD/VBS_WZ_PolarizLL_WLtoLEP_ZLtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_proc_card.dat new file mode 100644 index 000000000000..58a459d9eb87 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vt = ve vm vt + +generate p p > w+{0} z{T} j j QED=4 QCD=0, w+ > j j, z > l+ l- @1 +add process p p > w-{0} z{T} j j QED=4 QCD=0, w- > j j , z > l+ l- @2 +output VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP/VBS_WZ_PolarizTL_WLtoHAD_ZTtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_proc_card.dat new file mode 100644 index 000000000000..e2001397492f --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > w+{0} z{T} j j QED=4 QCD=0, w+ > l+ vl, z > j j @ 1 +add process p p > w-{0} z{T} j j QED=4 QCD=0, w- > l- vl~ , z > j j @ 2 +output VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD --nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD/VBS_WZ_PolarizTL_WLtoLEP_ZTtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_proc_card.dat new file mode 100644 index 000000000000..3a9da3b51aa7 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vt = ve vm vt + +generate p p > w+{T} z{0} j j QED=4 QCD=0, w+ > j j, z > l+ l- @1 +add process p p > w-{T} z{0} j j QED=4 QCD=0, w- > j j , z > l+ l- @2 +output VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP/VBS_WZ_PolarizTL_WTtoHAD_ZLtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_proc_card.dat new file mode 100644 index 000000000000..16675d2ffb07 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > w+{T} z{0} j j QED=4 QCD=0, w+ > l+ vl, z > j j @1 +add process p p > w-{T} z{0} j j QED=4 QCD=0, w- > l- vl~ , z > j j @2 +output VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTL_semilep/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD/VBS_WZ_PolarizTL_WTtoLEP_ZLtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_proc_card.dat new file mode 100644 index 000000000000..aa51fc846688 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > w+{T} z{T} j j QED=4 QCD=0, w+ > j j, z > l+ l- @1 +add process p p > w-{T} z{T} j j QED=4 QCD=0, w- > j j , z > l+ l- @2 +output VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP/VBS_WZ_PolarizTT_WTtoHAD_ZTtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_proc_card.dat new file mode 100644 index 000000000000..3b093d3c64bb --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > w+{T} z{T} j j QED=4 QCD=0, w+ > l+ vl, z > j j @1 +add process p p > w-{T} z{T} j j QED=4 QCD=0, w- > l- vl~ , z > j j @2 +output VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/WZ/VBS_WZ_PolarizTT_semilep/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD/VBS_WZ_PolarizTT_WTtoLEP_ZTtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_proc_card.dat new file mode 100644 index 000000000000..22929c2d29a9 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_proc_card.dat @@ -0,0 +1,8 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > z{0} z{0} j j QED=4 QCD=0, z > j j, z > l+ l- +output VBS_ZZ_PolarizLL_semilep --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_run_card.dat new file mode 100644 index 000000000000..28973d478cf8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizLL_semilep/VBS_ZZ_PolarizLL_semilep_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 1 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_proc_card.dat new file mode 100644 index 000000000000..7061ba8bc40b --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_proc_card.dat @@ -0,0 +1,8 @@ +import model sm-ckm_no_b_mass_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > z{T} z{0} j j QED=4 QCD=0, z > j j, z > l+ l- +output VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP/VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_proc_card.dat new file mode 100644 index 000000000000..994d9a6ae562 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_proc_card.dat @@ -0,0 +1,8 @@ +import model sm-ckm_no_b_mass_no_b_mass +define vl~ = ve~ vm~ vt~ +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ +define vl = ve vm vt + +generate p p > z{0} z{T} j j QED=4 QCD=0, z > j j, z > l+ l- +output VBS_ZZ_PolarizTL_ZTtoHAD_ZLtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_run_card.dat new file mode 100644 index 000000000000..0d4ea8e3e0ef --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTL_semilep/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD/VBS_ZZ_PolarizTL_ZTtoLEP_ZLtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_proc_card.dat new file mode 100644 index 000000000000..e3be6625c7dc --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_proc_card.dat @@ -0,0 +1,8 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ tm~ +define vl = ve vm tm +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > z{T} z{T} j j QED=4 QCD=0, z > j j, z > l+ l- +output VBS_ZZ_PolarizTT_semilep --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_run_card.dat new file mode 100644 index 000000000000..28973d478cf8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ZZ/VBS_ZZ_PolarizTT_semilep/VBS_ZZ_PolarizTT_semilep_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 1 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_proc_card.dat new file mode 100644 index 000000000000..c1a50b5abbdd --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > w+{0} w+{0} j j QED=4 QCD=0, w+ > l+ vl, w+ > j j @1 +add process p p > w-{0} w-{0} j j QED=4 QCD=0, w- > l- vl~, w- > j j @2 +output VBS_ssWW_PolarizLL_semilep --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_run_card.dat new file mode 100644 index 000000000000..28973d478cf8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizLL_semilep/VBS_ssWW_PolarizLL_semilep_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 1 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_proc_card.dat new file mode 100644 index 000000000000..09bc94511c37 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > w+{T} w+{0} j j QED=4 QCD=0, w+ > l+ vl, w+ > j j @1 +add process p p > w-{T} w-{0} j j QED=4 QCD=0, w- > l- vl~, w- > j j @2 +output VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_run_card.dat new file mode 100644 index 000000000000..ed17250b2abf --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP/VBS_ssWW_PolarizTL_WTtoHAD_WLtoLEP_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 1 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_proc_card.dat new file mode 100644 index 000000000000..bbf358dea2f5 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l- = e- mu- ta- +define l+ = e+ mu+ ta+ + +generate p p > w+{0} w+{T} j j QED=4 QCD=0, w+ > j j, w+ > l+ vl @ 1 +add process p p > w-{0} w-{T} j j QED=4 QCD=0, w- > j j, w- > l- vl~ @ 2 +output VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD --nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_run_card.dat new file mode 100644 index 000000000000..701d73fc88ad --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTL_semilep/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD/VBS_ssWW_PolarizTL_WTtoLEP_WLtoHAD_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 10000 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_proc_card.dat new file mode 100644 index 000000000000..051ffcb3ef3f --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_proc_card.dat @@ -0,0 +1,9 @@ +import model sm-ckm_no_b_mass +define vl~ = ve~ vm~ vt~ +define vl = ve vm vt +define l+ = e+ mu+ ta+ +define l- = e- mu- ta- + +generate p p > w+{T} w+{T} j j QED=4 QCD=0, w+ > l+ vl, w+ > j j @1 +add process p p > w-{T} w-{T} j j QED=4 QCD=0, w- > l- vl~, w- > j j @2 +output VBS_ssWW_PolarizTT_semilep --nojpeg \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_run_card.dat new file mode 100644 index 000000000000..ed17250b2abf --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/VBS/VVjj_semilep_Polarized/ssWW/VBS_ssWW_PolarizTT_semilep/VBS_ssWW_PolarizTT_semilep_run_card.dat @@ -0,0 +1,263 @@ +#********************************************************************* +# 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 * +#********************************************************************* + 1 = 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 * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type + 6800.0 = ebeam1 ! beam 1 total energy in GeV + 6800.0 = ebeam2 ! beam 2 total energy in GeV +# To see polarised beam options: type "update beam_pol" +#********************************************************************* +# Frame where to evaluate the matrix-element (not the cut!) for polarization +#********************************************************************* + 1, 2 = me_frame ! list of particles to sum-up to define the rest-frame + ! in which to evaluate the matrix-element + ! [1,2] means the partonic center of mass + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set + $DEFAULT_PDF_SETS = lhaid + $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 +#********************************************************************* + False = 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 + +#********************************************************************* +# +#********************************************************************* +# handling of the helicities: +# 0: sum over all helicities +# 1: importance sampling over helicities +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. +#********************************************************************* +# 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) + #********************************************************************* + # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not + # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) + #********************************************************************* + True = cut_decays ! Cut decay products +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 10.0 = ptj ! minimum pt for the jets + 0.0 = ptb ! minimum pt for the b + 0.0 = pta ! minimum pt for the photons + 0.0 = ptl ! minimum pt for the charged leptons + 0.0 = misset ! minimum missing Et (sum of neutrino's momenta) + -1.0 = ptjmax ! maximum pt for the jets + -1.0 = ptbmax ! maximum pt for the b + -1.0 = ptamax ! maximum pt for the photons + -1.0 = ptlmax ! maximum pt for the charged leptons + -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta) + {} = 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}) +#********************************************************************* +# Minimum and maximum E's (in the center of mass frame) * +#********************************************************************* + 0.0 = ej ! minimum E for the jets + 0.0 = eb ! minimum E for the b + 0.0 = ea ! minimum E for the photons + 0.0 = el ! minimum E for the charged leptons + -1.0 = ejmax ! maximum E for the jets + -1.0 = ebmax ! maximum E for the b + -1.0 = eamax ! maximum E for the photons + -1.0 = elmax ! maximum E for the charged leptons + {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle + {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50}) + +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 6.5 = etaj ! max rap for the jets + -1 = etab ! max rap for the b + 2.5 = etaa ! max rap for the photons + -1.0 = etal ! max rap for the charged leptons + 0.0 = etajmin ! min rap for the jets + 0.0 = etabmin ! min rap for the b + 0.0 = etaamin ! min rap for the photons + 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.4 = drjj ! min distance between jets + 0.0 = drbb ! min distance between b's + 0.4 = drll ! min distance between leptons + 0.0 = draa ! min distance between gammas + 0.0 = drbj ! min distance between b and jet + 0.0 = draj ! min distance between gamma and jet + 0.0 = drab ! min distance between gamma and b + 0.0 = drbl ! min distance between b and lepton + 0.0 = dral ! min distance between gamma and lepton + 0.4 = drjl ! min distance between jet and lepton + -1.0 = drjjmax ! max distance between jets + -1.0 = drbbmax ! max distance between b's + -1.0 = drllmax ! max distance between leptons + -1.0 = draamax ! max distance between gammas + -1.0 = drbjmax ! max distance between b and jet + -1.0 = drajmax ! max distance between gamma and jet + -1.0 = drjlmax ! max distance between jet and lepton + -1.0 = drabmax ! max distance between gamma and b + -1.0 = drblmax ! max distance between b and lepton + -1.0 = dralmax ! maxdistance between gamma and lepton +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + 100.0 = mmjj ! min invariant mass of a jet pair + 0.0 = mmbb ! min invariant mass of a b pair + 0.0 = mmaa ! min invariant mass of gamma gamma 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 = mmbbmax ! max invariant mass of a b pair + -1.0 = mmaamax ! max invariant mass of gamma gamma 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 * + #********************************************************************* + 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 = xptb ! minimum pt for at least one b + 0.0 = xpta ! minimum pt for at least one photon + 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 + 0.0 = ptl3min ! minimum pt for the third lepton in pt + 0.0 = ptl4min ! minimum pt for the fourth 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 + -1.0 = ptl3max ! maximum pt for the third lepton in pt + -1.0 = ptl4max ! maximum pt for the fourth 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 + #*********************************************************************** + # Photon-isolation cuts, according to hep-ph/9801442 * + # When ptgmin=0, all the other parameters are ignored * + # When ptgmin>0, pta and draj are not going to be used * + #*********************************************************************** + 0.0 = ptgmin ! Min photon transverse momentum + 0.4 = 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) + #********************************************************************* + # 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 +#********************************************************************* +# KT DURHAM CUT * +#********************************************************************* + -1.0 = ktdurham + 0.4 = dparameter +#********************************************************************* +# 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 \ No newline at end of file