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