ITS/AliITSDigits2RecPoints.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2
   3 #include <Riostream.h>
   4 #include <TClassTable.h>
   5 #include <TDatime.h>
   6 #include <TGeoManager.h>
   7 #include <TROOT.h>
   8 #include <TStopwatch.h>
   9 #include "AliRun.h"
  10 #include "AliRunDigitizer.h"
  11 #include "AliITSDetTypeRec.h"
  12 #include "AliITS.h"
  13 #include "AliITSresponseSDD.h"
  14 #include "AliITSreconstruction.h"
  15
  16 #endif
  17 #define DEBUG
  18
  19 Int_t AliITSDigits2RecPoints(TString filename="galice.root",TString fileRP=""){
  20   // Standard ITS Digits to RecPoints.
  21
  22   // Get geometry
  23   TGeoManager::Import("geometry.root");
  24
  25   // Dynamically link some shared libs
  26   if (gClassTable->GetID("AliRun") < 0) {
  27     gROOT->ProcessLine(".x $(ALICE_ROOT)/macros/loadlibs.C");
  28   }else if (gAlice){
  29     delete AliRunLoader::Instance();
  30     delete gAlice;
  31     gAlice=0;
  32   } // end if
  33
  34   TStopwatch timer;
  35 #ifdef DEBUG
  36   cout << "Creating reconstructed points from digits for the ITS..." << endl;
  37 #endif
  38   AliITSreconstruction *itsr = new AliITSreconstruction(filename);
  39
  40   timer.Start();
  41   if(!(fileRP.IsNull()))itsr->SetOutputFile(fileRP);
  42   itsr->Init();
  43   itsr->Exec();
  44   timer.Stop();
  45   timer.Print();
  46   delete itsr;
  47   return 0;
  48 }