From 9ae15fb41efb10d85d4c9dbf35cd6e66fc11c71b Mon Sep 17 00:00:00 2001 From: mverweij Date: Wed, 26 Feb 2014 11:57:08 +0100 Subject: [PATCH] add minimum seed energy, weight for cluster center of gravity, NxN setting, turn off all options by default we don't want to do here and add function with less params adjust comment --- PWG/EMCAL/AliEmcalClusterMaker.cxx | 2 +- PWG/EMCAL/macros/AddTaskClusterizerFast.C | 39 +++++++++++++++++------ 2 files changed, 30 insertions(+), 11 deletions(-) diff --git a/PWG/EMCAL/AliEmcalClusterMaker.cxx b/PWG/EMCAL/AliEmcalClusterMaker.cxx index e0e5d43e369..7f4d5a1a2a2 100644 --- 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 1b86cf5e0ba..221f1dda435 100644 --- 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); -- 2.39.3