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