- TLorentzVector p;
- gMC->TrackMomentum(p);
-
- TRACKR.spausr[0] = xsco; // x
- TRACKR.spausr[1] = ysco; // y
- TRACKR.spausr[2] = zsco; // z
- TRACKR.spausr[3] = gMC->TrackTime(); // t
- TRACKR.spausr[4] = p[0]; // px
- TRACKR.spausr[5] = p[1]; // py
- TRACKR.spausr[6] = p[2]; // pz
- TRACKR.spausr[7] = p[3]; // e
- TRACKR.spausr[8] = gMC->TrackLength(); // Length
- }
+ TLorentzVector p;
+ gMC->TrackMomentum(p);
+ EMFSTK.espark[npnw][0] = xsco; // x
+ EMFSTK.espark[npnw][1] = ysco; // y
+ EMFSTK.espark[npnw][2] = zsco; // z
+ EMFSTK.espark[npnw][3] = gMC->TrackTime(); // t
+ EMFSTK.espark[npnw][4] = p[0]; // px
+ EMFSTK.espark[npnw][5] = p[1]; // py
+ EMFSTK.espark[npnw][6] = p[2]; // pz
+ EMFSTK.espark[npnw][7] = p[3]; // e
+ EMFSTK.espark[npnw][8] = gMC->TrackLength(); // Length
+ } // Track found in stack
+ } // Loop over emf stack
+ } // Electromagnetic process
+
+