PWG1/RunAlignmentDataFilterITS.C

   1 void RunAlignmentDataFilterITS() {
   2   //
   3   // Macro to extract AliTrackPoints for ITS
   4   // A.Dainese, andrea.dainese@pd.infn.it
   5   //
   6
   7   Int_t nentries=1234567890;
   8   Int_t firstentry=0;
   9
  10   // Load geometry file
  11   AliGeomManager::LoadGeometry("geometry.root");
  12
  13   // Load PWG1 library
  14   gSystem->Load("libANALYSIS.so");
  15   gSystem->Load("libANALYSISalice.so");
  16   gSystem->Load("libPWG1.so");
  17
  18
  19   // Create the analysis manager
  20   AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
  21
  22   // Create the task
  23   AliAlignmentDataFilterITS *taskFilter = new AliAlignmentDataFilterITS("filterITS");
  24   AliLog::SetClassDebugLevel("AliAlignmentDataFilterITS",10);
  25   // configuration via AliITSRecoParam (should be taken from OCDB)
  26   AliITSReconstructor *itsReconstructor = new AliITSReconstructor();
  27   AliITSRecoParam *itsRecoParam = AliITSRecoParam::GetCosmicTestParam();
  28   itsReconstructor->SetRecoParam(itsRecoParam);
  29
  30
  31   // Add ESD handler
  32   AliESDInputHandler *esdH = new AliESDInputHandler();
  33   mgr->SetInputEventHandler(esdH);
  34
  35   TChain *chainESD = new TChain("esdTree");
  36   chainESD->Add("AliESDs.root");
  37
  38   // Attach input
  39   cInput = mgr->CreateContainer("cInput",TChain::Class(),AliAnalysisManager::kInputContainer);
  40   //mgr->ConnectInput(taskFilter, 0, cInput); // v4-16-Release
  41   mgr->ConnectInput(taskFilter,0,mgr->GetCommonInputContainer());
  42
  43   // Attach output
  44   cOutput0= mgr->CreateContainer("cOutput0",TTree::Class(),AliAnalysisManager::kOutputContainer,"AliTrackPoints.root");
  45   mgr->ConnectOutput(taskFilter,0,cOutput0);
  46   cOutput1= mgr->CreateContainer("cOutput1",TList::Class(),AliAnalysisManager::kOutputContainer,"AliTrackPoints.root");
  47   mgr->ConnectOutput(taskFilter,1,cOutput1);
  48
  49   // Enable debug printouts
  50   mgr->SetDebugLevel(10);
  51
  52   // Run analysis
  53   mgr->InitAnalysis();
  54   mgr->PrintStatus();
  55   mgr->StartAnalysis("local",chainESD,nentries,firstentry);
  56
  57  return;
  58 }
  59