fTimeMin(-1.),fTimeMax(1.),fTimeCut(1.),
fDefaultInit(kFALSE),fToUnfold(kFALSE),
fNumberOfECAClusters(0), fECAClusteringThreshold(0.),
- fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),
+ fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),fRejectBelowThreshold(0),
fClusterUnfolding(NULL)
{
// ctor
fTimeMin(-1.),fTimeMax(1.),fTimeCut(1.),
fDefaultInit(kFALSE),fToUnfold(kFALSE),
fNumberOfECAClusters(0), fECAClusteringThreshold(0.),
- fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),
+ fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),fRejectBelowThreshold(0),
fClusterUnfolding(NULL)
{
// Ctor with the indication of the file where header Tree and digits Tree are stored.
fTimeMin(-1.),fTimeMax(1.),fTimeCut(1.),
fDefaultInit(kFALSE),fToUnfold(kFALSE),
fNumberOfECAClusters(0), fECAClusteringThreshold(0.),
- fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),
+ fECALocMaxCut(0.),fECAW0(0.),fMinECut(0.),fRejectBelowThreshold(0),
fClusterUnfolding(NULL)
{
// ctor, geometry and calibration are initialized elsewhere.
fECAClusteringThreshold = recParam->GetClusteringThreshold();
fECAW0 = recParam->GetW0();
fMinECut = recParam->GetMinECut();
+ fRejectBelowThreshold = recParam->GetRejectBelowThreshold();
fToUnfold = recParam->GetUnfold();
fECALocMaxCut = recParam->GetLocMaxCut();
fTimeCut = recParam->GetTimeCut();
SetNColDiff(recParam->GetNColDiff());
AliDebug(1,Form("Reconstruction parameters: fECAClusteringThreshold=%.3f GeV, fECAW=%.3f, fMinECut=%.3f GeV, "
- "fToUnfold=%d, fECALocMaxCut=%.3f GeV, fTimeCut=%e s,fTimeMin=%e s,fTimeMax=%e s",
- fECAClusteringThreshold,fECAW0,fMinECut,fToUnfold,fECALocMaxCut,fTimeCut, fTimeMin, fTimeMax));
+ "fToUnfold=%d, fECALocMaxCut=%.3f GeV, fTimeCut=%e s,fTimeMin=%e s,fTimeMax=%e s,fRejectBelowThreshold=%d",
+ fECAClusteringThreshold,fECAW0,fMinECut,fToUnfold,fECALocMaxCut,fTimeCut, fTimeMin, fTimeMax,
+ fRejectBelowThreshold));
if (fToUnfold) {
Int_t i=0;
// Read the digits from the input tree
TBranch *branch = digitsTree->GetBranch("EMCAL");
- if (!branch) {
+ if (!branch)
+ {
AliError("can't get the branch with the EMCAL digits !");
return;
}
- if (!fDigitsArr)
- fDigitsArr = new TClonesArray("AliEMCALDigit",100);
+
+ if (!fDigitsArr) fDigitsArr = new TClonesArray("AliEMCALDigit",100);
+ else fDigitsArr->Clear("C"); // avoid leak
+
branch->SetAddress(&fDigitsArr);
branch->GetEntry(0);
}