PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C

   1 AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst,Bool_t MassWidthCut,Bool_t MassCal,Double_t asspTCut,Double_t angleCut){\r
   2   //\r
   3   // HFE standard task configuration\r
   4   //\r
   5 \r
   6   Bool_t kAnalyseTaggedTracks = kTRUE;\r
   7   \r
   8   AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");\r
   9   hfecuts->CreateStandardCuts();\r
  10   hfecuts->SetMinNClustersTPC(100);\r
  11   hfecuts->SetMinRatioTPCclusters(0.6);\r
  12   hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r
  13   hfecuts->SetMinNClustersITS(3);\r
  14   hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny);\r
  15   hfecuts->SetCheckITSLayerStatus(kFALSE);\r
  16   hfecuts->SetVertexRange(10.);\r
  17   hfecuts->SetTOFPIDStep(kFALSE);\r
  18   hfecuts->SetPtRange(2, 50);\r
  19   hfecuts->SetMaxImpactParam(3.,3.);\r
  20   \r
  21   AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");\r
  22   printf("task ------------------------ %p\n ", task);\r
  23   task->SetHFECuts(hfecuts);\r
  24   task->SetMassConstraint(MassConst);\r
  25   task->SetMassWidthCut(MassWidthCut);\r
  26   Double_t masscut = 0.05;\r
  27   if(!MassConst)masscut = 0.1;\r
  28   task->SetInvariantMassCut(masscut);\r
  29   task->SetOpeningAngleCut(angleCut);\r
  30   task->SetMimpTassCut(asspTCut);\r
  31   task->SetMassCalMethod(MassCal);\r
  32 \r
  33   // Define PID\r
  34   AliHFEpid *pid = task->GetPID();\r
  35   if(useMC) pid->SetHasMCData(kTRUE);\r
  36   pid->AddDetector("TPC", 0);\r
  37   pid->AddDetector("EMCAL", 1);\r
  38 \r
  39   Double_t params[4];\r
  40   char *cutmodel;\r
  41   if(useMC){\r
  42           // Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r
  43           // for high momentum it is consistent with a flat -0.94\r
  44           //cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
  45           //Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
  46           //for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
  47           cutmodel = "pol0";\r
  48           params[0] = -1.0; //sigma min\r
  49   } else {\r
  50           // Data is consistent with a flat 0.12\r
  51           cutmodel = "pol0";\r
  52           //params[0] = -0.0015;\r
  53           //params[0] = -3.0;\r
  54           //params[0] = -0.05; //sigma min\r
  55           params[0] = -1.0; //sigma min\r
  56   }\r
  57   //pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
  58   for(Int_t a=0;a<11;a++)pid->ConfigureTPCcentralityCut(a,cutmodel,params,3.0);\r
  59   \r
  60 \r
  61   // change E/p cuts\r
  62   AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r
  63   emcpid->SetEoPMax(1.3);\r
  64   emcpid->SetEoPMim(0.9);\r
  65 \r
  66   printf("*************************************\n");\r
  67   printf("Configuring standard Task:\n");\r
  68 //  task->PrintStatus();\r
  69   pid->PrintStatus();\r
  70   printf("*************************************\n"); \r
  71   return task;\r
  72 }\r