3 void ITSreadClustTestSPDDubna (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");
18 // Anyway, this macro needs to read a gAlice file, so it
19 // clears the gAlice object if there is already one in memory...
27 // Connect the Root Galice file containing Geometry, Kine and Hits
29 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
30 if (!file) file = new TFile("galice.root");
33 // Get AliRun object from file or create it if not on file
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");
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;
53 TTree *TH = gAlice->TreeH();
54 Int_t ntracks = TH->GetEntries();
55 cout<<"ntracks "<<ntracks<<endl;
59 //AliITSRawClusterSDD *ITSclust;
60 AliITSRawClusterSPD *ITSclust;
62 // Get pointers to Alice detectors and Digits containers
63 AliITS *ITS = (AliITS*)gAlice->GetModule("ITS");
65 TClonesArray *Particles = gAlice->Particles();
68 TTree *TC=ITS->TreeC();
69 Int_t nent=TC->GetEntries();
70 printf("Found %d entries in the tree (must be one per module per event!)\n",nent);
72 for (Int_t idet=0;idet<3;idet++) {
74 TClonesArray *ITSclu = ITS->ClustersAddress(idet);
76 if (idet != 0) continue;
77 for (Int_t mod=0; mod<nent; mod++) {
80 Int_t ncl = ITSclu->GetEntries();
81 if (ncl) printf("Found %d clusters for module %d in det type %d \n",ncl,mod,idet);
85 for (Int_t icl=0;icl<ncl;icl++) {
86 ITSclust = (AliITSRawClusterSPD*)ITSclu->UncheckedAt(icl);
87 printf("%d %d %f %f %f\n",ITSclust->NclZ(),ITSclust->NclX(),ITSclust->Q(),ITSclust->X(),ITSclust->Z());
89 //for (Int_t icl=0;icl<ncl;icl++) {
90 //ITSclust = (AliITSRawClusterSDD*)ITSclu->UncheckedAt(icl);
91 //printf("%d %d %f %f %f\n",ITSclust->Anodes(),ITSclust->Samples(),ITSclust->Q(),ITSclust->X(),ITSclust->Z());
98 cout<<"END test for clusters "<<endl;