add nested decay (i.e. allow daughters to decay to any requested depth)

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

diff --git a/EVGEN/AliGenReaderEcalHijing.cxx b/EVGEN/AliGenReaderEcalHijing.cxx
index c00cb9cc5d32838b2f4b6b25c6d5c843a8946eac..c9482220a610cb04f4cdcf9f183615544fd805ed 100644 (file)
--- a/EVGEN/AliGenReaderEcalHijing.cxx
+++ b/EVGEN/AliGenReaderEcalHijing.cxx
@@ -23,11 +23,12 @@
 //
 #include <TFile.h>
 #include <TParticle.h>
+#include <TDatabasePDG.h>
 #include <TTree.h>
 
 #include "AliGenReaderEcalHijing.h"
 
-ClassImp(AliGenReaderEcalHijing);
+ClassImp(AliGenReaderEcalHijing)
 
 AliGenReaderEcalHijing::AliGenReaderEcalHijing():
     fNcurrent(0),
@@ -37,7 +38,14 @@ AliGenReaderEcalHijing::AliGenReaderEcalHijing():
     fNahij(0),
     fNphij(0)
 {
-// Default constructor
+  // 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(const AliGenReaderEcalHijing &reader):
@@ -50,6 +58,13 @@ AliGenReaderEcalHijing::AliGenReaderEcalHijing(const AliGenReaderEcalHijing &rea
     fNphij(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);
 }