// $Id$
AliAnalysisTaskSE *AddTaskEMCALTender(
- const char *perstr = "LHC11h",//string defining period
Bool_t distBC = kTRUE, //distance to bad channel
Bool_t recalibClus = kTRUE, //recalibrate cluster energy
Bool_t recalcClusPos = kTRUE, //recalculate cluster position
Bool_t remBC = kTRUE, //remove bad channels
UInt_t nonLinFunct = AliEMCALRecoUtils::kBeamTestCorrected,
Bool_t reclusterize = kTRUE, //reclusterize
- Float_t seedthresh = 0.3, //seed threshold
- Float_t cellthresh = 0.05, //cell threshold
+ Float_t seedthresh = 0.100, //seed threshold
+ Float_t cellthresh = 0.050, //cell threshold
UInt_t clusterizer = AliEMCALRecParam::kClusterizerv2,
- Bool_t trackMatch = kFALSE, //track matching
+ Bool_t trackMatch = kTRUE, //track matching
Bool_t updateCellOnly = kFALSE, //only change if you run your own clusterizer task
Float_t timeMin = 100e-9, //minimum time of physical signal in a cell/digit (s)
Float_t timeMax = 900e-9, //maximum time of physical signal in a cell/digit (s)
- Float_t timeCut = 50e-9, //maximum time difference between the digits inside EMC cluster (s)
+ Float_t timeCut = 900e-9, //maximum time difference between the digits inside EMC cluster (s)
const char *pass = 0 //string defining pass (use none if figured out from path)
)
{
AliAnalysisTaskSE *ana = 0;
AliAnalysisDataContainer *coutput1 = 0;
- TString period(perstr);
- period.ToLower();
- if (period == "lhc12a15e")
- nonLinFunct = AliEMCALRecoUtils::kPi0MCv3;
- else if (period == "lhc12a15a")
- nonLinFunct = AliEMCALRecoUtils::kPi0MCv2;
gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/ConfigEmcalTenderSupply.C");
fidRegion, calibEnergy, calibTime, remBC, nonLinFunct, reclusterize, seedthresh,
cellthresh, clusterizer, trackMatch, updateCellOnly, timeMin, timeMax, timeCut);
- if (!period.Contains("lhc10h") && !period.Contains("lhc11h")) {
- ::Info("AddTaskEMCALTender.C","Switching on TrackMatching() since it is not PbPb");
- EMCALSupply->SwitchOnTrackMatch();
- }
-
if (pass)
EMCALSupply->SetPass(pass);
alitender->AddSupply(EMCALSupply);
alitender->SetDefaultCDBStorage("raw://");
ana = alitender;
-
coutput1 = mgr->CreateContainer("emcal_tender_event",
AliESDEvent::Class(),
AliAnalysisManager::kExchangeContainer,
"default_tender");
- }
- else if (evhand->InheritsFrom("AliAODInputHandler")) {
+ } else if (evhand->InheritsFrom("AliAODInputHandler")) {
AliEmcalTenderTask* emcaltender = mgr->GetTopTasks()->FindObject("AliEmcalTenderTask");
if (!emcaltender)
emcaltender = new AliEmcalTenderTask("AliEmcalTenderTask");
AliAODEvent::Class(),
AliAnalysisManager::kExchangeContainer,
"default_tender");
- }
- else {
+ } else {
::Error("AddTaskEMCALTender", "Input event handler not recognized, AOD/ESD expected. Returning...");
return NULL;
}
// $Id$
AliAnalysisTaskSE *AddTaskEMCALTenderUsingDatasetDef(
- const char *dstr = "LHC11h",
+ const char *perstr = "LHC11h",
const char *pass = 0 /*should not be needed*/
)
{
Bool_t remExotic = kTRUE; //remove exotic cells
Bool_t fidRegion = kFALSE; //apply fiducial cuts
Bool_t calibEnergy = kTRUE; //calibrate energy
- Bool_t calibTime = kTRUE, //calibrate timing
- Bool_t remBC = kTRUE, //remove bad channels
+ Bool_t calibTime = kTRUE; //calibrate timing
+ Bool_t remBC = kTRUE; //remove bad channels
UInt_t nonLinFunct = AliEMCALRecoUtils::kBeamTestCorrected;
- Bool_t reclusterize = kFALSE, //reclusterize
- Float_t seedthresh = 0.3, //seed threshold
- Float_t cellthresh = 0.05, //cell threshold
- UInt_t clusterizer = AliEMCALRecParam::kClusterizerv2,
- Bool_t trackMatch = kFALSE, //track matching
- Bool_t updateCellOnly = kFALSE, //only change if you run your own clusterizer task
- Float_t timeMin = 100e-9, //minimum time of physical signal in a cell/digit (s)
- Float_t timeMax = 900e-9, //maximum time of physical signal in a cell/digit (s)
- Float_t timeCut = 50e-9, //maximum time difference between the digits inside EMC cluster (s)
+ Bool_t reclusterize = kFALSE; //reclusterize
+ Float_t seedthresh = 0.100; //seed threshold
+ Float_t cellthresh = 0.050; //cell threshold
+ UInt_t clusterizer = AliEMCALRecParam::kClusterizerv2;
+ Bool_t trackMatch = kTRUE; //track matching
+ Bool_t updateCellOnly = kFALSE; //only change if you run your own clusterizer task
+ Float_t timeMin = 100e-9; //minimum time of physical signal in a cell/digit (s)
+ Float_t timeMax = 900e-9; //maximum time of physical signal in a cell/digit (s)
+ Float_t timeCut = 900e-9; //maximum time difference between the digits inside EMC cluster (s)
+ TString period(perstr);
+ period.ToLower();
+ if (period == "lhc12a15e")
+ nonLinFunct = AliEMCALRecoUtils::kPi0MCv3;
+ else if (period == "lhc12a15a")
+ nonLinFunct = AliEMCALRecoUtils::kPi0MCv2;
-AliAnalysisTaskSE *task = AddTaskEMCALTender(
- "tothink",distBC,recalibClus,recalcClusPos,nonLinearCorr,remExotic,fidRegion,calibEnergy,calibTime,
- remBC, nonLinFunct, reclusterize, seedthresh, cellthresh, clusterizer, trackMatch, updateCellOnly,
- timeMin, timeMax, timeCut, 0);
+ AliAnalysisTaskSE *task = AddTaskEMCALTender(
+ distBC, recalibClus, recalcClusPos, nonLinearCorr,
+ remExotic, fidRegion, calibEnergy, calibTime, remBC, nonLinFunct,
+ reclusterize, seedthresh, cellthresh, clusterizer, trackMatch, updateCellOnly,
+ timeMin, timeMax, timeCut, pass);
- return et;
+ return task;
}
Bool_t remBC = kTRUE, //remove bad channels
UInt_t nonLinFunct = AliEMCALRecoUtils::kBeamTestCorrected,
Bool_t reclusterize = kTRUE, //reclusterize
- Float_t seedthresh = 0.3, //seed threshold
- Float_t cellthresh = 0.05, //cell threshold
+ Float_t seedthresh = 0.100, //seed threshold
+ Float_t cellthresh = 0.050, //cell threshold
UInt_t clusterizer = AliEMCALRecParam::kClusterizerv2,
- Bool_t trackMatch = kFALSE, //track matching
+ Bool_t trackMatch = kTRUE, //track matching
Bool_t updateCellOnly = kFALSE, //only change if you run your own clusterizer task
Float_t timeMin = 100e-9, //minimum time of physical signal in a cell/digit (s)
Float_t timeMax = 900e-9, //maximum time of physical signal in a cell/digit (s)
- Float_t timeCut = 50e-9 //maximum time difference between the digits inside EMC cluster (s)
- )
+ Float_t timeCut = 900e-9 //maximum time difference between the digits inside EMC cluster (s)
+)
{
AliEMCALTenderSupply *EMCALSupply = new AliEMCALTenderSupply("EMCALtender");
EMCALSupply->SetDebugLevel(2);