Unused variable i removed
[u/mrichter/AliRoot.git] / PMD / AliPMDSDigitsRead.C
CommitLineData
0709b99f 1// ----------------------------------------------------//
2// //
3// This is a macro to read PMD.SDigits.root //
4// //
5// ----------------------------------------------------//
6#include <Riostream.h>
7void AliPMDSDigitsRead(Int_t nevt = 1)
8{
9 TStopwatch timer;
10 timer.Start();
11
12
13 TH2F *h2 = new TH2F("h2","Y vs. X",200,-100.,100.,200,-100.,100.);
14 AliRunLoader *fRunLoader = AliRunLoader::Open("galice.root");
15 if (!fRunLoader)
16 {
17 Error("Open","Can not open session for file %s.",file);
18 }
19
20 // fRunLoader->LoadgAlice();
21 fRunLoader->LoadHeader();
22 gAlice = fRunLoader->GetAliRun();
23
24 if (gAlice)
25 {
26 printf("AliRun object found on file.\n");
27 }
28 else
29 {
30 printf("Could not find AliRun object.\n");
31 }
32 fPMD = (AliPMD*)gAlice->GetDetector("PMD");
33 fPMDLoader = fRunLoader->GetLoader("PMDLoader");
34 if (fPMDLoader == 0x0)
35 {
36 cerr<<"OpengAlice : Can not find PMD or PMDLoader\n";
37 }
38
39 fPMDLoader->LoadSDigits("READ");
40 TClonesArray *fSDigits;
41
42 // -------------------------------------------------------------- //
43
44 Int_t det = 0,smn = 0;
45 Int_t xpos, ypos;
46 Int_t xpad, ypad;
47 Float_t edep;
48 Float_t xx,yy;
49
50 AliPMDUtility cc;
51
52 for (Int_t ievt = 0; ievt <nevt; ievt++)
53 {
54 fRunLoader->GetEvent(ievt);
55 fTreeS = fPMDLoader->TreeS();
56 if (fTreeS == 0x0)
57 {
58 cout << " Can not get TreeD" << endl;
59 }
60 AliPMDsdigit *pmdsdigit;
61 TBranch *branch = fTreeS->GetBranch("PMDSDigit");
62 branch->SetAddress(&fSDigits);
63
64 Int_t nmodules = (Int_t) fTreeS->GetEntries();
65
66 for (Int_t imodule = 0; imodule < nmodules; imodule++)
67 {
68 fTreeS->GetEntry(imodule);
69 Int_t nentries = fSDigits->GetLast();
70 for (Int_t ient = 0; ient < nentries+1; ient++)
71 {
72 pmdsdigit = (AliPMDsdigit*)fSDigits->UncheckedAt(ient);
73
74 det = pmdsdigit->GetDetector();
75 smn = pmdsdigit->GetSMNumber();
76 xpos = pmdsdigit->GetRow();
77 ypos = pmdsdigit->GetColumn();
78 edep = pmdsdigit->GetCellEdep();
79 Int_t trno = pmdsdigit->GetTrackNumber();
80
81 if(smn <12)
82 {
83 xpad = ypos;
84 ypad = xpos;
85 }
86 else if(smn >=12 && smn < 24)
87 {
88 xpad = xpos;
89 ypad = ypos;
90 }
91
92 if(det == 1)
93 {
94 cc.RectGeomCellPos(smn,xpad,ypad,xx,yy);
95 h2->Fill(xx,yy);
96 }
97 }
98 } // modules
99
100 }
101 h2->Draw();
102
103 timer.Stop();
104 timer.Print();
105}
106