evgenConfig.description = "Test of Pythia 8 + Photos++ in Z pT"
evgenConfig.keywords = ["QED", "Z"]
include("MC12JobOptions/Pythia8_AU2_CTEQ6L1_Common.py")
include("MC12JobOptions/PhotosPythia8_Fragment.py")
# Uncomment to turn exponentiation mode off
# Most physics channels should leave the default mode to on (i.e. leave this commented out)
#topAlg.Photospp.ExponentiationMode=False
topAlg.Pythia8.Commands += ["WeakSingleBoson:ffbar2gmZ=on"]
topAlg.Pythia8.Commands += ["23:onMode=off", "23:onIfAny=11"]
evgenConfig.minevents = 25000
from Rivet_i.Rivet_iConf import Rivet_i
rivet = Rivet_i("Rivet")
rivet.Analyses = ["ATLAS_2011_S9131140"]
topAlg += rivet
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from GaudiSvc.GaudiSvcConf import THistSvc
svcMgr += THistSvc()
svcMgr.THistSvc.Output = ["Rivet DATAFILE='Rivet.root' OPT='RECREATE'"]
## 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
cmds += """
## Use preweighting rather than jet slices to cover wide pT range
create ThePEG::ReweightMinPT /Herwig/Generators/preweight ReweightMinPT.so
set /Herwig/Generators/preweight:Power 6
insert /Herwig/Generators/LHCGenerator:EventHandler:SubProcessHandlers[0]:Preweights[0] /Herwig/Generators/preweight
set /Herwig/Generators/LHCGenerator:EventHandler:Weighted On
"""
topAlg.Herwigpp.Commands = cmds.splitlines()
## Add Rivet UE analysis
from Rivet_i.Rivet_iConf import Rivet_i
rivet = Rivet_i()
rivet.Analyses = ['MC_JETS', 'D0_2004_S5992206']
topAlg += rivet
