Commit 8c3270b3 authored by Marjorie Shapiro's avatar Marjorie Shapiro
Browse files

update Pythia8_ShowerWeights


git-svn-id: svn+ssh://svn.cern.ch/reps/atlasoff/Generators/MC15JobOptions/trunk@813045 4525493e-7705-40b1-a816-d608a930855b
parent 376253a2
2018-02-17 Marjorie Shapiro <mdshapiro@lbl.gov>
* Update Pythia8_ShowerWeights.py for release 21 and cases where pdf doesn't have variations
* tagging MC15JobOptions-00-09-15
2018-02-17 Christian Gutschow <chris.g@cern.ch>
* tagging MC15JobOptions-00-09-14
* csv update for 364250, 364253-5
......
......@@ -7,14 +7,22 @@
# 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.
ver = os.popen("cmt show versions External/Pythia8").read()
cmak_sys = os.getenv('CMAKE_PREFIX_PATH')
if cmak_sys:
ver = 'Pythia8-0230'
print "cmake system"
else:
ver = os.popen("cmt show versions External/Pythia8").read()
parsed = ver.split(" ")
pieces = parsed[1].split("-")
print "Pythia Version 8.",pieces[2]
if int(pieces[2])>=26:
if "ShowerWeightNames" in genSeq.Pythia8.__slots__.keys():
print "Initalizing Shower Weights in Pythia8_Powheg_Main31.py"
genSeq.Pythia8.Commands += [
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,\
......@@ -40,12 +48,12 @@ if int(pieces[2])>=26:
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\
hardLo fsr:cNS=-2.0 isr:cNS=-2.0,\
isr:PDF:plus=1,\
isr:PDF:minus=2,\
isr:PDF:minus=2\
}']
genSeq.Pythia8.ShowerWeightNames = [
genSeq.Pythia8.ShowerWeightNames = [
"Var3cUp",
"Var3cDown",
"isr:muRfac=2.0_fsr:muRfac=2.0",
......@@ -73,9 +81,64 @@ if int(pieces[2])>=26:
"isr:PDF:plus",
"isr:PDF:minus"
]
else:
print "Pythia8_i version too old to support isr/fsr shower weights"
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.2",pieces[2]
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment