Commit c37bfdba authored by Ewelina Maria Lobodzinska's avatar Ewelina Maria Lobodzinska Committed by Spyros Argyropoulos
Browse files

add common fragments [skip modfiles] [skip athena]

parent 9178b9a2
Pipeline #1003638 skipped
# MC16 Stop-Stau Simplified Model Grid, JIRA: ATLMCPROD-7163
include('MadGraphControl/MadGraphControl_SimplifiedModel_TT_stau.py')
include('MadGraphControl/SimplifiedModel_TT_stau.py')
# Warnings:
# UserHooks Warning: the calculation of uncertainty variations will only be consistent in the absence of any external modifications to the shower branching probabilities via the UserHooks framework. It is therefore strongly advised to avoid combining the automated uncertainty calculations with any such UserHooks modifications.
#
#Merging Warning: in multi-jet merging approaches, trial showers are used to generate missing Sudakov factor corrections to the hard matrix elements. Currently that framework is not consistently combined with the variations introduced here, so the two should not be used simultaneously. This means shower weights should NOT be used with CKKLW or FxFx merging!
#
#
# Pythia8 shower weights are only available in release 2.26 and later. The
# test below checks the Pythia8 release and also verifies that the Pythia_i
# tag is recent enought to store the shower weights.
if "PYTHIA8VER" in os.environ:
verph8str = str(os.environ['PYTHIA8VER'])[:3]
verph8 = int(verph8str)
print("verph8 ", verph8)
else:
print("PYTHIA8VER not available !!! assuming version 2.40)
verph8 = 240
if (verph8 >= 226):
if "ShowerWeightNames" in genSeq.Pythia8.__slots__.keys():
print "Initalizing Shower Weights from Pythia8_ShowerWeights.py"
if("NNPDF" in evgenConfig.tune):
genSeq.Pythia8.Commands += [
'UncertaintyBands:doVariations = on',
'UncertaintyBands:List = {\
Var3cUp isr:muRfac=0.549241,\
Var3cDown isr:muRfac=1.960832,\
isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
hardHi fsr:cNS=2.0 isr:cNS=2.0,\
hardLo fsr:cNS=-2.0 isr:cNS=-2.0,\
isr:PDF:plus=1,\
isr:PDF:minus=2\
}']
genSeq.Pythia8.ShowerWeightNames = [
"Var3cUp",
"Var3cDown",
"isr:muRfac=2.0_fsr:muRfac=2.0",
"isr:muRfac=2.0_fsr:muRfac=1.0",
"isr:muRfac=2.0_fsr:muRfac=0.5",
"isr:muRfac=1.0_fsr:muRfac=2.0",
"isr:muRfac=1.0_fsr:muRfac=0.5",
"isr:muRfac=0.5_fsr:muRfac=2.0",
"isr:muRfac=0.5_fsr:muRfac=1.0",
"isr:muRfac=0.5_fsr:muRfac=0.5",
"isr:muRfac=1.75_fsr:muRfac=1.0",
"isr:muRfac=1.5_fsr:muRfac=1.0",
"isr:muRfac=1.25_fsr:muRfac=1.0",
"isr:muRfac=0.625_fsr:muRfac=1.0",
"isr:muRfac=0.75_fsr:muRfac=1.0",
"isr:muRfac=0.875_fsr:muRfac=1.0",
"isr:muRfac=1.0_fsr:muRfac=1.75",
"isr:muRfac=1.0_fsr:muRfac=1.5",
"isr:muRfac=1.0_fsr:muRfac=1.25",
"isr:muRfac=1.0_fsr:muRfac=0.625",
"isr:muRfac=1.0_fsr:muRfac=0.75",
"isr:muRfac=1.0_fsr:muRfac=0.875",
"hardHi",
"hardLo",
"isr:PDF:plus",
"isr:PDF:minus"
]
else:
genSeq.Pythia8.Commands += [
'UncertaintyBands:doVariations = on',
'UncertaintyBands:List = {\
Var3cUp isr:muRfac=0.549241,\
Var3cDown isr:muRfac=1.960832,\
isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
hardHi fsr:cNS=2.0 isr:cNS=2.0,\
hardLo fsr:cNS=-2.0 isr:cNS=-2.0\
}']
genSeq.Pythia8.ShowerWeightNames = [
"Var3cUp",
"Var3cDown",
"isr:muRfac=2.0_fsr:muRfac=2.0",
"isr:muRfac=2.0_fsr:muRfac=1.0",
"isr:muRfac=2.0_fsr:muRfac=0.5",
"isr:muRfac=1.0_fsr:muRfac=2.0",
"isr:muRfac=1.0_fsr:muRfac=0.5",
"isr:muRfac=0.5_fsr:muRfac=2.0",
"isr:muRfac=0.5_fsr:muRfac=1.0",
"isr:muRfac=0.5_fsr:muRfac=0.5",
"isr:muRfac=1.75_fsr:muRfac=1.0",
"isr:muRfac=1.5_fsr:muRfac=1.0",
"isr:muRfac=1.25_fsr:muRfac=1.0",
"isr:muRfac=0.625_fsr:muRfac=1.0",
"isr:muRfac=0.75_fsr:muRfac=1.0",
"isr:muRfac=0.875_fsr:muRfac=1.0",
"isr:muRfac=1.0_fsr:muRfac=1.75",
"isr:muRfac=1.0_fsr:muRfac=1.5",
"isr:muRfac=1.0_fsr:muRfac=1.25",
"isr:muRfac=1.0_fsr:muRfac=0.625",
"isr:muRfac=1.0_fsr:muRfac=0.75",
"isr:muRfac=1.0_fsr:muRfac=0.875",
"hardHi",
"hardLo"
]
else:
print "Pythia8_i version too old for shower weights"
else:
print "No shower weights since not supported in Pythia8.",verph8
# Warnings:
# UserHooks Warning: the calculation of uncertainty variations will only be consistent in the absence of any external modifications to the shower branching probabilities via the UserHooks framework. It is therefore strongly advised to avoid combining the automated uncertainty calculations with any such UserHooks modifications.
#
#Merging Warning: in multi-jet merging approaches, trial showers are used to generate missing Sudakov factor corrections to the hard matrix elements. Currently that framework is not consistently combined with the variations introduced here, so the two should not be used simultaneously. This means shower weights should NOT be used with CKKLW or FxFx merging!
#
#
# Pythia8 shower weights are only available in release 2.26 and later. The
# test below checks the Pythia8 release and also verifies that the Pythia_i
# tag is recent enought to store the shower weights.
if "PYTHIA8VER" in os.environ:
verph8str = str(os.environ['PYTHIA8VER'])[:3]
verph8 = int(verph8str)
print("verph8 ", verph8)
else:
print("PYTHIA8VER not available !!! assuming version 2.40)
verph8 = 240
if (verph8 >= 226):
if "ShowerWeightNames" in genSeq.Pythia8.__slots__.keys():
print "Initalizing Shower Weights from Pythia8_ShowerWeights.py"
if("NNPDF" in evgenConfig.tune):
genSeq.Pythia8.Commands += [
'UncertaintyBands:doVariations = on',
'UncertaintyBands:List = {\
Var3cUp isr:muRfac=0.549241,\
Var3cDown isr:muRfac=1.960832,\
isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
hardHi fsr:cNS=2.0 isr:cNS=2.0,\
hardLo fsr:cNS=-2.0 isr:cNS=-2.0,\
isr:PDF:plus=1,\
isr:PDF:minus=2\
}']
genSeq.Pythia8.ShowerWeightNames = [
"Var3cUp",
"Var3cDown",
"isr:muRfac=2.0_fsr:muRfac=2.0",
"isr:muRfac=2.0_fsr:muRfac=1.0",
"isr:muRfac=2.0_fsr:muRfac=0.5",
"isr:muRfac=1.0_fsr:muRfac=2.0",
"isr:muRfac=1.0_fsr:muRfac=0.5",
"isr:muRfac=0.5_fsr:muRfac=2.0",
"isr:muRfac=0.5_fsr:muRfac=1.0",
"isr:muRfac=0.5_fsr:muRfac=0.5",
"isr:muRfac=1.75_fsr:muRfac=1.0",
"isr:muRfac=1.5_fsr:muRfac=1.0",
"isr:muRfac=1.25_fsr:muRfac=1.0",
"isr:muRfac=0.625_fsr:muRfac=1.0",
"isr:muRfac=0.75_fsr:muRfac=1.0",
"isr:muRfac=0.875_fsr:muRfac=1.0",
"isr:muRfac=1.0_fsr:muRfac=1.75",
"isr:muRfac=1.0_fsr:muRfac=1.5",
"isr:muRfac=1.0_fsr:muRfac=1.25",
"isr:muRfac=1.0_fsr:muRfac=0.625",
"isr:muRfac=1.0_fsr:muRfac=0.75",
"isr:muRfac=1.0_fsr:muRfac=0.875",
"hardHi",
"hardLo",
"isr:PDF:plus",
"isr:PDF:minus"
]
else:
genSeq.Pythia8.Commands += [
'UncertaintyBands:doVariations = on',
'UncertaintyBands:List = {\
Var3cUp isr:muRfac=0.549241,\
Var3cDown isr:muRfac=1.960832,\
isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
hardHi fsr:cNS=2.0 isr:cNS=2.0,\
hardLo fsr:cNS=-2.0 isr:cNS=-2.0\
}']
genSeq.Pythia8.ShowerWeightNames = [
"Var3cUp",
"Var3cDown",
"isr:muRfac=2.0_fsr:muRfac=2.0",
"isr:muRfac=2.0_fsr:muRfac=1.0",
"isr:muRfac=2.0_fsr:muRfac=0.5",
"isr:muRfac=1.0_fsr:muRfac=2.0",
"isr:muRfac=1.0_fsr:muRfac=0.5",
"isr:muRfac=0.5_fsr:muRfac=2.0",
"isr:muRfac=0.5_fsr:muRfac=1.0",
"isr:muRfac=0.5_fsr:muRfac=0.5",
"isr:muRfac=1.75_fsr:muRfac=1.0",
"isr:muRfac=1.5_fsr:muRfac=1.0",
"isr:muRfac=1.25_fsr:muRfac=1.0",
"isr:muRfac=0.625_fsr:muRfac=1.0",
"isr:muRfac=0.75_fsr:muRfac=1.0",
"isr:muRfac=0.875_fsr:muRfac=1.0",
"isr:muRfac=1.0_fsr:muRfac=1.75",
"isr:muRfac=1.0_fsr:muRfac=1.5",
"isr:muRfac=1.0_fsr:muRfac=1.25",
"isr:muRfac=1.0_fsr:muRfac=0.625",
"isr:muRfac=1.0_fsr:muRfac=0.75",
"isr:muRfac=1.0_fsr:muRfac=0.875",
"hardHi",
"hardLo"
]
else:
print "Pythia8_i version too old for shower weights"
else:
print "No shower weights since not supported in Pythia8.",verph8
from MadGraphControl.MadGraphUtils import *
# General settings
nevents=150000
mode=0
nJobs=1
gridpack_dir=None
gridpack_mode=False
cluster_type=None
cluster_queue=None
# MG Particle cuts
mllcut=40
# Merging settings
maxjetflavor=5
ickkw=0
nJetMax=4
ktdurham=30
dparameter=0.4
### DSID lists (extensions can include filters, systematics samples, etc.)
Zee_5fl_Np0=[361500]
Zee_5fl_Np1=[361501]
Zee_5fl_Np2=[361502]
Zee_5fl_Np3=[361503]
Zee_5fl_Np4=[361504]
Zmumu_5fl_Np0=[361505, 421107]
Zmumu_5fl_Np1=[361506]
Zmumu_5fl_Np2=[361507]
Zmumu_5fl_Np3=[361508]
Zmumu_5fl_Np4=[361509]
Ztautau_5fl_Np0=[361510]
Ztautau_5fl_Np1=[361511]
Ztautau_5fl_Np2=[361512]
Ztautau_5fl_Np3=[361513]
Ztautau_5fl_Np4=[361514]
Znunu_5fl_Np0=[361515]
Znunu_5fl_Np1=[361516]
Znunu_5fl_Np2=[361517]
Znunu_5fl_Np3=[361518]
Znunu_5fl_Np4=[361519]
### Electrons
if hasattr(runArgs,'runNumber'):
runNum = int(runArgs.runNumber)
else:
parDSID = str(runArgs.jobConfig[0])
if len(parDSID)==6 and parDSID.isdigit():
runNum = int(parDSID)
else:
runNum = int(parDSID[-6:])
if runNum in Zee_5fl_Np0:
mgproc="generate p p > e+ e- @0"
name='Zee_Np0'
process="pp>e+e-"
elif runNum in Zee_5fl_Np1:
mgproc="generate p p > e+ e- j @1"
name='Zee_Np1'
process="pp>e+e-"
elif runNum in Zee_5fl_Np2:
mgproc="generate p p > e+ e- j j @2"
name='Zee_Np2'
process="pp>e+e-"
elif runNum in Zee_5fl_Np3:
mgproc="generate p p > e+ e- j j j @3"
name='Zee_Np3'
process="pp>e+e-"
elif runNum in Zee_5fl_Np4:
mgproc="generate p p > e+ e- j j j j @4"
name='Zee_Np4'
process="pp>e+e-"
nevents=5000
gridpack_mode=True
gridpack_dir='madevent/'
mode=1
cluster_type='pbs'
cluster_queue='medium'
nJobs=20
### Muons
elif runNum in Zmumu_5fl_Np0:
mgproc="generate p p > mu+ mu- @0"
name='Zmumu_Np0'
process="pp>mu+mu-"
elif runNum in Zmumu_5fl_Np1:
mgproc="generate p p > mu+ mu- j @1"
name='Zmumu_Np1'
process="pp>mu+mu-"
elif runNum in Zmumu_5fl_Np2:
mgproc="generate p p > mu+ mu- j j @2"
name='Zmumu_Np2'
process="pp>mu+mu-"
elif runNum in Zmumu_5fl_Np3:
mgproc="generate p p > mu+ mu- j j j @3"
name='Zmumu_Np3'
process="pp>mu+mu-"
elif runNum in Zmumu_5fl_Np4:
mgproc="generate p p > mu+ mu- j j j j @4"
name='Zmumu_Np4'
process="pp>mu+mu-"
nevents=5000
gridpack_mode=True
gridpack_dir='madevent/'
mode=1
cluster_type='pbs'
cluster_queue='medium'
nJobs=20
### Taus
elif runNum in Ztautau_5fl_Np0:
mgproc="generate p p > ta+ ta- @0"
name='Ztautau_Np0'
process="pp>ta+ta-"
elif runNum in Ztautau_5fl_Np1:
mgproc="generate p p > ta+ ta- j @1"
name='Ztautau_Np1'
process="pp>ta+ta-"
elif runNum in Ztautau_5fl_Np2:
mgproc="generate p p > ta+ ta- j j @2"
name='Ztautau_Np2'
process="pp>ta+ta-"
elif runNum in Ztautau_5fl_Np3:
mgproc="generate p p > ta+ ta- j j j @3"
name='Ztautau_Np3'
process="pp>ta+ta-"
elif runNum in Ztautau_5fl_Np4:
mgproc="generate p p > ta+ ta- j j j j @4"
name='Ztautau_Np4'
process="pp>ta+ta-"
nevents=5000
gridpack_mode=True
gridpack_dir='madevent/'
mode=1
cluster_type='pbs'
cluster_queue='medium'
nJobs=20
### Neutrinos
elif runNum in Znunu_5fl_Np0:
mgproc="generate p p > vl vl~ @0"
name='Znunu_Np0'
process="pp>NEUTRINOS"
elif runNum in Znunu_5fl_Np1:
mgproc="generate p p > vl vl~ j @1"
name='Znunu_Np1'
process="pp>NEUTRINOS"
elif runNum in Znunu_5fl_Np2:
mgproc="generate p p > vl vl~ j j @2"
name='Znunu_Np2'
process="pp>NEUTRINOS"
elif runNum in Znunu_5fl_Np3:
mgproc="generate p p > vl vl~ j j j @3"
name='Znunu_Np3'
process="pp>NEUTRINOS"
elif runNum in Znunu_5fl_Np4:
mgproc="generate p p > vl vl~ j j j j @4"
name='Znunu_Np4'
process="pp>NEUTRINOS"
nevents=5000
gridpack_mode=True
gridpack_dir='madevent/'
mode=1
cluster_type='pbs'
cluster_queue='medium'
nJobs=20
else:
raise RuntimeError("runNumber %i not recognised in these jobOptions."%runNum)
stringy = 'madgraph.'+str(runNum)+'.MadGraph_'+str(name)
fcard = open('proc_card_mg5.dat','w')
fcard.write("""
import model sm-no_b_mass
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
"""+mgproc+"""
output -f
""")
fcard.close()
beamEnergy=-999
if hasattr(runArgs,'ecmEnergy'):
beamEnergy = runArgs.ecmEnergy / 2.
else:
raise RuntimeError("No center of mass energy found.")
process_dir = new_process(grid_pack=gridpack_dir)
#Fetch default LO run_card.dat and set parameters
extras = { 'lhe_version' : '2.0',
'cut_decays' : 'F',
'pdlabel' : "'nn23lo1'",
# 'pdlabel' : "'lhapdf'",
# 'lhaid' : 247000,
'maxjetflavor' : maxjetflavor,
'asrwgtflavor' : maxjetflavor,
'ickkw' : 0,
'ptj' : 20,
'ptb' : 20,
'mmll' : mllcut,
'mmjj' : 0,
'drjj' : 0,
'drll' : 0,
'drjl' : 0.4,
'ptl' : 0,
'etal' : 10,
'etab' : 6,
'etaj' : 6,
'ktdurham' : ktdurham,
'dparameter' : dparameter,
'use_syst' : 'F' }
build_run_card(run_card_old=get_default_runcard(proc_dir=process_dir),run_card_new='run_card.dat',
nevts=nevents,rand_seed=runArgs.randomSeed,beamEnergy=beamEnergy,xqcut=0.,
extras=extras)
print_cards()
generate(run_card_loc='run_card.dat',param_card_loc=None,mode=mode,njobs=nJobs,proc_dir=process_dir,
grid_pack=gridpack_mode,gridpack_dir=gridpack_dir,cluster_type=cluster_type,cluster_queue=cluster_queue,
nevents=nevents,random_seed=runArgs.randomSeed)
arrange_output(proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz')
#### Shower
evgenConfig.description = 'MadGraph_'+str(name)
evgenConfig.keywords+=['Z','electron','jets','drellYan']
evgenConfig.inputfilecheck = stringy
runArgs.inputGeneratorFile=stringy+'._00001.events.tar.gz'
include("Pythia8_i/Pythia8_A14_NNPDF23LO_EvtGen_Common.py")
include("Pythia8_i/Pythia8_MadGraph.py")
PYTHIA8_nJetMax=nJetMax
PYTHIA8_TMS=ktdurham
PYTHIA8_Dparameter=dparameter
PYTHIA8_Process=process
PYTHIA8_nQuarksMerge=maxjetflavor
include("Pythia8_i/Pythia8_CKKWL_kTMerge.py")
from MadGraphControl.MadGraphUtils import *
# General settings
minevents=5000
nevents=30000
mode=0
# MG Particle cuts
mllcut=-1
# MG merging settings
maxjetflavor=5
ickkw=0
# Pythia8 merging settings
nJetMax=2
ktdurham=30
dparameter=0.4
### DSID lists (extensions can include systematics samples)
ttZllonshell_Np0 =[