add minimum seed energy, weight for cluster center of gravity, NxN setting, turn...

adjust comment

diff --git a/PWG/EMCAL/AliEmcalClusterMaker.cxx b/PWG/EMCAL/AliEmcalClusterMaker.cxx
index e0e5d43e369ef0daefa31295b58f7d2ffe449142..7f4d5a1a2a2f83ca270a730479411dfea4129296 100644 (file)
--- a/PWG/EMCAL/AliEmcalClusterMaker.cxx
+++ b/PWG/EMCAL/AliEmcalClusterMaker.cxx
@@ -122,7 +122,7 @@ void AliEmcalClusterMaker::ExecOnce()
 //________________________________________________________________________
 Bool_t AliEmcalClusterMaker::Run() 
 {
-  // Run the hadronic correction
+  // Run the cluster maker
 
   // delete output
   fOutClusters->Delete();

diff --git a/PWG/EMCAL/macros/AddTaskClusterizerFast.C b/PWG/EMCAL/macros/AddTaskClusterizerFast.C
index 1b86cf5e0ba1e8b4a80815662b19c39b695a9313..221f1dda4355146e97424d1e203828f5b24aa034 100644 (file)
--- a/PWG/EMCAL/macros/AddTaskClusterizerFast.C
+++ b/PWG/EMCAL/macros/AddTaskClusterizerFast.C
@@ -4,15 +4,30 @@ AliAnalysisTaskEMCALClusterizeFast* AddTaskClusterizerFast(
   const char* clusName  = "",
   UInt_t inputCellType  = AliAnalysisTaskEMCALClusterizeFast::kFEEData,
   UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2,
-  Bool_t nonLinearCorr  = kFALSE, 
-  UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected,
-  Bool_t calcDistToBC   = kFALSE, 
-  Bool_t remBC          = kFALSE,
-  Bool_t remExotic      = kFALSE,
-  Bool_t fidRegion      = kFALSE,
-  Bool_t updateCells    = kFALSE,
-  Bool_t trackMatch     = kFALSE,
-  Double_t minE         = 0.05
+  Double_t seedE        = 0.1,
+  Double_t cellE        = 0.05,
+  Bool_t calcDistToBC   = kFALSE) {
+
+  return AddTaskClusterizerFast(taskname,cellsName,clusName,inputCellType,clusterizer,kFALSE,0,calcDistToBC,kFALSE,kFALSE,kFALSE,kFALSE,kFALSE,cellE,seedE);
+
+}
+
+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
+  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
 ) 
 {
   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
@@ -25,7 +40,11 @@ AliAnalysisTaskEMCALClusterizeFast* AddTaskClusterizerFast(
 
   AliEMCALRecParam *recparam = task->GetRecParam();
   recparam->SetClusterizerFlag(clusterizer);
-  recparam->SetMinECut(minE);
+  recparam->SetMinECut(cellE);
+  recparam->SetClusteringThreshold(seedE);
+  recparam->SetW0(4.5);
+  if (clusterizer == AliEMCALRecParam::kClusterizerNxN) //MV: copied from tender. please check
+    recparam->SetNxM(3,3);
 
   AliEMCALRecoUtils *recoUtils = new AliEMCALRecoUtils();
   recoUtils->SetNonLinearityFunction(nonLinFunct);