{
// Import particles from Pythia stack
fParticles->Clear();
+ Int_t nparts = 0;
Int_t ioff = 0;
- Int_t numpart = fPythia->event.size();
+ fNumberOfParticles = fPythia->event.size();
if (fPythia->event[0].id() == 90) {
- numpart--;
ioff = -1;
}
TClonesArray &a = *((TClonesArray*)fParticles);
- for (Int_t i = 1; i <= numpart; i++) {
- new(a[i]) TParticle(
+ for (Int_t i = 0; i < fNumberOfParticles; i++) {
+ if (fPythia->event[i].id() == 90) continue;
+ new(a[nparts]) TParticle(
fPythia->event[i].id(),
fPythia->event[i].isFinal(),
fPythia->event[i].mother1() + ioff,
fPythia->event[i].yProd(), // [mm]
fPythia->event[i].zProd(), // [mm]
fPythia->event[i].tProd()); // [mm/c]
+ nparts++;
}
+ if (ioff == -1) fNumberOfParticles--;
return fParticles;
}