Merge pull request #3753 from PaulVaucelle/addMyFile

Add DisplacedTop Cards

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
+

bin/MadGraph5_aMCatNLO/cards/RPV_2018/RPV_PROD_smuXXX_neuYYY_extramodels.dat

@@ -0,0 +1 @@
+DisplacedTopUDD.zip

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

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'

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