updates configuration ClusterizerFast

diff --git a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
index 489f1abb50fc3cdcaa4e745e2d20cc82c8758c98..7dfd45508f0ef0790ce6e180b1179b54265bb702 100644 (file)
--- a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
+++ b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx
@@ -89,7 +89,7 @@ AliAnalysisTaskEMCALClusterizeFast::AliAnalysisTaskEMCALClusterizeFast() :
   fCaloCellsName(),  
   fCaloClustersName("newCaloClusters"),
   fDoUpdateCells(kFALSE),
-  fDoClusterize(kFALSE),
+  fDoClusterize(kTRUE),
   fClusterBadChannelCheck(kFALSE),
   fRejectExoticClusters(kFALSE),
   fRejectExoticCells(kFALSE),
@@ -141,7 +141,7 @@ AliAnalysisTaskEMCALClusterizeFast::AliAnalysisTaskEMCALClusterizeFast(const cha
   fCaloCellsName(),  
   fCaloClustersName("newCaloClusters"),
   fDoUpdateCells(kFALSE),
-  fDoClusterize(kFALSE),
+  fDoClusterize(kTRUE),
   fClusterBadChannelCheck(kFALSE),
   fRejectExoticClusters(kFALSE),
   fRejectExoticCells(kFALSE),
@@ -338,7 +338,6 @@ void AliAnalysisTaskEMCALClusterizeFast::FillDigitsArray()
   // Fill digits array
 
   fDigitsArr->Clear("C");
-
   switch (fInputCellType) {
 
   case kFEEData :
@@ -516,7 +515,7 @@ Bool_t AliAnalysisTaskEMCALClusterizeFast::AcceptCell(Int_t cellNumber) {
     fRecoUtils->SwitchOnRejectExoticCell();//switch on and off
     Int_t bunchCrossNo = InputEvent()->GetBunchCrossNumber();
     Bool_t isEx = fRecoUtils->IsExoticCell(cellNumber, fCaloCells, bunchCrossNo);
-    accept = !isEx;
+    if(isEx) accept = kFALSE;
     fRecoUtils->SwitchOffRejectExoticCell();//switch on and off
   }
   return accept;

diff --git a/PWG/EMCAL/macros/AddTaskClusterizerFast.C b/PWG/EMCAL/macros/AddTaskClusterizerFast.C
index ee01d77be92591e042f0698e4c0c709eb252ed37..1e5d9b73a028186a9aecb6bcdcd772f675bcdc02 100644 (file)
--- a/PWG/EMCAL/macros/AddTaskClusterizerFast.C
+++ b/PWG/EMCAL/macros/AddTaskClusterizerFast.C
@@ -2,40 +2,16 @@ AliAnalysisTaskEMCALClusterizeFast* AddTaskClusterizerFast(
   const char* taskname  = "ClusterizerFast",
   const char* cellsName = "",
   const char* clusName  = "",
-  UInt_t inputCellType  = AliAnalysisTaskEMCALClusterizeFast::kFEEData,
   UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Double_t seedE        = 0.1,
   Double_t cellE        = 0.05,
+  Double_t seedE        = 0.1,
   const Float_t timeMin = -1,      //minimum time of physical signal in a cell/digit (s)
   const Float_t timeMax = +1,      //maximum time of physical signal in a cell/digit (s)
   const Float_t timeCut =  1,      //maximum time difference between the digits inside EMC cluster (s)
-  Bool_t calcDistToBC   = kFALSE) {
+  Bool_t remExoticCell  = kTRUE,
+  Bool_t calcDistToBC   = kFALSE,
+  UInt_t inputCellType  = AliAnalysisTaskEMCALClusterizeFast::kFEEData) {
 
-  return AddTaskClusterizerFast(taskname,cellsName,clusName,inputCellType,clusterizer,kFALSE,0,calcDistToBC,kFALSE,kFALSE,kFALSE,kFALSE,kFALSE,cellE,seedE,timeMin,timeMax,timeCut);
-
-}
-
-AliAnalysisTaskEMCALClusterizeFast* AddTaskClusterizerFast(
-  const char* taskname  = "ClusterizerFast",
-  const char* cellsName = "",
-  const char* clusName  = "",
-  UInt_t inputCellType  = AliAnalysisTaskEMCALClusterizeFast::kFEEData,
-  UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Bool_t nonLinearCorr  = kFALSE,                                   //MV: remove
-  UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected,    //MV: remove
-  Bool_t calcDistToBC   = kFALSE,                                   //MV: keep?
-  Bool_t remBC          = kFALSE,                                   //MV: remove
-  Bool_t remExotic      = kFALSE,                                   //MV: remove but replace with exotic cell removal
-  Bool_t fidRegion      = kFALSE,                                   //MV: keep?
-  Bool_t updateCells    = kFALSE,                                   //MV: remove
-  Bool_t trackMatch     = kFALSE,                                   //MV: remove
-  Double_t cellE        = 0.05,                                    //minimum energy for cells to be clustered
-  Double_t seedE        = 0.1,                                      //minimum energy for seed
-  const Float_t timeMin = -1,      //minimum time of physical signal in a cell/digit (s)
-  const Float_t timeMax = +1,      //maximum time of physical signal in a cell/digit (s)
-  const Float_t timeCut =  1      //maximum time difference between the digits inside EMC cluster (s)
-) 
-{
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
   if (!mgr) {
     ::Error("AddTaskClusterizerFast", "No analysis manager found.");
@@ -57,19 +33,19 @@ AliAnalysisTaskEMCALClusterizeFast* AddTaskClusterizerFast(
     recparam->SetNxM(3,3);
 
   AliEMCALRecoUtils *recoUtils = new AliEMCALRecoUtils();
-  recoUtils->SetNonLinearityFunction(nonLinFunct);
+  recoUtils->SetNonLinearityFunction(0);
   task->SetEMCALRecoUtils(recoUtils);
 
   task->SetAttachClusters(kTRUE);
   task->SetCaloClustersName(clusName);
   task->SetCaloCellsName(cellsName);
   task->SetInputCellType(inputCellType);
+  Printf("inputCellType: %d",inputCellType);
+
+  task->SetClusterize(kTRUE);
 
-  task->SetUpdateCells(updateCells);
-  task->SetClusterBadChannelCheck(remBC);
-  task->SetRejectExoticClusters(remExotic);
-  task->SetFiducial(fidRegion);
-  task->SetDoNonLinearity(nonLinearCorr);
+  task->SetClusterBadChannelCheck(kTRUE);
+  task->SetRejectExoticCells(remExoticCell);
   task->SetRecalDistToBadChannels(calcDistToBC);
 
   mgr->AddTask(task);