//=========================================================================//
//                                                                         //
//           Analysis AddTask for Particle Ratio Fluctuation Study         //
//              Author: Deepika Rathee  || Satyajit Jenara                 //
//                      drathee@cern.ch || sjena@cern.ch 
//                        Thu Jun 19 11:44:51 CEST 2014
//                                                                         //
//=========================================================================//

void AddAliEbyEPidBayesTTask(Double_t ptl=0.5, 
			Int_t AODfilterBit = 768, 
			Int_t pidtype = 3, 
			Int_t requestTofPid = 1, 
			Double_t nSigmaCut = 3., 
			Bool_t isMC = 0, 
			TString centralityEstimator = "V0M", 
			TString ctaskname = "2011") {

 
  TString taskname = "EbyECF_";
  taskname += ctaskname;
  taskname += "_";
  taskname += Form("%d",isMC);
  taskname += "_";
  taskname += centralityEstimator;
  taskname += "_";
  taskname += Form("%d",AODfilterBit);
  
 
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  if (!mgr) {
    ::Error("AddTaskFluctuations", "No analysis manager to connect to.");
    return NULL;
  }

  if (!mgr->GetInputEventHandler()) {
    ::Error("AddTaskFluctuations", "This task requires an input event handler");
    return NULL;
  }
  TString type = mgr->GetInputEventHandler()->GetDataType(); 

 
  TString basefilename = AliAnalysisManager::GetCommonFileName();
  

  AliHelperPID* help = new AliHelperPID();
  help->SetNSigmaCut(nSigmaCut);
  help->SetPIDType(pidtype);
  if (requestTofPid) help->SetfRequestTOFPID(requestTofPid);
  if (ptl != 0 ) help->SetfPtTOFPID(ptl);
  
  if (pidtype == 3){
    AliPIDCombined *pidc=new AliPIDCombined();
    pidc->SetDefaultTPCPriors();
    help->SetPIDCombined(pidc);
  }
  
  AliEbyEPidTTask *taskqa;
  taskqa = new AliEbyEPidTTask(taskname.Data());
  taskqa->SetAODtrackCutBit(AODfilterBit);
  taskqa->SetCentralityEstimator(centralityEstimator.Data());
  taskqa->SetAnalysisType(isMC);
  taskqa->SetHelperPID(help);
  mgr->AddTask(taskqa);
  
  AliAnalysisDataContainer *coutqa = mgr->CreateContainer(Form("%s_QA",taskname.Data()),TList::Class(),AliAnalysisManager::kOutputContainer, Form("%s",basefilename.Data()));
  AliAnalysisDataContainer *coutt = mgr->CreateContainer("fEventTree",TTree::Class(),AliAnalysisManager::kOutputContainer, Form("%s",basefilename.Data()));

  mgr->ConnectInput(taskqa, 0, mgr->GetCommonInputContainer());
  mgr->ConnectOutput(taskqa, 1, coutqa);
 mgr->ConnectOutput(taskqa, 2, coutt);
 
 
  return;
}