PWG/EMCAL/macros/AddTaskEMCALTenderUsingDatasetDef.C

   1 // $Id$
   2
   3 AliAnalysisTaskSE *AddTaskEMCALTenderUsingDatasetDef(
   4   const char *perstr  = "LHC11h",
   5   const char *pass    = 0 /*should not be needed*/
   6 )
   7 {
   8   gROOT->LoadMacro("$ALICE_ROOT/PWG/EMCAL/macros/AddTaskEMCALTender.C");
   9
  10   Bool_t distBC         = kTRUE;   //distance to bad channel
  11   Bool_t recalibClus    = kTRUE;   //recalibrate cluster energy
  12   Bool_t recalcClusPos  = kTRUE;   //recalculate cluster position
  13   Bool_t nonLinearCorr  = kTRUE;   //apply non-linearity
  14   Bool_t remExotic      = kTRUE;   //remove exotic cells
  15   Bool_t fidRegion      = kFALSE;  //apply fiducial cuts
  16   Bool_t calibEnergy    = kTRUE;   //calibrate energy
  17   Bool_t calibTime      = kTRUE;   //calibrate timing
  18   Bool_t remBC          = kTRUE;   //remove bad channels
  19   UInt_t nonLinFunct    = AliEMCALRecoUtils::kBeamTestCorrected;
  20   Bool_t reclusterize   = kFALSE;  //reclusterize
  21   Float_t seedthresh    = 0.100;   //seed threshold
  22   Float_t cellthresh    = 0.050;   //cell threshold
  23   UInt_t clusterizer    = AliEMCALRecParam::kClusterizerv2;
  24   Bool_t trackMatch     = kTRUE;   //track matching
  25   Bool_t updateCellOnly = kFALSE;  //only change if you run your own clusterizer task
  26   Float_t timeMin       = 100e-9;  //minimum time of physical signal in a cell/digit (s)
  27   Float_t timeMax       = 900e-9;  //maximum time of physical signal in a cell/digit (s)
  28   Float_t timeCut       = 900e-9;  //maximum time difference between the digits inside EMC cluster (s)
  29
  30   TString period(perstr);
  31   period.ToLower();
  32   if (period == "lhc12a15e")
  33     nonLinFunct = AliEMCALRecoUtils::kPi0MCv3;
  34   else if (period == "lhc12a15a")
  35     nonLinFunct = AliEMCALRecoUtils::kPi0MCv2;
  36   else if(period == "lhc13b" || period == "lhc13c" || period == "lhc13d" || period == "lhc13e" || period == "lhc13f" || period == "lhc13g" ||
  37           period == "lhc13b4" || period == "lhc13b4_fix" || period == "lhc13b4_plus") {
  38     reclusterize = kTRUE;
  39     seedthresh = 0.3;
  40     cellthresh = 0.05;
  41     if(period == "lhc13b" || period == "lhc13c" || period == "lhc13d" || period == "lhc13e" || period == "lhc13f" || period == "lhc13g") {
  42       timeMin = -50e-9;
  43       timeMax =  50e-9;
  44       timeCut =  1e6;
  45     }
  46     else if(period == "lhc13b4" || period == "lhc13b4_fix" || period == "lhc13b4_plus") {
  47       timeMin = -1;
  48       timeMax =  1e6;
  49       timeCut =  1e6;
  50     }
  51   }
  52
  53   AliAnalysisTaskSE *task = AddTaskEMCALTender(
  54     distBC, recalibClus, recalcClusPos, nonLinearCorr,
  55     remExotic, fidRegion, calibEnergy, calibTime, remBC, nonLinFunct,
  56     reclusterize, seedthresh, cellthresh, clusterizer, trackMatch, updateCellOnly,
  57     timeMin, timeMax, timeCut, pass);
  58
  59   return task;
  60 }