ITS/ITSHitsToFastPoints.C

   1 #include "iostream.h"
   2
   3 void ITSHitsToFastPoints (Int_t evNumber1=0,Int_t evNumber2=0,Int_t nsignal=25, Int_t size=-1)
   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    } else {
  18       delete gAlice;
  19       gAlice=0;
  20    }
  21
  22 // Connect the Root Galice file containing Geometry, Kine and Hits
  23
  24    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
  25    if (file) file->Close();
  26    file = new TFile("galice.root","UPDATE");
  27    file->ls();
  28
  29    printf ("I'm after Map \n");
  30
  31 // Get AliRun object from file or create it if not on file
  32
  33    if (!gAlice) {
  34        gAlice = (AliRun*)file->Get("gAlice");
  35        if (gAlice) printf("AliRun object found on file\n");
  36        if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  37    }
  38    printf ("I'm after gAlice \n");
  39
  40    AliITS *ITS  = (AliITS*) gAlice->GetModule("ITS");
  41    if (!ITS) return;
  42
  43    // Set the simulation model
  44    AliITSsimulationFastPoints *sim = new AliITSsimulationFastPoints();
  45
  46    for (Int_t i=0;i<3;i++) {
  47        ITS->SetSimulationModel(i,sim);
  48    }
  49
  50 //
  51 // Event Loop
  52 //
  53
  54    Int_t nbgr_ev=0;
  55   TStopwatch timer;
  56
  57    for (int ev=evNumber1; ev<= evNumber2; ev++) {
  58        Int_t nparticles = gAlice->GetEvent(ev);
  59        cout << "event         " <<ev<<endl;
  60        cout << "nparticles  " <<nparticles<<endl;
  61        gAlice->SetEvent(ev);
  62        if(!gAlice->TreeR()) gAlice-> MakeTree("R");
  63        ITS->MakeBranch("R");
  64        if (ev < evNumber1) continue;
  65        if (nparticles <= 0) return;
  66
  67        Int_t bgr_ev=Int_t(ev/nsignal);
  68        //printf("bgr_ev %d\n",bgr_ev);
  69        timer.Start();
  70        ITS->HitsToFastRecPoints(ev,bgr_ev,size," ","All"," ");
  71        timer.Stop(); timer.Print();
  72    } // event loop
  73
  74    //delete gAlice;
  75    //gAlice=0;
  76    file->Close();
  77 }
  78
  79
  80
  81
  82
  83
  84
  85
  86
  87
  88
  89
  90
  91