1 void AliTPCDigitsDisplay(int sec, int row,
2 int max_t_chan=500, float min_t=0., float max_t=500.,
3 int max_p_chan=200, float min_p=0., float max_p=200.)
5 // Dynamically link some shared libs
6 if (gClassTable->GetID("AliRun") < 0) {
7 gROOT->LoadMacro("loadlibs.C");
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("TreeD0_Param1");
21 t->GetBranch("Digits")->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,max_p_chan,min_p,max_p);
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 < 0) 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);