ROOT logo
AliAnalysisTaskElecHadronCorrel* ConfigHFEElecHadronCorrelPbPb(Bool_t useMC,
                                                               Bool_t EventTrigSelMB=kTRUE,
                                                               Bool_t TrigSelCen = kTRUE,
                                                               Double_t CentMin = 0,
                                                               Double_t CentMax = 7,
                                                               Double_t TPCNsigMinE = -2,
                                                               Double_t TPCNsigMaxE = 2,
                                                               Double_t TPCNsigMinH = -10,
                                                               Double_t TPCNsigMaxH = -3.5,
                                                               Double_t SSM02Min = 0.03,  
                                                               Double_t SSM02Max = 0.5,
                                                               Double_t SSM20Min = 0.03,  
                                                               Double_t SSM20Max = 0.3,
                                                               Double_t Disp = 1,
                                                               Double_t EovPMin = 0.8,    
                                                               Double_t EovPMax = 1.2,    
                                                               Double_t InvM = 0.1,
                                                               const char* CentralityMet = "VOM",
                                                               Int_t TPCNClsEle = 100,
                                                               Int_t TPCNClsHad = 80,
                                                               Bool_t AssoEleITSref=kTRUE,
                                                               Int_t AssoElecTPCNCls=80,
                                                               TString ContNameExt = "Central",
                                                               TString TaskName="hfeCorrl"){

//
  // HFE standard task configuration
  //

  Bool_t kAnalyseTaggedTracks = kTRUE;
  
  AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");
//  hfecuts->CreateStandardCuts();
  hfecuts->SetMinNClustersTPC(TPCNClsEle);
  hfecuts->SetMinNClustersITS(3);
  hfecuts->SetMinNTrackletsTRD(0);
  hfecuts->SetMinRatioTPCclusters(0.6);
   // hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
  hfecuts->SetRequireITSPixel();
  hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); 
  hfecuts->SetMaxChi2perClusterITS(-1);
  hfecuts->SetMaxChi2perClusterTPC(3.5);
  hfecuts->SetCheckITSLayerStatus(kFALSE); // shud be put back
//  hfecuts->UnsetVertexRequirement();
  hfecuts->SetVertexRange(10.);
  hfecuts->SetRequireSigmaToVertex();
  //hfecuts->SetSigmaToVertex(10);
  hfecuts->SetTOFPIDStep(kFALSE);
//  hfecuts->SetQAOn();
 hfecuts->SetPtRange(0, 30);

 TString taskName = TaskName;

  AliAnalysisTaskElecHadronCorrel *task = new AliAnalysisTaskElecHadronCorrel(taskName);
  printf("task ------------------------ %p\n ", task);
  task->SetHFECuts(hfecuts);
//  task->SetRemovePileUp(kTRUE);
//  task->SetInvariantMassCut(0.1);

  task->SetEventTriggerSelectionMB(EventTrigSelMB);
  task->SetTriggerSelection(TrigSelCen);
  task->SetCentralityParameters(CentMin, CentMax, CentralityMet);
  task->SetInvariantMassCut(InvM);
  task->SetAssoElecTPCNCls(AssoElecTPCNCls);
  task->SetTPCnsigmaCutsElecSelection(TPCNsigMinE,TPCNsigMaxE);
  task->SetTPCnsigmaCutsHadSelection(TPCNsigMinH,TPCNsigMaxH);
  task->SetShowerShapeCutsM02(SSM02Min,SSM02Max);
  task->SetShowerShapeCutsM20(SSM20Min,SSM20Max);
  task->SetShowerShapeCutsDisp(0,Disp);
  task->SetEovPCuts(EovPMin,EovPMax);
  task->SetRejectKinkMother(kTRUE);
  task->SetTPCNClsHad(TPCNClsHad);
  task->SetAssoElecITSrefit(AssoEleITSref);

  // 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.8);

  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] = -1.0; //sigma min
  }
//  pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); 
  pid->ConfigureTPCasymmetric(0,30,-1,3.0); 
*/
  printf("*************************************\n");
  printf("Configuring standard Task:\n");
//  task->PrintStatus();
  pid->PrintStatus();
  printf("*************************************\n"); 
  return task;
}
 ConfigHFEElecHadronCorrelPbPb.C:1
 ConfigHFEElecHadronCorrelPbPb.C:2
 ConfigHFEElecHadronCorrelPbPb.C:3
 ConfigHFEElecHadronCorrelPbPb.C:4
 ConfigHFEElecHadronCorrelPbPb.C:5
 ConfigHFEElecHadronCorrelPbPb.C:6
 ConfigHFEElecHadronCorrelPbPb.C:7
 ConfigHFEElecHadronCorrelPbPb.C:8
 ConfigHFEElecHadronCorrelPbPb.C:9
 ConfigHFEElecHadronCorrelPbPb.C:10
 ConfigHFEElecHadronCorrelPbPb.C:11
 ConfigHFEElecHadronCorrelPbPb.C:12
 ConfigHFEElecHadronCorrelPbPb.C:13
 ConfigHFEElecHadronCorrelPbPb.C:14
 ConfigHFEElecHadronCorrelPbPb.C:15
 ConfigHFEElecHadronCorrelPbPb.C:16
 ConfigHFEElecHadronCorrelPbPb.C:17
 ConfigHFEElecHadronCorrelPbPb.C:18
 ConfigHFEElecHadronCorrelPbPb.C:19
 ConfigHFEElecHadronCorrelPbPb.C:20
 ConfigHFEElecHadronCorrelPbPb.C:21
 ConfigHFEElecHadronCorrelPbPb.C:22
 ConfigHFEElecHadronCorrelPbPb.C:23
 ConfigHFEElecHadronCorrelPbPb.C:24
 ConfigHFEElecHadronCorrelPbPb.C:25
 ConfigHFEElecHadronCorrelPbPb.C:26
 ConfigHFEElecHadronCorrelPbPb.C:27
 ConfigHFEElecHadronCorrelPbPb.C:28
 ConfigHFEElecHadronCorrelPbPb.C:29
 ConfigHFEElecHadronCorrelPbPb.C:30
 ConfigHFEElecHadronCorrelPbPb.C:31
 ConfigHFEElecHadronCorrelPbPb.C:32
 ConfigHFEElecHadronCorrelPbPb.C:33
 ConfigHFEElecHadronCorrelPbPb.C:34
 ConfigHFEElecHadronCorrelPbPb.C:35
 ConfigHFEElecHadronCorrelPbPb.C:36
 ConfigHFEElecHadronCorrelPbPb.C:37
 ConfigHFEElecHadronCorrelPbPb.C:38
 ConfigHFEElecHadronCorrelPbPb.C:39
 ConfigHFEElecHadronCorrelPbPb.C:40
 ConfigHFEElecHadronCorrelPbPb.C:41
 ConfigHFEElecHadronCorrelPbPb.C:42
 ConfigHFEElecHadronCorrelPbPb.C:43
 ConfigHFEElecHadronCorrelPbPb.C:44
 ConfigHFEElecHadronCorrelPbPb.C:45
 ConfigHFEElecHadronCorrelPbPb.C:46
 ConfigHFEElecHadronCorrelPbPb.C:47
 ConfigHFEElecHadronCorrelPbPb.C:48
 ConfigHFEElecHadronCorrelPbPb.C:49
 ConfigHFEElecHadronCorrelPbPb.C:50
 ConfigHFEElecHadronCorrelPbPb.C:51
 ConfigHFEElecHadronCorrelPbPb.C:52
 ConfigHFEElecHadronCorrelPbPb.C:53
 ConfigHFEElecHadronCorrelPbPb.C:54
 ConfigHFEElecHadronCorrelPbPb.C:55
 ConfigHFEElecHadronCorrelPbPb.C:56
 ConfigHFEElecHadronCorrelPbPb.C:57
 ConfigHFEElecHadronCorrelPbPb.C:58
 ConfigHFEElecHadronCorrelPbPb.C:59
 ConfigHFEElecHadronCorrelPbPb.C:60
 ConfigHFEElecHadronCorrelPbPb.C:61
 ConfigHFEElecHadronCorrelPbPb.C:62
 ConfigHFEElecHadronCorrelPbPb.C:63
 ConfigHFEElecHadronCorrelPbPb.C:64
 ConfigHFEElecHadronCorrelPbPb.C:65
 ConfigHFEElecHadronCorrelPbPb.C:66
 ConfigHFEElecHadronCorrelPbPb.C:67
 ConfigHFEElecHadronCorrelPbPb.C:68
 ConfigHFEElecHadronCorrelPbPb.C:69
 ConfigHFEElecHadronCorrelPbPb.C:70
 ConfigHFEElecHadronCorrelPbPb.C:71
 ConfigHFEElecHadronCorrelPbPb.C:72
 ConfigHFEElecHadronCorrelPbPb.C:73
 ConfigHFEElecHadronCorrelPbPb.C:74
 ConfigHFEElecHadronCorrelPbPb.C:75
 ConfigHFEElecHadronCorrelPbPb.C:76
 ConfigHFEElecHadronCorrelPbPb.C:77
 ConfigHFEElecHadronCorrelPbPb.C:78
 ConfigHFEElecHadronCorrelPbPb.C:79
 ConfigHFEElecHadronCorrelPbPb.C:80
 ConfigHFEElecHadronCorrelPbPb.C:81
 ConfigHFEElecHadronCorrelPbPb.C:82
 ConfigHFEElecHadronCorrelPbPb.C:83
 ConfigHFEElecHadronCorrelPbPb.C:84
 ConfigHFEElecHadronCorrelPbPb.C:85
 ConfigHFEElecHadronCorrelPbPb.C:86
 ConfigHFEElecHadronCorrelPbPb.C:87
 ConfigHFEElecHadronCorrelPbPb.C:88
 ConfigHFEElecHadronCorrelPbPb.C:89
 ConfigHFEElecHadronCorrelPbPb.C:90
 ConfigHFEElecHadronCorrelPbPb.C:91
 ConfigHFEElecHadronCorrelPbPb.C:92
 ConfigHFEElecHadronCorrelPbPb.C:93
 ConfigHFEElecHadronCorrelPbPb.C:94
 ConfigHFEElecHadronCorrelPbPb.C:95
 ConfigHFEElecHadronCorrelPbPb.C:96
 ConfigHFEElecHadronCorrelPbPb.C:97
 ConfigHFEElecHadronCorrelPbPb.C:98
 ConfigHFEElecHadronCorrelPbPb.C:99
 ConfigHFEElecHadronCorrelPbPb.C:100
 ConfigHFEElecHadronCorrelPbPb.C:101
 ConfigHFEElecHadronCorrelPbPb.C:102
 ConfigHFEElecHadronCorrelPbPb.C:103
 ConfigHFEElecHadronCorrelPbPb.C:104
 ConfigHFEElecHadronCorrelPbPb.C:105
 ConfigHFEElecHadronCorrelPbPb.C:106
 ConfigHFEElecHadronCorrelPbPb.C:107
 ConfigHFEElecHadronCorrelPbPb.C:108
 ConfigHFEElecHadronCorrelPbPb.C:109
 ConfigHFEElecHadronCorrelPbPb.C:110
 ConfigHFEElecHadronCorrelPbPb.C:111