| 42 |
|
for (UInt_t i=0; i<fParticles->GetEntries(); ++i) { |
| 43 |
|
GenParticle* p = fParticles->At(i); |
| 44 |
|
int I = i+1; // Particle index (starts at 1) |
| 45 |
< |
int KF = p->GetPdgId(); // Pdg code |
| 45 |
> |
int KF = p->PdgId(); // Pdg code |
| 46 |
|
double p_x = p->Px(); if (fabs(p_x)<0.0005) p_x = 0.; // Momenta. We only compare the |
| 47 |
|
double p_y = p->Py(); if (fabs(p_y)<0.0005) p_y = 0.; // first 3 digits after the |
| 48 |
|
double p_z = p->Pz(); if (fabs(p_z)<0.0005) p_z = 0.; // decimal. |
| 49 |
|
double E = p->E(); if (fabs(E )<0.0005) E = 0.; // Energy |
| 50 |
|
int mind=0; |
| 51 |
|
if(p->HasMother()) { |
| 52 |
< |
const GenParticle *mother = p->GetMother(); |
| 52 |
> |
const GenParticle *mother = p->Mother(); |
| 53 |
|
if(mother) { |
| 54 |
|
for (UInt_t j=0; j<fParticles->GetEntries(); ++j) { |
| 55 |
|
const GenParticle *test = fParticles->At(j); |
