ITS/ITSreadRecPointsTest.C

   1 #include "iostream.h"
   2
   3 void ITSreadRecPointsTest (Int_t evNumber1=0,Int_t evNumber2=0)
   4 {
   5 /////////////////////////////////////////////////////////////////////////
   6 //   This macro is a small example of a ROOT macro
   7 //   illustrating how to read the output of GALICE
   8 //   and do some analysis.
   9 //
  10 /////////////////////////////////////////////////////////////////////////
  11
  12 // Dynamically link some shared libs
  13
  14    if (gClassTable->GetID("AliRun") < 0) {
  15       gROOT->LoadMacro("loadlibs.C");
  16       loadlibs();
  17    }
  18      // Anyway, this macro needs to read a gAlice file, so it
  19   // clears the gAlice object if there is already one in memory...
  20   else {
  21                 if(gAlice){
  22                         delete gAlice;
  23                         gAlice = 0;
  24                 }
  25         }
  26
  27 // Connect the Root Galice file containing Geometry, Kine and Hits
  28
  29    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
  30    if (!file) file = new TFile("galice.root");
  31    file->ls();
  32
  33 // Get AliRun object from file or create it if not on file
  34
  35    if (!gAlice) {
  36       gAlice = (AliRun*)file->Get("gAlice");
  37       if (gAlice) printf("AliRun object found on file\n");
  38       if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  39    }
  40
  41 //
  42 //   Loop over events
  43 //
  44    Int_t Nh=0;
  45    Int_t Nh1=0;
  46    for (int nev=0; nev<= evNumber2; nev++) {
  47      Int_t nparticles = gAlice->GetEvent(nev);
  48      cout << "nev         " << nev <<endl;
  49      cout << "nparticles  " << nparticles <<endl;
  50      if (nev < evNumber1) continue;
  51      if (nparticles <= 0) return;
  52
  53      TTree *TH = gAlice->TreeH();
  54      Int_t ntracks = TH->GetEntries();
  55      cout<<"ntracks "<<ntracks<<endl;
  56
  57      Int_t nbytes = 0;
  58
  59      AliITSRecPoint  *recp;
  60
  61      // Get pointers to Alice detectors and Digits containers
  62      AliITS *ITS  = (AliITS*)gAlice->GetModule("ITS");
  63      if (!ITS) return;
  64      TClonesArray *Particles = gAlice->Particles();
  65      TTree *TR = gAlice->TreeR();
  66      Int_t nent=TR->GetEntries();
  67      printf("Found %d entries in the tree (must be one per module per event!)\n",nent);
  68
  69      TClonesArray *ITSrec  = ITS->RecPoints();
  70      for (Int_t mod=0; mod<nent; mod++) {
  71        ITS->ResetRecPoints();
  72        TR->GetEvent(mod);
  73        Int_t nrecp = ITSrec->GetEntries();
  74        if (nrecp) printf("Found %d rec points for module %d \n",nrecp,mod);
  75        if (!nrecp) continue;
  76
  77        for (Int_t irec=0;irec<nrecp;irec++) {
  78                 recp   = (AliITSRecPoint*)ITSrec->UncheckedAt(irec);
  79                 printf("%d %f %f %d %d %d\n",irec,recp->GetX(),recp->GetZ(),recp->GetLabel(0),recp->GetLabel(1),recp->GetLabel(2));
  80
  81        }
  82      }
  83    }   // event loop
  84
  85      cout<<"END  test for rec points "<<endl;
  86
  87      file->Close();
  88 }
  89
  90
  91