Updated fragmentation model

[u/mrichter/AliRoot.git] / EVGEN / AliGenReaderEcalHijing.cxx

diff --git a/EVGEN/AliGenReaderEcalHijing.cxx b/EVGEN/AliGenReaderEcalHijing.cxx
index 7864818fb52f6077e6e7228d8276aa1dad8d2986..c9482220a610cb04f4cdcf9f183615544fd805ed 100644 (file)
--- a/EVGEN/AliGenReaderEcalHijing.cxx
+++ b/EVGEN/AliGenReaderEcalHijing.cxx
@@ -14,22 +14,58 @@
  **************************************************************************/
 
 /* $Id$ */
-
+//
+// Realisation of AliGenReader to be used with AliGenExtFile
+// It reads Hijing events from a ntuple like event structure.
+// The event format is defined in Init()
+// NextEvent() is used to loop over events and NextParticle() to loop over particles.  
+// Author: andreas.morsch@cern.ch
+//
 #include <TFile.h>
 #include <TParticle.h>
+#include <TDatabasePDG.h>
 #include <TTree.h>
-#include <TVirtualMC.h>
 
 #include "AliGenReaderEcalHijing.h"
 
 ClassImp(AliGenReaderEcalHijing)
 
+AliGenReaderEcalHijing::AliGenReaderEcalHijing():
+    fNcurrent(0),
+    fNparticle(0),
+    fTreeNtuple(0),
+    fNjatt(0),
+    fNahij(0),
+    fNphij(0)
+{
+  // Default constructor
+  for (Int_t i = 0; i < 10000; i++) {
+    fKhij[i]  = 0;     
+    fPxhij[i] = 0.;    
+    fPyhij[i] = 0.;    
+    fPzhij[i] = 0.;
+    fEhij[i]  = 0.;     
+  }
+}
 
-AliGenReaderEcalHijing::AliGenReaderEcalHijing() 
+AliGenReaderEcalHijing::AliGenReaderEcalHijing(const AliGenReaderEcalHijing &reader):
+    AliGenReader(reader),
+    fNcurrent(0),
+    fNparticle(0),
+    fTreeNtuple(0),
+    fNjatt(0),
+    fNahij(0),
+    fNphij(0)
 {
-// Default constructor
-    fNcurrent   = 0;
-    fTreeNtuple = 0;
+    // Copy constructor
+  for (Int_t i = 0; i < 10000; i++) {
+    fKhij[i]  = 0;     
+    fPxhij[i] = 0.;    
+    fPyhij[i] = 0.;    
+    fPzhij[i] = 0.;
+    fEhij[i]  = 0.;     
+  }
+    reader.Copy(*this);
 }
 
 void AliGenReaderEcalHijing::Init() 
@@ -82,8 +118,9 @@ Int_t AliGenReaderEcalHijing::NextEvent()
 
 TParticle* AliGenReaderEcalHijing::NextParticle() 
 {
-    Float_t p[4];
 // Read the next particle
+
+    Float_t p[4];
     Int_t ipart = fKhij[fNparticle];
     p[0] = fPxhij[fNparticle];
     p[1] = fPyhij[fNparticle];      
@@ -112,7 +149,7 @@ AliGenReaderEcalHijing& AliGenReaderEcalHijing::operator=(const  AliGenReaderEca
     return (*this);
 }
 
-void AliGenReaderEcalHijing::Copy(AliGenReaderEcalHijing&) const
+void AliGenReaderEcalHijing::Copy(TObject&) const
 {
     //
     // Copy