ProcessLevel:all = off
pythia.next()
PartonLevel
pythia.forceTimeShower(...)
HadronLevel
pythia.event
pythia.init()
Pythia pythia; // Declare generator. Event& event = pythia.event // Convenient shorthand. pythia.readString("ProcessLevel:all = off"); // The trick! pythia.init(); // Initialization. for (int iEvent = 0; iEvent < nEvent; ++iEvent) { // Insert filling of event here! pythia.next(); // Do the hadron level. }
event
event.reset()
event.append(...)
append
Vec4
append( id, status, col, acol, p, m) append( id, status, col, acol, px, py, pz, e, m) append( id, status, mother1, mother2, daughter1, daughter2, col, acol, p, m) append( id, status, mother1, mother2, daughter1, daughter2, col, acol, px, py, pz, e, m)
scale
pythia.forceTimeShower()
pol
id
col
acol
Check:event = on
pythia.forceTimeShower( int iBeg, int iEnd, double pTmax, int nBranchMax = 0)
iBeg
iEnd
pTmax
nBranchMax
pythia.event[i].scale(...)
main21.cc
<init>
</init>
<LesHouchesEvents version="1.0">
<event>....</event>
Beams:frameType = 4
Beams:LHEF
LHAup
Beams:frameType = 5
ProcessLevel:resonanceDecays = off
PartonLevel:FSR = off
PartonLevel:FSRinResonances = off
pythia.forceHadronLevel()
HadronLevel:all
HadronLevel:all = off
Event savedEvent = pythia.event
pythia.event = savedEvent
pythia.moreDecays()
pythia.particleData.mayDecay( id, false/true)
event.saveSize()
event.restoreSize()
event[i].statusPos()
main15.cc