PWGLF/SPECTRA/PiKaPr/TestAOD/runAODProof.C

   1 void runAODProof(const char * proofMode = "full")
   2 {
   3
   4    gEnv->SetValue("XSec.GSI.DelegProxy", "2");
   5    //  TProof::Open("rbertens@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("libOADB.so");
  16    gSystem->Load("libANALYSISalice.so");
  17    gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  18
  19
  20    AliAnalysisAlien * handler = new AliAnalysisAlien("test");
  21    handler->SetOverwriteMode();
  22    handler->SetRunMode(proofMode);
  23    handler->SetProofReset(0);
  24    handler->SetAliROOTVersion("v4-21-29-AN");
  25    handler->SetProofCluster(Form("%s@alice-caf.cern.ch", gSystem->Getenv("CAFUSER")));
  26    //handler->SetProofCluster("rbertens@skaf.saske.sk");
  27    //   handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree");
  28    handler->SetProofDataSet("/alice/sim/LHC11a10a_000138653_AOD048#aodTree");
  29    handler->SetNproofWorkersPerSlave(1);
  30    handler->SetAliRootMode("default");
  31    handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
  32    handler->SetAnalysisSource("AliSpectraAODHistoManager.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODTrackCuts.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
  33    handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
  34
  35    //  handler->SetAliRootMode("");
  36    handler->SetClearPackages();
  37
  38
  39    AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
  40    mgr->SetGridHandler(handler);
  41    AliAODInputHandler* aodH = new AliAODInputHandler();
  42    mgr->SetInputEventHandler(aodH);
  43
  44    gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g");
  45    gROOT->LoadMacro("AliSpectraAODEventCuts.cxx+g");
  46    gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
  47    gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
  48
  49    // Add PID task
  50    gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
  51    AliAnalysisTask * taskPID = AddTaskPIDResponse();
  52    mgr->AddTask(taskPID);
  53
  54    AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise");
  55    task->SetIsMC(1);
  56    mgr->AddTask(task);
  57
  58    // Set the cuts
  59    AliSpectraAODEventCuts * vcuts = new AliSpectraAODEventCuts("Event Cuts");
  60    AliSpectraAODTrackCuts  * tcuts = new AliSpectraAODTrackCuts("Track Cuts");
  61    tcuts->SetTrackType(6);
  62    tcuts->SetEta(1.);
  63    tcuts->SetP(.7);
  64    // vcuts->SetCentralityCutMin(0.0) // default
  65    // vcuts->SetCentralityCutMax(0.0) // default
  66    task->SetEventCuts(vcuts);
  67    task->SetTrackCuts(tcuts);
  68    task->SetNSigmaForIdentification(3.);
  69
  70    AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
  71    AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer("chistpt", AliSpectraAODHistoManager::Class(),  AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  72    AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer("cvcutpt", AliSpectraAODEventCuts::Class(),    AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  73    AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer("ctcutpt", AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root");
  74
  75    mgr->ConnectInput(task, 0, cinput);
  76    mgr->ConnectOutput(task, 1, coutputpt1);
  77    mgr->ConnectOutput(task, 2, coutputpt2);
  78    mgr->ConnectOutput(task, 3, coutputpt3);
  79    mgr->SetDebugLevel(2);
  80
  81    if (!mgr->InitAnalysis()) return;
  82    mgr->PrintStatus();
  83    mgr->StartAnalysis("proof");
  84 }