ROOT logo
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