]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TEPEMGEN/AliGenEpEmv1.cxx
Completely reengineered version of CMake build system (Johny)
[u/mrichter/AliRoot.git] / TEPEMGEN / AliGenEpEmv1.cxx
index f4a7f6addb6a90a3f5c2542a99e81e85b249bd95..3aa0f17c204832424062b11b389da8de04fc8165 100644 (file)
  *                                                                        *
  **************************************************************************/
 
-/*
-$Log$
-Revision 1.1  2002/11/06 10:26:45  hristov
-Event generator for e+e- pair production (Yu.Kharlov)
-
-*/
+/* $Id$ */
 
 // Event generator of single e+e- pair production in ultraperipheral PbPb collisions
 // at 5.5 TeV/nucleon.
@@ -84,20 +79,18 @@ ClassImp(AliGenEpEmv1)
 
 //------------------------------------------------------------
 
-AliGenEpEmv1::AliGenEpEmv1()
+AliGenEpEmv1::AliGenEpEmv1():
+  fMass(0),
+  fEpEmGen(0),
+  fDebug(0),
+  fEvent(0),
+  fHeader(AliGenEventHeader())
 {
   // Default constructor
   // Avoid zero pt
   if (fPtMin == 0) fPtMin = 1.E-04;
 }
 
-//____________________________________________________________
-AliGenEpEmv1::AliGenEpEmv1(const AliGenEpEmv1 & gen)
-{
-  // copy constructor
-  gen.Copy(*this);
-}
-
 //____________________________________________________________
 AliGenEpEmv1::~AliGenEpEmv1()
 {
@@ -114,7 +107,7 @@ void AliGenEpEmv1::Init()
   fMass = TDatabasePDG::Instance()->GetParticle(11)->Mass();
 
   SetMC(new TEpEmGen());
-  fEpEmGen = (TEpEmGen*) fgMCEvGen;
+  fEpEmGen = (TEpEmGen*) fMCEvGen;
   fEpEmGen ->Initialize(fYMin,fYMax,fPtMin,fPtMax);
   fEvent = 0;
 }
@@ -175,7 +168,7 @@ void AliGenEpEmv1::Generate()
   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);
@@ -185,12 +178,14 @@ void AliGenEpEmv1::Generate()
   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) {
     printf("=====> AliGenEpEmv1::Generate(): \n   Event %d, sigma=%f +- %f kb\n",
           fEvent,fEpEmGen->GetXsection(),fEpEmGen->GetDsection());
   }
+  fHeader.SetEventWeight(weight);
+  AddHeader(&fHeader);
 }