Added Albericas code

[u/mrichter/AliRoot.git] / PWG2 / EVCHAR / runCentrality.C

void runCentrality(const char * type = "a")
{
  // Load common libraries
  gSystem->Load("libCore.so");  
  gSystem->Load("libTree.so");
  gSystem->Load("libGeom.so");
  gSystem->Load("libVMC.so");
  gSystem->Load("libPhysics.so");
  gSystem->Load("libSTEERBase");
  gSystem->Load("libESD");
  gSystem->Load("libAOD");
  gSystem->Load("libANALYSIS");
  gSystem->Load("libANALYSISalice");   
  // Use AliRoot includes to compile our task
  gROOT->ProcessLine(".include $ALICE_ROOT/include");
  // gROOT->ProcessLine(gSystem->ExpandPathName(".include $ALICE_ROOT/PWG2/FORWARD/analysis"));

  // form filename 
  TString filenameStr = Form("%s",type);
  filenameStr = TString("LHC10f8")+filenameStr;
  const char * filename = filenameStr.Data();
  
  // Create and configure the alien handler plugin
  gROOT->LoadMacro("CreateAlienHandler.C");
  AliAnalysisGrid *alienHandler = CreateAlienHandler(filename);  
  if(!alienHandler) return;
 
 // Create the analysis manager
  AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisManager");

  // Connect plug-in to the analysis manager
  mgr->SetGridHandler(alienHandler);
  
  // My task
  gROOT->LoadMacro("AliAnalysisTaskCentralityTreeMaker.cxx++g");   
  AliAnalysisTaskCentralityTreeMaker *task = new AliAnalysisTaskCentralityTreeMaker("CentralityTask");  
  // Writing (or not) output tree
  //task->SetTreeFilling(writeTree);
  task->SetMCInput();
  mgr->AddTask(task);

  AliMCEventHandler*  mcHandler = new AliMCEventHandler();
  mgr->SetMCtruthEventHandler(mcHandler);

  AliESDInputHandler* esdH = new AliESDInputHandler();
  mgr->SetInputEventHandler(esdH);
  
  // Physics selection
  //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
  //AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(kFALSE);
  // Selecting collision candidates
  //task->SelectCollisionCandidates();

  // Create containers for input/output
  mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
  
  AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("coutput1",TList::Class(),
                                            AliAnalysisManager::kOutputContainer, 
  					    "cenHistos.root");
  mgr->ConnectOutput(task, 1, coutput1);
  
  AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("coutput2",TTree::Class(),
  					    AliAnalysisManager::kOutputContainer,
        				    "cenTree.root");
  coutput2->SetSpecialOutput();  
  mgr->ConnectOutput(task, 2, coutput2);

  // Enable debug printouts
  mgr->SetDebugLevel(2);

  if(!mgr->InitAnalysis()) return;

  mgr->PrintStatus();
  // Start analysis in grid.
  mgr->StartAnalysis("grid");

};