diff --git a/.gitignore b/.gitignore index 7478d65..b0ad2c0 100644 --- 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 new file mode 100644 index 0000000..7de0bbe --- /dev/null +++ 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 new file mode 100644 index 0000000..635c0b6 --- /dev/null +++ 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 new file mode 100644 index 0000000..afd237e --- /dev/null +++ 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 new file mode 100644 index 0000000..510d9f8 --- /dev/null +++ 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 new file mode 100644 index 0000000..56264e8 --- /dev/null +++ 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 new file mode 100644 index 0000000..7de0bbe --- /dev/null +++ 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 new file mode 100644 index 0000000..635c0b6 --- /dev/null +++ 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 new file mode 100644 index 0000000..afd237e --- /dev/null +++ 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 new file mode 100644 index 0000000..ec1db6f --- /dev/null +++ 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 diff --git a/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_run_card.dat b/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_ckm_NLO_run_card.dat new file mode 100644 index 0000000..90b1784 --- /dev/null +++ b/Cards/SM_process/st_tch_5f_ckm_NLO/st_tch_5f_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 + 292200 = 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 + 292201 = PDF_set_min ! First of the error PDF sets + 292302 = 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. * +#*********************************************************************** + 5 = maxjetflavor +#*********************************************************************** diff --git a/Cards/mg5_commands.txt b/Cards/mg5_commands.txt new file mode 100644 index 0000000..7de0bbe --- /dev/null +++ b/Cards/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/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LmNuWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat b/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LmNuWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 603c443..0000000 --- a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LmNuWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,31 +0,0 @@ -set group_subprocesses Auto -set ignore_six_quark_processes False -set loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ -define l- = e- mu- -define vl = ve vm -define vl~ = ve~ vm~ -#generate p p > l+ vl w- j j QED=5 QCD=0 -generate p p > l- vl~ w+ j j QED=5 QCD=0, w+ > j j -output test --nojpeg - -launch - set nevents 5000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat b/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 5ea801c..0000000 --- a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,31 +0,0 @@ -set group_subprocesses Auto -set ignore_six_quark_processes False -set loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ -define l- = e- mu- -define vl = ve vm -define vl~ = ve~ vm~ -#generate p p > l+ vl w- j j QED=5 QCD=0 -generate p p > l+ vl w- j j QED=5 QCD=0, w- > j j -output LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10 --nojpeg - -launch - set nevents 5000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WMlepWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat b/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WMlepWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat deleted file mode 100644 index bb69453..0000000 --- a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WMlepWPhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,31 +0,0 @@ -set group_subprocesses Auto -set ignore_six_quark_processes False -set loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ -define l- = e- mu- -define vl = ve vm -define vl~ = ve~ vm~ -#generate p p > l+ vl w- j j QED=5 QCD=0 -generate p p > w+ w- j j QED=4 QCD=0, w+ > j j, w- > l- vl~ -output test --nojpeg - -launch - set nevents 5000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WPlepWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat b/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WPlepWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 063fa8c..0000000 --- a/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/WPlepWMhadJJ_EWK_LO_SM_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,31 +0,0 @@ -set group_subprocesses Auto -set ignore_six_quark_processes False -set loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ -define l- = e- mu- -define vl = ve vm -define vl~ = ve~ vm~ -#generate p p > l+ vl w- j j QED=5 QCD=0 -generate p p > w+ w- j j QED=4 QCD=0, w- > j j, w+ > l+ vl -output test --nojpeg - -launch - set nevents 5000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/Condor_submit.sh b/Older/Condor_submit.sh deleted file mode 100644 index 74e4475..0000000 --- a/Older/Condor_submit.sh +++ /dev/null @@ -1,20 +0,0 @@ -sed -i 's/Ext6/Ext7/g' RunMadGraph_EWK.jdl -sed -i 's/Ext6/Ext7/g' RunMadGraph_EWK.sh -sed -i 's/Ext6/Ext7/g' RunMadGraph_QCD.jdl -sed -i 's/Ext6/Ext7/g' RunMadGraph_QCD.sh -sed -i 's/Ext6/Ext7/g' RunMadGraph_EWKaQCD.jdl -sed -i 's/Ext6/Ext7/g' RunMadGraph_EWKaQCD.sh -sed -i 's/set iseed 29/set iseed 39/g' EWK_cards/*.dat -sed -i 's/set iseed 29/set iseed 39/g' QCD_Cards/*.dat -sed -i 's/set iseed 29/set iseed 39/g' EWKaQCD_cards/*.dat -mkdir condor_logs/Ext7/ -mkdir condor_logs/Ext7/EWK_cards -mkdir condor_logs/Ext7/QCD_Cards -mkdir condor_logs/Ext7/EWKaQCD_cards -eosmkdir /eos/uscms/store/user/rasharma/MonteCarlo_Samples/Ext7 -eosmkdir /eos/uscms/store/user/rasharma/MonteCarlo_Samples/Ext7/EWK -eosmkdir /eos/uscms/store/user/rasharma/MonteCarlo_Samples/Ext7/QCD -eosmkdir /eos/uscms/store/user/rasharma/MonteCarlo_Samples/Ext7/EWKaQCD -condor_submit RunMadGraph_EWK.jdl -condor_submit RunMadGraph_QCD.jdl -condor_submit RunMadGraph_EWKaQCD.jdl diff --git a/Older/EWK_cards/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index a36a235..0000000 --- a/Older/EWK_cards/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=4 QCD=0, w- > j j , z > l+ l- -output WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 778e9f6..0000000 --- a/Older/EWK_cards/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,55 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w- -define w2 = w- -generate p p > w1 w2 j j QED=4 QCD=0, w1 > l- vl~, w2 > j j -output WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 3e92234..0000000 --- a/Older/EWK_cards/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=4 QCD=0, w- > l- vl~, z > j j -output WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index eaeffb4..0000000 --- a/Older/EWK_cards/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=4 QCD=0, w+ > j j , w- > l- vl~ -output WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index a4c7e49..0000000 --- a/Older/EWK_cards/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=4 QCD=0, w+ > j j, z > l+ l- -output WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index c8fc2a6..0000000 --- a/Older/EWK_cards/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=4 QCD=0, w+ > l+ vl, w- > j j -output WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 8849338..0000000 --- a/Older/EWK_cards/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w+ -define w2 = w+ -generate p p > w1 w2 j j QED=4 QCD=0, w1 > l+ vl, w2 > j j -output WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index dce9cb8..0000000 --- a/Older/EWK_cards/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=4 QCD=0, w+ > l+ vl, z > j j -output WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWK_cards/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat b/Older/EWK_cards/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 83141c8..0000000 --- a/Older/EWK_cards/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define z1 = z -define z2 = z -generate p p > z1 z2 j j QED=4 QCD=0 , z1 > l+ l- , z2 > j j -output ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 1e5b954..0000000 --- a/Older/EWKaQCD_cards/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=4 QCD=99, w- > j j , z > l+ l- -output WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index cb651c7..0000000 --- a/Older/EWKaQCD_cards/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,55 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w- -define w2 = w- -generate p p > w1 w2 j j QED=4 QCD=99, w1 > l- vl~, w2 > j j -output WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 8085490..0000000 --- a/Older/EWKaQCD_cards/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=4 QCD=99, w- > l- vl~, z > j j -output WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 6313d3f..0000000 --- a/Older/EWKaQCD_cards/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=4 QCD=99, w+ > j j , w- > l- vl~ -output WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index b6752a3..0000000 --- a/Older/EWKaQCD_cards/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=4 QCD=99, w+ > j j, z > l+ l- -output WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index f33051d..0000000 --- a/Older/EWKaQCD_cards/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=4 QCD=99, w+ > l+ vl, w- > j j -output WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index baad297..0000000 --- a/Older/EWKaQCD_cards/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w+ -define w2 = w+ -generate p p > w1 w2 j j QED=4 QCD=99, w1 > l+ vl, w2 > j j -output WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 8223b46..0000000 --- a/Older/EWKaQCD_cards/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=4 QCD=99, w+ > l+ vl, z > j j -output WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/EWKaQCD_cards/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat b/Older/EWKaQCD_cards/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 23cfd80..0000000 --- a/Older/EWKaQCD_cards/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define z1 = z -define z2 = z -generate p p > z1 z2 j j QED=4 QCD=99 , z1 > l+ l- , z2 > j j -output ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/List_all_LHEfile.py b/Older/List_all_LHEfile.py deleted file mode 100644 index 7bc0ef4..0000000 --- a/Older/List_all_LHEfile.py +++ /dev/null @@ -1,18 +0,0 @@ -import sys -import os - -source = "/eos/uscms/store/user/rasharma/MonteCarlo_Samples/" - -print source - -neglect1="VBFNLO_outputs" -neglect2="tWch" - -for root, dirs, filenames in os.walk(source): - for f in filenames: - filepath = root + os.sep + f - if filepath.find(neglect1) == -1: - if filepath.find(neglect2) == -1: - if filepath.endswith(".lhe"): - #print root - print filepath diff --git a/Older/List_all_LHEfile_For_HADD.py b/Older/List_all_LHEfile_For_HADD.py deleted file mode 100644 index 1da2d3c..0000000 --- a/Older/List_all_LHEfile_For_HADD.py +++ /dev/null @@ -1,34 +0,0 @@ -import sys -import os - -source = "/eos/uscms/store/user/rasharma/MonteCarlo_Samples/" - -print "\n=======================================" -print "This is the name of the script: ", sys.argv[0] -print "Number of arguments: ", len(sys.argv) -print "The arguments are: " , str(sys.argv) - -print source - -neglect1="VBFNLO_outputs" -neglect2="tWch" - -#include="WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10" -include=str(sys.argv[1]) -#print "Include = ",include - -#print "hadd "+include+".root ", -command="hadd -f "+include+".root " - -for root, dirs, filenames in os.walk(source): - for f in filenames: - filepath = root + os.sep + f - if filepath.find(neglect1) == -1: - if filepath.find(neglect2) == -1: - if filepath.find(include) != -1: - if filepath.endswith(".root"): - #print root - print filepath - command+=" "+filepath+" " -print command -os.system(command) diff --git a/Older/List_all_LHEfile_For_HADD_VBFNLO.py b/Older/List_all_LHEfile_For_HADD_VBFNLO.py deleted file mode 100644 index e003464..0000000 --- a/Older/List_all_LHEfile_For_HADD_VBFNLO.py +++ /dev/null @@ -1,34 +0,0 @@ -import sys -import os - -source = "/eos/uscms/store/user/rasharma/MonteCarlo_Samples/" - -print "\n=======================================" -print "This is the name of the script: ", sys.argv[0] -print "Number of arguments: ", len(sys.argv) -print "The arguments are: " , str(sys.argv) - -print source - -neglect1="Events" -neglect2="tWch" - -#include="WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10" -include=str(sys.argv[1]) -#print "Include = ",include - -#print "hadd "+include+".root ", -command="hadd -f "+include+"_VBFNLO.root " - -for root, dirs, filenames in os.walk(source): - for f in filenames: - filepath = root + os.sep + f - if filepath.find(neglect1) == -1: - if filepath.find(neglect2) == -1: - if filepath.find(include) != -1: - if filepath.endswith("event.root"): - #print root - print filepath - command+=" "+filepath+" " -print command -os.system(command) diff --git a/Older/List_for_LHEanalyzer.py b/Older/List_for_LHEanalyzer.py deleted file mode 100644 index 9424f7d..0000000 --- a/Older/List_for_LHEanalyzer.py +++ /dev/null @@ -1,24 +0,0 @@ -import sys -import os - -source = "/eos/uscms/store/user/rasharma/MonteCarlo_Samples/" - -print source - -#neglect1="VBFNLO_outputs" -neglect1="Events" -neglect2="tWch" - -#include="202_WPWMhadjj" -include="201_WPhadWMjj" - -for root, dirs, filenames in os.walk(source): - for f in filenames: - filepath = root + os.sep + f - if filepath.find(neglect1) == -1: - if filepath.find(neglect2) == -1: - if filepath.find(include) != -1: - if filepath.endswith(".lhe"): - #print root - #print filepath - print "./LHEanalyzer "+filepath+" "+filepath.replace("lhe","root") diff --git a/Older/QCD_Cards/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index d0a397f..0000000 --- a/Older/QCD_Cards/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=2 QCD=99, w- > j j , z > l+ l- -output WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index b5be897..0000000 --- a/Older/QCD_Cards/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,55 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w- -define w2 = w- -generate p p > w1 w2 j j QED=2 QCD=99, w1 > l- vl~, w2 > j j -output WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 52feaa5..0000000 --- a/Older/QCD_Cards/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w- z j j QED=2 QCD=99, w- > l- vl~, z > j j -output WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 54656a6..0000000 --- a/Older/QCD_Cards/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=2 QCD=99, w+ > j j , w- > l- vl~ -output WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 01ec9dc..0000000 --- a/Older/QCD_Cards/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=2 QCD=99, w+ > j j, z > l+ l- -output WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 09dab92..0000000 --- a/Older/QCD_Cards/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ w- j j QED=2 QCD=99, w+ > l+ vl, w- > j j -output WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index bddf2b6..0000000 --- a/Older/QCD_Cards/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define w1 = w+ -define w2 = w+ -generate p p > w1 w2 j j QED=2 QCD=99, w1 > l+ vl, w2 > j j -output WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index b3058cd..0000000 --- a/Older/QCD_Cards/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,54 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -generate p p > w+ z j j QED=2 QCD=99, w+ > l+ vl, z > j j -output WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/QCD_Cards/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat b/Older/QCD_Cards/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat deleted file mode 100644 index 602d86a..0000000 --- a/Older/QCD_Cards/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat +++ /dev/null @@ -1,56 +0,0 @@ -#************************************************************ -#* MadGraph5_aMC@NLO * -#* * -#* * * * -#* * * * * * -#* * * * * 5 * * * * * -#* * * * * * -#* * * * -#* * -#* * -#* VERSION 2.2.1 2014-09-25 * -#* * -#* 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 loop_optimized_output True -set complex_mass_scheme False -import model sm-no_b_mass -define p = g u c d s b u~ c~ d~ s~ b~ -define j = p -define l+ = e+ mu+ ta+ -define l- = e- mu- ta- -define vl = ve vm vt -define vl~ = ve~ vm~ vt~ -define z1 = z -define z2 = z -generate p p > z1 z2 j j QED=2 QCD=99 , z1 > l+ l- , z2 > j j -output ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10 --nojpeg - -launch - - set nevents 50000 - set ptj 10 - set ptb 10 - set pta 0 - set ptl 0 - set etaj 6.5 - set etab 6.5 - set etaa -1 - set etal -1 - set drbb 0.4 - set drbj 0.4 - set mmjj 100 - set mmbb 100 - set maxjetflavor 5 - set iseed 39 diff --git a/Older/README.md b/Older/README.md deleted file mode 100644 index aa51194..0000000 --- a/Older/README.md +++ /dev/null @@ -1,49 +0,0 @@ -# MadGraph_With_Condor - -* This repo contains script for submitting the madgraph jobs. -* Just we have to prepare the proc card and give its name in the \*.jdl file -* Also, if you need to generate a lot of events then you can use the **Queue** feature of condor job. - * The script sets "iSeed" for every run so the seed should be unique and hence you are not going to replicate it. - -# How to generate events - - cmsrel CMSSW_8_0_11 - cd CMSSW_8_0_11/src - cmsenv - wget https://launchpad.net/mg5amcnlo/2.0/2.5.x/+download/MG5_aMC_v2.5.5.tar.gz - -* Modify files RunMadGraph_condor.sh and RunMadGraph_condor.jdl, then submit it using - - condor_submit RunMadGraph_condor.jdl - -* Note that all input \*_proc_card.dat's prefix should be exactly same as name of output file that madgraph will create. - -* Before submitting the condor you have to - - voms-proxy-init --voms cms --valid 168:00 - -Extract cross-section - - grep "Total:\|Cross-section" - -Check the random seed from stdout file - - grep "Using random number seed offset" */*.stdout - -Change the seed - - sed -i 's/set iseed 15/set iseed 17/g' EWK_cards/*.dat - sed -i 's/set iseed 15/set iseed 17/g' QCD_Cards/*.dat - sed -i 's/set iseed 15/set iseed 17/g' EWKaQCD_cards/*.dat - -change directory name - - sed -i 's/Ext3/Ext4/g' RunMadGraph_EWK.jdl - sed -i 's/Ext3/Ext4/g' RunMadGraph_EWK.sh - sed -i 's/Ext3/Ext4/g' RunMadGraph_QCD.jdl - sed -i 's/Ext3/Ext4/g' RunMadGraph_QCD.sh - sed -i 's/Ext3/Ext4/g' RunMadGraph_EWKaQCD.jdl - sed -i 's/Ext3/Ext4/g' RunMadGraph_EWKaQCD.sho - -condor_q -submitter rasharma | grep "running,\|idle,\|.fnal.gov" | awk -F "," '{print $3,$4}' -condor_q -submitter rasharma | grep "running,\|idle,\|.fnal.gov" diff --git a/Older/RunMadGraph_EWK.jdl b/Older/RunMadGraph_EWK.jdl deleted file mode 100644 index e0ef54a..0000000 --- a/Older/RunMadGraph_EWK.jdl +++ /dev/null @@ -1,55 +0,0 @@ -Executable = RunMadGraph_EWK.sh -Universe = vanilla -Requirements =FileSystemDomain=="fnal.gov" && Arch=="X86_64" -Notification = ERROR -Should_Transfer_Files = YES -WhenToTransferOutput = ON_EXIT -request_memory = 3000 -request_cpus = 8 -transfer_input_files = RunMadGraph_EWK.sh,MG5_aMC_v2.5.4.tar.gz,EWK_cards/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat,EWK_cards/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -x509userproxy = $ENV(X509_USER_PROXY) -Output = condor_logs/July262017/EWK/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WMhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WMlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WMlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WPhadWMlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WPhadZlepJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WPlepWMhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WPlepWPhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = WPlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWK/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWK/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWK/ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.log -Arguments = ZlepZhadJJ_SM_LO_EWK_mjj100_pTj10_proc_card.dat -Queue 50 diff --git a/Older/RunMadGraph_EWK.sh b/Older/RunMadGraph_EWK.sh deleted file mode 100755 index a11037c..0000000 --- a/Older/RunMadGraph_EWK.sh +++ /dev/null @@ -1,52 +0,0 @@ -#!/bin/bash -date -cd /uscms_data/d3/rasharma/aQGC_analysis/CMS_FulllSimulation_April2017/LHE_GEN/CMSSW_8_0_21/src -cmsenv -echo $PWD -eval `scram runtime -sh` -cd - -echo $PWD -a=$PWD -echo $a - -RandNum=`python -c "import random; print random.randint(0,10000000)"` - -echo "Random Seed = ${RandNum}" - -tar -xf MG5_aMC_v2.5.4.tar.gz -cd MG5_aMC_v2_5_4 -cp ../*.dat . -sed "s/39/${RandNum}/g" $* -./bin/mg5_aMC $* - -echo -e "\n\n=====================================================================\n\n" -test=`echo $* | sed 's/_proc_card.dat//g'` -echo "=====> Variable test = "$test -echo "List all root files = " -ls -echo "List all files" -ls $test/Events -ls $test/Events/* - -# copy output to eos -OUTDIR=root://cmseos.fnal.gov//store/user/rasharma/MonteCarlo_Samples/MadGraph/EWK - -eos root://cmseos.fnal.gov mkdir $OUTDIR/$RandNum - -OUTDIR=$OUTDIR/$RandNum -# -# -echo "xrdcp output for condor" -for FILE in $test/Events/*/* -do - echo "xrdcp -f ${FILE} ${OUTDIR}/${FILE}" - xrdcp -f ${FILE} ${OUTDIR}/${FILE} 2>&1 - XRDEXIT=$? - if [[ $XRDEXIT -ne 0 ]]; then - rm *.root - echo "exit code $XRDEXIT, failure in xrdcp" - exit $XRDEXIT - fi - rm ${FILE} -done -date diff --git a/Older/RunMadGraph_EWKaQCD.jdl b/Older/RunMadGraph_EWKaQCD.jdl deleted file mode 100644 index 2a38b40..0000000 --- a/Older/RunMadGraph_EWKaQCD.jdl +++ /dev/null @@ -1,55 +0,0 @@ -Executable = RunMadGraph_EWKaQCD.sh -Universe = vanilla -Requirements =FileSystemDomain=="fnal.gov" && Arch=="X86_64" -Notification = ERROR -Should_Transfer_Files = YES -WhenToTransferOutput = ON_EXIT -request_memory = 3000 -request_cpus = 8 -transfer_input_files = RunMadGraph_EWKaQCD.sh,MG5_aMC_v2.5.4.tar.gz,EWKaQCD_cards/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat,EWKaQCD_cards/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -x509userproxy = $ENV(X509_USER_PROXY) -Output = condor_logs/July262017/EWKaQCD/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WMhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WMlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WMlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WPhadWMlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WPhadZlepJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWMhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWPhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = WPlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/EWKaQCD/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/EWKaQCD/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/EWKaQCD/ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.log -Arguments = ZlepZhadJJ_SM_LO_EWKaQCD_mjj100_pTj10_proc_card.dat -Queue 50 diff --git a/Older/RunMadGraph_EWKaQCD.sh b/Older/RunMadGraph_EWKaQCD.sh deleted file mode 100755 index 5b0747a..0000000 --- a/Older/RunMadGraph_EWKaQCD.sh +++ /dev/null @@ -1,52 +0,0 @@ -#!/bin/bash -date -cd /uscms_data/d3/rasharma/aQGC_analysis/CMS_FulllSimulation_April2017/LHE_GEN/CMSSW_8_0_21/src -cmsenv -echo $PWD -eval `scram runtime -sh` -cd - -echo $PWD -a=$PWD -echo $a - -RandNum=`python -c "import random; print random.randint(0,10000000)"` - -echo "Random Seed = ${RandNum}" - -tar -xf MG5_aMC_v2.5.4.tar.gz -cd MG5_aMC_v2_5_4 -cp ../*.dat . -sed "s/39/${RandNum}/g" $* -./bin/mg5_aMC $* - -echo -e "\n\n=====================================================================\n\n" -test=`echo $* | sed 's/_proc_card.dat//g'` -echo "=====> Variable test = "$test -echo "List all root files = " -ls -echo "List all files" -ls $test/Events -ls $test/Events/* - -# copy output to eos -OUTDIR=root://cmseos.fnal.gov//store/user/rasharma/MonteCarlo_Samples/MadGraph/EWKaQCD - -eos root://cmseos.fnal.gov mkdir $OUTDIR/$RandNum - -OUTDIR=$OUTDIR/$RandNum -# -# -echo "xrdcp output for condor" -for FILE in $test/Events/*/* -do - echo "xrdcp -f ${FILE} ${OUTDIR}/${FILE}" - xrdcp -f ${FILE} ${OUTDIR}/${FILE} 2>&1 - XRDEXIT=$? - if [[ $XRDEXIT -ne 0 ]]; then - rm *.root - echo "exit code $XRDEXIT, failure in xrdcp" - exit $XRDEXIT - fi - rm ${FILE} -done -date diff --git a/Older/RunMadGraph_QCD.jdl b/Older/RunMadGraph_QCD.jdl deleted file mode 100644 index 0338c92..0000000 --- a/Older/RunMadGraph_QCD.jdl +++ /dev/null @@ -1,55 +0,0 @@ -Executable = RunMadGraph_QCD.sh -Universe = vanilla -Requirements =FileSystemDomain=="fnal.gov" && Arch=="X86_64" -Notification = ERROR -Should_Transfer_Files = YES -WhenToTransferOutput = ON_EXIT -request_memory = 3000 -request_cpus = 8 -transfer_input_files = RunMadGraph_QCD.sh,MG5_aMC_v2.5.4.tar.gz,QCD_Cards/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat,QCD_Cards/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -x509userproxy = $ENV(X509_USER_PROXY) -Output = condor_logs/July262017/QCD/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPhadWMlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPhadZlepJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWMhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWPhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 -Output = condor_logs/July262017/QCD/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/July262017/QCD/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/July262017/QCD/ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = ZlepZhadJJ_SM_LO_QCD_mjj100_pTj10_proc_card.dat -Queue 50 diff --git a/Older/RunMadGraph_QCD.sh b/Older/RunMadGraph_QCD.sh deleted file mode 100755 index 147c19a..0000000 --- a/Older/RunMadGraph_QCD.sh +++ /dev/null @@ -1,52 +0,0 @@ -#!/bin/bash -date -cd /uscms_data/d3/rasharma/aQGC_analysis/CMS_FulllSimulation_April2017/LHE_GEN/CMSSW_8_0_21/src -cmsenv -echo $PWD -eval `scram runtime -sh` -cd - -echo $PWD -a=$PWD -echo $a - -RandNum=`python -c "import random; print random.randint(0,10000000)"` - -echo "Random Seed = ${RandNum}" - -tar -xf MG5_aMC_v2.5.4.tar.gz -cd MG5_aMC_v2_5_4 -cp ../*.dat . -sed "s/39/${RandNum}/g" $* -./bin/mg5_aMC $* - -echo -e "\n\n=====================================================================\n\n" -test=`echo $* | sed 's/_proc_card.dat//g'` -echo "=====> Variable test = "$test -echo "List all root files = " -ls -echo "List all files" -ls $test/Events -ls $test/Events/* - -# copy output to eos -OUTDIR=root://cmseos.fnal.gov//store/user/rasharma/MonteCarlo_Samples/MadGraph/QCD - -eos root://cmseos.fnal.gov mkdir $OUTDIR/$RandNum - -OUTDIR=$OUTDIR/$RandNum -# -# -echo "xrdcp output for condor" -for FILE in $test/Events/*/* -do - echo "xrdcp -f ${FILE} ${OUTDIR}/${FILE}" - xrdcp -f ${FILE} ${OUTDIR}/${FILE} 2>&1 - XRDEXIT=$? - if [[ $XRDEXIT -ne 0 ]]; then - rm *.root - echo "exit code $XRDEXIT, failure in xrdcp" - exit $XRDEXIT - fi - rm ${FILE} -done -date diff --git a/Older/RunMadGraph_condor.jdl b/Older/RunMadGraph_condor.jdl deleted file mode 100644 index fcd6755..0000000 --- a/Older/RunMadGraph_condor.jdl +++ /dev/null @@ -1,53 +0,0 @@ -Executable = RunGENSIM_condor.sh -Universe = vanilla -Requirements =FileSystemDomain=="fnal.gov" && Arch=="X86_64" -Notification = ERROR -Should_Transfer_Files = YES -WhenToTransferOutput = ON_EXIT -transfer_input_files = RunGENSIM_condor.sh,MG5_aMC_v2.5.4.tar.gz,WMhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WMlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WMlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WPhadWMlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WPhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WPlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WPlepWPhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,WPlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat,ZlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -x509userproxy = $ENV(X509_USER_PROXY) -Output = condor_logs/WMhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WMhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WMhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WMlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WMlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WMlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WMlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WMlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WMlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WMlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WPhadWMlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WPhadWMlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WPhadWMlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPhadWMlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WPhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WPhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WPhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPhadZlepJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WPlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WPlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WPlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWMhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WPlepWPhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WPlepWPhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WPlepWPhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepWPhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/WPlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/WPlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/WPlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = WPlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue -Output = condor_logs/ZlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Error = condor_logs/ZlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.stdout -Log = condor_logs/ZlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.log -Arguments = ZlepZhadJJ_EWK_LO_QCD_mjj100_pTj10_proc_card.dat -Queue diff --git a/Older/RunMadGraph_condor.sh b/Older/RunMadGraph_condor.sh deleted file mode 100755 index ce61c42..0000000 --- a/Older/RunMadGraph_condor.sh +++ /dev/null @@ -1,43 +0,0 @@ -#!/bin/bash -date -cd /uscms_data/d3/rasharma/aQGC_analysis/CMS_FulllSimulation_April2017/LHE_GEN/CMSSW_8_0_21/src -cmsenv -echo $PWD -eval `scram runtime -sh` -cd - -echo $PWD -a=$PWD -echo $a - -tar -xf MG5_aMC_v2.5.4.tar.gz -cd MG5_aMC_v2_5_4 -cp ../*.dat . -./bin/mg5_aMC $* - -echo -e "\n\n=====================================================================\n\n" -test=`echo $* | sed 's/_proc_card.dat//g'` -echo "=====> Variable test = "$test -echo "List all root files = " -ls -echo "List all files" -ls $test/Events -ls $test/Events/* - -# copy output to eos -OUTDIR=root://cmseos.fnal.gov//store/user/rasharma/MadGraph_OutPuts/ -# -# -echo "xrdcp output for condor" -for FILE in $test/Events/*/* -do - echo "xrdcp -f ${FILE} ${OUTDIR}/${FILE}" - xrdcp -f ${FILE} ${OUTDIR}/${FILE} 2>&1 - XRDEXIT=$? - if [[ $XRDEXIT -ne 0 ]]; then - rm *.root - echo "exit code $XRDEXIT, failure in xrdcp" - exit $XRDEXIT - fi - rm ${FILE} -done -date diff --git a/Older/Unzip_all_LHEfile.py b/Older/Unzip_all_LHEfile.py deleted file mode 100644 index be944f9..0000000 --- a/Older/Unzip_all_LHEfile.py +++ /dev/null @@ -1,19 +0,0 @@ -import sys -import os - -source = "/eos/uscms/store/user/rasharma/MonteCarlo_Samples/" - -print source - -for root, dirs, filenames in os.walk(source): - for f in filenames: - filepath = root + os.sep + f - if filepath.endswith(".gz"): - print root - print filepath - os.system("mv "+filepath+" /tmp/rasharma/") - print("listing all file in $PWD:") - os.system("ls") - os.system("gunzip /tmp/rasharma/"+f) - os.system("mv /tmp/rasharma/unweighted_events.lhe "+root) - print("\n\n#####################################\n\n") diff --git a/Older/condor_logs/.gitignore b/Older/condor_logs/.gitignore deleted file mode 100644 index 29cd89b..0000000 --- a/Older/condor_logs/.gitignore +++ /dev/null @@ -1,2 +0,0 @@ -*.stdout -*.log diff --git a/Older/general_condor_functions.py b/Older/general_condor_functions.py deleted file mode 100644 index a58b7c7..0000000 --- a/Older/general_condor_functions.py +++ /dev/null @@ -1,129 +0,0 @@ -#! /usr/bin/env python - -import os -import argparse -import datetime - -def getbasic_parser(): - parser = argparse.ArgumentParser(description='User inputs') - - parser.add_argument('-i', '--inputpath', - default='/store/user/lnujj/VVjj_aQGC/Gridpacks/', - help='tar file path' - ) - parser.add_argument('-f', '--tarfile', - required=True, - help='input tar file name' - ) - parser.add_argument('-o', '--outputpath', - default='/store/user/lnujj/VVjj_aQGC/LHEFiles/', - help='outputpath path of LHEFiles' - ) - parser.add_argument('-od', '--outputdir', - default='LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10', - help='Name of output directory' - ) - parser.add_argument('-t', '--testrun', - default=True, - help='is it a test run for check? True/False' - ) - parser.add_argument('-c', '--cmsswversion', - default='CMSSW_9_3_8', - help='cmssw version to be used' - ) - parser.add_argument('-j', '--jdlfilename', - default='run_mg5_condor', - help='name of jdl file and its sh file' - ) - return parser - -def create_output_directory(args): - """ - This function will create two directories. - - First directory will be created in PWD for storing the log files - - Another directory will be created in eos for storing the outputs - """ - if args.testrun: - output_folder = args.outputpath+datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M')+"_TEST/" - output_log_path = "output_logs/" + datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') + "_TEST" - print "Name of output dir: ", output_folder - else: - os.system('xrdfs root://cmseos.fnal.gov/ mkdir ' + args.outputpath + args.outputdir) - output_folder = args.outputpath+args.outputdir+"/"+datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') - output_log_path = "output_logs/" + args.outputdir+"/"+ datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') - print "Name of output directory : ", output_folder - print "Name of output log directory: ",output_log_path - - # create directory in eos for output files - os.system('xrdfs root://cmseos.fnal.gov/ mkdir ' + output_folder) - # create directory in pwd for log files - os.system('mkdir -p ' + output_log_path) - return output_folder, output_log_path - -def exclude_function(filename): - """ - Function that returns True/False if the filename ends with - *.log or *.stdout - Input arguments: - filename Name of file - """ - return bool(filename.endswith('.log') or filename.endswith('.stdout')) - -def make_tarfile(output_filename, source_dir): - """ - Function to create the tar file. - Input arguments: - output_filename Name of output file - source_dir path of source directory to make tar file - """ - import tarfile - with tarfile.open(output_filename, "w:gz") as tar: - tar.add(source_dir, arcname=os.path.basename(source_dir)) - -def create_jdl_file_for_condor(args, inputlist, output_log_path): - outjdl_file = open(args.jdlfilename+".jdl", "w") - outjdl_file.write("Executable = "+args.jdlfilename+".sh\n") - outjdl_file.write("Universe = vanilla\n") - #outjdl_file.write("Requirements =FileSystemDomain==\"fnal.gov\" && Arch==\"X86_64\"") - outjdl_file.write("Notification = ERROR\n") - outjdl_file.write("Should_Transfer_Files = YES\n") - outjdl_file.write("WhenToTransferOutput = ON_EXIT\n") - #outjdl_file.write("include : list-infiles.sh |\n") - if inputlist != "": - outjdl_file.write("Transfer_Input_Files = "+inputlist+"\n") - outjdl_file.write("x509userproxy = $ENV(X509_USER_PROXY)\n") - - outjdl_file.write("Output = "+output_log_path+"/"+args.outputdir+".stdout\n") - outjdl_file.write("Error = "+output_log_path+"/"+args.outputdir+".stdout\n") - outjdl_file.write("Log = "+output_log_path+"/"+args.outputdir+".log\n") - # outjdl_file.write("Arguments = -n "+args.outputdir+" + \n") - outjdl_file.write("Queue\n") - -def create_sh_file_for_condor(args, command, output_folder): - outscript = open(args.jdlfilename+".sh", "w") - outscript.write('#!/bin/bash') - outscript.write("\n"+'echo "Starting job on " `date`') - outscript.write("\n"+'echo "Running on: `uname -a`"') - outscript.write("\n"+'echo "System software: `cat /etc/redhat-release`"') - outscript.write("\n"+'source /cvmfs/cms.cern.ch/cmsset_default.sh') - outscript.write("\n"+'export SCRAM_ARCH=slc6_amd64_gcc530') - outscript.write("\n"+'eval `scramv1 project CMSSW '+args.cmsswversion+'`') - outscript.write("\n"+'cd '+ args.cmsswversion + '/src/') - outscript.write("\n"+'eval `scram runtime -sh`') - outscript.write("\n"+'xrdcp -s root://cmseos.fnal.gov/'+args.inputpath+args.tarfile+' .') - outscript.write("\n"+'tar -xf '+args.tarfile) - outscript.write("\n"+'echo "====> List files : " ') - outscript.write("\n"+'ls -alh') - outscript.write("\n"+command) - outscript.write("\n"+'echo "====> List files : " ') - outscript.write("\n"+'ls -alh') - outscript.write("\n"+'echo "====> List only LHE files : " ') - outscript.write("\n"+'ls *.lhe') - outscript.write("\n"+'echo "====> copying *.lhe file to stores area..." ') - outscript.write("\n"+'xrdcp -f *.lhe root://cmseos.fnal.gov/' + output_folder) - outscript.write("\n"+'rm *.lhe') - outscript.write("\n"+'cd ${_CONDOR_SCRATCH_DIR}') - outscript.write("\n"+'rm -rf ' + args.cmsswversion) - outscript.write("\n") - outscript.close() - os.system("chmod 777 "+args.jdlfilename+".sh") diff --git a/Older/generate_lhe_condor.py b/Older/generate_lhe_condor.py deleted file mode 100644 index c3912f5..0000000 --- a/Older/generate_lhe_condor.py +++ /dev/null @@ -1,51 +0,0 @@ -#! /usr/bin/env python - -from general_condor_functions import getbasic_parser -from general_condor_functions import create_output_directory -from general_condor_functions import create_jdl_file_for_condor -from general_condor_functions import create_sh_file_for_condor - -def getargs(): - parser = getbasic_parser() - parser.add_argument('-n', '--nevents', - default=1100, - help='Total number of events to generate.' - ) - parser.add_argument('-cpu', '--ncpu', - default=8, - help='number of cpu to run' - ) - parser.add_argument('-r', '--randomnumber', - default=11, - help='random seed' - ) - return parser.parse_args() - - -def main(): - # get input arguments - args = getargs() - - # list of input files to be added in jdl file - inputlist = args.jdlfilename+".sh" - - # command to run - command = './runcmsgrid.sh '+str(args.nevents)+' '+str(args.randomnumber)+' '+str(args.ncpu) - - # Get Output directory name - output_folder, output_log_path = create_output_directory(args) - - # create the jdl file for condor - create_jdl_file_for_condor(args, inputlist, output_log_path) - - # create the sh file for condor - create_sh_file_for_condor(args, command, output_folder) - - # running help - print "===> Set Proxy Using:" - print "\tvoms-proxy-init --voms cms --valid 168:00" - print "\"condor_submit "+args.jdlfilename+".jdl\" to submit" - - -if __name__ == "__main__": - main() diff --git a/README.md b/README.md index f84f015..8a63ed5 100644 --- a/README.md +++ b/README.md @@ -1,89 +1,23 @@ -# Run Standalone Madgraph Using Condor +# Introduction +- It is a simple script which takes the input cards and run the MadGraph5_aMC@NLO. It also helps to run the MadGraph5_aMC@NLO on condor. +- We just need to provide the MG5 tarball url, cards and output directory. It will take care of the rest. +- For condor jobs, once the jobs are done, it will transfer the full working directory to the output directory. The reason for copying full directory is because of feynman diagrams and the web page which is generated by MadGraph5_aMC@NLO. -The main script is `generate_lhe_condor.py`. Its help argument goes like: +# Setup -```bash -python generate_lhe_condor.py -h -usage: generate_lhe_condor.py [-h] [-i INPUTPATH] -f TARFILE [-o OUTPUTPATH] - [-od OUTPUTDIR] [-t TESTRUN] [-c CMSSWVERSION] - [-j JDLFILENAME] [-m MEMORY] [-n NEVENTS] - [-cpu NCPU] [-r RANDOMNUMBER] +1. Setup environment + ```bash + source /cvmfs/sft.cern.ch/lcg/views/dev4cuda/latest/x86_64-el9-gcc11-opt/setup.sh + ``` -User inputs +2. Run it using command: + 2.1 Local -optional arguments: - -h, --help show this help message and exit - -i INPUTPATH, --inputpath INPUTPATH - tar file path - -f TARFILE, --tarfile TARFILE - input tar file name - -o OUTPUTPATH, --outputpath OUTPUTPATH - outputpath path of LHEFiles - -od OUTPUTDIR, --outputdir OUTPUTDIR - Name of output directory - -t TESTRUN, --testrun TESTRUN - is it a test run for check? True/False - -c CMSSWVERSION, --cmsswversion CMSSWVERSION - cmssw version to be used - -j JDLFILENAME, --jdlfilename JDLFILENAME - name of jdl file and its sh file - -m MEMORY, --memory MEMORY - memory for condor jobs - -n NEVENTS, --nevents NEVENTS - Total number of events to generate. - -cpu NCPU, --ncpu NCPU - number of cpu to run - -r RANDOMNUMBER, --randomnumber RANDOMNUMBER - random seed -``` + ```bash + python mainScript.py --cards Cards/SM_process/st_tch_5f_ckm_NLO --local + ``` + 2.2 Using condor -# How to use - -1. STEP: 1: - ```bash - python generate_lhe_condor.py -f - ``` - This will create the `jdl` and `sh` file for the condor submission. - -2. STEP: 2: Submit the condor job using the file created by STEP:1. - - -# Few points noted/learned - -1. If the seed value given to madgraph is too long then it throws error like: - ```bash - set maxjetflavor 5 - INFO: modify parameter maxjetflavor of the run_card.dat to 5 - set iseed 6132039411 - INFO: modify parameter iseed of the run_card.dat to 6132039411 - Do you want to edit a card (press enter to bypass editing)? - /------------------------------------------------------------\ - | 1. param : param_card.dat | - | 2. run : run_card.dat | - \------------------------------------------------------------/ - you can also - - enter the path to a valid card or banner. - - use the 'set' command to modify a parameter directly. - The set option works only for param_card and run_card. - Type 'help set' for more information on this command. - - call an external program (ASperGE/MadWidth/...). - Type 'help' for the list of available command - [0, done, 1, param, 2, run, enter path] - The answer to the previous question is not set in your input file - Use 0 value - INFO: Update the dependent parameter of the param_card.dat - Generating 5000 events with run name run_01 - survey run_01 - INFO: compile directory - Not able to open file /storage/local/data1/condor/execute/dir_41567/CMSSW_9_3_8/src/MG5_aMC_v2_6_5/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt - compile Source Directory - Using random number seed offset = 6132039411 - Error detected in "generate_events run_01" - write debug file /storage/local/data1/condor/execute/dir_41567/CMSSW_9_3_8/src/MG5_aMC_v2_6_5/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10/run_01_tag_1_debug.log - If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo - MadGraph5Error : Random seed too large 6132039414 > 30081*30081 - INFO: - quit - INFO: - quit - ``` + ```bash + python mainScript.py --url https://launchpad.net/mg5amcnlo/3.0/3.5.x/+download/MG5_aMC_v3.5.4.tar.gz --cards Cards/SM_process/ st_tch_5f_ckm_NLO --output /eos/user/r/rasharma/www/MadGraph_Study/HwH + ``` diff --git a/general_condor_functions.py b/general_condor_functions.py deleted file mode 100644 index de09a67..0000000 --- a/general_condor_functions.py +++ /dev/null @@ -1,163 +0,0 @@ -#! /usr/bin/env python - -import os -import argparse -import datetime -import random - -def getbasic_parser(): - parser = argparse.ArgumentParser(description='User inputs') - - parser.add_argument('-i', '--inputpath', - default='/store/user/lnujj/VVjj_aQGC/Gridpacks/', - help='tar file path' - ) - parser.add_argument('-il', '--InProcCardPath', - default='/uscms/home/rasharma/nobackup/aQGC_analysis/Run2/MadGraph_With_Condor/LpNuWMhadJJ_EWK_LO_SM_mjj100_pTj10', - help='path of input LHE file' - ) - parser.add_argument('-n', '--nevents', - default=5000, # if number of jobs are more than 5k then MG5 splits them - help='Total number of events to generate.' - ) - parser.add_argument('-f', '--tarfile', - required=True, - help='input tar file name' - ) - parser.add_argument('-o', '--outputpath', - default='/store/user/lnujj/VVjj_aQGC/LHEFiles/', - help='outputpath path of LHEFiles' - ) - parser.add_argument('-od', '--outputdir', - default='test', - help='Name of output directory' - ) - parser.add_argument('-t', '--testrun', - default=False, - help='is it a test run for check? True/False' - ) - parser.add_argument('-c', '--cmsswversion', - default='CMSSW_9_3_8', - help='cmssw version to be used' - ) - parser.add_argument('-j', '--jdlfilename', - default='run_mg5_condor', - help='name of jdl file and its sh file' - ) - parser.add_argument('-m', '--memory', - default=5000, - help='memory for condor jobs' - ) - parser.add_argument('-tj', '--totaljobs', - default=1, - help='number of condor jobs copy' - ) - return parser - -def create_output_directory(args): - """ - This function will create two directories. - - First directory will be created in PWD for storing the log files - - Another directory will be created in eos for storing the outputs - """ - if args.testrun: - output_folder = args.outputpath+datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M')+"_TEST/" - output_log_path = "output_logs/" + datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') + "_TEST" - print "Name of output dir: ", output_folder - else: - os.system('xrdfs root://cmseos.fnal.gov/ mkdir ' + args.outputpath + args.outputdir) - output_folder = args.outputpath+args.outputdir+"/"+datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') - output_log_path = "output_logs/" + args.outputdir+"/"+ datetime.datetime.now().strftime('%Y_%m_%d_%Hh%M') - print "Name of output directory : ", output_folder - print "Name of output log directory: ",output_log_path - - # create directory in eos for output files - os.system('xrdfs root://cmseos.fnal.gov/ mkdir ' + output_folder) - # create directory in pwd for log files - os.system('mkdir -p ' + output_log_path) - return output_folder, output_log_path - -def exclude_function(filename): - """ - Function that returns True/False if the filename ends with - *.log or *.stdout - Input arguments: - filename Name of file - """ - return bool(filename.endswith('.log') or filename.endswith('.stdout')) - -def make_tarfile(output_filename, source_dir): - """ - Function to create the tar file. - Input arguments: - output_filename Name of output file - source_dir path of source directory to make tar file - """ - import tarfile - with tarfile.open(output_filename, "w:gz") as tar: - tar.add(source_dir, arcname=os.path.basename(source_dir)) - -def create_jdl_file_for_condor(args, inputlist, output_log_path): - outjdl_file = open(args.jdlfilename+".jdl", "w") - outjdl_file.write("Executable = "+args.jdlfilename+".sh\n") - outjdl_file.write("Universe = vanilla\n") - #outjdl_file.write("Requirements =FileSystemDomain==\"fnal.gov\" && Arch==\"X86_64\"") - outjdl_file.write("Notification = ERROR\n") - outjdl_file.write("Should_Transfer_Files = YES\n") - outjdl_file.write("WhenToTransferOutput = ON_EXIT\n") - outjdl_file.write("request_memory = "+str(args.memory)+"\n") - outjdl_file.write("request_cpus = 8\n") - #outjdl_file.write("include : list-infiles.sh |\n") - if inputlist != "": - outjdl_file.write("Transfer_Input_Files = "+inputlist+"\n") - outjdl_file.write("x509userproxy = $ENV(X509_USER_PROXY)\n") - - outjdl_file.write("Output = "+output_log_path+"/"+args.outputdir+"_$(Process).stdout\n") - outjdl_file.write("Error = "+output_log_path+"/"+args.outputdir+"_$(Process).stdout\n") - outjdl_file.write("Log = "+output_log_path+"/"+args.outputdir+"_$(Process).log\n") - outjdl_file.write("Arguments = $(Process) $(Cluster)\n") - outjdl_file.write("Queue "+args.totaljobs+"\n") - -def create_sh_file_for_condor(args, command, output_folder): - outscript = open(args.jdlfilename+".sh", "w") - outscript.write('#!/bin/bash') - outscript.write("\n"+'echo "Starting job on " `date`') - outscript.write("\n"+'echo "Running on: `uname -a`"') - outscript.write("\n"+'echo "System software: `cat /etc/redhat-release`"') - outscript.write("\n"+'source /cvmfs/cms.cern.ch/cmsset_default.sh') - outscript.write("\n"+'export SCRAM_ARCH=slc6_amd64_gcc530') - outscript.write("\n"+'eval `scramv1 project CMSSW '+args.cmsswversion+'`') - outscript.write("\n"+'cd '+ args.cmsswversion + '/src/') - outscript.write("\n"+'eval `scram runtime -sh`') - outscript.write("\n"+'wget https://launchpad.net/mg5amcnlo/2.0/2.6.x/+download/MG5_aMC_v2.6.5.tar.gz') - outscript.write("\n"+'tar -xf MG5_aMC_v2.6.5.tar.gz') - outscript.write("\n"+'cd MG5_aMC_v2_6_5') - outscript.write("\n"+'cp ../../../'+args.tarfile+' .') - outscript.write("\n"+'echo "====> List files : " ') - outscript.write("\n"+'ls -alh') - outscript.write("\n"+'sed -i "s/39/'+str(random.randint(1,10))+'$1/g" '+args.tarfile) - outscript.write("\n"+'sed -i "14s/test/'+args.outputdir+'/g" '+args.tarfile) - outscript.write("\n"+'sed -i "17s/5000/'+str(args.nevents)+'/g" '+args.tarfile) - outscript.write("\n"+"=================================") - outscript.write("\n"+"cat "+args.tarfile) - outscript.write("\n"+"=================================") - outscript.write("\n"+command) - outscript.write("\n"+'echo "====> List files : " ') - outscript.write("\n"+'ls -alh') - outscript.write("\n"+'echo "====> List files in directory '+args.outputdir+': " ') - outscript.write("\n"+'ls '+args.outputdir) - outscript.write("\n"+'echo "====> List file in '+args.outputdir+'/Events : " ') - outscript.write("\n"+'ls '+args.outputdir+'/Events/') - outscript.write("\n"+'echo "====> List file in '+args.outputdir+'/Events/run_01/ : " ') - outscript.write("\n"+'ls '+args.outputdir+'/Events/run_01/') - outscript.write("\n"+'echo "====> copying *.lhe file to stores area..." ') - outscript.write("\n"+'echo "xrdcp -f '+args.outputdir+'/Events/run_01/*.lhe.gz root://cmseos.fnal.gov/' + output_folder+'/'+args.tarfile+'_'+'${1}'+'.lhe.gz"') - outscript.write("\n"+'xrdcp -f '+args.outputdir+'/Events/run_01/*.lhe.gz root://cmseos.fnal.gov/' + output_folder+'/'+args.tarfile+'_'+'${1}'+'.lhe.gz') - outscript.write("\n"+'echo "xrdcp -f '+args.outputdir+'/*.log root://cmseos.fnal.gov/' + output_folder+'/'+args.tarfile+'_'+'${1}'+'.log"') - outscript.write("\n"+'xrdcp -f '+args.outputdir+'/*.log root://cmseos.fnal.gov/' + output_folder+'/'+args.tarfile+'_'+'${1}'+'.log') - outscript.write("\n"+'echo "End job on " `date`') - outscript.write("\n"+'cd ${_CONDOR_SCRATCH_DIR}') - outscript.write("\n"+'rm -rf ' + args.cmsswversion) - outscript.write("\n") - outscript.close() - os.system("chmod 777 "+args.jdlfilename+".sh") diff --git a/generate_lhe_condor.py b/generate_lhe_condor.py deleted file mode 100644 index 684564b..0000000 --- a/generate_lhe_condor.py +++ /dev/null @@ -1,47 +0,0 @@ -#! /usr/bin/env python - -from general_condor_functions import getbasic_parser -from general_condor_functions import create_output_directory -from general_condor_functions import create_jdl_file_for_condor -from general_condor_functions import create_sh_file_for_condor - -def getargs(): - parser = getbasic_parser() - parser.add_argument('-cpu', '--ncpu', - default=8, - help='number of cpu to run' - ) - parser.add_argument('-r', '--randomnumber', - default=11, - help='random seed' - ) - return parser.parse_args() - - -def main(): - # get input arguments - args = getargs() - - # list of input files to be added in jdl file - inputlist = args.jdlfilename+".sh, "+args.InProcCardPath+"/"+args.tarfile - - # command to run - command = './bin/mg5_aMC '+args.tarfile - - # Get Output directory name - output_folder, output_log_path = create_output_directory(args) - - # create the jdl file for condor - create_jdl_file_for_condor(args, inputlist, output_log_path) - - # create the sh file for condor - create_sh_file_for_condor(args, command, output_folder) - - # running help - print "===> Set Proxy Using:" - print "\tvoms-proxy-init --voms cms --valid 168:00" - print "\"condor_submit "+args.jdlfilename+".jdl\" to submit" - - -if __name__ == "__main__": - main() diff --git a/mainScript.py b/mainScript.py new file mode 100644 index 0000000..a320a4d --- /dev/null +++ b/mainScript.py @@ -0,0 +1,188 @@ +import argparse +import os +import subprocess +import tarfile +import urllib.request +from contextlib import contextmanager + +from scripts.update_madgraph import update_cards +from utils.condor_job_template import CondorJobManager + + +@contextmanager +def change_dir(target_directory): + """A context manager to temporarily change the working directory.""" + original_directory = os.getcwd() # Save the current directory + try: + os.chdir(target_directory) # Change to the target directory + yield + finally: + os.chdir(original_directory) # Revert to the original directory + + +class MadGraphManager: + def __init__(self, condor_job_manager, url, nevents, iseed, cards_path, output_dir="output"): + self.condor_job_manager = condor_job_manager + self.url = url + self.nevents = nevents + self.iseed = iseed + self.cards_path = os.path.join(os.getcwd(), cards_path) + self.cards_dir_name = os.path.basename(cards_path) + self.mg5_directory = None + self.output_dir = output_dir + + def get_mg5_directory(self): + # Extract the expected directory name from the URL, formatted correctly + tar_file_name = self.url.split("/")[-1] + expected_dir_name = tar_file_name.replace(".tar.gz", "").replace(".", "_") + self.mg5_directory = expected_dir_name + + def confirm_directory(self, directory_path): + """Ask user to confirm if the found directory is the correct one.""" + answer = input( + f"Found existing directory '{directory_path}'. Use this directory? (y/n): " + ) + return answer.lower() == "y" + + def download_and_extract(self, extract_to="."): + """Download and extract tar.gz file, returning the path of the main directory. + If the directory already exists, skip the download and return the existing directory path. + """ + # if no URL is provided, check the expected MG5 directory in the current directory and confirm from user if this is what we want + # to check look for the directory name similar to "MG5_aMC_v3_5_4" + if not self.url: + for item in os.listdir(extract_to): + if os.path.isdir(os.path.join(extract_to, item)) and "MG5_aMC_v" in item: + full_path = os.path.join(extract_to, item) + if self.confirm_directory(full_path): + print("Using existing directory:", full_path) + self.mg5_directory = full_path + return + print( + "No suitable directory found or confirmed. Please provide a valid URL for download." + ) + + # Extract the expected directory name from the URL, formatted correctly + tar_file_name = self.url.split("/")[-1] + expected_dir_name = tar_file_name.replace(".tar.gz", "").replace(".", "_") + + # Check if the directory already exists + expected_path = os.path.join(extract_to, expected_dir_name) + if os.path.exists(expected_path) and os.path.isdir(expected_path): + print("Directory already exists:", expected_path) + self.mg5_directory = expected_path + + # Directory does not exist, proceed with download and extraction + file_path = os.path.join(extract_to, tar_file_name) + urllib.request.urlretrieve(self.url, file_path) + + with tarfile.open(file_path, "r:gz") as tar: + tar.extractall(path=extract_to) + directory_name = [member for member in tar.getmembers() if member.isdir()][0].name + + os.remove(file_path) + self.mg5_directory = os.path.join(extract_to, directory_name) + + def run_madgraph(self): + """Run MadGraph simulation locally.""" + # change the directory then run the command finally go back to original directory + with change_dir(self.mg5_directory): + self.mg5_command = 'bin/mg5_aMC' + # proc_card should be with the same name as the cards directory appended with "_proc_card.dat" + self.proc_card = os.path.join(self.cards_path, self.cards_dir_name + "_proc_card.dat") + if not os.path.exists(self.proc_card): + # check if proc card exists, if not terminate the program + print(f"Proc card not found: {self.proc_card}") + print("Please check the path and try again.") + return + + try: + # Start the MadGraph process + process = subprocess.Popen( + [self.mg5_command, self.proc_card], + stdin=subprocess.PIPE, + stdout=subprocess.PIPE, + stderr=subprocess.PIPE, + text=True, + ) + + # Interact with the process + while True: + # Read a line of output from MadGraph + output = process.stdout.readline() + if output == '' and process.poll() is not None: + break + if output: + print(output.strip()) + # Automatically answering prompts + if "Enter your command" in output: + process.stdin.write('command_to_send\n') + process.stdin.flush() + + # Check for any remaining errors + stderr_output = process.stderr.read() + if stderr_output: + print("Errors:", stderr_output.strip()) + + except Exception as e: + print(f"Failed to run MadGraph: {str(e)}") + + def setup_simulation(self): + self.download_and_extract() + # update_cards(self.nevents, self.iseed) + # create_condor_job() + + def create_condor_job(self): + job_details = { + "download_link": self.url, + "tarball_name": self.url.split("/")[-1], + "directory": self.mg5_directory, + } + + log_directory = os.path.join(base_directory, "logs") + + # Prepare and manage HTCondor jobs + self.condor_job_manager.prepare_job_files( + "my_simulation", + log_directory, + self.output_dir, + self.cards_path, + job_details, + num_jobs=1, + ) + + +if __name__ == "__main__": + parser = argparse.ArgumentParser(description="Setup and submit MG5 simulations.") + parser.add_argument( + "--url", type=str, required=False, help="URL of the MG5 tar file." + ) + parser.add_argument( + "--cards", type=str, required=True, help="Path to the MadGraph cards." + ) + parser.add_argument( + "--nevents", type=int, default=1000, help="Number of events to simulate." + ) + parser.add_argument( + "--iseed", type=int, default=11, help="Random seed for the simulation." + ) + parser.add_argument( + "--local", action='store_true', help="Run MadGraph locally." + ) + parser.add_argument( "--output", type=str, default="output", help="Output directory for the simulation.") + args = parser.parse_args() + + base_directory = os.getcwd() + scripts_directory = os.path.join(base_directory, "scripts") + condor_manager = CondorJobManager(base_directory, scripts_directory) + mg_manager = MadGraphManager( + condor_manager, args.url, args.nevents, args.iseed, args.cards, args.output + ) + if args.local: + # mg_manager.setup_simulation() + mg_manager.download_and_extract() + mg_manager.run_madgraph() + else: + mg_manager.get_mg5_directory() + mg_manager.create_condor_job() + pass diff --git a/scripts/job_description.jdl b/scripts/job_description.jdl new file mode 100644 index 0000000..d31c3cc --- /dev/null +++ b/scripts/job_description.jdl @@ -0,0 +1,17 @@ +Executable = {executable} +Universe = vanilla +Notification = ERROR +Should_Transfer_Files = YES +Transfer_Input_Files = {transfer_input_files} +Transfer_Output_Files = "" +Requirements = (TARGET.OpSysAndVer =?= "AlmaLinux9") +x509userproxy = $ENV(X509_USER_PROXY) +getenv = True ++JobFlavour = "longlunch" +request_memory = 12000 +request_cpus = 8 +Output = {log_path}/job_$(ClusterId)_$(ProcId).out +Error = {log_path}/job_$(ClusterId)_$(ProcId).err +Log = {log_path}/job_$(ClusterId)_$(ProcId).log +Arguments = $(ClusterId) $(ProcId) {output_dir} {card_dir} {card_dir_name} +Queue {n_jobs} diff --git a/scripts/run_madgraph.sh b/scripts/run_madgraph.sh new file mode 100644 index 0000000..8ceeaae --- /dev/null +++ b/scripts/run_madgraph.sh @@ -0,0 +1,46 @@ +#!/bin/bash +# Bash script to set up and execute a MadGraph simulation + +# Stop execution on any error +set -e + +echo "Job started on: $(date)" +echo "Running on: $(uname -a)" +echo "System software: $(cat /etc/redhat-release)" +source /cvmfs/cms.cern.ch/cmsset_default.sh +echo "###################################################" + +# Print input arguments +echo "Cluster ID: $1" +echo "Process ID: $2" +echo "Output Directory: $3" +echo "Input Cards Directory: $4" +echo "Cards Directory Name: $5" +echo "" + +# Setup environment +source /cvmfs/sft.cern.ch/lcg/views/dev4cuda/latest/x86_64-el9-gcc11-opt/setup.sh + +# Directory and file management +OUTPUT_DIR="$3/$5/$(date +%Y%m%d_%H%M%S)/$1_$2" +mkdir -p ${{OUTPUT_DIR}} + +# Download and setup MadGraph +wget -N "{mg_download_link}" -O "{mg_tarball}" +tar -zxf "{mg_tarball}" +cd {mg_dir} + +# Execute MadGraph using provided card +echo "Running MadGraph with: ../${{5}}/${{5}}_proc_card.dat" +./bin/mg5_aMC ../${{5}}/${{5}}_proc_card.dat + +# output directory, where the LHE file and other log files are stored +SOURCE_DIR="`pwd`" # copy full process directory + +# Copy the entire directory structure and files from SOURCE_DIR to OUTPUT_DIR +echo "Sync started on " $(date) +rsync -a "$SOURCE_DIR" "${{OUTPUT_DIR}}/" +echo "Sync finished on " $(date) + +echo "Job completed on: $(date)" +echo "==================================" diff --git a/scripts/update_madgraph.py b/scripts/update_madgraph.py new file mode 100644 index 0000000..c501189 --- /dev/null +++ b/scripts/update_madgraph.py @@ -0,0 +1,3 @@ +def update_cards(nevents, iseed): + # Read the card file, update parameters, and save. + print(f"Updating cards with nevents={nevents} and iseed={iseed}") diff --git a/utils/condor_job_template.py b/utils/condor_job_template.py new file mode 100644 index 0000000..fd95de4 --- /dev/null +++ b/utils/condor_job_template.py @@ -0,0 +1,75 @@ +import os + + +class CondorJobManager: + def __init__(self, base_dir, scripts_dir): + self.base_dir = base_dir + self.scripts_dir = scripts_dir + + def load_template(self, filename): + """Load a template file from the scripts directory.""" + path = os.path.join(self.scripts_dir, filename) + with open(path, "r") as file: + return file.read() + + def write_file(self, path, content): + """Write content to a file.""" + with open(path, "w") as file: + file.write(content) + print(f"File written successfully: {path}") + + def create_directory(self, path): + """Ensure directory exists.""" + os.makedirs(path, exist_ok=True) + print(f"Directory created or verified: {path}") + + def prepare_job_files( + self, job_name, log_path, output_dir, card_dir, mg_details, num_jobs=1 + ): + """Prepare and write HTCondor job files based on templates.""" + self.create_directory(log_path) + self.create_directory(output_dir) + + bash_script_content = self.load_template("run_madgraph.sh").format( + mg_download_link=mg_details["download_link"], + mg_tarball=mg_details["tarball_name"], + mg_dir=mg_details["directory"], + ) + + jdl_content = self.load_template("job_description.jdl").format( + executable=f"{job_name}.sh", + transfer_input_files=f"{job_name}.sh, {card_dir}", + log_path=log_path, + output_dir=output_dir, + card_dir=card_dir.split("/")[-1], + card_dir_name=os.path.basename(card_dir), + n_jobs=num_jobs, + ) + + self.write_file(f"{job_name}.sh", bash_script_content) + self.write_file(f"{job_name}.jdl", jdl_content) + print("HTCondor job files prepared successfully.") + + +if __name__ == "__main__": + job_details = { + "download_link": "https://launchpad.net/mg5amcnlo/3.0/3.5.x/+download/MG5_aMC_v3.5.4.tar.gz", + "tarball_name": "MG5_aMC_v3.5.4.tar.gz", + "directory": "MG5_aMC_v3_5_4", + } + + base_directory = "./" + scripts_directory = "./scripts" + log_directory = os.path.join(base_directory, "logs") + output_directory = os.path.join(base_directory, "output") + cards_directory = os.path.join(base_directory, "Cards/SM_process/st_tch_5f_ckm_NLO") + + manager = CondorJobManager(base_directory, scripts_directory) + manager.prepare_job_files( + "my_simulation", + log_directory, + output_directory, + cards_directory, + job_details, + num_jobs=1, + )