1 void TPCDigitsDisplay(int sec, int row,
2 int max_t_chan=500, float min_t=0., float max_t=500.)
4 // Dynamically link some shared libs
5 if (gClassTable->GetID("AliRun") < 0) {
6 gSystem->Load("libGeant3Dummy.so"); // a dummy version of Geant3
7 gSystem->Load("PHOS/libPHOSdummy.so"); // the standard Alice classes
8 gSystem->Load("libgalice.so"); // the standard Alice classes
14 // Connect the Root Galice file containing Geometry, Kine and Hits
15 TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
17 f = new TFile("galice.root");
19 TClonesArray *fDigits=new TClonesArray("AliTPCdigit",10000);
20 TTree *t=(TTree*)f->Get("TreeD1");
21 t->GetBranch("TPC")->SetAddress(&fDigits);
22 Int_t sectors_by_rows=(Int_t)t->GetEntries();
23 for (Int_t n=0; n<sectors_by_rows; n++) {
24 if (!t->GetEvent(n)) continue;
25 AliTPCdigit *dig=(AliTPCdigit*)fDigits->UncheckedAt(0);
27 if (sec < dig->fSector) break;
28 if (sec != dig->fSector) continue;
29 if (row != dig->fPadRow) continue;
32 sprintf(s,"Sector %d Row %d\n",sec,row);
33 TH2F *h = new TH2F("h",s,max_t_chan,min_t,max_t,200,0,200);
34 Int_t ndigits=fDigits->GetEntriesFast();
35 for (Int_t ndig=0; ndig<ndigits; ndig++) {
36 dig=(AliTPCdigit*)fDigits->UncheckedAt(ndig);
37 if (dig->fSignal < 10) continue; //cluster finder threshold
38 h->Fill(dig->fTime,dig->fPad,dig->fSignal);
41 gStyle->SetOptStat(0);
42 TCanvas *c1=new TCanvas("c1","TPC digits display",0,0,1110,680);
43 TPad *p1=new TPad("p1","",0,0,1,0.5);
45 TPad *p2=new TPad("p2","",0,0.5,1,1);