// --- ROOT system ---
-
-
+#include "TBranch.h"
+#include "TClonesArray.h"
+#include "TTree.h"
// --- Standard library ---
-#include <iostream.h>
-#include <stdlib.h>
-
// --- AliRoot header files ---
-
+#include "AliConfig.h"
+#include "AliLog.h"
#include "AliPHOSPID.h"
ClassImp(AliPHOSPID)
//____________________________________________________________________________
- AliPHOSPID::AliPHOSPID():TTask("","")
+AliPHOSPID::AliPHOSPID():
+ TObject(),
+ fGeom(NULL),
+ fESD(0x0),
+ fEMCRecPoints(NULL),
+ fCPVRecPoints(NULL),
+ fTrackSegments(NULL),
+ fRecParticles(NULL),
+ fEnergyCorrectionOn(kTRUE)
{
// ctor
}
+
//____________________________________________________________________________
-AliPHOSPID::AliPHOSPID(const char* headerFile, const char * name ):TTask(name, headerFile)
+AliPHOSPID::AliPHOSPID(AliPHOSGeometry *geom):
+ TObject(),
+ fGeom(geom),
+ fESD(0x0),
+ fEMCRecPoints(NULL),
+ fCPVRecPoints(NULL),
+ fTrackSegments(NULL),
+ fRecParticles(NULL),
+ fEnergyCorrectionOn(kTRUE)
{
// ctor
+ fEMCRecPoints = new TObjArray(100) ;
+ fCPVRecPoints = new TObjArray(100) ;
+ fRecParticles = new TClonesArray("AliPHOSRecParticle",100) ;
+ fRecParticles->SetName("RECPARTICLES");
+
+}
+
+//____________________________________________________________________________
+AliPHOSPID::AliPHOSPID(const AliPHOSPID & pid) :
+ TObject(pid),
+ fGeom(pid.fGeom),
+ fESD(pid.fESD),
+ fEMCRecPoints(pid.fEMCRecPoints),
+ fCPVRecPoints(pid.fCPVRecPoints),
+ fTrackSegments(pid.fTrackSegments),
+ fRecParticles(pid.fRecParticles),
+ fEnergyCorrectionOn(pid.fEnergyCorrectionOn)
+{
+ // Copy constructor
}
//____________________________________________________________________________
AliPHOSPID::~AliPHOSPID()
{
// dtor
+ if (fEMCRecPoints) {
+ fEMCRecPoints->Delete();
+ delete fEMCRecPoints;
+ }
+ if (fCPVRecPoints) {
+ fCPVRecPoints->Delete();
+ delete fCPVRecPoints;
+ }
+ if (fRecParticles) {
+ fRecParticles->Delete();
+ delete fRecParticles;
+ }
+}
+
+//____________________________________________________________________________
+void AliPHOSPID::SetInput(TTree *clustersTree, TClonesArray *trackSegments)
+{
+ // Read the clusters tree and creates the
+ // arrays with the EMC and CPV
+ // clusters.
+ // and set the corresponding branch addresses
+
+ fTrackSegments = trackSegments;
+
+ TBranch *emcbranch = clustersTree->GetBranch("PHOSEmcRP");
+ if (!emcbranch) {
+ AliError("can't get the branch with the PHOS EMC clusters !");
+ return;
+ }
+ emcbranch->SetAddress(&fEMCRecPoints);
+ fEMCRecPoints->Delete();
+ emcbranch->GetEntry(0);
+
+ TBranch *cpvbranch = clustersTree->GetBranch("PHOSCpvRP");
+ if (!cpvbranch) {
+ AliError("can't get the branch with the PHOS CPV clusters !");
+ return;
+ }
+ cpvbranch->SetAddress(&fCPVRecPoints);
+ fCPVRecPoints->Delete();
+ cpvbranch->GetEntry(0);
}