7b443030de73979f83bf9600380715b2df42e7cb
[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / configs / PbPb / ConfigHFEElecV2.C
1 AliAnalysisTaskElecV2* ConfigHFEElecV2(Bool_t useMC){
2   //
3   // HFE standard task configuration
4   //
5
6   Bool_t kAnalyseTaggedTracks = kTRUE;
7   
8   AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");
9   hfecuts->CreateStandardCuts();
10   hfecuts->SetMinNClustersTPC(100);
11   hfecuts->SetMinRatioTPCclusters(0.6);
12   hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
13   hfecuts->SetMinNClustersITS(3);
14   hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);
15   hfecuts->SetCheckITSLayerStatus(kFALSE);
16   hfecuts->SetVertexRange(10.);
17   hfecuts->SetTOFPIDStep(kFALSE);
18   hfecuts->SetPtRange(2, 50);
19   hfecuts->SetMaxImpactParam(1,2);
20   
21   AliAnalysisTaskElecV2 *task = new AliAnalysisTaskElecV2("HFE v2");
22   printf("task ------------------------ %p\n ", task);
23   task->SetHFECuts(hfecuts);
24   task->SetInvariantMassCut(0.05);
25
26   // Define PID
27   AliHFEpid *pid = task->GetPID();
28   if(useMC) pid->SetHasMCData(kTRUE);
29   pid->AddDetector("TPC", 0);
30   pid->AddDetector("EMCAL", 1);
31   // change E/p cuts
32   AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
33   emcpid->SetEoPMax(1.2);
34   emcpid->SetEoPMim(0.9);
35
36   Double_t params[4];
37   char *cutmodel;
38   if(useMC){
39           // Monte-Carlo needs modelling of the falling mean with momentum at low momentum
40           // for high momentum it is consistent with a flat -0.94
41           cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";
42           Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};
43           for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];
44   } else {
45           // Data is consistent with a flat 0.12
46           cutmodel = "pol0";
47           //params[0] = -0.0015;
48           //params[0] = -3.0;
49           //params[0] = -0.05; //sigma min
50           params[0] = -1.0; //sigma min
51   }
52   pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); 
53
54   printf("*************************************\n");
55   printf("Configuring standard Task:\n");
56 //  task->PrintStatus();
57   pid->PrintStatus();
58   printf("*************************************\n"); 
59   return task;
60 }