new config for systematic study TPC+EMCal

[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / configs / PbPb / ConfigHFECalSys.C

AliAnalysisTaskHFECal* ConfigHFECalSys(Bool_t useMC, int TPCclust, int TPCclustPID, int Nits, int ITSstat){\r
  //\r
  // HFE standard task configuration\r
  //\r
\r
  Bool_t kAnalyseTaggedTracks = kTRUE;\r
  \r
  AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsEMCAL","HFE Standard Cuts");\r
  hfecuts->CreateStandardCuts();\r
  hfecuts->SetMinNClustersTPC(TPCclust);\r
  hfecuts->SetMinNClustersTPCPID(TPCclustPID);\r
  hfecuts->SetMinRatioTPCclusters(0.6);\r
  hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);\r
  hfecuts->SetMinNClustersITS(Nits);\r
  if(ITSstat==0)hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny);\r
  if(ITSstat==1)hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);\r
  hfecuts->SetCheckITSLayerStatus(kFALSE);\r
  hfecuts->SetVertexRange(10.);\r
  hfecuts->SetTOFPIDStep(kFALSE);\r
  hfecuts->SetPtRange(2, 50);\r
  hfecuts->SetMaxImpactParam(3.,3.);\r
  \r
  AliAnalysisTaskHFECal *task = new AliAnalysisTaskHFECal("HFEanalysisEMCal");\r
  printf("task ------------------------ %p\n ", task);\r
  task->SetHFECuts(hfecuts);\r
  task->SetInvariantMassCut(0.05);\r
\r
  // Define PID\r
  AliHFEpid *pid = task->GetPID();\r
  if(useMC) pid->SetHasMCData(kTRUE);\r
  pid->AddDetector("TPC", 0);\r
  pid->AddDetector("EMCAL", 1);\r
\r
  Double_t params[4];\r
  char *cutmodel;\r
  if(useMC){\r
	  // Monte-Carlo needs modelling of the falling mean with momentum at low momentum\r
	  // for high momentum it is consistent with a flat -0.94\r
	  //cutmodel = "[0]*TMath::Exp([1]*x) + [2] + [3]*x";\r
	  //Double_t paramsMC[4] = {0.7174, -1.588, -0.9395, 0.0246};\r
	  //for(int ipar = 0; ipar < 4; ipar++) params[ipar] = paramsMC[ipar];\r
	  cutmodel = "pol0";\r
	  params[0] = -1.0; //sigma min\r
  } else {\r
	  // Data is consistent with a flat 0.12\r
	  cutmodel = "pol0";\r
	  //params[0] = -0.0015;\r
	  //params[0] = -3.0;\r
	  //params[0] = -0.05; //sigma min\r
	  params[0] = -1.0; //sigma min\r
  }\r
  //pid->ConfigureTPCdefaultCut(cutmodel, params,3.0); \r
  for(Int_t a=0;a<11;a++)pid->ConfigureTPCcentralityCut(a,cutmodel,params,3.0);\r
  \r
\r
  // change E/p cuts\r
  AliHFEpidEMCAL *emcpid = pid->AliHFEpid::GetDetPID(AliHFEpid::kEMCALpid);\r
  emcpid->SetEoPMax(1.3);\r
  emcpid->SetEoPMim(0.9);\r
\r
  printf("*************************************\n");\r
  printf("Configuring standard Task:\n");\r
//  task->PrintStatus();\r
  pid->PrintStatus();\r
  printf("*************************************\n"); \r
  return task;\r
}\r