//____________________________________________________________
AliGenEpEmv1::AliGenEpEmv1(const AliGenEpEmv1 & gen)
+ :AliGenMC(gen)
{
// copy constructor
gen.Copy(*this);
fMass = TDatabasePDG::Instance()->GetParticle(11)->Mass();
SetMC(new TEpEmGen());
- fEpEmGen = (TEpEmGen*) fgMCEvGen;
+ fEpEmGen = (TEpEmGen*) fMCEvGen;
fEpEmGen ->Initialize(fYMin,fYMax,fPtMin,fPtMax);
fEvent = 0;
}
id = 11;
if (fDebug == 2)
printf("id=%+3d, p = (%+11.4e,%+11.4e,%+11.4e) GeV\n",id,p[0],p[1],p[2]);
- SetTrack(fTrackIt,-1, id,p,origin,polar,0,kPPrimary,nt,weight);
+ PushTrack(fTrackIt,-1, id,p,origin,polar,0,kPPrimary,nt,weight);
// Produce positron
mt = TMath::Sqrt(ptPositron*ptPositron + fMass*fMass);
id = -11;
if (fDebug == 2)
printf("id=%+3d, p = (%+11.4e,%+11.4e,%+11.4e) GeV\n",id,p[0],p[1],p[2]);
- SetTrack(fTrackIt,-1, id,p,origin,polar,0,kPPrimary,nt,weight);
+ PushTrack(fTrackIt,-1, id,p,origin,polar,0,kPPrimary,nt,weight);
fEvent++;
if (fEvent%1000 == 0) {