updated parametrisation for pi0, eta, omega, phi for different energies, allow select...

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

diff --git a/EVGEN/AliGenReaderEcalJets.cxx b/EVGEN/AliGenReaderEcalJets.cxx
index 8646e5a2f5ad4f90a50a8adce30e3b4046442da1..a6875ed2beb852cbf71618857f5fee2d05d33e56 100644 (file)
--- a/EVGEN/AliGenReaderEcalJets.cxx
+++ b/EVGEN/AliGenReaderEcalJets.cxx
@@ -23,19 +23,83 @@
 //
 #include <TFile.h>
 #include <TParticle.h>
+#include <TDatabasePDG.h>
 #include <TTree.h>
-#include <TVirtualMC.h>
 
 #include "AliGenReaderEcalJets.h"
 
 ClassImp(AliGenReaderEcalJets)
 
 
-AliGenReaderEcalJets::AliGenReaderEcalJets() 
+AliGenReaderEcalJets::AliGenReaderEcalJets():
+    fNcurrent(0),
+    fNparticle(0),
+    fNev(0),
+    fNpart(0),
+    fNjet(0),     
+    fNsjet(0),
+    fNpjet(0),
+    fTreeNtuple(0) 
 {
 // Default constructor
-    fNcurrent   = 0;
-    fTreeNtuple = 0;
+    for (Int_t i = 0; i < 200; i++) {
+       if (i < 2) {
+           fX[i]     = 0.;    
+           fXtyp[i]  = 0; 
+       } else if (i < 10) {
+           fJet[i]   = 0.;  
+           fJeta[i]  = 0.; 
+           fJphi[i]  = 0.; 
+           fJset[i]  = 0.; 
+           fJseta[i] = 0.;
+           fJsphi[i] = 0.;
+           fJpet[i]  = 0.; 
+           fJpeta[i] = 0.;
+           fJpphi[i] = 0.;
+       } else {
+           fXpt[i]  = 0.;
+           fXeta[i] = 0.;
+           fXphi[i] = 0.;
+           fXid[i]  = 0; 
+       }
+    }
+}
+
+ AliGenReaderEcalJets::AliGenReaderEcalJets(const AliGenReaderEcalJets &reader)
+     :AliGenReader(reader),
+      fNcurrent(0),
+      fNparticle(0),
+      fNev(0),
+      fNpart(0),
+      fNjet(0),     
+      fNsjet(0),
+      fNpjet(0),
+      fTreeNtuple(0) 
+{
+    for (Int_t i = 0; i < 200; i++) {
+       if (i < 2) {
+           fX[i]     = 0.;    
+           fXtyp[i]  = 0; 
+       } else if (i < 10) {
+           fJet[i]   = 0.;  
+           fJeta[i]  = 0.; 
+           fJphi[i]  = 0.; 
+           fJset[i]  = 0.; 
+           fJseta[i] = 0.;
+           fJsphi[i] = 0.;
+           fJpet[i]  = 0.; 
+           fJpeta[i] = 0.;
+           fJpphi[i] = 0.;
+       } else {
+           fXpt[i]  = 0.;
+           fXeta[i] = 0.;
+           fXphi[i] = 0.;
+           fXid[i]  = 0; 
+       }
+    }
+
+    // Copy Constructor
+    reader.Copy(*this);
 }
 
 void AliGenReaderEcalJets::Init()