- else if (kp == EMFSTK.np) {
- mech = kPDeltaRay;
- cppstack->SetTrack(done, parent, pdg,
- px, py, pz, e, vx, vy, vz, tof,
- polx, poly, polz, mech, ntr, weight, is);
-cout << endl << " !!! stupre: ntr=" << ntr << endl;
- EMFSTK.iespak[kp][mkbmx2-1] = ntr;
- } // end of Delta ray
+ else if (kp == EMFSTK.np) {
+ mech = kPDeltaRay;
+ cppstack->SetTrack(done, parent, pdg,
+ px, py, pz, e, vx, vy, vz, tof,
+ polx, poly, polz, mech, ntr, weight, is);
+ cout << endl << " !!! stupre: ntr=" << ntr << " parent=" << parent << endl;
+ EMFSTK.iespak[kp][mkbmx2-1] = ntr;
+ } // end of Delta ray