from Herwigpp_i.Herwigpp_iConf import Herwigpp genSeq += Herwigpp() evgenConfig.generators += ["Herwigpp"] genSeq.Herwigpp.Commands += ["set /Herwig/Model:EW/Sin2ThetaW 0.23113"]
## Job option script for POWHEG with Herwig++ ## 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
#-------------------------------------------------------------- # Number of events to be processed (default is 10) theApp.EvtMax = 1000 #-------------------------------------------------------------- # Algorithms Private Options #-------------------------------------------------------------- from AthenaServices.AthenaServicesConf import AtRndmGenSvc ServiceMgr += AtRndmGenSvc() ## ServiceMgr.AtRndmGenSvc.Seeds = ["PYTHIA 4789899 989240512", ## "PYTHIA_INIT 820021 2347532"] # ServiceManager.AtRndmGenSvc.ReadFromFile = true; from AthenaCommon.AlgSequence import AlgSequence job = AlgSequence() from Herwigpp_i.Herwigpp_iConf import Herwigpp job += Herwigpp() #herwigpp = job.Herwigpp cmds = [] include("herwigpp.common.py") include("herwigpp.DiJet.py") job.Herwigpp.Commands += cmds #sj#from TruthExamples.TruthExamplesConf import DumpMC #sj#job += DumpMC() # start my own #from HepMCAnalysis_i.HepMCAnalysis_iConf import HepMCAnalysis from HepMCAnalysis_i.HepMCAnalysis_iConfig import HepMCAnalysis_i