1 AliAnalysisTaskFlowTPCEMCalEP* ConfigHFE_FLOW_TPCEMCal_EP(Bool_t useMC, Double_t openingAngle, Double_t invMass,TString nonHFEalgorithm){
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(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(1.5, 50);
19 hfecuts->SetMaxImpactParam(1,2);
21 AliAnalysisTaskFlowTPCEMCalEP *task = new AliAnalysisTaskFlowTPCEMCalEP("HFE v2");
22 printf("task ------------------------ %p\n ", task);
23 task->SetHFECuts(hfecuts);
24 task->SetOpeningAngleCut(openingAngle);
25 task->SetInvariantMassCut(invMass);
26 task->SetNonHFEalgorithm(nonHFEalgorithm);
29 AliHFEpid *pid = task->GetPID();
30 if(useMC) pid->SetHasMCData(kTRUE);
31 pid->AddDetector("TPC", 0);
32 pid->AddDetector("EMCAL", 1);
34 AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);
35 emcpid->SetEoPMax(1.2);
36 emcpid->SetEoPMim(0.9);
41 // Monte-Carlo needs modelling of the falling mean with momentum at low momentum
42 // for high momentum it is consistent with a flat -0.94
43 cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";
44 Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};
45 for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];
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);
56 printf("*************************************\n");
57 printf("Configuring standard Task:\n");
58 // task->PrintStatus();
60 printf("*************************************\n");