X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EVGEN%2FAliGenReaderEcalJets.cxx;h=a6875ed2beb852cbf71618857f5fee2d05d33e56;hb=d78f01ff1adf84d01f94ac918af31932d610c773;hp=9ae491b0b3329f100ecefb99b7ebe61cdf821c17;hpb=88cb7938ca21d4a80991d4e7aa564008c29340f7;p=u%2Fmrichter%2FAliRoot.git

diff --git a/EVGEN/AliGenReaderEcalJets.cxx b/EVGEN/AliGenReaderEcalJets.cxx
index 9ae491b0b33..a6875ed2beb 100644
--- a/EVGEN/AliGenReaderEcalJets.cxx
+++ b/EVGEN/AliGenReaderEcalJets.cxx
@@ -14,22 +14,92 @@
  **************************************************************************/
 
 /* $Id$ */
-
+//
+// Realisation of AliGenReader to be used with AliGenExtFile
+// It reads Pythia Jet 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 "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() 
@@ -94,8 +164,9 @@ Int_t AliGenReaderEcalJets::NextEvent()
 
 TParticle* AliGenReaderEcalJets::NextParticle() 
 {
-    Float_t p[4];
 // Read the next particle
+
+    Float_t p[4];
     Int_t    ipart  = fXid[fNparticle];
     Float_t  pt     = fXpt[fNparticle];
     Float_t  eta    = fXeta[fNparticle];
@@ -121,9 +192,19 @@ TParticle* AliGenReaderEcalJets::NextParticle()
 AliGenReaderEcalJets& AliGenReaderEcalJets::operator=(const  AliGenReaderEcalJets& rhs)
 {
 // Assignment operator
-    return *this;
+    rhs.Copy(*this);
+    return (*this);
 }
 
+void AliGenReaderEcalJets::Copy(TObject&) const
+{
+    //
+    // Copy 
+    //
+    Fatal("Copy","Not implemented!\n");
+}
+
+