updated for Non-HFE systemaitc study

[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / configs / PbPb / ConfigHFECal.C

diff --git a/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C b/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C
index da80daeb4edd06bad52faf4b710756bb6958c98b..b5d5a22879725965e4125969128a4d87b021d370 100644 (file)
--- a/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C
+++ b/PWGHF/hfe/macros/configs/PbPb/ConfigHFECal.C
@@ -1,4 +1,4 @@
-AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC){\r
+AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC,Bool_t MassConst,Bool_t MassWidthCut){\r
   //\r
   // HFE standard task configuration\r
   //\r
@@ -21,7 +21,11 @@ AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC){
   AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");\r
   printf("task ------------------------ %p\n ", task);\r
   task->SetHFECuts(hfecuts);\r
-  task->SetInvariantMassCut(0.05);\r
+  task->SetMassConstraint(MassConst);\r
+  task->SetMassWidthCut(MassWidthCut);\r
+  Double_t masscut = 0.05;\r
+  if(!MassConst)masscut = 0.1;\r
+  task->SetInvariantMassCut(masscut);\r
 \r
   // Define PID\r
   AliHFEpid *pid = task->GetPID();\r
@@ -47,7 +51,9 @@ AliAnalysisTaskHFECal* ConfigHFECal(Bool_t useMC){
          //params[0] = -0.05; //sigma min\r
          params[0] = -1.0; //sigma min\r
   }\r
-  pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
+  //pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
+  for(Int_t a=0;a<11;a++)pid->ConfigureTPCcentralityCut(a,cutmodel,params,3.0);\r
+  \r
 \r
   // change E/p cuts\r
   AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r