1 AliAnalysisTaskHFECal* ConfigHFECalSys(Bool_t useMC, int TPCclust, int TPCclustPID, int Nits, int ITSstat, int QAhist){
3 // HFE standard task configuration
6 Bool_t kAnalyseTaggedTracks = kTRUE;
8 AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");
9 hfecuts->CreateStandardCuts();
10 hfecuts->SetMinNClustersTPC(TPCclust);
11 hfecuts->SetMinNClustersTPCPID(TPCclustPID);
12 hfecuts->SetMinRatioTPCclusters(0.6);
13 hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
14 hfecuts->SetMinNClustersITS(Nits);
15 if(ITSstat==0)hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny);
16 if(ITSstat==1)hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);
17 hfecuts->SetCheckITSLayerStatus(kFALSE);
18 hfecuts->SetVertexRange(10.);
19 hfecuts->SetTOFPIDStep(kFALSE);
20 hfecuts->SetPtRange(2, 50);
21 hfecuts->SetMaxImpactParam(3.,3.);
23 AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");
24 printf("task ------------------------ %p\n ", task);
25 task->SetHFECuts(hfecuts);
26 task->SetInvariantMassCut(0.05);
27 printf("<<<<<< ------------------ QA status %d",QAhist);
28 task->SetQAHist(QAhist);
31 AliHFEpid *pid = task->GetPID();
32 if(useMC) pid->SetHasMCData(kTRUE);
33 pid->AddDetector("TPC", 0);
34 pid->AddDetector("EMCAL", 1);
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];
45 params[0] = -1.0; //sigma min
47 // Data is consistent with a flat 0.12
49 //params[0] = -0.0015;
51 //params[0] = -0.05; //sigma min
52 params[0] = -1.0; //sigma min
54 //pid->ConfigureTPCdefaultCut(cutmodel, params,3.0);
55 for(Int_t a=0;a<11;a++)pid->ConfigureTPCcentralityCut(a,cutmodel,params,3.0);
59 AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
60 emcpid->SetEoPMax(1.3);
61 emcpid->SetEoPMim(0.9);
63 printf("*************************************\n");
64 printf("Configuring standard Task:\n");
65 // task->PrintStatus();
67 printf("*************************************\n");