PWG/muon/RunAODGeneration.C

   1 //
   2 //  Macro for AOD generation
   3 //  Gines Martinez, Subatech, October. 2007
   4 //  Generated from Andrea Morsch macro JetAnalysisManagerLoc.C
   5 //
   6 //  In this example the libraries PWG0base and PWG3base (not really needed in this version)
   7 //  are supposed to be loaded from a par file via RunAnalysis
   8 //
   9 void RunAODGeneration(char* filein = "AliESDs.root", char* fileout = "AliAOD.root" )
  10 {
  11     gSystem->Load("libTree");
  12     gSystem->Load("libGeom");
  13     gSystem->Load("libVMC");
  14     gSystem->Load("libANALYSIS");
  15     gSystem->Load("libSTEERBase");
  16     gSystem->Load("libAOD");
  17     gSystem->Load("libESD");
  18
  19     // Input ESD files
  20     TChain* chain = new TChain("esdTree");
  21     chain->Add(filein);
  22
  23     // Make aod handler
  24     AliAODHandler* aodHandler = new AliAODHandler();
  25     aodHandler->SetOutputFileName(fileout);
  26
  27     // Make the analysis manager
  28     AliAnalysisManager *mgr  = new AliAnalysisManager("AOD Manager", "AOD Manager");
  29     mgr->SetOutputEventHandler(aodHandler);
  30     mgr-> SetDebugLevel(10);
  31
  32     // Task for AOD generation from PWG0base directory
  33     AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");
  34     esdfilter->SetDebugLevel(10);
  35     mgr->AddTask(esdfilter);
  36
  37     //
  38     // Create containers for input/output
  39     AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
  40                                                              AliAnalysisManager::kInputContainer);
  41     AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
  42                                                               AliAnalysisManager::kOutputContainer, "default");
  43
  44     mgr->ConnectInput  (esdfilter,  0, cinput1  );
  45     mgr->ConnectOutput (esdfilter,  0, coutput1 );
  46
  47
  48     // Run the analysis
  49     mgr->InitAnalysis();
  50     mgr->PrintStatus();
  51     mgr->StartAnalysis("local",chain);
  52 }