if (pdgD == 91 || pdgD == 92) {
Int_t jmin = daughter->GetFirstDaughter() - 1;
Int_t jmax = daughter->GetLastDaughter() - 1;
- for (Int_t j = jmin; j <= jmax; j++)
- ((TParticle *) fParticles.At(j))->SetFirstMother(i+1);
+ for (Int_t jp = jmin; jp <= jmax; jp++)
+ ((TParticle *) fParticles.At(jp))->SetFirstMother(i+1);
} // is string or cluster
} // has daughter
} // heavy quark
Int_t jmin = (fPythia->Version() == 6) ? (daughter->GetFirstDaughter() - 1) : (daughter->GetFirstDaughter());
Int_t jmax = (fPythia->Version() == 6) ? (daughter->GetLastDaughter() - 1) : (daughter->GetLastDaughter());
- for (Int_t j = jmin; j <= jmax; j++)
- ((TParticle *) fParticles.At(j))->SetFirstMother(i+1);
+ for (Int_t jp = jmin; jp <= jmax; jp++)
+ ((TParticle *) fParticles.At(jp))->SetFirstMother(i+1);
} // is string or cluster
} // has daughter
} // heavy quark
// Isotropic decay ????
Double_t cost = 2. * gRandom->Rndm() - 1.;
Double_t sint = TMath::Sqrt(1. - cost * cost);
- Double_t phi = 2. * TMath::Pi() * gRandom->Rndm();
+ Double_t phis = 2. * TMath::Pi() * gRandom->Rndm();
Double_t pz1 = pst * cost;
Double_t pz2 = -pst * cost;
Double_t pt1 = pst * sint;
Double_t pt2 = -pst * sint;
- Double_t px1 = pt1 * TMath::Cos(phi);
- Double_t py1 = pt1 * TMath::Sin(phi);
- Double_t px2 = pt2 * TMath::Cos(phi);
- Double_t py2 = pt2 * TMath::Sin(phi);
+ Double_t px1 = pt1 * TMath::Cos(phis);
+ Double_t py1 = pt1 * TMath::Sin(phis);
+ Double_t px2 = pt2 * TMath::Cos(phis);
+ Double_t py2 = pt2 * TMath::Sin(phis);
fPyjets->P[0][iGlu] = px1;
fPyjets->P[1][iGlu] = py1;
// Isotropic decay ????
Double_t cost = 2. * gRandom->Rndm() - 1.;
Double_t sint = TMath::Sqrt(1. - cost * cost);
- Double_t phi = 2. * TMath::Pi() * gRandom->Rndm();
+ Double_t phis = 2. * TMath::Pi() * gRandom->Rndm();
Double_t pz1 = pst * cost;
Double_t pz2 = -pst * cost;
Double_t pt1 = pst * sint;
Double_t pt2 = -pst * sint;
- Double_t px1 = pt1 * TMath::Cos(phi);
- Double_t py1 = pt1 * TMath::Sin(phi);
- Double_t px2 = pt2 * TMath::Cos(phi);
- Double_t py2 = pt2 * TMath::Sin(phi);
+ Double_t px1 = pt1 * TMath::Cos(phis);
+ Double_t py1 = pt1 * TMath::Sin(phis);
+ Double_t px2 = pt2 * TMath::Cos(phis);
+ Double_t py2 = pt2 * TMath::Sin(phis);
fPyjets->P[0][iGlu] = px1;
fPyjets->P[1][iGlu] = py1;