//////////////////////////////////////////////////////////////////////////////
// --- 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
- fEventFolderName = "" ;
}
//____________________________________________________________________________
-AliEMCALClusterizer::AliEMCALClusterizer(const TString alirunFileName, const TString eventFolderName):
- TTask("EMCAL"+AliConfig::fgkReconstructionerTaskName, alirunFileName), fEventFolderName(eventFolderName)
+AliEMCALClusterizer::~AliEMCALClusterizer()
{
- // ctor
-
+ // dtor
+ if (fDigitsArr) {
+ fDigitsArr->Delete();
+ delete fDigitsArr;
+ }
+ if (fRecPoints) {
+ fRecPoints->Delete();
+ delete fRecPoints;
+ }
}
//____________________________________________________________________________
-AliEMCALClusterizer::~AliEMCALClusterizer()
+void AliEMCALClusterizer::SetInput(TTree *digitsTree)
{
- // dtor
+ // 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;
+ }
+ if (!fDigitsArr)
+ fDigitsArr = new TClonesArray("AliEMCALDigit",100);
+ branch->SetAddress(&fDigitsArr);
+ branch->GetEntry(0);
}
+//____________________________________________________________________________
+void AliEMCALClusterizer::SetOutput(TTree *clustersTree)
+{
+ // Read the digits from the input tree
+ fTreeR = clustersTree;
+
+ 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);
+}