## 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"
## Min bias 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(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 min bias matrix element insert /Herwig/MatrixElements/SimpleQCD:MatrixElements[0] /Herwig/MatrixElements/MEMinBias ## Remove all hard process cuts set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV ## x_min values for the "LO**_UE7000" tune set /Herwig/Cuts/QCDCuts:X1Min 0.055 set /Herwig/Cuts/QCDCuts:X2Min 0.055 ## Turn on the soft MPI modelling, and treat the "hard" process as equivalent to the secondary scatters set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
## 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()
## Herwig++ config for the CTEQ6L1 UE-EE-5 tune series with an 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("CT10.LHgrid", "cteq6ll.LHpdf") \ + hw.ue_tune_cmds("UE-EE-5-CTEQ6L1") genSeq.Herwigpp.Commands += cmds.splitlines() del cmds include("MC15JobOptions/Herwigpp_EvtGen.py") evgenConfig.tune = "CTEQ6L1-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 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"
## 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