## Herwig++ config for the CTEQ6L1 UE-EE-5 tune series with a NNPDF3.0 NLO ME PDF
include("MC15JobOptions/Herwigpp_Base_Fragment.py")
## Construct command set
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(runArgs.ecmEnergy) + hw.base_cmds() \
+ hw.nlo_pdf_cmds("NNPDF30_nlo_as_0118.LHgrid", "cteq6ll.LHpdf") \
+ hw.ue_tune_cmds("UE-EE-5-CTEQ6L1")
genSeq.Herwigpp.Commands += cmds.splitlines()
del cmds
evgenConfig.tune = "CTEQ6L1-UE-EE-5"
## Job options file for Herwig++, W -> tau nu_tau production
include("GeneratorUtils/StdEvgenSetup.py")
# ## Control log levels
# svcMgr.MessageSvc.OutputLevel = DEBUG
## Add H++ to the job list
from Herwigpp_i.Herwigpp_iConf import Herwigpp
topAlg += Herwigpp()
## Get basic Herwig++ Atlas tune params
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(7000) + hw.base_cmds() + hw.lo_pdf_cmds("cteq6ll.LHpdf") + hw.ue_tune_cmds("CTEQ6L1-UE-EE-7000-3")
## Add to commands
cmds += """
## Set up qq -> W -> tau nu_tau process
insert /Herwig/MatrixElements/SimpleQCD:MatrixElements[0] /Herwig/MatrixElements/MEqq2W2ff
set /Herwig/MatrixElements/SimpleQCD:MatrixElements[0]:Process Tau
"""
## Set commands
topAlg.Herwigpp.Commands = cmds.splitlines()
## Riccardo Di Sipio
##
## See bottom of this file for an example POWHEG config
include("GeneratorUtils/StdEvgenSetup.py")
# ## Control log levels
# svcMgr.MessageSvc.OutputLevel = DEBUG
## Add H++ to the job list
from Herwigpp_i.Herwigpp_iConf import Herwigpp
topAlg += Herwigpp()
from Herwigpp_i import config as hw
## According to Paolo Nason, using NLO PDFs for ME and LO for PS/UE is okay
cmds = hw.energy_cmds(7000) + hw.base_cmds() + hw.lo_pdf_cmds(
"cteq6ll.LHpdf") + hw.ue_tune_cmds("CTEQ6L1-UE-EE-7000-3")
# TODO: Check that mass resets (to match POWHEG) are approved.
cmds += """
setup /Herwig/Particles/t 6 t 172.5 1.4 14 0 2 3 2 0
#setup /Herwig/Particles/tbar -6 tbar 174.2 1.4 14 0 -2 -3 2 0
set /Herwig/Particles/u:NominalMass 0.320
set /Herwig/Particles/ubar:NominalMass 0.320
set /Herwig/Particles/d:NominalMass 0.320
set /Herwig/Particles/dbar:NominalMass 0.320
set /Herwig/Particles/s:NominalMass 0.5
set /Herwig/Particles/sbar:NominalMass 0.5
## Herwig++ config for the MRSTMCal UE-EE-5 tune series with a NNPDF3.0 NLO ME PDF
include("MC15JobOptions/Herwigpp_Base_Fragment.py")
## Construct command set
from Herwigpp_i import config as hw
cmds = (
hw.energy_cmds(runArgs.ecmEnergy)
+ hw.base_cmds()
+ hw.nlo_pdf_cmds("NNPDF30_nlo_as_0118.LHgrid", "MRSTMCal.LHgrid")
+ hw.ue_tune_cmds("UE-EE-5-LO**")
)
genSeq.Herwigpp.Commands += cmds.splitlines()
del cmds
evgenConfig.tune = "LO**-UE-EE-5"
## Generic job options file for Herwig++
include("GeneratorUtils/StdEvgenSetup.py")
# ## Control log levels
# svcMgr.MessageSvc.OutputLevel = DEBUG
## Add H++ to the job list
from Herwigpp_i.Herwigpp_iConf import Herwigpp
topAlg += Herwigpp()
## Get basic Herwig++ Atlas tune params
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(1960) + hw.base_cmds() + hw.lo_pdf_cmds(
"cteq6ll.LHpdf") + hw.ue_tune_cmds("CTEQ6L1-UE-EE-1800-1")
## Set jet min pT cut: override on Athena command line with e.g. "-c JETPTCUT=30"
JET_PT_CUT = 1
if 'JETPTCUT' in dir():
JET_PT_CUT = JETPTCUT
## Add to commands
cmds += """
set /Herwig/Generators/LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
insert /Herwig/MatrixElements/SimpleQCD:MatrixElements[0] /Herwig/MatrixElements/MEQCD2to2
set /Herwig/Cuts/JetKtCut:MinKT %f*GeV
set /Herwig/Cuts/QCDCuts:MHatMin %f*GeV
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
""" % (JET_PT_CUT, 2 * JET_PT_CUT)
## Use weighted events with an enhanced pT spectrum
## Herwig++ config for the CTEQ6L1 UE-EE-4 tune series with an LO ME PDF
include("MC15JobOptions/Herwigpp_Base_Fragment.py")
## Construct command set
from Herwigpp_i import config as hw
cmds = (
hw.energy_cmds(int(runArgs.ecmEnergy))
+ hw.base_cmds()
+ hw.lo_pdf_cmds("cteq6ll.LHpdf")
+ hw.ue_tune_cmds("UE-EE-4-CTEQ6L1")
)
genSeq.Herwigpp.Commands += cmds.splitlines()
del cmds
evgenConfig.tune = "CTEQ6L1-UE-EE-4"
## Job options file for Herwig++, LEP1 Z -> mu+ mu- production
include("GeneratorUtils/StdEvgenSetup.py")
# ## Control log levels
# svcMgr.MessageSvc.OutputLevel = DEBUG
## Add H++ to the job list
from Herwigpp_i.Herwigpp_iConf import Herwigpp
topAlg += Herwigpp()
## Get basic Herwig++ Atlas tune params
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(91.2) + hw.base_cmds()
## Add to commands
cmds += """
## Set up LEP1 e+ e- initial state
set /Herwig/Generators/LHCGenerator:EventHandler:CascadeHandler:MPIHandler NULL
set /Herwig/Generators/LHCGenerator:EventHandler:BeamA /Herwig/Particles/e+
set /Herwig/Generators/LHCGenerator:EventHandler:BeamB /Herwig/Particles/e-
## e+ e- -> Z/gamma* -> mu+ mu- process
insert /Herwig/MatrixElements/SimpleQCD:MatrixElements[0] /Herwig/MatrixElements/MEee2gZ2ll
set /Herwig/MatrixElements/SimpleQCD:MatrixElements[0]:Allowed Muon
## Tight windowing around Z mass (for MC-TESTER comparison to PHOTOS)
set /Herwig/Cuts/MassCut:MinM 90.2*GeV
set /Herwig/Cuts/MassCut:MaxM 92.2*GeV
"""
## Job options file for Herwig++, LEP1 Z -> qq production
include("GeneratorUtils/StdEvgenSetup.py")
evgenConfig.description = "Herwigpp e+e-"
evgenConfig.generators += ["Herwigpp"]
include("MC12JobOptions/Herwigpp_Base_Fragment.py")
## Construct command set
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(runArgs.ecmEnergy) + hw.base_cmds() + hw.lo_pdf_cmds("cteq6ll.LHpdf")
topAlg.Herwigpp.Commands = cmds.splitlines()
evgenConfig.minevents = 5000
## Add to commands
cmds += """
cd /Herwig/MatrixElements
insert SimpleEE:MatrixElements 0 MEee2gZ2qq
cd /Herwig/Generators
set LEPGenerator:EventHandler:LuminosityFunction:Energy 91.2
set LEPGenerator:EventHandler:CascadeHandler:MPIHandler NULL
set /Herwig/Particles/e-:PDF /Herwig/Partons/LeptonPDF
set /Herwig/Particles/e+:PDF /Herwig/Partons/LeptonPDF
# These next lines are put here just to avoid a meaningless error that p and pbar do not have PDF's
##set /Herwig/Particles/p+:PDF /Herwig/Partons/AtlasPDFsetLO
##set /Herwig/Particles/pbar-:PDF /Herwig/Partons/AtlasPDFsetLO
set /Herwig/Generators/LEPGenerator:EventHandler:BeamA /Herwig/Particles/e+
## Job options file for Herwig++, reading an NLO Les Houches Accord event file
include("GeneratorUtils/StdEvgenSetup.py")
# ## Control log levels
# svcMgr.MessageSvc.OutputLevel = DEBUG
## Add H++ to the job list
from Herwigpp_i.Herwigpp_iConf import Herwigpp
topAlg += Herwigpp()
## Get basic Herwig++ Atlas tune params
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(7000) + hw.base_cmds() + \
hw.nlo_pdf_cmds("CT10.LHgrid", "cteq6ll.LHpdf") + hw.ue_tune_cmds("CTEQ6L1-UE-EE-7000-3") + \
hw.lhef_cmds(nlo=True)
## Add to commands
cmds += """
## Some particle decay modes
# cd /Herwig/Particles
# set b:Stable Stable
# set bbar:Stable Stable
# set tau+:Stable Stable
# set tau-:Stable Stable
# set t:Synchronized Not_synchronized
# set tbar:Synchronized Not_synchronized
# setup t 6 t 174.2 1.5759 14 0 2 3 2 0
# setup tbar -6 tbar 174.2 1.5759 14 0 -2 -3 2 0
## Common job option for gamma gamma processes in Herwig++
## MPI tune is not relevant as the pdf_gammagamma_cmds() function disables MPI
## Contact: Oldrich Kepka
include("MC15JobOptions/Herwigpp_Base_Fragment.py")
from Herwigpp_i import config as hw
cmds = hw.energy_cmds(
runArgs.ecmEnergy) + hw.base_cmds() + hw.pdf_gammagamma_cmds()
cmds += "create ThePEG::O1AlphaS /Herwig/AlphaQCD_O1 O1AlphaS.so\n"
cmds += "set /Herwig/Generators/LHCGenerator:StandardModelParameters:QCD/RunningAlphaS /Herwig/AlphaQCD_O1\n"
genSeq.Herwigpp.Commands += cmds.splitlines()
del cmds
