evgenConfig.description = "bjet particle gun"
evgenConfig.contact = ['*****@*****.**']
evgenConfig.keywords = ["bottom", "jets"]
#### include Common file
###include Pythia8_Base_Fragment.py
## Base config for Pythia8
from Pythia8_i.Pythia8_iConf import Pythia8_i
pythia8_i = Pythia8_i("Pythia8")
genSeq += pythia8_i
evgenConfig.generators += ["Pythia8"]
genSeq.Pythia8.Commands += [
"Main:timesAllowErrors = 500",
"6:m0 = 172.5",
"23:m0 = 91.1876",
"23:mWidth = 2.4952",
"24:m0 = 80.399",
"24:mWidth = 2.085",
"StandardModel:sin2thetaW = 0.23113",
"StandardModel:sin2thetaWbar = 0.23146",
"ParticleDecays:limitTau0 = on",
"ParticleDecays:tau0Max = 10.0",
"PartonLevel:ISR = off",
"PartonLevel:MPI = off"]
## Base config for Pythia8
from Pythia8_i.Pythia8_iConf import Pythia8_i
genSeq += Pythia8_i("Pythia8")
evgenConfig.generators += ["Pythia8"]
genSeq.Pythia8.Commands += [
"Main:timesAllowErrors = 500", "6:m0 = 172.5", "23:m0 = 91.1876",
"23:mWidth = 2.4952", "24:m0 = 80.399", "24:mWidth = 2.085",
"StandardModel:sin2thetaW = 0.23113",
"StandardModel:sin2thetaWbar = 0.23146", "ParticleDecays:limitTau0 = on",
"ParticleDecays:tau0Max = 10.0"
]
#--------------------------------------------------------------
# Event related parameters
#--------------------------------------------------------------
# Number of events to be processed (default is 10)
theApp.EvtMax = 1000
#--------------------------------------------------------------
# Algorithms Private Options
#--------------------------------------------------------------
from AthenaServices.AthenaServicesConf import AtRndmGenSvc
ServiceMgr += AtRndmGenSvc()
from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()
from Pythia8_i.Pythia8_iConf import Pythia8_i
job += Pythia8_i()
include("pythia8.common.py")
include("pythia8.top.py")
#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
myHepMCAnalysis = HepMCAnalysis_i("HepMCAnalysis_i",
file='Pythia8_Process_ttbar_test.root')
#myHepMCAnalysis.JetAnalysis=True
#myHepMCAnalysis.WplusJetAnalysis=True
# Number of events to be processed (default is 10)
theApp.EvtMax = 10
#--------------------------------------------------------------
# Algorithms Private Options
#--------------------------------------------------------------
from AthenaServices.AthenaServicesConf import AtRndmGenSvc
ServiceMgr += AtRndmGenSvc()
ServiceMgr.AtRndmGenSvc.Seeds = ["PYTHIA8 4789899 989240512",
"PYTHIA8_INIT 820021 2347532"]
svcMgr.MessageSvc.OutputLevel = INFO
from AthenaCommon.AlgSequence import AlgSequence
topSequence = AlgSequence()
from Pythia8_i.Pythia8_iConf import Pythia8_i
topSequence+=Pythia8_i()
topSequence.Pythia8_i.CollisionEnergy = 7000
#topSequence.Pythia8_i.Commands += ['WeakSingleBoson:ffbar2gmZ = on']
topSequence.Pythia8_i.LHEFile = "Powheg.ZMu.MC11.events"
topSequence.Pythia8_i.Commands += ['SpaceShower:pTmaxMatch = 2']
topSequence.Pythia8_i.Commands += ['TimeShower:pTmaxMatch = 2']
topSequence.Pythia8_i.UserHook = "QCDVetoedShower"
from TruthExamples.TruthExamplesConf import DumpMC
topSequence += DumpMC()
from G4AtlasApps.SimFlags import simFlags
simFlags.load_atlas_flags()
simFlags.EventFilter.set_Off()
simFlags.MagneticField.set_Off()
simFlags.ForwardDetectors.set_On()
simFlags.ForwardDetectors = 2
from AthenaCommon.AlgSequence import AlgSequence
topSeq = AlgSequence()
from Pythia8_i.Pythia8_iConf import Pythia8_i
topSeq += Pythia8_i()
include("ForwardTransport/Pythia8Config.ALFA.py")
include("G4AtlasApps/G4Atlas.flat.configuration.py")
from AthenaCommon.CfgGetter import getAlgorithm
topSeq += getAlgorithm("G4AtlasAlg", tryDefaultConfigurable=True)
include(
"ForwardTransportSvc/preInclude.ForwardTransportFlags_4.0TeV_0090.00m_nominal_v01.py"
)
include("ForwardTransportSvc/ForwardTransportSvcConfig.ALFA.py")
include(
"ForwardTransportSvc/postInclude.ForwardTransportSvcConfig.FillRootTree.py"
)
include("GeneratorUtils/StdEvgenSetup.py")
svcMgr.MessageSvc.OutputLevel = INFO
svcMgr.AtRndmGenSvc.Seeds = [
"PYTHIA8 4789899 989240512", "PYTHIA8_INIT 820021 2347532"
]
from Pythia8_i.Pythia8_iConf import Pythia8_i
topAlg += Pythia8_i()
topAlg.Pythia8_i.CollisionEnergy = 7000
topAlg.Pythia8_i.Commands += ['HardQCD:all = on']
include("GeneratorUtils/postJO.DumpMC.py")
generate(run_card_loc='run_card.dat',param_card_loc='param_card.dat',mode=2,njobs=8,run_name='Test',proc_dir=process_dir)
stringy = 'madgraph.'+str(runArgs.runNumber)+'.MadGraph_Pythia8_RPV_Gluino_'+Mgo[0:3]
if hasattr(runArgs,'skipEvents'): skip_events=runArgs.skipEvents
skip_events=0
arrange_output(run_name='Test',proc_dir=process_dir,outputDS=stringy+'._00001.events.tar.gz',skip_events=skip_events)
####PYTHIA 8 SHOWERING TO EVGEN############################
from AthenaCommon.AlgSequence import AlgSequence
topAlg = AlgSequence("TopAlg")
from Pythia8_i.Pythia8_iConf import Pythia8_i
Pythia8 = Pythia8_i()
topAlg += Pythia8
Pythia8.Commands += ["Tune:pp=5"]
Pythia8.Commands += [ "Main:timesAllowErrors = 1000"]
Pythia8.Commands += ["6:m0 = 172.5"]
Pythia8.Commands += ["23:m0 = 91.1876"]
Pythia8.Commands += ["23:mWidth = 2.4952"]
Pythia8.Commands += ["24:m0 = 80.399"]
Pythia8.Commands += ["24:mWidth = 2.085"]
Pythia8.Commands += ["ParticleDecays:limitTau0 = on"] # set K_S, Lambda stable
Pythia8.Commands += ["PDF:useLHAPDF=on"]
Pythia8.Commands += ["PDF:LHAPDFset = cteq6ll.LHpdf"]
Pythia8.LHEFile = stringy+'._00001.events'
Pythia8.MaxFailures = 1000000
## Example job option script to run an event generator with FixHepMC filtering
## Author: Andy Buckley <*****@*****.**>
include("GeneratorUtils/StdEvgenSetup.py")
theApp.EvtMax = 100
## Configure and add an event generator to the alg seq
from Pythia8_i.Pythia8_iConf import Pythia8_i
topAlg += Pythia8_i("Pythia8")
topAlg.Pythia8.CollisionEnergy = 7000.0
topAlg.Pythia8.Commands += ['HardQCD:all = on', 'PhaseSpace:pTHatMin = 30.0']
## Now include FixHepMC
include("EvgenProdTools/EnableFixHepMC.py")
#topAlg.FixHepMC.OutputLevel = DEBUG
