/*
$Log$
+Revision 1.14 2000/04/05 08:36:13 morsch
+Check status code of particles in Pythia event
+to avoid double counting as partonic state and final state particle.
+
Revision 1.13 1999/11/09 07:38:48 fca
Changes for compatibility with version 2.23 of ROOT
while(1)
{
fPythia->Pyevnt();
+ fPythia->Lulist(1);
fTrials++;
fPythia->ImportParticles(particles,"All");
Int_t np = particles->GetEntriesFast();
Int_t nc=0;
if (np == 0 ) continue;
if (fProcess != mb) {
- for (Int_t i = 0; i<np; i++) {
+ for (Int_t i = 0; i<np-1; i++) {
TParticle * iparticle = (TParticle *) particles->At(i);
Int_t ks = iparticle->GetStatusCode();
- if (ks==21) continue;
kf = CheckPDGCode(iparticle->GetPdgCode());
+ if (ks==21) continue;
+
fChildWeight=(fPythia->GetBraPart(kf))*fParentWeight;
//
// Parent
} // select particle
} // particle loop
} else {
- for (Int_t i = 0; i<np; i++) {
+ for (Int_t i = 0; i<np-1; i++) {
TParticle * iparticle = (TParticle *) particles->At(i);
kf = CheckPDGCode(iparticle->GetPdgCode());
Int_t ks = iparticle->GetStatusCode();
+
if (ks==1 && kf!=0 && KinematicSelection(iparticle)) {
nc++;
//