1 void bcmhits(Int_t events = 10)
3 TH1::AddDirectory(kFALSE);
4 zH = new TH1F("zH", "z coordinate of hit", 210., -2100., 2100.);
5 rH = new TH1F("rH", "r coordinate of hit", 200., 0., 20.);
6 tH = new TH1F("tH", "time of hits", 100., 0., 100.);
8 // Creating Run Loader and openning file containing Hits
9 TClonesArray* hits = new TClonesArray ("AliBCMHit", 1000);
10 AliRunLoader* rl = AliRunLoader::Open("galice.root");
12 printf(">>> Error : Error Opening file \n");
18 gAlice = rl->GetAliRun();
19 AliBCM* pBCM = (AliBCM*) gAlice->GetDetector("BCM");
22 for (Int_t nev = 0; nev < events; nev++) {
25 AliLoader* bcmL = rl->GetLoader("BCMLoader");
26 // bcmL->InitDefaults();
28 TTree* treeH = bcmL->TreeH();
29 Int_t ntracks = (Int_t) treeH->GetEntries();
30 printf("Number of particles %6d %6d \n", nev, ntracks);
31 for (Int_t itr = 0; itr < ntracks; itr++) {
35 Int_t nhits = hits->GetEntriesFast();
36 for (Int_t hit = 0; hit < nhits; hit++)
39 AliBCMHit* mHit = (AliBCMHit *) hits->UncheckedAt(hit);
40 Float_t x = mHit->X();
41 Float_t y = mHit->Y();
42 Float_t z = mHit->Z();
43 Float_t r = TMath::Sqrt(x * x + y * y);
46 tH->Fill(mHit->Time() * 1.e9);