Commit 33dc547b authored by Ewelina Maria Lobodzinska's avatar Ewelina Maria Lobodzinska
Browse files

Merge branch '21.6-sherpa-common-jos' into '21.6'

Include JobOption base fragments in Sherpa_i.

See merge request atlas/athena!16707

Former-commit-id: ff7b1be64c86944ccaebe1014f259d29e872fce2
parents 0dbe0ad3 feb4cb35
6c7bc0a41c74f6e6e1734523748ec3f948bc63b7
\ No newline at end of file
......@@ -37,3 +37,5 @@ atlas_add_component( Sherpa_i
GaudiKernel GeneratorModulesLib TruthUtils AthenaBaseComps AthenaKernel )
add_definitions(-DSHERPA_ROOT="${SHERPA_ROOT}")
atlas_install_joboptions( share/common/*.py )
## Base config for Sherpa
from Sherpa_i.Sherpa_iConf import Sherpa_i
genSeq += Sherpa_i()
evgenConfig.generators = ["Sherpa"]
# Unfortunately this is broken in TestHepMC,
# cf. https://its.cern.ch/jira/browse/AGENE-1412
# -> disabled again
# ## Exclude status=11 and 20 particles from no-decay-vertex check
# testSeq.TestHepMC.NoDecayVertexStatuses = [1, 3, 4, 11, 20]
if hasattr(testSeq, "TestHepMC"):
testSeq.remove(TestHepMC())
## Tell Sherpa to read its run card sections from the jO
## TODO: write out Run.dat from genSeq.Sherpa_i.RunCard and read from it
## instead of the JO
genSeq.Sherpa_i.Parameters += [ 'RUNDATA=%s' % runArgs.jobConfig[0] ]
## Tell Sherpa to write logs into a separate file
## (need for production, looping job detection, Wolfgang Ehrenfeld)
genSeq.Sherpa_i.Parameters += [ 'LOG_FILE=sherpa.log' ]
## General ATLAS parameters
genSeq.Sherpa_i.Parameters += [
"BEAM_1=2212",
"BEAM_2=2212",
"MAX_PROPER_LIFETIME=10.0",
"HEPMC_TREE_LIKE=1",
"PRETTY_PRINT=Off",
]
## Cap event weights at |w|<10 in unweighted evgen; set H-event shower mode
## in the soft region to avoid instable weights
genSeq.Sherpa_i.Parameters += [
"OVERWEIGHT_THRESHOLD=10",
"PP_HPSMODE=0"
]
## Enable scale/PDF variations by default
import os
if os.environ.has_key('LHAPATH'):
pdfgridpath = os.environ['LHAPATH'].split(":")[-1]
if not os.path.isdir(pdfgridpath):
pdfgridpath = '/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current/'
genSeq.Sherpa_i.Parameters += [
"HEPMC_USE_NAMED_WEIGHTS=1",
"SCALE_VARIATIONS=0.25,0.25 0.25,1. 1.,0.25 1.,1. 1.,4. 4.,1. 4.,4.",
"PDF_VARIATIONS=NNPDF30NNLO[all] NNPDF30_nnlo_as_0117 NNPDF30_nnlo_as_0119 MMHT2014nnlo68cl CT14nnlo",
"NNPDF_GRID_PATH=%s" % pdfgridpath,
]
## Particle masses/widths
genSeq.Sherpa_i.Parameters += [
"MASS[6]=172.5",
"WIDTH[6]=1.32",
"MASS[15]=1.777",
"WIDTH[15]=2.26735e-12",
"MASS[23]=91.1876",
"WIDTH[23]=2.4952",
"MASS[24]=80.399",
"WIDTH[24]=2.085",
]
## Switch to EW_SCHEME=0 to be able to set PDG value of thetaW
genSeq.Sherpa_i.Parameters += [
"EW_SCHEME=0",
"SIN2THETAW=0.23113",
]
## improve b-fragmentation (https://its.cern.ch/jira/browse/AGENE-1426)
genSeq.Sherpa_i.Parameters += [
"SPLIT_LEADEXPONENT=1.0",
]
## set/add partial widths for H, W, Z to PDG values
## cf. https://sherpa.hepforge.org/doc/SHERPA-MC-2.2.4.html#HDH_005fWIDTH
genSeq.Sherpa_i.Parameters += [
"HDH_WIDTH[6,24,5]=1.32", #updated to comply latest updated top results
"HDH_WIDTH[-6,-24,-5]=1.32", #updated to comply latest updated top results
"HDH_WIDTH[25,5,-5]=2.35e-3",
"HDH_WIDTH[25,15,-15]=2.57e-4",
"HDH_WIDTH[25,13,-13]=8.91e-7",
"HDH_WIDTH[25,4,-4]=1.18e-4",
"HDH_WIDTH[25,3,-3]=1.00e-6",
"HDH_WIDTH[25,21,21]=3.49e-4",
"HDH_WIDTH[25,22,22]=9.28e-6",
"HDH_WIDTH[24,2,-1]=0.7041",
"HDH_WIDTH[24,4,-3]=0.7041",
"HDH_WIDTH[24,12,-11]=0.2256",
"HDH_WIDTH[24,14,-13]=0.2256",
"HDH_WIDTH[24,16,-15]=0.2256",
"HDH_WIDTH[-24,-2,1]=0.7041",
"HDH_WIDTH[-24,-4,3]=0.7041",
"HDH_WIDTH[-24,-12,11]=0.2256",
"HDH_WIDTH[-24,-14,13]=0.2256",
"HDH_WIDTH[-24,-16,15]=0.2256",
"HDH_WIDTH[23,1,-1]=0.3828",
"HDH_WIDTH[23,2,-2]=0.2980",
"HDH_WIDTH[23,3,-3]=0.3828",
"HDH_WIDTH[23,4,-4]=0.2980",
"HDH_WIDTH[23,5,-5]=0.3828",
"HDH_WIDTH[23,11,-11]=0.0840",
"HDH_WIDTH[23,12,-12]=0.1663",
"HDH_WIDTH[23,13,-13]=0.0840",
"HDH_WIDTH[23,14,-14]=0.1663",
"HDH_WIDTH[23,15,-15]=0.0840",
"HDH_WIDTH[23,16,-16]=0.1663",
]
## OpenLoops parameters
genSeq.Sherpa_i.Parameters += [
"OL_PARAMETERS=preset=2 write_parameters=1",
]
## Sherpa config with NNPDF3.0 NNLO PDF
include("Sherpa_i/Common_2.2.4_Base_Fragment.py")
## NNPDF3.0 NNLO is Sherpa's default PDF/tune, thus no need to set anything up
evgenConfig.tune = "NNPDF3.0 NNLO"
## Base config for Sherpa
from Sherpa_i.Sherpa_iConf import Sherpa_i
genSeq += Sherpa_i()
evgenConfig.generators = ["Sherpa"]
## Tell Sherpa to read its run card sections from the jO
## TODO: write out Run.dat from genSeq.Sherpa_i.RunCard and read from it
## instead of the JO
genSeq.Sherpa_i.Parameters += [ 'RUNDATA=%s' % runArgs.jobConfig[0] ]
## Tell Sherpa to write logs into a separate file
## (need for production, looping job detection, Wolfgang Ehrenfeld)
genSeq.Sherpa_i.Parameters += [ 'LOG_FILE=sherpa.log' ]
## General ATLAS parameters
genSeq.Sherpa_i.Parameters += [
"BEAM_1=2212",
"BEAM_2=2212",
"MAX_PROPER_LIFETIME=10.0",
"HEPMC_TREE_LIKE=1",
"PRETTY_PRINT=Off",
]
## Cap event weights at |w|<10 in unweighted evgen; set H-event shower mode
## in the soft region to avoid instable weights
genSeq.Sherpa_i.Parameters += [
"OVERWEIGHT_THRESHOLD=10",
"PP_HPSMODE=0"
]
## Enable scale variations by default
genSeq.Sherpa_i.Parameters += [
"HEPMC_USE_NAMED_WEIGHTS=1",
"SCALE_VARIATIONS=0.25,0.25 0.25,1. 1.,0.25 1.,1. 1.,4. 4.,1. 4.,4.",
]
## Particle masses/widths
genSeq.Sherpa_i.Parameters += [
"MASS[6]=172.5",
"WIDTH[6]=1.32",
"MASS[15]=1.777",
"WIDTH[15]=2.26735e-12",
"MASS[23]=91.1876",
"WIDTH[23]=2.4952",
"MASS[24]=80.399",
"WIDTH[24]=2.085",
]
## Switch to EW_SCHEME=0 to be able to set PDG value of thetaW
genSeq.Sherpa_i.Parameters += [
"EW_SCHEME=0",
"SIN2THETAW=0.23113",
]
## set/add partial widths for H, W, Z to PDG values
## cf. https://sherpa.hepforge.org/doc/SHERPA-MC-2.2.4.html#HDH_005fWIDTH
genSeq.Sherpa_i.Parameters += [
"HDH_WIDTH[6,24,5]=1.32", #updated to comply latest updated top results
"HDH_WIDTH[-6,-24,-5]=1.32", #updated to comply latest updated top results
"HDH_WIDTH[25,5,-5]=2.35e-3",
"HDH_WIDTH[25,15,-15]=2.57e-4",
"HDH_WIDTH[25,13,-13]=8.91e-7",
"HDH_WIDTH[25,4,-4]=1.18e-4",
"HDH_WIDTH[25,3,-3]=1.00e-6",
"HDH_WIDTH[25,21,21]=3.49e-4",
"HDH_WIDTH[25,22,22]=9.28e-6",
"HDH_WIDTH[24,2,-1]=0.7041",
"HDH_WIDTH[24,4,-3]=0.7041",
"HDH_WIDTH[24,12,-11]=0.2256",
"HDH_WIDTH[24,14,-13]=0.2256",
"HDH_WIDTH[24,16,-15]=0.2256",
"HDH_WIDTH[-24,-2,1]=0.7041",
"HDH_WIDTH[-24,-4,3]=0.7041",
"HDH_WIDTH[-24,-12,11]=0.2256",
"HDH_WIDTH[-24,-14,13]=0.2256",
"HDH_WIDTH[-24,-16,15]=0.2256",
"HDH_WIDTH[23,1,-1]=0.3828",
"HDH_WIDTH[23,2,-2]=0.2980",
"HDH_WIDTH[23,3,-3]=0.3828",
"HDH_WIDTH[23,4,-4]=0.2980",
"HDH_WIDTH[23,5,-5]=0.3828",
"HDH_WIDTH[23,11,-11]=0.0840",
"HDH_WIDTH[23,12,-12]=0.1663",
"HDH_WIDTH[23,13,-13]=0.0840",
"HDH_WIDTH[23,14,-14]=0.1663",
"HDH_WIDTH[23,15,-15]=0.0840",
"HDH_WIDTH[23,16,-16]=0.1663",
]
## OpenLoops parameters
genSeq.Sherpa_i.Parameters += [
"OL_PARAMETERS=preset=2 write_parameters=1",
]
## Sherpa config with NNPDF3.0 NNLO PDF
include("Sherpa_i/Common_2.2.5_Base_Fragment.py")
evgenConfig.tune = "NNPDF3.0 NNLO"
genSeq.Sherpa_i.Parameters += [
"PDF_LIBRARY=LHAPDFSherpa",
"USE_PDF_ALPHAS=1",
"PDF_SET=NNPDF30_nnlo_as_0118",
"PDF_VARIATIONS=NNPDF30_nnlo_as_0118[all] NNPDF30_nnlo_as_0117 NNPDF30_nnlo_as_0119 MMHT2014nnlo68cl CT14nnlo PDF4LHC15_nnlo_30_pdfas[all]",
]
Supports Markdown
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