--- /dev/null
+AliAnalysisTaskElecHadronCorrel* ConfigHFEElecHadronCorrelPbPb(Bool_t useMC){\r
+ //\r
+ // HFE standard task configuration\r
+ //\r
+\r
+ Bool_t kAnalyseTaggedTracks = kTRUE;\r
+ \r
+ AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");\r
+// hfecuts->CreateStandardCuts();\r
+ hfecuts->SetMinNClustersTPC(120);\r
+ hfecuts->SetMinNClustersITS(3);\r
+ hfecuts->SetMinNTrackletsTRD(0);\r
+ hfecuts->SetMinRatioTPCclusters(0.6);\r
+ // hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r
+ hfecuts->SetRequireITSPixel();\r
+ hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); \r
+ hfecuts->SetMaxChi2perClusterITS(-1);\r
+ hfecuts->SetMaxChi2perClusterTPC(3.5);\r
+ hfecuts->SetCheckITSLayerStatus(kFALSE); // shud be put back\r
+// hfecuts->UnsetVertexRequirement();\r
+ hfecuts->SetVertexRange(10.);\r
+ hfecuts->SetRequireSigmaToVertex();\r
+ //hfecuts->SetSigmaToVertex(10);\r
+ hfecuts->SetTOFPIDStep(kFALSE);\r
+// hfecuts->SetQAOn();\r
+ hfecuts->SetPtRange(0, 30);\r
+\r
+ AliAnalysisTaskElecHadronCorrel *task = new AliAnalysisTaskElecHadronCorrel("HFE-hadron correlations");\r
+ printf("task ------------------------ %p\n ", task);\r
+ task->SetHFECuts(hfecuts);\r
+ task->SetInvariantMassCut(0.07);\r
+ task->SetRejectKinkMother(kTRUE);\r
+// task->SetRemovePileUp(kTRUE);\r
+\r
+ // Define PID\r
+ AliHFEpid *pid = task->GetPID();\r
+ if(useMC) pid->SetHasMCData(kTRUE);\r
+ pid->AddDetector("TPC", 0);\r
+ pid->AddDetector("EMCAL", 1);\r
+ // change E/p cuts\r
+ AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r
+ emcpid->SetEoPMax(1.2);\r
+ emcpid->SetEoPMim(0.8);\r
+\r
+ Double_t params[4];\r
+ char *cutmodel;\r
+ if(useMC){\r
+ // Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r
+ // for high momentum it is consistent with a flat -0.94\r
+ cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
+ Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
+ for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
+ } else {\r
+ // Data is consistent with a flat 0.12\r
+ cutmodel = "pol0";\r
+ //params[0] = -0.0015;\r
+ //params[0] = -3.0;\r
+ params[0] = -1.0; //sigma min\r
+ }\r
+// pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
+ pid->ConfigureTPCasymmetric(0,30,-1,3.0); \r
+\r
+ printf("*************************************\n");\r
+ printf("Configuring standard Task:\n");\r
+// task->PrintStatus();\r
+ pid->PrintStatus();\r
+ printf("*************************************\n"); \r
+ return task;\r
+}\r