diff --git a/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_customizecards.dat new file mode 100644 index 000000000000..7cee25b4605d --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_customizecards.dat @@ -0,0 +1,18 @@ +set param_card mass 1000006 1.000000e+03 +set param_card mass 2000006 1.000000e+03 +set param_card mass 1000013 XXX +set param_card mass 2000013 XXX +set param_card mass 1000023 YYY +set param_card SELMIX 1 1 7.071070e-01 # rrl1x1 +set param_card SELMIX 1 4 7.071070e-01 # rrl1x4 +set param_card SELMIX 2 2 7.071070e-01 # rrl2x2 +set param_card SELMIX 2 5 7.071070e-01 # rrl2x5 +set param_card SELMIX 3 3 7.071070e-01 # rrl3x3 +set param_card SELMIX 3 6 7.071070e-01 # rrl3x6 +set param_card SELMIX 4 1 7.071070e-01 # rrl4x1 +set param_card SELMIX 4 4 7.071070e-01 # rrl4x4 +set param_card SELMIX 5 2 7.071070e-01 # rrl5x2 +set param_card SELMIX 5 5 7.071070e-01 # rrl5x5 +set param_card SELMIX 6 3 7.071070e-01 # rrl6x3 +set param_card SELMIX 6 6 7.071070e-01 # rrl6x6 + diff --git a/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_extramodels.dat new file mode 100644 index 000000000000..6a96ff4d5ff8 --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_extramodels.dat @@ -0,0 +1 @@ +DisplacedTopUDD.zip \ No newline at end of file diff --git a/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_proc_card.dat new file mode 100644 index 000000000000..b7ac2e66eeee --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_proc_card.dat @@ -0,0 +1,74 @@ +#************************************************************ +#* MadGraph5_aMC@NLO * +#* * +#* * * * +#* * * * * * +#* * * * * 5 * * * * * +#* * * * * * +#* * * * +#* * +#* * +#* VERSION 2.9.15 2023-05-12 * +#* * +#* The MadGraph5_aMC@NLO Development Team - Find us at * +#* https://server06.fynu.ucl.ac.be/projects/madgraph * +#* * +#************************************************************ +#* * +#* Command File for MadGraph5_aMC@NLO * +#* * +#* run as ./bin/mg5_aMC filename * +#* * +#************************************************************ +set group_subprocesses Auto +set ignore_six_quark_processes False +set max_t_for_channel 99 +set loop_optimized_output True +set low_mem_multicore_nlo_generation False +set default_unset_couplings 99 +set zerowidth_tchannel True +set loop_color_flows False +set gauge unitary +set complex_mass_scheme False +set max_npoint_for_channel 0 +import model sm +define p = g u c d s u~ c~ d~ s~ +define j = g u c d s u~ c~ d~ s~ +define l+ = e+ mu+ +define l- = e- mu- +define vl = ve vm vt +define vl~ = ve~ vm~ vt~ +import model DisplacedTopUDD --modelname +generate p p > sl2+ sl2- / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @1 +add process p p > sl2+ sl2- / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @2 +add process p p > sl2+ sl2- / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @3 +add process p p > sl2+ sl2- / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @4 +add process p p > sl2+ sl2- j / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @5 +add process p p > sl2+ sl2- j / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @6 +add process p p > sl2+ sl2- j / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t d s / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @7 +add process p p > sl2+ sl2- j / h01 h02 a0 n1 n2 n3 n4, (sl2+ > n2 mu\ ++, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd4~ sd5 sd5~ sd6 sd\ +6~), (sl2- > n2 mu-, n2 > t~ d~ s~ / sd1 sd1~ sd2 sd2~ sd3 sd3~ sd4 sd\ +4~ sd5 sd5~ sd6 sd6~) @8 +output RPV_PROD_smu200_neu180 -f diff --git a/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_run_card.dat new file mode 100644 index 000000000000..c6e5b2a17eee --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_run_card.dat @@ -0,0 +1,191 @@ +#********************************************************************* +# 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 to_full * +#********************************************************************* +# +#********************************************************************* +# 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 * +#********************************************************************* +#$$ + 20000 = nevents ! Number of unweighted events requested +#$$ + 0 = iseed ! rnd seed (0=assigned automatically=default)) +#$$ +#********************************************************************* +# Collider type and energy * +# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * +# 3=photon from electron, 4=photon from muon * +#********************************************************************* + 1 = lpp1 ! beam 1 type + 1 = lpp2 ! beam 2 type +#$$ + 6500.0 = ebeam1 ! beam 1 total energy in GeV + 6500.0 = ebeam2 ! beam 2 total energy in GeV +#$$ +# To see polarised beam options: type "update beam_pol" + +#********************************************************************* +# PDF CHOICE: this automatically fixes also alpha_s and its evol. * +#********************************************************************* + lhapdf = pdlabel ! PDF set +#$$ + $DEFAULT_PDF_SETS = lhaid ! if pdlabel=lhapdf, this is the lhapdf number + $DEFAULT_PDF_MEMBERS = reweight_PDF +#$$ +# To see heavy ion options: type "update ion_pdf" +#********************************************************************* +# Renormalization and factorization scales * +#********************************************************************* + False = fixed_ren_scale ! if .true. use fixed ren scale + False = fixed_fac_scale ! if .true. use fixed fac scale + 91.188 = scale ! fixed ren scale + 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 + 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 + -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices + 1.0 = scalefact ! scale factor for event-by-event scales +#********************************************************************* +# Type and output format +#********************************************************************* + 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 +#********************************************************************* +# Matching parameter (MLM only) +#********************************************************************* +#$$ + 1 = ickkw ! 0 no matching, 1 MLM +#$$ + 1.0 = alpsfact ! scale factor for QCD emission vx + False = chcluster ! cluster only according to channel diag + 4 = asrwgtflavor ! highest quark flavor for a_s reweight + True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 + ! (turn off for VBF and single top processes) +#$$ + 30.0 = xqcut ! minimum kt jet measure between partons +#$$ + +#*********************************************************************** +# Turn on either the ktdurham or ptlund cut to activate * +# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] * +#*********************************************************************** + -1.0 = ktdurham + 0.4 = dparameter + -1.0 = ptlund + 1, 2, 3, 4, 5, 6, 21, 1000001, 1000002, 1000003, 1000004, 1000005, 1000006, 1000021, 2000001, 2000002, 2000003, 2000004, 2000005, 2000006 = pdgs_for_merging_cut ! PDGs for two cuts above + +#********************************************************************* +# +#********************************************************************* +# Phase-Space Optimization strategy (basic options) +#********************************************************************* + 0 = nhel ! using helicities importance sampling or not. + ! 0: sum over helicity, 1: importance sampling +#$$ 2 = sde_strategy ! default integration strategy (hep-ph/2021.00773) + ! 1 is old strategy (using amp square) + ! 2 is new strategy (using only the denominator) +#$$ + 2 = hard_survey +#$$ +# To see advanced option forPhase-Space optimization: type "update psoptim" +#********************************************************************* +# Generation bias, check the wiki page below for more information: * +# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' * +#********************************************************************* + None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-] + {} = bias_parameters ! Specifies the parameters of the module. +# +#******************************* +# Parton level cuts definition * +#******************************* +# +# +#********************************************************************* +# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay +#********************************************************************* + 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) +#********************************************************************* +# Standard Cuts * +#********************************************************************* +# Minimum and maximum pt's (for max, -1 means no cut) * +#********************************************************************* + 20.0 = ptj ! minimum pt for the jets + -1.0 = ptjmax ! maximum pt for the jets + {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle + {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) +# +# For display option for energy cut in the partonic center of mass frame type 'update ecut' +# +#********************************************************************* +# Maximum and minimum absolute rapidity (for max, -1 means no cut) * +#********************************************************************* + 5.0 = etaj ! max rap for the jets + {} = 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 * +#********************************************************************* +#********************************************************************* +# Minimum and maximum invariant mass for pairs * +#********************************************************************* + {} = 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. +#********************************************************************* +# Inclusive cuts * +#********************************************************************* + 0.0 = ptheavy ! minimum pt for at least one heavy final state +#********************************************************************* +# maximal pdg code for quark to be considered as a light jet * +# (otherwise b cuts are applied) * +#********************************************************************* + 4 = 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 +#$$ +# + systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++] + ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset', '--alps=0.5,1,2'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule +# Syscalc is deprecated but to see the associate options type'update syscalc' + +#*********************************************************************** +# Turn on either the ktdurham or ptlund cut to activate * +# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] * +#*********************************************************************** + +#********************************************************************* +# Phase-Space Optim (advanced) +#********************************************************************* + 0 = job_strategy ! see appendix of 1507.00020 (page 26) + -1.0 = tmin_for_channel ! limit the non-singular reach of --some-- channel of integration related to T-channel diagram (value between -1 and 0), -1 is no impact + -1 = survey_splitting ! for loop-induced control how many core are used at survey for the computation of a single iteration. + 2 = survey_nchannel_per_job ! control how many Channel are integrated inside a single job on cluster/multicore + -1 = refine_evt_by_job ! control the maximal number of events for the first iteration of the refine (larger means less jobs) +#********************************************************************* +# Compilation flag. No automatic re-compilation (need manual "make clean" in Source) +#********************************************************************* + -O = global_flag ! fortran optimization flag use for the all code. + = aloha_flag ! fortran optimization flag for aloha function. Suggestions: '-ffast-math' + = matrix_flag ! fortran optimization flag for matrix.f function. Suggestions: '-O3' diff --git a/bin/MadGraph5_aMCatNLO/cards/RPV_2018/set_cards.sh b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/set_cards.sh new file mode 100644 index 000000000000..3e379468a11b --- /dev/null +++ b/bin/MadGraph5_aMCatNLO/cards/RPV_2018/set_cards.sh @@ -0,0 +1,38 @@ +#!/bin/bash + +# Define of Smuon Mass (1) and Neutralino Mass (2) +list1=(200 250 300 350 400 450 500) +list2=(180 200 230 250 280 300 330 350 380 400 430 450 480) + +echo "Pull Request : Add DisplacedTop cards" + +# Loop over both lists +for smuon in "${list1[@]}"; do + for neu in "${list2[@]}"; do + # Check if neutralino mass is smaller than smuon mass (with specific values) + dm=$((smuon - neu)) + if [[ $neu -lt $smuon ]] && [[ !(($dm -eq 70 && $neu -ne 180) || \ + ($dm -eq 120 && $neu -ne 180) || \ + ($dm -eq 170 && $neu -ne 180) || \ + ($dm -eq 220 && $neu -ne 180) || \ + ($dm -eq 270 && $neu -ne 180)) ]]; then + + # Copy the reference param_card to create the one of the new sample + # cp ./RPV_PROD_smuXXX_neuYYY_param_card.dat ./RPV_PROD_smu${smuon}_neu${neu}_param_card.dat + cp ./RPV_PROD_smuXXX_neuYYY_run_card.dat ./RPV_PROD_smu${smuon}_neu${neu}_run_card.dat + cp ./RPV_PROD_smuXXX_neuYYY_proc_card.dat ./RPV_PROD_smu${smuon}_neu${neu}_proc_card.dat + cp ./RPV_PROD_smuXXX_neuYYY_customizecards.dat ./RPV_PROD_smu${smuon}_neu${neu}_customizecards.dat + cp ./RPV_PROD_smuXXX_neuYYY_extramodels.dat ./RPV_PROD_smu${smuon}_neu${neu}_extramodels.dat + + # Edit the param_card.dat based on smuon and neutralino masses + + sed -i "s/XXX/${smuon}/g" ./RPV_PROD_smu${smuon}_neu${neu}_customizecards.dat + sed -i "s/YYY/${neu}/g" ./RPV_PROD_smu${smuon}_neu${neu}_customizecards.dat + + # sed -i "s/XXX/${smuon}/g" ./RPV_PROD_smu${smuon}_neu${neu}_param_card.dat + # sed -i "s/YYY/${neu}/g" ./RPV_PROD_smu${smuon}_neu${neu}_param_card.dat + + echo "Created and modified files for smuon mass $smuon and neutralino mass $neu" + fi + done +done diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_customizecards.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_customizecards.dat deleted file mode 100644 index a93dc506881f..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_customizecards.dat +++ /dev/null @@ -1,12 +0,0 @@ -set param_card frblock 1 1.00e+00 -set param_card frblock 2 9.9e-01 -set param_card higgs 3 0.00e+00 -set param_card mass 25 125.0 -set param_card mass 35 HMASS -set param_card mass 36 AMASS -set param_card mass 37 AMASS -set param_card width 35 HWIDTH -set param_card width 36 AWIDTH -set param_card width 37 AWIDTH -set param_card higgs 1 lambda2 -set param_card higgs 2 lambda3 diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_extramodels.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_extramodels.dat deleted file mode 100644 index 36ee330b5a18..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_extramodels.dat +++ /dev/null @@ -1,2 +0,0 @@ -#http://feynrules.irmp.ucl.ac.be/attachment/wiki/2HDM/2HDMtII_NLO.tar.gz -2HDMtII_NLO.tar.gz diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_madspin_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_madspin_card.dat deleted file mode 100644 index dcff333a980c..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_madspin_card.dat +++ /dev/null @@ -1,19 +0,0 @@ -set ms_dir ./madspingrid - -#initialization parameters -#set Nevents_for_max_weight 250 # number of events for the estimate of the max. weight -#set BW_cut 15 # cut on how far the particle can be off-shell -set max_weight_ps_point 400 # number of PS to estimate the maximum for each event -set max_running_process 1 - -define v = ve vm vt -define v~ = ve~ vm~ vt~ -define f = u c d s u~ c~ d~ s~ ve vm vt e- mu- ta- ve~ vm~ vt~ e+ mu+ ta+ - -# specify the decay for the final state particles -set spinmode none -decay h2 > t t~, ( t > w+ b, w+ > f f ), ( t~ > w- b~, w- > f f ) -decay z > v v~ - -# running the actual code -launch diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_proc_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_proc_card.dat deleted file mode 100644 index 490ef18534bd..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_proc_card.dat +++ /dev/null @@ -1,9 +0,0 @@ -set group_subprocesses Auto -set ignore_six_quark_processes False -#set loop_color_flows False -set gauge unitary -set complex_mass_scheme False - -import model 2HDMtII_NLO -generate p p > h3 > h2 z [QCD] -output AToZHToNuNuTTbar -nojpeg diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_run_card.dat b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_run_card.dat deleted file mode 100644 index 3ad154aa1916..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/AToZHToNuNuTTbar_run_card.dat +++ /dev/null @@ -1,269 +0,0 @@ -#********************************************************************* -# 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 * -#********************************************************************* -# -#******************* -# Running parameters -#******************* -# -#********************************************************************* -# Tag name for the run (one word) * -#********************************************************************* - tag_1 = run_tag ! name of the run -#********************************************************************* -# Run to generate the grid pack * -#********************************************************************* - .true. = gridpack !True = setting up the grid pack -#********************************************************************* -# Number of events and rnd seed * -# Warning: Do not generate more than 1M events in a single run * -# If you want to run Pythia, avoid more than 50k events in a run. * -#********************************************************************* - 2500 = 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 = ebeam1 ! beam 1 total energy in GeV - 6800 = ebeam2 ! beam 2 total energy in GeV -#********************************************************************* -# Beam polarization from -100 (left-handed) to 100 (right-handed) * -#********************************************************************* - 0 = polbeam1 ! beam polarization for beam 1 - 0 = polbeam2 ! beam polarization for beam 2 -#********************************************************************* -# PDF CHOICE: this automatically fixes also alpha_s and its evol. * -#********************************************************************* - 'lhapdf' = pdlabel ! PDF set - $DEFAULT_PDF_SETS = lhaid - $DEFAULT_PDF_MEMBERS = reweight_PDF ! if pdlabel=lhapdf, this is the lhapdf number -#********************************************************************* -# Renormalization and factorization scales * -#********************************************************************* - F = fixed_ren_scale ! if .true. use fixed ren scale - F = fixed_fac_scale ! if .true. use fixed fac scale - 91.1880 = scale ! fixed ren scale - 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1 - 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2 - 1 = scalefact ! scale factor for event-by-event scales -#********************************************************************* -# Matching - Warning! ickkw > 1 is still beta -#********************************************************************* - 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching - 1 = highestmult ! for ickkw=2, highest mult group - 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE - 1 = alpsfact ! scale factor for QCD emission vx - F = chcluster ! cluster only according to channel diag - F = pdfwgt ! for ickkw=1, perform pdf reweighting - 4 = asrwgtflavor ! highest quark flavor for a_s reweight - T = clusinfo ! include clustering tag in output - 3.0 = lhe_version ! Change the way clustering information pass to shower. -#********************************************************************* -#********************************************************** -# -#********************************************************** -# Automatic ptj and mjj cuts if xqcut > 0 -# (turn off for VBF and single top processes) -#********************************************************** - T = auto_ptj_mjj ! Automatic setting of ptj and mjj -#********************************************************** -# -#********************************** -# BW cutoff (M+/-bwcutoff*Gamma) -#********************************** - 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) -#********************************************************** -# Apply pt/E/eta/dr/mij cuts on decay products or not -# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) -#********************************************************** - T = cut_decays ! Cut decay products -#************************************************************* -# Number of helicities to sum per event (0 = all helicities) -# 0 gives more stable result, but longer run time (needed for -# long decay chains e.g.). -# Use >=2 if most helicities contribute, e.g. pure QCD. -#************************************************************* - 0 = nhel ! Number of helicities used per event -#******************* -# Standard Cuts -#******************* -# -#********************************************************************* -# Minimum and maximum pt's (for max, -1 means no cut) * -#********************************************************************* - 10 = ptj ! minimum pt for the jets - 0 = ptb ! minimum pt for the b - 10 = pta ! minimum pt for the photons - 0 = ptl ! minimum pt for the charged leptons - 0 = misset ! minimum missing Et (sum of neutrino's momenta) - 0 = ptheavy ! minimum pt for one heavy final state - -1 = ptjmax ! maximum pt for the jets - -1 = ptbmax ! maximum pt for the b - -1 = ptamax ! maximum pt for the photons - -1 = ptlmax ! maximum pt for the charged leptons - -1 = missetmax ! maximum missing Et (sum of neutrino's momenta) -#********************************************************************* -# Minimum and maximum E's (in the center of mass frame) * -#********************************************************************* - 0 = ej ! minimum E for the jets - 0 = eb ! minimum E for the b - 0 = ea ! minimum E for the photons - 0 = el ! minimum E for the charged leptons - -1 = ejmax ! maximum E for the jets - -1 = ebmax ! maximum E for the b - -1 = eamax ! maximum E for the photons - -1 = elmax ! maximum E for the charged leptons -#********************************************************************* -# Maximum and minimum absolute rapidity (for max, -1 means no cut) * -#********************************************************************* - 5 = etaj ! max rap for the jets - -1 = etab ! max rap for the b - -1 = etaa ! max rap for the photons - -1 = etal ! max rap for the charged leptons - 0 = etajmin ! min rap for the jets - 0 = etabmin ! min rap for the b - 0 = etaamin ! min rap for the photons - 0 = etalmin ! main rap for the charged leptons -#********************************************************************* -# Minimum and maximum DeltaR distance * -#********************************************************************* - 0.1 = drjj ! min distance between jets - 0 = drbb ! min distance between b's - 0 = drll ! min distance between leptons - 0 = draa ! min distance between gammas - 0 = drbj ! min distance between b and jet - 0.1 = draj ! min distance between gamma and jet - 0 = drjl ! min distance between jet and lepton - 0 = drab ! min distance between gamma and b - 0 = drbl ! min distance between b and lepton - 0.1 = dral ! min distance between gamma and lepton - -1 = drjjmax ! max distance between jets - -1 = drbbmax ! max distance between b's - -1 = drllmax ! max distance between leptons - -1 = draamax ! max distance between gammas - -1 = drbjmax ! max distance between b and jet - -1 = drajmax ! max distance between gamma and jet - -1 = drjlmax ! max distance between jet and lepton - -1 = drabmax ! max distance between gamma and b - -1 = drblmax ! max distance between b and lepton - -1 = dralmax ! maxdistance between gamma and lepton -#********************************************************************* -# Minimum and maximum invariant mass for pairs * -# WARNING: for four lepton final state mmll cut require to have * -# different lepton masses for each flavor! * -#********************************************************************* - 0 = mmjj ! min invariant mass of a jet pair - 0 = mmbb ! min invariant mass of a b pair - 0 = mmaa ! min invariant mass of gamma gamma pair - 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair - -1 = mmjjmax ! max invariant mass of a jet pair - -1 = mmbbmax ! max invariant mass of a b pair - -1 = mmaamax ! max invariant mass of gamma gamma pair - -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair -#********************************************************************* -# Minimum and maximum invariant mass for all letpons * -#********************************************************************* - 0 = mmnl ! min invariant mass for all letpons (l+- and vl) - -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl) -#********************************************************************* -# Minimum and maximum pt for 4-momenta sum of leptons * -#********************************************************************* - 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) - -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) -#********************************************************************* -# Inclusive cuts * -#********************************************************************* - 0 = xptj ! minimum pt for at least one jet - 0 = xptb ! minimum pt for at least one b - 0 = xpta ! minimum pt for at least one photon - 0 = xptl ! minimum pt for at least one charged lepton -#********************************************************************* -# Control the pt's of the jets sorted by pt * -#********************************************************************* - 0 = ptj1min ! minimum pt for the leading jet in pt - 0 = ptj2min ! minimum pt for the second jet in pt - 0 = ptj3min ! minimum pt for the third jet in pt - 0 = ptj4min ! minimum pt for the fourth jet in pt - -1 = ptj1max ! maximum pt for the leading jet in pt - -1 = ptj2max ! maximum pt for the second jet in pt - -1 = ptj3max ! maximum pt for the third jet in pt - -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt - 0 = ptl2min ! minimum pt for the second lepton in pt - 0 = ptl3min ! minimum pt for the third lepton in pt - 0 = ptl4min ! minimum pt for the fourth lepton in pt - -1 = ptl1max ! maximum pt for the leading lepton in pt - -1 = ptl2max ! maximum pt for the second lepton in pt - -1 = ptl3max ! maximum pt for the third lepton in pt - -1 = ptl4max ! maximum pt for the fourth lepton in pt -#********************************************************************* -# Control the Ht(k)=Sum of k leading jets * -#********************************************************************* - 0 = htjmin ! minimum jet HT=Sum(jet pt) - -1 = htjmax ! maximum jet HT=Sum(jet pt) - 0 = ihtmin !inclusive Ht for all partons (including b) - -1 = ihtmax !inclusive Ht for all partons (including b) - 0 = ht2min ! minimum Ht for the two leading jets - 0 = ht3min ! minimum Ht for the three leading jets - 0 = ht4min ! minimum Ht for the four leading jets - -1 = ht2max ! maximum Ht for the two leading jets - -1 = ht3max ! maximum Ht for the three leading jets - -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case - 0 = deltaeta ! minimum rapidity for two jets in the WBF case -#********************************************************************* -# KT DURHAM CUT * -#********************************************************************* - -1 = ktdurham - 0.4 = dparameter -#********************************************************************* -# maximal pdg code for quark to be considered as a light jet * -# (otherwise b cuts are applied) * -#********************************************************************* -4 = maxjetflavor ! Maximum jet pdg code -#********************************************************************* -# Jet measure cuts * -#********************************************************************* -0 = xqcut ! minimum kt jet measure between partons -#********************************************************************* -# -#********************************************************************* -# Store info for systematics studies * -# WARNING: If use_syst is T, matched Pythia output is * -# meaningful ONLY if plotted taking matchscale * -# reweighting into account! * -#********************************************************************* - T = use_syst ! Enable systematics studies -# -#************************************** diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/copy.sh b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/copy.sh deleted file mode 100755 index a731bf2c02aa..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/copy.sh +++ /dev/null @@ -1,35 +0,0 @@ -#!/bin/bash - -default=$(basename $(pwd)) - -skip_headers=1 -while IFS=, read -r MA MH lambda2 lambda3 -do - if ((skip_headers)) - then - ((skip_headers--)) - else - echo "Creating cards for: MA=$MA MH=$MH" - - # specify the natural line width ks=width/mass - ks=0.03 - - newdir="$default"_MA-"$MA"_MH-"$MH" - mkdir $newdir - cp "$default"_customizecards.dat $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - cp "$default"_extramodels.dat $newdir/"$default"_MA-"$MA"_MH-"$MH"_extramodels.dat - cp "$default"_madspin_card.dat $newdir/"$default"_MA-"$MA"_MH-"$MH"_madspin_card.dat - cp "$default"_proc_card.dat $newdir/"$default"_MA-"$MA"_MH-"$MH"_proc_card.dat - cp "$default"_run_card.dat $newdir/"$default"_MA-"$MA"_MH-"$MH"_run_card.dat - # modify output name - sed -i "s/"$default"/"$default"_MA-"$MA"_MH-"$MH"/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_proc_card.dat - # Modify mass parameter - sed -i "s/AMASS/"$MA".0/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - sed -i "s/HMASS/"$MH".0/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - sed -i "s/HWIDTH/"$(bc <<< "$ks * $MH")"/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - sed -i "s/AWIDTH/"$(bc <<< "$ks * $MA")"/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - sed -i "s/lambda2/"$lambda2"/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - sed -i "s/lambda3/"$lambda3"/g" $newdir/"$default"_MA-"$MA"_MH-"$MH"_customizecards.dat - - fi -done < points.csv diff --git a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/points.csv b/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/points.csv deleted file mode 100644 index cc123b3453e8..000000000000 --- a/bin/MadGraph5_aMCatNLO/cards/production/13p6TeV/AToZHToNuNuTTbar/points.csv +++ /dev/null @@ -1,247 +0,0 @@ -mA,mH,lambda2,lambda3 -500,370,0.257734,3.98889 -550,350,0.257734,6.19592 -550,400,0.257734,4.9588 -600,350,0.257734,8.09284 -600,400,0.257734,6.85572 -600,450,0.257734,5.45365 -650,350,0.257734,10.1547 -650,400,0.257734,8.91759 -650,450,0.257734,7.51552 -650,500,0.257734,5.9485 -700,450,0.257734,9.74234 -700,500,0.257734,8.17532 -700,550,0.257734,6.44334 -750,400,0.257734,13.5362 -750,350,0.257734,14.7733 -750,450,0.257734,12.1341 -750,500,0.257734,10.5671 -750,550,0.257734,8.83511 -750,600,0.257734,6.93819 -800,350,0.257734,17.33 -800,400,0.257734,16.0929 -800,450,0.257734,14.6908 -800,500,0.257734,13.1238 -800,550,0.257734,11.3918 -850,400,0.257734,18.8146 -850,350,0.257734,20.0517 -850,450,0.257734,17.4125 -850,500,0.257734,15.8455 -850,550,0.257734,14.1135 -850,600,0.257734,12.2166 -850,650,0.257734,10.1547 -850,700,0.257734,7.92789 -900,450,0.257734,20.2991 -900,500,0.257734,18.7321 -900,550,0.257734,17.0001 -900,600,0.257734,15.1032 -900,650,0.257734,13.0413 -900,700,0.257734,10.8145 -900,750,0.257734,8.42274 -950,400,0.257734,24.7528 -950,350,0.257734,25.9899 -950,450,0.257734,23.3507 -950,500,0.257734,21.7837 -950,550,0.257734,20.0517 -950,600,0.257734,18.1548 -950,650,0.257734,16.0929 -950,700,0.257734,13.8661 -950,750,0.257734,11.4743 -950,800,0.257734,8.91759 -1000,350,0.257734,29.2064 -1000,400,0.257734,27.9693 -1000,450,0.257734,26.5672 -1000,500,0.257734,25.0002 -1000,550,0.257734,23.2682 -1000,650,0.257734,19.3094 -1000,750,0.257734,14.6908 -1000,850,0.257734,9.41244 -1050,350,0.257734,32.5879 -1050,400,0.257734,31.3507 -1050,450,0.257734,29.9487 -1050,500,0.257734,28.3816 -1050,550,0.257734,26.6497 -1050,600,0.257734,24.7528 -1050,650,0.257734,22.6909 -1050,700,0.257734,20.4641 -1050,750,0.257734,18.0723 -1050,800,0.257734,15.5156 -1050,850,0.257734,12.7939 -1050,900,0.257734,9.90729 -1100,350,0.257734,36.1343 -1100,400,0.257734,34.8972 -1100,450,0.257734,33.4951 -1100,500,0.257734,31.9281 -1100,550,0.257734,30.1961 -1100,600,0.257734,28.2992 -1100,650,0.257734,26.2373 -1100,700,0.257734,24.0105 -1100,750,0.257734,21.6187 -1100,800,0.257734,19.062 -1100,850,0.257734,16.3403 -1100,900,0.257734,13.4537 -1100,950,0.257734,10.4021 -1150,350,0.257734,39.8456 -1150,450,0.257734,37.2065 -1150,550,0.257734,33.9075 -1150,650,0.257734,29.9487 -1150,750,0.257734,25.3301 -1150,850,0.257734,20.0517 -1150,950,0.257734,14.1135 -1200,350,0.257734,43.722 -1200,400,0.257734,42.48484 -1200,500,0.257734,39.5157 -1200,600,0.257734,35.8869 -1200,700,0.257734,31.5982 -1200,800,0.257734,26.6497 -1200,900,0.257734,21.0414 -450,350,0.257734,2.89693 -500,400,0.257734,3.22683 -550,450,0.257734,3.55673 -600,500,0.257734,3.88663 -650,550,0.257734,4.21652 -700,600,0.257734,4.54642 -750,650,0.257734,4.87632 -800,700,0.257734,5.20622 -850,750,0.257734,5.53612 -900,800,0.257734,5.86602 -950,850,0.257734,6.19592 -1000,900,0.257734,6.52582 -1050,950,0.257734,6.85572 -1100,1000,0.257734,7.18562 -1150,1050,0.257734,7.51552 -1200,1100,0.257734,7.84542 -450,330,0.257734,3.34559 -500,330,0.257734,4.91261 -550,330,0.257734,6.64458 -600,330,0.257734,8.5415 -650,330,0.257734,10.6034 -700,330,0.257734,12.8302 -750,330,0.257734,15.222 -800,330,0.257734,17.7787 -850,330,0.257734,20.5004 -900,330,0.257734,23.387 -950,330,0.257734,26.4385 -1000,330,0.257734,29.6551 -1050,330,0.257734,33.0365 -1100,330,0.257734,36.5829 -1150,330,0.257734,40.2943 -1200,330,0.257734,44.1706 -430,330,0.257734,2.76497 -1300,400,0.257734,50.7323 -1300,500,0.257734,47.7632 -1300,600,0.257734,44.1343 -1300,700,0.257734,39.8456 -1300,800,0.257734,34.8972 -1300,900,0.257734,29.2889 -1300,1000,0.257734,23.0208 -1300,1100,0.257734,16.0929 -1300,1200,0.257734,8.50522 -1400,400,0.257734,59.6396 -1400,500,0.257734,56.6705 -1400,600,0.257734,53.0416 -1400,700,0.257734,48.7529 -1400,800,0.257734,43.8044 -1400,900,0.257734,38.1961 -1400,1000,0.257734,31.9281 -1400,1100,0.257734,25.0002 -1400,1200,0.257734,17.4125 -1400,1300,0.257734,9.16501 -1500,500,0.257734,66.2376 -1500,600,0.257734,62.6087 -1500,700,0.257734,58.32 -1500,800,0.257734,53.3715 -1500,900,0.257734,47.7632 -1500,1100,0.257734,34.5673 -1500,1200,0.257734,26.9796 -1500,1300,0.257734,18.7321 -1600,400,0.257734,79.4336 -1600,500,0.257734,76.4645 -1600,600,0.257734,72.8356 -1600,700,0.257734,68.5469 -1600,800,0.257734,63.5984 -1600,900,0.257734,57.9901 -1600,1000,0.257734,51.722 -1600,1100,0.257734,44.7941 -1600,1200,0.257734,37.2065 -1600,1300,0.257734,28.959 -1600,1400,0.257734,20.0517 -1600,1500,0.257734,10.4846 -1700,400,0.257734,90.3202 -1700,500,0.257734,87.3511 -1700,600,0.257734,83.7222 -1700,700,0.257734,79.4336 -1700,800,0.257734,74.4851 -1700,900,0.257734,68.8768 -1700,1000,0.257734,62.6087 -1700,1100,0.257734,55.6808 -1700,1200,0.257734,48.0931 -1700,1300,0.257734,39.8456 -1700,1400,0.257734,30.9384 -1700,1500,0.257734,21.3713 -1700,1600,0.257734,11.1444 -1800,500,0.257734,98.8976 -1800,600,0.257734,95.2687 -1800,700,0.257734,90.98 -1800,800,0.257734,86.0315 -1800,900,0.257734,80.4233 -1800,1000,0.257734,74.1552 -1800,1100,0.257734,67.2273 -1800,1200,0.257734,59.6396 -1800,1300,0.257734,51.3921 -1800,1400,0.257734,42.4848 -1800,1500,0.257734,32.9178 -1800,1700,0.257734,11.8042 -1900,400,0.257734,114.073 -1900,500,0.257734,111.104 -1900,600,0.257734,107.475 -1900,700,0.257734,103.186 -1900,800,0.257734,98.2378 -1900,900,0.257734,92.6295 -1900,1000,0.257734,86.3614 -1900,1100,0.257734,79.4336 -1900,1200,0.257734,71.8459 -1900,1300,0.257734,63.5984 -1900,1400,0.257734,54.6911 -1900,1500,0.257734,45.124 -1900,1600,0.257734,34.8972 -1900,1700,0.257734,24.0105 -1900,1800,0.257734,12.464 -2000,400,0.257734,126.939 -2000,500,0.257734,123.97 -2000,600,0.257734,120.341 -2000,700,0.257734,116.052 -2000,800,0.257734,111.104 -2000,900,0.257734,105.496 -2000,1000,0.257734,99.2275 -2000,1100,0.257734,92.2996 -2000,1200,0.257734,84.7119 -2000,1300,0.257734,76.4645 -2000,1400,0.257734,67.5572 -2000,1500,0.257734,57.9901 -2000,1600,0.257734,47.7632 -2000,1700,0.257734,36.8766 -2000,1800,0.257734,25.3301 -2000,1900,0.257734,13.1238 -2100,500,0.257734,137.496 -2100,600,0.257734,133.867 -2100,700,0.257734,129.578 -2100,800,0.257734,124.63 -2100,900,0.257734,119.021 -2100,1100,0.257734,105.825 -2100,1200,0.257734,98.2378 -2100,1300,0.257734,89.9903 -2100,1400,0.257734,81.0831 -2100,1500,0.257734,71.516 -2100,1700,0.257734,50.4024 -2100,1800,0.257734,38.8559 -2100,1900,0.257734,26.6497 -1300,350,0.257734,51.9694 -1400,350,0.257734,60.8767 -1500,350,0.257734,70.4438 -1600,350,0.257734,80.6707 -1700,350,0.257734,91.5574 -1800,350,0.257734,103.104 -1900,350,0.257734,115.31 -2000,350,0.257734,128.176 -2100,350,0.257734,141.702 diff --git a/bin/utils/run_req_new.sh b/bin/utils/run_req_new.sh deleted file mode 100644 index 7defa47e16f3..000000000000 --- a/bin/utils/run_req_new.sh +++ /dev/null @@ -1,8 +0,0 @@ -python3.9 -m venv venv && source ./venv/bin/activate -echo "@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" -echo "If necessary, run:" -echo "pip install git+https://github.com/cms-PdmV/mcm_scripts.git" -echo "@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@" -echo -#pip install git+https://github.com/cms-PdmV/mcm_scripts.git -python3 request_fragment_check.py --bypass_status --develop --prepid $1