Skip to content
Snippets Groups Projects
Commit 51f33639 authored by Frank Siegert's avatar Frank Siegert
Browse files

Update to Sherpa 2.2.7

parent 95827201
No related branches found
No related tags found
No related merge requests found
## 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",
]
## Add Lund hadronisation model switches from arXiv:1905.09127
genSeq.Sherpa_i.Parameters += [
"PARJ(21)=0.36",
"PARJ(41)=0.3",
"PARJ(42)=0.6",
]
## 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/2.2.7_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]",
]
include("Sherpa_i/2.2.5_NNPDF30NNLO.py")
include("Sherpa_i/2.2.7_NNPDF30NNLO.py")
evgenConfig.description = "Sherpa 2.2.x example JO, Z+0,1-jet production."
evgenConfig.keywords = [ "2lepton" ]
......
......@@ -2,6 +2,6 @@
# File specifying the location of Sherpa to use.
#
set( SHERPA_VERSION 2.2.6p1 )
set( SHERPA_VERSION 2.2.7p1 )
set( SHERPA_ROOT
${LCG_RELEASE_DIR}/MCGenerators/sherpa/${SHERPA_VERSION}/${LCG_PLATFORM} )
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment