centrality selection

[u/mrichter/AliRoot.git] / PWGCF / FLOW / macros / AddTaskFlowEPCascade.C

AliAnalysisTask *AddTaskFlowEPCascade(Float_t centrMin=0, Float_t centrMax=100,\r
				      TString fileName, TString suffixName="",\r
				      Double_t vtxCut = 10., \r
				      Double_t etaCut = 0.9, \r
				      Int_t nTPCcls = 70,\r
				      Double_t mSigma = 0.0024)\r
{\r
  /*\r
  //-E-V-E-N-T- -c-u-t-s-----------------------------------------------------\r
  AliFlowEventCuts* cutsEvent \r
    = new AliFlowEventCuts(Form("event_cuts_%s", suffixName.Data()));\r
  cutsEvent->SetCentralityPercentileRange(centrMin,centrMax);\r
  cutsEvent->SetCentralityPercentileMethod(AliFlowEventCuts::kV0);\r
  cutsEvent->SetRefMultMethod(AliFlowEventCuts::kV0);\r
  cutsEvent->SetNContributorsRange(2);\r
  cutsEvent->SetPrimaryVertexZrange( -10., 10. );\r
  cutsEvent->SetCutSPDvertexerAnomaly();\r
  cutsEvent->SetCutZDCtiming();\r
  cutsEvent->SetCutTPCmultiplicityOutliers();\r
\r
  AliFlowTrackCuts* cutsRP \r
    = (AliFlowTrackCuts*) AliFlowTrackCuts::GetStandardGlobalTrackCuts2010();\r
  cutsRP->SetPtRange(0.15, 10.);  //added to extend pt range\r
  */\r
\r
  AliFlowTrackCuts * cutsDaughter \r
    = new AliFlowTrackCuts(Form("daughter_cuts_%s",suffixName.Data()));\r
  cutsDaughter->SetPtRange(0.15,10.0);\r
  cutsDaughter->SetEtaRange(-etaCut,etaCut);\r
  cutsDaughter->SetMinNClustersTPC(nTPCcls);\r
  cutsDaughter->SetMaxChi2PerClusterTPC(4.0);\r
  cutsDaughter->SetAODfilterBit(1); //TPC track only\r
  cutsDaughter->SetRequireITSRefit(kFALSE);\r
  cutsDaughter->SetRequireTPCRefit(kTRUE);\r
  cutsDaughter->SetMinNClustersITS(0);\r
  cutsDaughter->SetAcceptKinkDaughters(kFALSE);\r
\r
  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
  AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();\r
\r
  //____________________________________________//                           \r
  double pdg = 1.32171; double hwid = mSigma;\r
  double sep = 6*hwid; double XiStart = pdg-sep; \r
  double XiBands[3][2];\r
  XiBands[0][0]=(XiStart+6*hwid*0)-2*hwid; \r
  XiBands[0][1]=(XiStart+6*hwid*0)+2*hwid;\r
  XiBands[1][0]=(XiStart+6*hwid*1)-2*hwid;   \r
  XiBands[1][1]=(XiStart+6*hwid*1)+2*hwid;\r
  XiBands[2][0]=(XiStart+6*hwid*2)-2*hwid; \r
  XiBands[2][1]=(XiStart+6*hwid*2)+2*hwid;\r
\r
  double OmegaPDG = 1.67245;   \r
  double Omegasg = mSigma; // 1 sigma?\r
  double OmegaStart = OmegaPDG-6*Omegasg; \r
  double OmegaBands[3][2];\r
  OmegaBands[0][0]=(OmegaStart+6*Omegasg*0)-2*Omegasg;   \r
  OmegaBands[0][1]=(OmegaStart+6*Omegasg*0)+2*Omegasg;\r
  OmegaBands[1][0]=(OmegaStart+6*Omegasg*1)-2*Omegasg;     \r
  OmegaBands[1][1]=(OmegaStart+6*Omegasg*1)+2*Omegasg;\r
  OmegaBands[2][0]=(OmegaStart+6*Omegasg*2)-2*Omegasg;   \r
  OmegaBands[2][1]=(OmegaStart+6*Omegasg*2)+2*Omegasg;\r
\r
  AliAnalysisTaskFlowEPCascade * task\r
    = new AliAnalysisTaskFlowEPCascade(Form("v2EP_%s", \r
					    suffixName.Data()), \r
				       centrMin, centrMax,\r
				       XiBands, OmegaBands);\r
  // task->SetFlowEventCuts(cutsEvent);\r
  //task->SetFlowTrackCuts(cutsRP);\r
  task->SetVertexCut(vtxCut);\r
  task->SetFlowDauTrackCuts(cutsDaughter);\r
  task->SelectCollisionCandidates(AliVEvent::kMB);\r
  mgr->AddTask(task);\r
\r
  AliAnalysisDataContainer *coutput1 =\r
    mgr->CreateContainer(Form("hist_%s", suffixName.Data()), \r
			 TList::Class(),\r
                         AliAnalysisManager::kOutputContainer,\r
                         Form("%s.root",fileName.Data()));\r
  // Connect input/output                                                      \r
  mgr->ConnectInput(task, 0, cinput);\r
  //mgr->ConnectInput(task, 1, cinput1);\r
  mgr->ConnectOutput(task, 1, coutput1);\r
  \r
  return task;\r
}\r
\r