PWGLF/SPECTRA/PiKaPr/TestAOD/runAOD.C

   1 void runAOD(Bool_t isMC = 0)
   2 {
   3
   4 //gEnv->SetValue("XSec.GSI.DelegProxy", "2");
   5 //TProof::Open("pverstee@alice-caf.cern.ch");
   6
   7    gSystem->Load("libTree.so");
   8    gSystem->Load("libGeom.so");
   9    gSystem->Load("libVMC.so");
  10    gSystem->Load("libPhysics.so");
  11    gSystem->Load("libSTEERBase.so");
  12    gSystem->Load("libESD.so");
  13    gSystem->Load("libAOD.so");
  14    gSystem->Load("libANALYSIS.so");
  15    gSystem->Load("libANALYSISalice.so");
  16    gSystem->Load("libANALYSIS");
  17    gSystem->Load("libANALYSISalice");
  18    gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
  19
  20    ifstream flist;
  21    flist.open("filelist.txt");
  22    TChain * chain = new TChain("aodTree");
  23    TString line;
  24    while (line.ReadLine(flist))
  25    {
  26       gSystem->ExpandPathName(line);
  27       chain->AddFile(line.Data());
  28       cout << "Adding file " << line.Data() << endl;
  29    }
  30 //  gProof->UploadPackage("AF-v4-19-04-AN");
  31 //  gProof->EnablePackage("AF-v4-19-04-AN");
  32    gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  33    gStyle->SetPalette(1);
  34
  35
  36
  37    AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
  38    AliAODInputHandler* aodH = new AliAODInputHandler();
  39    mgr->SetInputEventHandler(aodH);
  40
  41    gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
  42    gROOT->LoadMacro("AliSpectraAODVertexCuts.cxx+g");
  43    gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
  44    gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
  45
  46    gROOT->LoadMacro("$(ALICE_ROOT)/ANALYSIS/macros/AddTaskPIDResponse.C");
  47    AliAnalysisTask * taskPID = AddTaskPIDResponse();
  48    mgr->AddTask(taskPID);
  49
  50
  51    AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
  52    task->SetIsMC(isMC);
  53    mgr->AddTask(task);
  54
  55    // Set the cuts
  56    AliSpectraAODVertexCuts * vcuts = new AliSpectraAODVertexCuts("Vertex Cuts");
  57    AliSpectraAODTrackCuts  * tcuts = new AliSpectraAODTrackCuts("Tracks Cuts");
  58    tcuts->SetTrackType(6);
  59    tcuts->SetEta(1.);
  60    task->SetVertexCuts(vcuts);
  61    task->SetTrackCuts(tcuts);
  62
  63
  64    AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  65    AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer("chistpt", AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  66    AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer("cvcutpt", AliSpectraAODVertexCuts::Class(),    AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  67    AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer("ctcutpt", AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  68
  69    mgr->ConnectInput(task, 0, cinput);
  70    mgr->ConnectOutput(task, 1, coutputpt1);
  71    mgr->ConnectOutput(task, 2, coutputpt2);
  72    mgr->ConnectOutput(task, 3, coutputpt3);
  73    mgr->SetDebugLevel(2);
  74
  75    if (!mgr->InitAnalysis()) return;
  76    mgr->PrintStatus();
  77    mgr->StartAnalysis("local", chain, 100);
  78 }