Changed default mutliplicity estimator to V0A.

[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / PiKaPr / TPCTOFpA / AddTaskTPCTOFpA.C

\r
\r
AliAnalysisTask *AddTaskTPCTOFpA(Int_t identifier = 0, Bool_t isMC = kFALSE, Bool_t isTPConly = kFALSE, Bool_t writeOwnFile = kFALSE, Bool_t saveMotherPDG = kFALSE, Bool_t useEvenDcaBinning = kFALSE, Bool_t smallTHnSparse = kFALSE, Double_t nSigmaTPCLow= -3., Double_t nSigmaTPCHigh = 3., Double_t rapidityLow = -0.5, Double_t rapidityHigh = 0.5, Bool_t ispA=kTRUE, Bool_t rapCMS = kFALSE, TString multEst = "V0A",  Bool_t setTrackCuts = kFALSE, AliESDtrackCuts *ESDtrackCuts = 0){\r
\r
\r
  //get the current analysis manager\r
  AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
  if (!mgr) {\r
    Error("AddTask_janielsk_TPCTOFpA", "No analysis manager found.");\r
    return 0;\r
  }\r
  //============= Set Task Name ===================\r
  //TString taskName=("AliAnalysisTPCTOFpA.cxx+g");\r
  //===============================================\r
  //            Load the task\r
  //gROOT->LoadMacro(taskName.Data());\r
\r
\r
  \r
  //========= Add task to the ANALYSIS manager =====\r
\r
  //normal tracks\r
  AliAnalysisTPCTOFpA *task = new AliAnalysisTPCTOFpA("janielskTaskTPCTOFpA");\r
  task->SelectCollisionCandidates(AliVEvent::kCINT5);\r
\r
  //switches\r
  if (isMC) task->SetIsMCtrue(isMC);\r
  if (isTPConly)task->SetUseTPConlyTracks(isTPConly);\r
  if (saveMotherPDG) task->SetSaveMotherPDG(saveMotherPDG);\r
  if (useEvenDcaBinning) task->SetEvenDCAbinning(kTRUE);\r
  if (smallTHnSparse){\r
    task->SetSmallTHnSparse(kTRUE);\r
    task->SetTPCnSigmaCuts(nSigmaTPCLow,nSigmaTPCHigh);\r
    task->SetRapidityCuts(rapidityLow,rapidityHigh);\r
  }\r
  if (ispA) task->SetIspA(kTRUE);\r
  if (rapCMS) task->SetRapCMS(kTRUE);\r
  task->SetCentEst(multEst.Data());\r
\r
  //initialize task\r
  task->Initialize();\r
\r
  //esd cuts need to be set after initialize or cuts will be replaced by standard cuts in initialize\r
  if (setTrackCuts) task->SetESDtrackCuts(ESDtrackCuts);\r
\r
  //add task to manager\r
  mgr->AddTask(task);\r
\r
\r
  \r
\r
  //================================================\r
  //              data containers\r
  //================================================\r
  //            find input container\r
  //below the trunk version\r
  AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();\r
\r
/*\r
  //dumm output container\r
  AliAnalysisDataContainer *coutput0 =\r
      mgr->CreateContainer(Form("akalweit_tree%i",identifier),\r
                           TTree::Class(),\r
                           AliAnalysisManager::kExchangeContainer,\r
                           Form("akalweit_default%i",identifier));\r
\r
\r
  //define output containers, please use 'username'_'somename'\r
  AliAnalysisDataContainer *coutput1 = \r
      mgr->CreateContainer(Form("akalweit_TPCTOFpA%i",identifier), TList::Class(),\r
                           AliAnalysisManager::kOutputContainer,Form("akalweit_TPCTOFpA%i.root",identifier));\r
*/\r
\r
  AliAnalysisDataContainer *coutput1;\r
\r
  if (!writeOwnFile) {\r
  	coutput1 =  mgr->CreateContainer(Form("janielsk_TPCTOFpA%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("%s:janielsk_TPCTOFpA", AliAnalysisManager::GetCommonFileName())); \r
	}\r
  else {\r
	coutput1 =  mgr->CreateContainer(Form("janielsk_TPCTOFpA%i",identifier), TList::Class(), AliAnalysisManager::kOutputContainer, Form("janielsk_TPCTOFpA.root"));\r
	}\r
\r
\r
\r
  //connect containers\r
\r
  //\r
  mgr->ConnectInput  (task,  0, cinput );\r
  //mgr->ConnectOutput (task,  0, coutput0);\r
  mgr->ConnectOutput (task,  1, coutput1);\r
\r
  return task;\r
}\r