// --- AliRoot header files ---
#include "AliEMCALReconstructor.h"
+#include "AliCodeTimer.h"
#include "AliESDEvent.h"
#include "AliESDCaloCluster.h"
#include "AliESDCaloCells.h"
AliEMCALReconstructor::~AliEMCALReconstructor()
{
// dtor
+
+ AliCodeTimer::Instance()->Print();
}
//____________________________________________________________________________
// the global tracking.
// Works on the current event.
+ AliCodeTimerAuto("")
+
InitRecParam();
AliEMCALClusterizerv1 clu;
clu.SetInput(digitsTree);
// EMCAL digits.
// Works on a single-event basis
+
+ AliCodeTimerAuto("")
+
rawReader->Reset() ;
TClonesArray *digitsArr = new TClonesArray("AliEMCALDigit",100);
//######################################################
//#########Calculate trigger and set trigger info###########
//######################################################
-
+
+ AliCodeTimerStart(Form("JLK trigger info"));
AliEMCALTrigger tr ;
// tr.SetPatchSize(1);//create 4x4 patches
tr.Trigger();
esd->AddEMCALTriggerPosition(triggerPosition);
esd->AddEMCALTriggerAmplitudes(triggerAmplitudes);
-
+
+ AliCodeTimerStop(Form("JLK trigger info"));
+
//########################################
//##############Fill CaloCells###############
//########################################
-
+
+ AliCodeTimerStart(Form("JLK fill calocells"));
+
TClonesArray *digits = new TClonesArray("AliEMCALDigit",1000);
TBranch *branchdig = digitsTree->GetBranch("EMCAL");
if (!branchdig) {
emcCells.SetNumberOfCells(idignew);
emcCells.Sort();
+ AliCodeTimerStop(Form("JLK fill calocells"));
+
//------------------------------------------------------------
//-----------------CLUSTERS-----------------------------
//------------------------------------------------------------
//######################################################
//Fill list of integers, each one is index of track to which the cluster belongs.
+ AliCodeTimerStart(Form("JLK track matching"));
+
// step 1 - initialize array of matched track indexes
Int_t *matchedTrack = new Int_t[nClusters];
for (Int_t iclus = 0; iclus < nClusters; iclus++)
if(iemcalMatch >= 0) matchedTrack[iemcalMatch] = itrack;
}
+ AliCodeTimerStop(Form("JLK track matching"));
+
//########################################
//##############Fill CaloClusters############
//########################################
+ AliCodeTimerStart(Form("JLK fill caloclusters"));
+
esd->SetNumberOfEMCALClusters(nClusters);
for (Int_t iClust = 0 ; iClust < nClusters ; iClust++) {
const AliEMCALRecPoint * clust = (const AliEMCALRecPoint*)clusters->At(iClust);
pid->SetReconstructor(kTRUE);
pid->RunPID(esd);
delete pid;
+
+ AliCodeTimerStop(Form("JLK fill caloclusters"));
}