3 void ITSreadRecPointsTest (Int_t evNumber1=0,Int_t evNumber2=0)
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.
10 /////////////////////////////////////////////////////////////////////////
12 // Dynamically link some shared libs
14 if (gClassTable->GetID("AliRun") < 0) {
15 gROOT->LoadMacro("loadlibs.C");
19 // Connect the Root Galice file containing Geometry, Kine and Hits
21 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
22 if (!file) file = new TFile("galice.root");
25 // Get AliRun object from file or create it if not on file
28 gAlice = (AliRun*)file->Get("gAlice");
29 if (gAlice) printf("AliRun object found on file\n");
30 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
38 for (int nev=0; nev<= evNumber2; nev++) {
39 Int_t nparticles = gAlice->GetEvent(nev);
40 cout << "nev " << nev <<endl;
41 cout << "nparticles " << nparticles <<endl;
42 if (nev < evNumber1) continue;
43 if (nparticles <= 0) return;
45 TTree *TH = gAlice->TreeH();
46 Int_t ntracks = TH->GetEntries();
47 cout<<"ntracks "<<ntracks<<endl;
53 // Get pointers to Alice detectors and Digits containers
54 AliITS *ITS = (AliITS*)gAlice->GetModule("ITS");
56 TClonesArray *Particles = gAlice->Particles();
57 TTree *TR = gAlice->TreeR();
58 Int_t nent=TR->GetEntries();
59 printf("Found %d entries in the tree (must be one per module per event!)\n",nent);
61 TClonesArray *ITSrec = ITS->RecPoints();
62 for (Int_t mod=0; mod<nent; mod++) {
63 ITS->ResetRecPoints();
65 Int_t nrecp = ITSrec->GetEntries();
66 if (nrecp) printf("Found %d rec points for module %d \n",nrecp,mod);
69 for (Int_t irec=0;irec<nrecp;irec++) {
70 recp = (AliITSRecPoint*)ITSrec->UncheckedAt(irec);
71 printf("%d %f %f %d %d %d\n",irec,recp->GetX(),recp->GetZ(),recp->fTracks[0],recp->fTracks[1],recp->fTracks[2]);
77 cout<<"END test for rec points "<<endl;