ram1123 · ram1123 · Feb 16, 2021 · Jun 25, 2022 · Nov 20, 2023 · May 3, 2024
diff --git a/.gitignore b/.gitignore
@@ -1,4 +1,12 @@
-# User defined 
+# User defined
+*.err
+*.out
+MG5_aMC_v*/
+py.py
+MG5_debug
+CMSSW_*/
+rasharma.cc
+*.tar.gz
 docker_stderror
 screen*
 output_logs

diff --git a/Cards/SM_process/st_sch_4f_ckm_NLO/mg5_commands.txt b/Cards/SM_process/st_sch_4f_ckm_NLO/mg5_commands.txt
@@ -0,0 +1,9 @@
+launch
+   0                  # Turn off showering/hadronization
+   set run_tag "test"
+   set use_syst False
+   set nevents 50
+   set iseed 11
+   # set reweight_pdf False
+   # set store_rwgt_info True
+   # reweight
diff --git a/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_customizecards.dat b/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_customizecards.dat
@@ -0,0 +1,3 @@
+#put card customizations here (change top and higgs mass for example)
+set param_card mass 6 172.5
+set param_card yukawa 6 172.5
diff --git a/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_madspin_card.dat b/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_madspin_card.dat
@@ -0,0 +1,40 @@
+#************************************************************
+#*                        MadSpin                           *
+#*                                                          *
+#*    P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer * 
+#*                                                          *
+#*    Part of the MadGraph5_aMC@NLO Framework:              *
+#*    The MadGraph5_aMC@NLO Development Team - Find us at   *
+#*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
+#*                                                          *
+#************************************************************
+#Some options (uncomment to apply)
+#
+#directory for gridpack mode
+set ms_dir ./madspingrid
+
+#initialization parameters
+ set Nevents_for_max_weigth 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
+# 
+#to properly limit the number of concurrent processes for grid running
+set max_running_process 1
+#
+
+define ell+ = e+ mu+ ta+
+define ell- = e- mu- ta-
+# specify the decay for the final state particles
+decay t > w+ b, w+ > ell+ vl
+decay t~ > w- b~, w- > ell- vl~
+decay w+ > ell+ vl
+decay w- > ell- vl~
+
+#define ell+ = e+ mu+ ta+
+#define ell- = e- mu- ta-
+
+# specify the decay for the final state particles
+#decay w+ > ell+ vl
+#decay w- > ell- vl~
+# running the actual code
+launch
diff --git a/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_proc_card.dat b/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_proc_card.dat
@@ -0,0 +1,18 @@
+import model loop_sm-ckm
+#switch to non-diagonal ckm matrix for speed
+#import model loop_sm
+
+generate p p > w+ > t b~ [QCD] @0
+add process p p > w- > t~ b  [QCD] @1
+
+output st_sch_4f_ckm_NLO -nojpeg
+
+launch
+   0                  # Turn off showering/hadronization
+   set run_tag "test"
+   set use_syst False
+   set nevents 50
+   set iseed 11
+   # set reweight_pdf False
+   # set store_rwgt_info True
+   # reweight
diff --git a/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_run_card.dat b/Cards/SM_process/st_sch_4f_ckm_NLO/st_sch_4f_ckm_NLO_run_card.dat
@@ -0,0 +1,146 @@
+#***********************************************************************
+#                        MadGraph5_aMC@NLO                             *
+#                                                                      *
+#                      run_card.dat aMC@NLO                            *
+#                                                                      *
+#  This file is used to set the parameters of the run.                 *
+#                                                                      *
+#  Some notation/conventions:                                          *
+#                                                                      *
+#   Lines starting with a hash (#) are info or comments                *
+#                                                                      *
+#   mind the format:   value    = variable     ! comment               *
+#***********************************************************************
+#
+#*******************                                                 
+# Running parameters
+#*******************                                                 
+#
+#***********************************************************************
+# Tag name for the run (one word)                                      *
+#***********************************************************************
+  tag_1     = run_tag ! name of the run 
+#***********************************************************************
+# Number of events (and their normalization) and the required          *
+# (relative) accuracy on the Xsec.                                     *
+# These values are ignored for fixed order runs                        *
+#***********************************************************************
+  1500 = nevents ! Number of unweighted events requested 
+ 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents)
+   -1  = nevt_job! Max number of events per job in event generation. 
+                 !  (-1= no split).
+average = event_norm ! Normalize events to sum or average to the X sect.
+#***********************************************************************
+# Number of points per itegration channel (ignored for aMC@NLO runs)   *
+#***********************************************************************
+ 0.01   = req_acc_FO       ! Required accuracy (-1=ignored, and use the 
+ 	                   ! number of points and iter. below)
+# These numbers are ignored except if req_acc_FO is equal to -1
+ 5000   = npoints_FO_grid  ! number of points to setup grids
+ 4      = niters_FO_grid   ! number of iter. to setup grids
+ 10000  = npoints_FO       ! number of points to compute Xsec
+ 6      = niters_FO        ! number of iter. to compute Xsec
+#***********************************************************************
+# Random number seed                                                   *
+#***********************************************************************
+     0    = iseed       ! rnd seed (0=assigned automatically=default))
+#***********************************************************************
+# Collider type and energy                                             *
+#***********************************************************************
+    1   = lpp1    ! beam 1 type (0 = no PDF)
+    1   = lpp2    ! beam 2 type (0 = no PDF)
+ 6500   = ebeam1  ! beam 1 energy in GeV
+ 6500   = ebeam2  ! beam 2 energy in GeV
+#***********************************************************************
+# PDF choice: this automatically fixes also alpha_s(MZ) and its evol.  *
+#***********************************************************************
+ lhapdf    = pdlabel   ! PDF set                                     
+ 292000    = lhaid     ! if pdlabel=lhapdf, this is the lhapdf number
+#***********************************************************************
+# Include the NLO Monte Carlo subtr. terms for the following parton    *
+# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8)         *
+# WARNING: PYTHIA6PT works only for processes without FSR!!!!          *
+#***********************************************************************
+  PYTHIA8   = parton_shower
+#***********************************************************************
+# Renormalization and factorization scales                             *
+# (Default functional form for the non-fixed scales is the sum of      *
+# the transverse masses of all final state particles and partons. This *
+# can be changed in SubProcesses/set_scales.f)                         *
+#***********************************************************************
+ T        = fixed_ren_scale  ! if .true. use fixed ren scale
+ T        = fixed_fac_scale  ! if .true. use fixed fac scale
+ 172.5   = muR_ref_fixed    ! fixed ren reference scale 
+ 172.5   = muF1_ref_fixed   ! fixed fact reference scale for pdf1
+ 172.5   = muF2_ref_fixed   ! fixed fact reference scale for pdf2
+#***********************************************************************
+# Renormalization and factorization scales (advanced and NLO options)  *
+#***********************************************************************
+ T        = fixed_QES_scale  ! if .true. use fixed Ellis-Sexton scale
+ 172.5   = QES_ref_fixed    ! fixed Ellis-Sexton reference scale
+ 1        = muR_over_ref     ! ratio of current muR over reference muR
+ 1        = muF1_over_ref    ! ratio of current muF1 over reference muF1
+ 1        = muF2_over_ref    ! ratio of current muF2 over reference muF2
+ 1        = QES_over_ref     ! ratio of current QES over reference QES
+#*********************************************************************** 
+# Reweight flags to get scale dependence and PDF uncertainty           *
+# For scale dependence: factor rw_scale_up/down around central scale   *
+# For PDF uncertainty: use LHAPDF with supported set                   *
+#***********************************************************************
+ .true.   = reweight_scale   ! reweight to get scale dependence
+  0.5     = rw_Rscale_down   ! lower bound for ren scale variations
+  2.0     = rw_Rscale_up     ! upper bound for ren scale variations
+  0.5     = rw_Fscale_down   ! lower bound for fact scale variations
+  2.0     = rw_Fscale_up     ! upper bound for fact scale variations
+ .true.  = reweight_PDF     ! reweight to get PDF uncertainty
+ 292001   = PDF_set_min      ! First of the error PDF sets 
+ 292102   = PDF_set_max      ! Last of the error PDF sets
+#***********************************************************************
+# Merging - WARNING! Applies merging only at the hard-event level.     *
+# After showering an MLM-type merging should be applied as well.       *
+# See http://amcatnlo.cern.ch/FxFx_merging.htm for more details.       *
+#***********************************************************************
+ 0        = ickkw            ! 0 no merging, 3 FxFx merging
+#***********************************************************************
+#
+#***********************************************************************
+# BW cutoff (M+/-bwcutoff*Gamma)                                       *
+#***********************************************************************
+ 15  = bwcutoff
+#***********************************************************************
+# Cuts on the jets                                                     *
+# Jet clustering is performed by FastJet.
+# When matching to a parton shower, these generation cuts should be    *
+# considerably softer than the analysis cuts.                          *
+# (more specific cuts can be specified in SubProcesses/cuts.f)         *
+#***********************************************************************
+   1  = jetalgo   ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
+ 0.7  = jetradius ! The radius parameter for the jet algorithm
+   0  = ptj       ! Min jet transverse momentum
+  -1  = etaj      ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
+#***********************************************************************
+# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-)        *
+# (more specific gen cuts can be specified in SubProcesses/cuts.f)     *
+#***********************************************************************
+   0  = ptl     ! Min lepton transverse momentum
+  -1  = etal    ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
+   0  = drll    ! Min distance between opposite sign lepton pairs
+   0  = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
+   0  = mll     ! Min inv. mass of all opposite sign lepton pairs
+  30  = mll_sf  ! Min inv. mass of all opp. sign same-flavor lepton pairs
+#***********************************************************************
+# Photon-isolation cuts, according to hep-ph/9801442                   *
+# When ptgmin=0, all the other parameters are ignored                  *
+#***********************************************************************
+  20  = ptgmin    ! Min photon transverse momentum
+  -1  = etagamma  ! Max photon abs(pseudo-rap)
+ 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)
+#***********************************************************************
+# Maximal PDG code for quark to be considered a jet when applying cuts.*
+# At least all massless quarks of the model should be included here.   *
+#***********************************************************************
+ 4 = maxjetflavor
+#***********************************************************************
diff --git a/Cards/SM_process/st_tch_5f_ckm_NLO/mg5_commands.txt b/Cards/SM_process/st_tch_5f_ckm_NLO/mg5_commands.txt
@@ -0,0 +1,9 @@
+launch
+   0                  # Turn off showering/hadronization
+   set run_tag "test"
+   set use_syst False
+   set nevents 50
+   set iseed 11
+   # set reweight_pdf False
+   # set store_rwgt_info True
+   # reweight
diff --git a/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_customizecards.dat b/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_customizecards.dat
@@ -0,0 +1,3 @@
+#put card customizations here (change top and higgs mass for example)
+set param_card mass 6 172.5
+set param_card yukawa 6 172.5
diff --git a/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_madspin_card.dat b/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_madspin_card.dat
@@ -0,0 +1,40 @@
+#************************************************************
+#*                        MadSpin                           *
+#*                                                          *
+#*    P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer * 
+#*                                                          *
+#*    Part of the MadGraph5_aMC@NLO Framework:              *
+#*    The MadGraph5_aMC@NLO Development Team - Find us at   *
+#*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
+#*                                                          *
+#************************************************************
+#Some options (uncomment to apply)
+#
+#directory for gridpack mode
+set ms_dir ./madspingrid
+
+#initialization parameters
+ set Nevents_for_max_weigth 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
+# 
+#to properly limit the number of concurrent processes for grid running
+set max_running_process 1
+#
+
+define ell+ = e+ mu+ ta+
+define ell- = e- mu- ta-
+# specify the decay for the final state particles
+decay t > w+ b, w+ > ell+ vl
+decay t~ > w- b~, w- > ell- vl~
+decay w+ > ell+ vl
+decay w- > ell- vl~
+
+#define ell+ = e+ mu+ ta+
+#define ell- = e- mu- ta-
+
+# specify the decay for the final state particles
+#decay w+ > ell+ vl
+#decay w- > ell- vl~
+# running the actual code
+launch
diff --git a/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_proc_card.dat b/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_proc_card.dat
@@ -0,0 +1,21 @@
+import model loop_sm-ckm
+#switch to non-diagonal ckm matrix for speed
+#import model loop_sm-no_b_mass
+
+define p = p b b~
+define j = j b b~
+
+generate p p > t j $$ w+ w- [QCD] @0
+add process p p > t~ j $$ w+ w- [QCD] @1
+
+output st_tch_5f_ckm_NLO -nojpeg
+
+launch
+   0                  # Turn off showering/hadronization
+   set run_tag "test"
+   set use_syst False
+   set nevents 50
+   set iseed 11
+   # set reweight_pdf False
+   # set store_rwgt_info True
+   # reweight