+ Float_t originP[3] = {0., 0., 0.};
+ originP[2] = zPrimary;
+
+ Float_t pP[3] = {0., 0., 0.};
+ Int_t ntP;
+
+ if (fSide == -1) {
+ originP[2] = -zPrimary;
+ origin[2] = -origin[2];
+ p[2] = -p[2];
+ }
+
+ gAlice->SetTrack(0,-1,kProton,pP,originP,polar,0,kPNoProcess,ntP);
+ gAlice->KeepTrack(ntP);
+
+