* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-*/
-
+/* $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 <TTree.h>
#include <TParticle.h>
+#include <TDatabasePDG.h>
+#include <TTree.h>
#include "AliGenReaderEcalJets.h"
-#include "AliMC.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;
+}
+
+ AliGenReaderEcalJets::AliGenReaderEcalJets(const AliGenReaderEcalJets &reader)
+ :AliGenReader(reader),
+ fNcurrent(0),
+ fNparticle(0),
+ fNev(0),
+ fNpart(0),
+ fNjet(0),
+ fNsjet(0),
+ fNpjet(0),
+ fTreeNtuple(0)
+{
+ // Copy Constructor
+ reader.Copy(*this);
}
void AliGenReaderEcalJets::Init()
Int_t AliGenReaderEcalJets::NextEvent()
{
// Read the next event
- Int_t nTracks, nread;
+ Int_t nTracks=0, nread=0;
TFile* pFile = fTreeNtuple->GetCurrentFile();
pFile->cd();
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];
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");
+}
+
+