* provided "as is" without express or implied warranty. *
**************************************************************************/
-/* $Id$ */
-
//_________________________________________________________________________
// Base class for the clusterization algorithm (pure abstract)
//*--
//*-- Author: Yves Schutz SUBATECH
//////////////////////////////////////////////////////////////////////////////
-// --- ROOT system ---
-#include "TGeometry.h"
-#include "TDirectory.h"
-#include "TFile.h"
-#include "TTree.h"
-
-// --- Standard library ---
-#include <stdlib.h>
+#include <TClonesArray.h>
+#include <TTree.h>
-// --- AliRoot header files ---
-#include "AliRun.h"
#include "AliPHOSClusterizer.h"
-#include "AliHeader.h"
-#include "AliPHOSGetter.h"
-#include "AliPHOSSDigitizer.h"
-#include "AliPHOSDigitizer.h"
+#include "AliPHOSDigit.h"
+#include "AliLog.h"
ClassImp(AliPHOSClusterizer)
+AliPHOSCalibData * AliPHOSClusterizer::fgCalibData = 0 ;
+
//____________________________________________________________________________
- AliPHOSClusterizer::AliPHOSClusterizer():TTask("","")
+AliPHOSClusterizer::AliPHOSClusterizer():
+ fGeom(NULL),
+ fDigitsArr(0),
+ fTreeR(0),
+ fEMCRecPoints(0),
+ fCPVRecPoints(0)
{
// ctor
-
- fEventFolderName = "" ;
+ fDigitsArr = new TClonesArray("AliPHOSDigit",100);
+ fEMCRecPoints = new TObjArray(100) ;
+ fEMCRecPoints ->SetName("EMCRECPOINTS") ;
+ fCPVRecPoints = new TObjArray(100) ;
+ fCPVRecPoints ->SetName("CPVRECPOINTS") ;
}
//____________________________________________________________________________
-AliPHOSClusterizer::AliPHOSClusterizer(const TString alirunFileName, const TString eventFolderName):
- TTask("PHOS"+AliConfig::fgkReconstructionerTaskName, alirunFileName), fEventFolderName(eventFolderName)
+AliPHOSClusterizer::AliPHOSClusterizer(AliPHOSGeometry *geom):
+ fGeom(geom),
+ fDigitsArr(0),
+ fTreeR(0),
+ fEMCRecPoints(0),
+ fCPVRecPoints(0)
{
// ctor
-
+ fDigitsArr = new TClonesArray("AliPHOSDigit",100);
+ fEMCRecPoints = new TObjArray(100) ;
+ fEMCRecPoints ->SetName("EMCRECPOINTS") ;
+ fCPVRecPoints = new TObjArray(100) ;
+ fCPVRecPoints ->SetName("CPVRECPOINTS") ;
}
//____________________________________________________________________________
AliPHOSClusterizer::~AliPHOSClusterizer()
{
// dtor
-
+ if (fDigitsArr) {
+ fDigitsArr->Delete();
+ delete fDigitsArr;
+ }
+ if (fEMCRecPoints) {
+ fEMCRecPoints->Delete();
+ delete fEMCRecPoints;
+ }
+ if (fCPVRecPoints) {
+ fCPVRecPoints->Delete();
+ delete fCPVRecPoints;
+ }
+}
+
+//____________________________________________________________________________
+void AliPHOSClusterizer::SetInput(TTree * digitsTree)
+{
+ // Get the tree with digits and sets
+ // the input array with digits for PHOS
+ TBranch *branch = digitsTree->GetBranch("PHOS");
+ if (!branch) {
+ AliError("can't get the branch with the PHOS digits !");
+ return;
+ }
+ fDigitsArr->Clear();
+ branch->SetAddress(&fDigitsArr);
+ branch->GetEntry(0);
}
+//____________________________________________________________________________
+void AliPHOSClusterizer::SetOutput(TTree * clustersTree)
+{
+ // Set the output clusters tree,
+ // creates the arrays for EMC and CPV,
+ // and set the corresponding branch addresses
+ fTreeR = clustersTree;
+
+ AliDebug(9, "Making array for EMC clusters");
+ Int_t split = 0;
+ Int_t bufsize = 32000;
+ fTreeR->Branch("PHOSEmcRP", "TObjArray", &fEMCRecPoints, bufsize, split);
+ AliDebug(9, "Making array for CPV clusters");
+ fTreeR->Branch("PHOSCpvRP", "TObjArray", &fCPVRecPoints, bufsize, split);
+}