//////////////////////////////////////////////////////////////////////////////
// --- ROOT system ---
-#include "TGeometry.h"
-#include "TDirectory.h"
-#include "TFile.h"
-#include "TTree.h"
+#include "TClonesArray.h"
+#include "TTree.h"
// --- Standard library ---
-#include <Riostream.h>
-#include <stdlib.h>
+
// --- AliRoot header files ---
-#include "AliRun.h"
#include "AliEMCALClusterizer.h"
-#include "AliHeader.h"
-#include "AliEMCALGetter.h"
-#include "AliEMCALSDigitizer.h"
-#include "AliEMCALDigitizer.h"
+#include "AliLog.h"
ClassImp(AliEMCALClusterizer)
//____________________________________________________________________________
- AliEMCALClusterizer::AliEMCALClusterizer():TTask("","")
+AliEMCALClusterizer::AliEMCALClusterizer():
+ fDigitsArr(NULL),
+ fTreeR(NULL),
+ fRecPoints(NULL)
{
// ctor
- fSplitFile = 0 ;
- fToSplit = kFALSE ;
-
}
//____________________________________________________________________________
-AliEMCALClusterizer::AliEMCALClusterizer(const char* headerFile, const char* name, const Bool_t toSplit):
-TTask(name, headerFile)
+AliEMCALClusterizer::~AliEMCALClusterizer()
{
- // ctor
- fToSplit = toSplit ;
- fSplitFile = 0 ;
+ // dtor
+ if (fDigitsArr) {
+ fDigitsArr->Delete();
+ delete fDigitsArr;
+ }
+ if (fRecPoints) {
+ fRecPoints->Delete();
+ delete fRecPoints;
+ }
+}
+//____________________________________________________________________________
+void AliEMCALClusterizer::SetInput(TTree *digitsTree)
+{
+ // Read the digits from the input tree
+ TBranch *branch = digitsTree->GetBranch("EMCAL");
+ if (!branch) {
+ AliError("can't get the branch with the EMCAL digits !");
+ return;
+ }
+ fDigitsArr = new TClonesArray("AliEMCALDigit",100);
+ branch->SetAddress(&fDigitsArr);
+ branch->GetEntry(0);
}
//____________________________________________________________________________
-AliEMCALClusterizer::~AliEMCALClusterizer()
+void AliEMCALClusterizer::SetOutput(TTree *clustersTree)
{
- // dtor
+ // Read the digits from the input tree
+ fTreeR = clustersTree;
- fSplitFile = 0 ;
- }
-
+ AliDebug(9, "Making array for EMCAL clusters");
+ fRecPoints = new TObjArray(100) ;
+ Int_t split = 0;
+ Int_t bufsize = 32000;
+ fTreeR->Branch("EMCALECARP", "TObjArray", &fRecPoints, bufsize, split);
+}