1 void AliITSPrintHits(TString hfn="galice.root",Int_t mod=-1,
3 // Macro to print out the recpoints for all or a specific module
5 // Dynamically link some shared libs
6 if (gClassTable->GetID("AliRun") < 0) {
7 gROOT->LoadMacro("loadlibs.C");
12 delete gAlice->GetRunLoader();
17 gROOT->LoadMacro("$(ALICE_ROOT)/ITS/AliITSstandard.C");
19 AliRunLoader *rl = AccessFile(hfn); // Set up to read in Data
20 Int_t retval = rl->LoadHeader();
22 cerr<<"AliITSPrintHits.C : LoadHeader returned error"<<endl;
26 AliITSLoader* ITSloader = (AliITSLoader*) rl->GetLoader("ITSLoader");
29 cerr<<"AliITSPrintHits.C : ITS loader not found"<<endl;
33 ITSloader->LoadHits("read");
34 AliITS *ITS = (AliITS*)gAlice->GetDetector("ITS");
36 cout << "Error: no ITS found. Aborting"<<endl;
41 Int_t evNumber2 = gAlice->GetEventsPerRun();
47 Int_t mod2 = ITS->GetITSgeom()->GetIndexMax();
54 Int_t nmodules,size=-1;
55 Int_t event,m,i,i2,hit,trk;
56 for(event = evNumber1; event < evNumber2; event++){
57 cout<<"Processing event "<<event<<endl;
59 ITS->InitModules(size,nmodules);
60 ITS->FillModules(event,0,-1," "," ");
61 for(m=mod1;m<mod2;m++){
62 i2 = (ITS->GetModule(m))->GetNhits();
63 cout << "Event=" << event << " module=" << m <<
64 " Number of Hits=" << i2 <<endl;
66 trk = (ITS->GetModule(m))->GetHitTrackIndex(i);
67 hit = (ITS->GetModule(m))->GetHitHitIndex(i);
68 hp = (ITS->GetModule(m))->GetHit(i);
69 cout << i << " trk#="<<trk<<" hit#="<< hit << " ";
70 hp->Print((ostream*)cout);