StRoot  1
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
plotPadOcc.C
1 // Plot pad occupancy for all clusters in all sectors which pass the
2 // cuts (3<=padlength<=6 and 4<=timelength<=8) for a "good" cluster
3 
4 // structures definitions
5 struct HIT
6 {
7  Float_t x,y,z;
8  Float_t rad,phi;
9  Float_t raderror,phierror;
10 };
11 
12 struct CLUSTER
13 {
14  Float_t timepos,padpos,timesigma,padsigma;
15  Float_t peakheight, charge;
16  Int_t timebin,pad;
17  Int_t padlength,timelength;
18  Int_t row,sec;
19  Int_t flag;
20  Int_t numpeaks;
21 };
22 
23 struct EVENT
24 {
25  Float_t run;
26  Int_t nevent;
27  Float_t temperature,pressure;
28 };
29 
30 // =====================================
31 
32 void plotPadOcc(const Char_t *fileName="filelist.txt")
33 {
34 
35  CLUSTER cluster;
36  HIT hit;
37  EVENT event;
38 
39  // plain Style
40  gROOT->Reset();
41  gStyle->SetTitleOffset(1.25);
42  gStyle->SetCanvasBorderMode(0);
43  gStyle->SetPadBorderMode(0);
44  gStyle->SetPadColor(0);
45  gStyle->SetCanvasColor(0);
46  gStyle->SetTitleColor(0);
47  gStyle->SetStatColor(0);
48  gStyle->SetPalette(1);
49  gStyle->SetOptStat(1111);
50  gStyle->SetOptFit();
51 
52  Int_t maxentries;
53 
54  Int_t secrow[61][3][161];
55 
56  for (Int_t isec=1;isec<61;isec++) {
57  for (Int_t irow=1;irow<3;irow++) {
58  for (Int_t ipad=1;ipad<161;ipad++) {
59  secrow[isec][irow][ipad]=0;
60  }
61  }
62  }
63 
64 // Ftpc West
65  sec1row1 = new TH1I("hsec1_hrow1","hsec1_hrow1",160,0,160);
66  sec2row1 = new TH1I("hsec2_hrow1","hsec2_hrow1",160,0,160);
67  sec3row1 = new TH1I("hsec3_hrow1","hsec3_hrow1",160,0,160);
68  sec4row1 = new TH1I("hsec4_hrow1","hsec4_hrow1",160,0,160);
69  sec5row1 = new TH1I("hsec5_hrow1","hsec5_hrow1",160,0,160);
70  sec6row1 = new TH1I("hsec6_hrow1","hsec6_hrow1",160,0,160);
71  sec1row2 = new TH1I("hsec1_hrow2","hsec1_hrow2",160,0,160);
72  sec2row2 = new TH1I("hsec2_hrow2","hsec2_hrow2",160,0,160);
73  sec3row2 = new TH1I("hsec3_hrow2","hsec3_hrow2",160,0,160);
74  sec4row2 = new TH1I("hsec4_hrow2","hsec4_hrow2",160,0,160);
75  sec5row2 = new TH1I("hsec5_hrow2","hsec5_hrow2",160,0,160);
76  sec6row2 = new TH1I("hsec6_hrow2","hsec6_hrow2",160,0,160);
77 
78  sec7row1 = new TH1I("hsec7_hrow1","hsec7_hrow1",160,0,160);
79  sec8row1 = new TH1I("hsec8_hrow1","hsec8_hrow1",160,0,160);
80  sec9row1 = new TH1I("hsec9_hrow1","hsec9_hrow1",160,0,160);
81  sec10row1 = new TH1I("hsec10_hrow1","hsec10_hrow1",160,0,160);
82  sec11row1 = new TH1I("hsec11_hrow1","hsec11_hrow1",160,0,160);
83  sec12row1 = new TH1I("hsec12_hrow1","hsec12_hrow1",160,0,160);
84  sec7row2 = new TH1I("hsec7_hrow2","hsec7_hrow2",160,0,160);
85  sec8row2 = new TH1I("hsec8_hrow2","hsec8_hrow2",160,0,160);
86  sec9row2 = new TH1I("hsec9_hrow2","hsec9_hrow2",160,0,160);
87  sec10row2 = new TH1I("hsec10_hrow2","hsec10_hrow2",160,0,160);
88  sec11row2 = new TH1I("hsec11_hrow2","hsec11_hrow2",160,0,160);
89  sec12row2 = new TH1I("hsec12_hrow2","hsec12_hrow2",160,0,160);
90 
91  sec13row1 = new TH1I("hsec13_hrow1","hsec13_hrow1",160,0,160);
92  sec14row1 = new TH1I("hsec14_hrow1","hsec14_hrow1",160,0,160);
93  sec15row1 = new TH1I("hsec15_hrow1","hsec15_hrow1",160,0,160);
94  sec16row1 = new TH1I("hsec16_hrow1","hsec16_hrow1",160,0,160);
95  sec17row1 = new TH1I("hsec17_hrow1","hsec17_hrow1",160,0,160);
96  sec18row1 = new TH1I("hsec18_hrow1","hsec18_hrow1",160,0,160);
97  sec13row2 = new TH1I("hsec13_hrow2","hsec13_hrow2",160,0,160);
98  sec14row2 = new TH1I("hsec14_hrow2","hsec14_hrow2",160,0,160);
99  sec15row2 = new TH1I("hsec15_hrow2","hsec15_hrow2",160,0,160);
100  sec16row2 = new TH1I("hsec16_hrow2","hsec16_hrow2",160,0,160);
101  sec17row2 = new TH1I("hsec17_hrow2","hsec17_hrow2",160,0,160);
102  sec18row2 = new TH1I("hsec18_hrow2","hsec18_hrow2",160,0,160);
103 
104  sec19row1 = new TH1I("hsec19_hrow1","hsec19_hrow1",160,0,160);
105  sec20row1 = new TH1I("hsec20_hrow1","hsec20_hrow1",160,0,160);
106  sec21row1 = new TH1I("hsec21_hrow1","hsec21_hrow1",160,0,160);
107  sec22row1 = new TH1I("hsec22_hrow1","hsec22_hrow1",160,0,160);
108  sec23row1 = new TH1I("hsec23_hrow1","hsec23_hrow1",160,0,160);
109  sec24row1 = new TH1I("hsec24_hrow1","hsec24_hrow1",160,0,160);
110  sec19row2 = new TH1I("hsec19_hrow2","hsec19_hrow2",160,0,160);
111  sec20row2 = new TH1I("hsec20_hrow2","hsec20_hrow2",160,0,160);
112  sec21row2 = new TH1I("hsec21_hrow2","hsec21_hrow2",160,0,160);
113  sec22row2 = new TH1I("hsec22_hrow2","hsec22_hrow2",160,0,160);
114  sec23row2 = new TH1I("hsec23_hrow2","hsec23_hrow2",160,0,160);
115  sec24row2 = new TH1I("hsec24_hrow2","hsec24_hrow2",160,0,160);
116 
117  sec25row1 = new TH1I("hsec25_hrow1","hsec25_hrow1",160,0,160);
118  sec26row1 = new TH1I("hsec26_hrow1","hsec26_hrow1",160,0,160);
119  sec27row1 = new TH1I("hsec27_hrow1","hsec27_hrow1",160,0,160);
120  sec28row1 = new TH1I("hsec28_hrow1","hsec28_hrow1",160,0,160);
121  sec29row1 = new TH1I("hsec29_hrow1","hsec29_hrow1",160,0,160);
122  sec30row1 = new TH1I("hsec30_hrow1","hsec30_hrow1",160,0,160);
123  sec25row2 = new TH1I("hsec25_hrow2","hsec25_hrow2",160,0,160);
124  sec26row2 = new TH1I("hsec26_hrow2","hsec26_hrow2",160,0,160);
125  sec27row2 = new TH1I("hsec27_hrow2","hsec27_hrow2",160,0,160);
126  sec28row2 = new TH1I("hsec28_hrow2","hsec28_hrow2",160,0,160);
127  sec29row2 = new TH1I("hsec29_hrow2","hsec29_hrow2",160,0,160);
128  sec30row2 = new TH1I("hsec30_hrow2","hsec30_hrow2",160,0,160);
129 
130 // Ftpc East
131  sec31row1 = new TH1I("hsec31_hrow1","hsec31_hrow1",160,0,160);
132  sec32row1 = new TH1I("hsec32_hrow1","hsec32_hrow1",160,0,160);
133  sec33row1 = new TH1I("hsec33_hrow1","hsec33_hrow1",160,0,160);
134  sec34row1 = new TH1I("hsec34_hrow1","hsec34_hrow1",160,0,160);
135  sec35row1 = new TH1I("hsec35_hrow1","hsec35_hrow1",160,0,160);
136  sec36row1 = new TH1I("hsec36_hrow1","hsec36_hrow1",160,0,160);
137  sec31row2 = new TH1I("hsec31_hrow2","hsec31_hrow2",160,0,160);
138  sec32row2 = new TH1I("hsec32_hrow2","hsec32_hrow2",160,0,160);
139  sec33row2 = new TH1I("hsec33_hrow2","hsec33_hrow2",160,0,160);
140  sec34row2 = new TH1I("hsec34_hrow2","hsec34_hrow2",160,0,160);
141  sec35row2 = new TH1I("hsec35_hrow2","hsec35_hrow2",160,0,160);
142  sec36row2 = new TH1I("hsec36_hrow2","hsec36_hrow2",160,0,160);
143 
144  sec37row1 = new TH1I("hsec37_hrow1","hsec37_hrow1",160,0,160);
145  sec38row1 = new TH1I("hsec38_hrow1","hsec38_hrow1",160,0,160);
146  sec39row1 = new TH1I("hsec39_hrow1","hsec39_hrow1",160,0,160);
147  sec40row1 = new TH1I("hsec40_hrow1","hsec40_hrow1",160,0,160);
148  sec41row1 = new TH1I("hsec41_hrow1","hsec41_hrow1",160,0,160);
149  sec42row1 = new TH1I("hsec42_hrow1","hsec42_hrow1",160,0,160);
150  sec37row2 = new TH1I("hsec37_hrow2","hsec37_hrow2",160,0,160);
151  sec38row2 = new TH1I("hsec38_hrow2","hsec38_hrow2",160,0,160);
152  sec39row2 = new TH1I("hsec39_hrow2","hsec39_hrow2",160,0,160);
153  sec40row2 = new TH1I("hsec40_hrow2","hsec40_hrow2",160,0,160);
154  sec41row2 = new TH1I("hsec41_hrow2","hsec41_hrow2",160,0,160);
155  sec42row2 = new TH1I("hsec42_hrow2","hsec42_hrow2",160,0,160);
156 
157  sec43row1 = new TH1I("hsec43_hrow1","hsec43_hrow1",160,0,160);
158  sec44row1 = new TH1I("hsec44_hrow1","hsec44_hrow1",160,0,160);
159  sec45row1 = new TH1I("hsec45_hrow1","hsec45_hrow1",160,0,160);
160  sec46row1 = new TH1I("hsec46_hrow1","hsec46_hrow1",160,0,160);
161  sec47row1 = new TH1I("hsec47_hrow1","hsec47_hrow1",160,0,160);
162  sec48row1 = new TH1I("hsec48_hrow1","hsec48_hrow1",160,0,160);
163  sec43row2 = new TH1I("hsec43_hrow2","hsec43_hrow2",160,0,160);
164  sec44row2 = new TH1I("hsec44_hrow2","hsec44_hrow2",160,0,160);
165  sec45row2 = new TH1I("hsec45_hrow2","hsec45_hrow2",160,0,160);
166  sec46row2 = new TH1I("hsec46_hrow2","hsec46_hrow2",160,0,160);
167  sec47row2 = new TH1I("hsec47_hrow2","hsec47_hrow2",160,0,160);
168  sec48row2 = new TH1I("hsec48_hrow2","hsec48_hrow2",160,0,160);
169 
170  sec49row1 = new TH1I("hsec49_hrow1","hsec49_hrow1",160,0,160);
171  sec50row1 = new TH1I("hsec50_hrow1","hsec50_hrow1",160,0,160);
172  sec51row1 = new TH1I("hsec51_hrow1","hsec51_hrow1",160,0,160);
173  sec52row1 = new TH1I("hsec52_hrow1","hsec52_hrow1",160,0,160);
174  sec53row1 = new TH1I("hsec53_hrow1","hsec53_hrow1",160,0,160);
175  sec54row1 = new TH1I("hsec54_hrow1","hsec54_hrow1",160,0,160);
176  sec49row2 = new TH1I("hsec49_hrow2","hsec49_hrow2",160,0,160);
177  sec50row2 = new TH1I("hsec50_hrow2","hsec50_hrow2",160,0,160);
178  sec51row2 = new TH1I("hsec51_hrow2","hsec51_hrow2",160,0,160);
179  sec52row2 = new TH1I("hsec52_hrow2","hsec52_hrow2",160,0,160);
180  sec53row2 = new TH1I("hsec53_hrow2","hsec53_hrow2",160,0,160);
181  sec54row2 = new TH1I("hsec54_hrow2","hsec54_hrow2",160,0,160);
182 
183  sec55row1 = new TH1I("hsec55_hrow1","hsec55_hrow1",160,0,160);
184  sec56row1 = new TH1I("hsec56_hrow1","hsec56_hrow1",160,0,160);
185  sec57row1 = new TH1I("hsec57_hrow1","hsec57_hrow1",160,0,160);
186  sec58row1 = new TH1I("hsec58_hrow1","hsec58_hrow1",160,0,160);
187  sec59row1 = new TH1I("hsec59_hrow1","hsec59_hrow1",160,0,160);
188  sec60row1 = new TH1I("hsec60_hrow1","hsec60_hrow1",160,0,160);
189  sec55row2 = new TH1I("hsec55_hrow2","hsec55_hrow2",160,0,160);
190  sec56row2 = new TH1I("hsec56_hrow2","hsec56_hrow2",160,0,160);
191  sec57row2 = new TH1I("hsec57_hrow2","hsec57_hrow2",160,0,160);
192  sec58row2 = new TH1I("hsec58_hrow2","hsec58_hrow2",160,0,160);
193  sec59row2 = new TH1I("hsec59_hrow2","hsec59_hrow2",160,0,160);
194  sec60row2 = new TH1I("hsec60_hrow2","hsec60_hrow2",160,0,160);
195 
196 
197  // input
198  cout<<"Histogram the number of good clusters per pad for each sector"<<endl;
199 
200 //++++++++++++++++++++++++++++++++++++++++++++++++++
201 // read file list from text file (maximum of 200 files)
202  char **fList = new char[200][256];
203  int nfile=0;
204  ifstream fin(fileName);
205  if (fin) {
206  int j=0;
207  cout << "Input file " << fileName << " found." << endl;
208  for (Int_t i=0;i<200;i++){
209  char *line = new char[256];
210  fin >> line;
211  if (fin.eof()) break;
212  fList[j] = line;
213  cout << " Added file : " << fList[j] << endl;
214  j++;
215  }
216  nfile = j;
217  }
218  fin.close();
219 
220 //++++++++++++++++++++++++++++++++++++++++++++++++++
221 
222  for (Int_t ifile=0;ifile<nfile;ifile++) {
223  TString eingabe = fList[ifile];
224  TFile f(eingabe);
225 
226  // read in ntuple
227  // ==============
228 
229  TTree *dtree=(TTree*) f.Get("cl");
230  //TBranch *bhit=dtree->GetBranch("hit");
231  //bhit->SetAddress(&hit);
232  TBranch *bcluster=dtree->GetBranch("cluster");
233  bcluster->SetAddress(&cluster);
234  TBranch *bevent=dtree->GetBranch("event");
235  bevent->SetAddress(&event);
236 
237  bevent->GetEntry(0);
238  //cout<<"Event number "<<event.nevent<<endl;
239  Int_t CurrentEvent = event.nevent;
240  maxentries=bcluster->GetEntries();
241 
242  for (int i=0;i<=maxentries;i++) {
243  bevent->GetEntry(i);
244  if (event.nevent != CurrentEvent){
245  //cout<<"Event number "<<event.nevent<<endl;
246  CurrentEvent = event.nevent;
247  }
248  bcluster->GetEntry(i);
249 
250  // Cut on #pad and #timebins for "good" clusters
251  if ((cluster.padlength>=3&&cluster.padlength<=6) &&
252  (cluster.timelength>=4&&cluster.timelength<=8)){
253  secrow[cluster.sec][cluster.row][cluster.pad]++;
254  } // End of cut loop
255  } // End of event loop
256  } // End of input file loop
257 
258  for (Int_t i=0;i<160;i++) {
259  //FTPC West
260  sec1row1->SetBinContent(i,secrow[1][1][i]);
261  sec2row1->SetBinContent(i,secrow[2][1][i]);
262  sec3row1->SetBinContent(i,secrow[3][1][i]);
263  sec4row1->SetBinContent(i,secrow[4][1][i]);
264  sec5row1->SetBinContent(i,secrow[5][1][i]);
265  sec6row1->SetBinContent(i,secrow[6][1][i]);
266  sec1row2->SetBinContent(i,secrow[1][2][i]);
267  sec2row2->SetBinContent(i,secrow[2][2][i]);
268  sec3row2->SetBinContent(i,secrow[3][2][i]);
269  sec4row2->SetBinContent(i,secrow[4][2][i]);
270  sec5row2->SetBinContent(i,secrow[5][2][i]);
271  sec6row2->SetBinContent(i,secrow[6][2][i]);
272 
273  sec7row1->SetBinContent(i,secrow[7][1][i]);
274  sec8row1->SetBinContent(i,secrow[8][1][i]);
275  sec9row1->SetBinContent(i,secrow[9][1][i]);
276  sec10row1->SetBinContent(i,secrow[10][1][i]);
277  sec11row1->SetBinContent(i,secrow[11][1][i]);
278  sec12row1->SetBinContent(i,secrow[12][1][i]);
279  sec7row2->SetBinContent(i,secrow[7][2][i]);
280  sec8row2->SetBinContent(i,secrow[8][2][i]);
281  sec9row2->SetBinContent(i,secrow[8][2][i]);
282  sec10row2->SetBinContent(i,secrow[10][2][i]);
283  sec11row2->SetBinContent(i,secrow[11][2][i]);
284  sec12row2->SetBinContent(i,secrow[12][2][i]);
285 
286  sec13row1->SetBinContent(i,secrow[13][1][i]);
287  sec14row1->SetBinContent(i,secrow[14][1][i]);
288  sec15row1->SetBinContent(i,secrow[15][1][i]);
289  sec16row1->SetBinContent(i,secrow[16][1][i]);
290  sec17row1->SetBinContent(i,secrow[17][1][i]);
291  sec18row1->SetBinContent(i,secrow[18][1][i]);
292  sec13row2->SetBinContent(i,secrow[13][2][i]);
293  sec14row2->SetBinContent(i,secrow[14][2][i]);
294  sec15row2->SetBinContent(i,secrow[15][2][i]);
295  sec16row2->SetBinContent(i,secrow[16][2][i]);
296  sec17row2->SetBinContent(i,secrow[17][2][i]);
297  sec18row2->SetBinContent(i,secrow[18][2][i]);
298 
299  sec19row1->SetBinContent(i,secrow[19][1][i]);
300  sec20row1->SetBinContent(i,secrow[20][1][i]);
301  sec21row1->SetBinContent(i,secrow[21][1][i]);
302  sec22row1->SetBinContent(i,secrow[22][1][i]);
303  sec23row1->SetBinContent(i,secrow[23][1][i]);
304  sec24row1->SetBinContent(i,secrow[24][1][i]);
305  sec19row2->SetBinContent(i,secrow[19][2][i]);
306  sec20row2->SetBinContent(i,secrow[20][2][i]);
307  sec21row2->SetBinContent(i,secrow[21][2][i]);
308  sec22row2->SetBinContent(i,secrow[22][2][i]);
309  sec23row2->SetBinContent(i,secrow[23][2][i]);
310  sec24row2->SetBinContent(i,secrow[24][2][i]);
311 
312  sec25row1->SetBinContent(i,secrow[25][1][i]);
313  sec26row1->SetBinContent(i,secrow[26][1][i]);
314  sec27row1->SetBinContent(i,secrow[27][1][i]);
315  sec28row1->SetBinContent(i,secrow[28][1][i]);
316  sec29row1->SetBinContent(i,secrow[29][1][i]);
317  sec30row1->SetBinContent(i,secrow[30][1][i]);
318  sec25row2->SetBinContent(i,secrow[25][2][i]);
319  sec26row2->SetBinContent(i,secrow[26][2][i]);
320  sec27row2->SetBinContent(i,secrow[27][2][i]);
321  sec28row2->SetBinContent(i,secrow[28][2][i]);
322  sec29row2->SetBinContent(i,secrow[29][2][i]);
323  sec30row2->SetBinContent(i,secrow[30][2][i]);
324  //FTPC East
325  sec31row1->SetBinContent(i,secrow[31][1][i]);
326  sec32row1->SetBinContent(i,secrow[32][1][i]);
327  sec33row1->SetBinContent(i,secrow[33][1][i]);
328  sec34row1->SetBinContent(i,secrow[34][1][i]);
329  sec35row1->SetBinContent(i,secrow[35][1][i]);
330  sec36row1->SetBinContent(i,secrow[36][1][i]);
331  sec31row2->SetBinContent(i,secrow[31][2][i]);
332  sec32row2->SetBinContent(i,secrow[32][2][i]);
333  sec33row2->SetBinContent(i,secrow[33][2][i]);
334  sec34row2->SetBinContent(i,secrow[34][2][i]);
335  sec35row2->SetBinContent(i,secrow[35][2][i]);
336  sec36row2->SetBinContent(i,secrow[36][2][i]);
337 
338  sec37row1->SetBinContent(i,secrow[37][1][i]);
339  sec38row1->SetBinContent(i,secrow[38][1][i]);
340  sec39row1->SetBinContent(i,secrow[39][1][i]);
341  sec40row1->SetBinContent(i,secrow[40][1][i]);
342  sec41row1->SetBinContent(i,secrow[41][1][i]);
343  sec42row1->SetBinContent(i,secrow[42][1][i]);
344  sec37row2->SetBinContent(i,secrow[37][2][i]);
345  sec38row2->SetBinContent(i,secrow[38][2][i]);
346  sec39row2->SetBinContent(i,secrow[39][2][i]);
347  sec40row2->SetBinContent(i,secrow[40][2][i]);
348  sec41row2->SetBinContent(i,secrow[41][2][i]);
349  sec42row2->SetBinContent(i,secrow[42][2][i]);
350 
351  sec43row1->SetBinContent(i,secrow[43][1][i]);
352  sec44row1->SetBinContent(i,secrow[44][1][i]);
353  sec45row1->SetBinContent(i,secrow[45][1][i]);
354  sec46row1->SetBinContent(i,secrow[46][1][i]);
355  sec47row1->SetBinContent(i,secrow[47][1][i]);
356  sec48row1->SetBinContent(i,secrow[48][1][i]);
357  sec43row2->SetBinContent(i,secrow[43][2][i]);
358  sec44row2->SetBinContent(i,secrow[44][2][i]);
359  sec45row2->SetBinContent(i,secrow[45][2][i]);
360  sec46row2->SetBinContent(i,secrow[46][2][i]);
361  sec47row2->SetBinContent(i,secrow[47][2][i]);
362  sec48row2->SetBinContent(i,secrow[48][2][i]);
363 
364  sec49row1->SetBinContent(i,secrow[49][1][i]);
365  sec50row1->SetBinContent(i,secrow[50][1][i]);
366  sec51row1->SetBinContent(i,secrow[51][1][i]);
367  sec52row1->SetBinContent(i,secrow[52][1][i]);
368  sec53row1->SetBinContent(i,secrow[53][1][i]);
369  sec54row1->SetBinContent(i,secrow[54][1][i]);
370  sec49row2->SetBinContent(i,secrow[49][2][i]);
371  sec50row2->SetBinContent(i,secrow[50][2][i]);
372  sec51row2->SetBinContent(i,secrow[51][2][i]);
373  sec52row2->SetBinContent(i,secrow[52][2][i]);
374  sec53row2->SetBinContent(i,secrow[53][2][i]);
375  sec54row2->SetBinContent(i,secrow[54][2][i]);
376 
377  sec55row1->SetBinContent(i,secrow[55][1][i]);
378  sec56row1->SetBinContent(i,secrow[56][1][i]);
379  sec57row1->SetBinContent(i,secrow[57][1][i]);
380  sec58row1->SetBinContent(i,secrow[58][1][i]);
381  sec59row1->SetBinContent(i,secrow[59][1][i]);
382  sec60row1->SetBinContent(i,secrow[60][1][i]);
383  sec55row2->SetBinContent(i,secrow[55][2][i]);
384  sec56row2->SetBinContent(i,secrow[56][2][i]);
385  sec57row2->SetBinContent(i,secrow[57][2][i]);
386  sec58row2->SetBinContent(i,secrow[58][2][i]);
387  sec59row2->SetBinContent(i,secrow[59][2][i]);
388  sec60row2->SetBinContent(i,secrow[60][2][i]);
389  } //End of pad loop
390 
391  TCanvas *c1 = new TCanvas("c1","ps",200,10,700,500);
392  c1->Divide(3,2);
393  psFile ="plotPadOcc.ps";
394  TPostScript *fps=new TPostScript(psFile,112);
395 
396  // PadRow 1
397  fps->NewPage();
398  c1->cd(1);
399  sec1row1->Draw();
400  c1->cd(2);
401  sec2row1->Draw();
402  c1->cd(3);
403  sec3row1->Draw();
404  c1->cd(4);
405  sec4row1->Draw();
406  c1->cd(5);
407  sec5row1->Draw();
408  c1->cd(6);
409  sec6row1->Draw();
410  c1->Update();
411  // PadRow 2
412  fps->NewPage();
413  c1->cd(1);
414  sec1row2->Draw();
415  c1->cd(2);
416  sec2row2->Draw();
417  c1->cd(3);
418  sec3row2->Draw();
419  c1->cd(4);
420  sec4row2->Draw();
421  c1->cd(5);
422  sec5row2->Draw();
423  c1->cd(6);
424  sec6row2->Draw();
425  c1->Update();
426  // PadRow 3
427  fps->NewPage();
428  c1->cd(1);
429  sec7row1->Draw();
430  c1->cd(2);
431  sec8row1->Draw();
432  c1->cd(3);
433  sec9row1->Draw();
434  c1->cd(4);
435  sec10row1->Draw();
436  c1->cd(5);
437  sec11row1->Draw();
438  c1->cd(6);
439  sec12row1->Draw();
440  c1->Update();
441  // PadRow 4
442  fps->NewPage();
443  c1->cd(1);
444  sec7row2->Draw();
445  c1->cd(2);
446  sec8row2->Draw();
447  c1->cd(3);
448  sec9row2->Draw();
449  c1->cd(4);
450  sec10row2->Draw();
451  c1->cd(5);
452  sec11row2->Draw();
453  c1->cd(6);
454  sec12row2->Draw();
455  c1->Update();
456  // PadRow 5
457  fps->NewPage();
458  c1->cd(1);
459  sec13row1->Draw();
460  c1->cd(2);
461  sec14row1->Draw();
462  c1->cd(3);
463  sec15row1->Draw();
464  c1->cd(4);
465  sec16row1->Draw();
466  c1->cd(5);
467  sec17row1->Draw();
468  c1->cd(6);
469  sec18row1->Draw();
470  c1->Update();
471  // PadRow 6
472  fps->NewPage();
473  c1->cd(1);
474  sec13row2->Draw();
475  c1->cd(2);
476  sec14row2->Draw();
477  c1->cd(3);
478  sec15row2->Draw();
479  c1->cd(4);
480  sec16row2->Draw();
481  c1->cd(5);
482  sec17row2->Draw();
483  c1->cd(6);
484  sec18row2->Draw();
485  c1->Update();
486  // PadRow 7
487  fps->NewPage();
488  c1->cd(1);
489  sec19row1->Draw();
490  c1->cd(2);
491  sec20row1->Draw();
492  c1->cd(3);
493  sec21row1->Draw();
494  c1->cd(4);
495  sec22row1->Draw();
496  c1->cd(5);
497  sec23row1->Draw();
498  c1->cd(6);
499  sec24row1->Draw();
500  c1->Update();
501  // PadRow 8
502  fps->NewPage();
503  c1->cd(1);
504  sec19row2->Draw();
505  c1->cd(2);
506  sec20row2->Draw();
507  c1->cd(3);
508  sec21row2->Draw();
509  c1->cd(4);
510  sec22row2->Draw();
511  c1->cd(5);
512  sec23row2->Draw();
513  c1->cd(6);
514  sec24row2->Draw();
515  c1->Update();
516  // PadRow 9
517  fps->NewPage();
518  c1->cd(1);
519  sec25row1->Draw();
520  c1->cd(2);
521  sec26row1->Draw();
522  c1->cd(3);
523  sec27row1->Draw();
524  c1->cd(4);
525  sec28row1->Draw();
526  c1->cd(5);
527  sec29row1->Draw();
528  c1->cd(6);
529  sec30row1->Draw();
530  c1->Update();
531  // PadRow 10
532  fps->NewPage();
533  c1->cd(1);
534  sec25row2->Draw();
535  c1->cd(2);
536  sec26row2->Draw();
537  c1->cd(3);
538  sec27row2->Draw();
539  c1->cd(4);
540  sec28row2->Draw();
541  c1->cd(5);
542  sec29row2->Draw();
543  c1->cd(6);
544  sec30row2->Draw();
545  c1->Update();
546 
547  // PadRow 11
548  fps->NewPage();
549  c1->cd(1);
550  sec31row1->Draw();
551  c1->cd(2);
552  sec32row1->Draw();
553  c1->cd(3);
554  sec33row1->Draw();
555  c1->cd(4);
556  sec34row1->Draw();
557  c1->cd(5);
558  sec35row1->Draw();
559  c1->cd(6);
560  sec36row1->Draw();
561  c1->Update();
562  // PadRow 12
563  fps->NewPage();
564  c1->cd(1);
565  sec31row2->Draw();
566  c1->cd(2);
567  sec32row2->Draw();
568  c1->cd(3);
569  sec33row2->Draw();
570  c1->cd(4);
571  sec34row2->Draw();
572  c1->cd(5);
573  sec35row2->Draw();
574  c1->cd(6);
575  sec36row2->Draw();
576  c1->Update();
577  // PadRow 13
578  fps->NewPage();
579  c1->cd(1);
580  sec37row1->Draw();
581  c1->cd(2);
582  sec38row1->Draw();
583  c1->cd(3);
584  sec39row1->Draw();
585  c1->cd(4);
586  sec40row1->Draw();
587  c1->cd(5);
588  sec41row1->Draw();
589  c1->cd(6);
590  sec42row1->Draw();
591  c1->Update();
592  // PadRow 14
593  fps->NewPage();
594  c1->cd(1);
595  sec37row2->Draw();
596  c1->cd(2);
597  sec38row2->Draw();
598  c1->cd(3);
599  sec39row2->Draw();
600  c1->cd(4);
601  sec40row2->Draw();
602  c1->cd(5);
603  sec41row2->Draw();
604  c1->cd(6);
605  sec42row2->Draw();
606  c1->Update();
607  // PadRow 15
608  fps->NewPage();
609  c1->cd(1);
610  sec43row1->Draw();
611  c1->cd(2);
612  sec44row1->Draw();
613  c1->cd(3);
614  sec45row1->Draw();
615  c1->cd(4);
616  sec46row1->Draw();
617  c1->cd(5);
618  sec47row1->Draw();
619  c1->cd(6);
620  sec48row1->Draw();
621  c1->Update();
622  // PadRow 16
623  fps->NewPage();
624  c1->cd(1);
625  sec43row2->Draw();
626  c1->cd(2);
627  sec44row2->Draw();
628  c1->cd(3);
629  sec45row2->Draw();
630  c1->cd(4);
631  sec46row2->Draw();
632  c1->cd(5);
633  sec47row2->Draw();
634  c1->cd(6);
635  sec48row2->Draw();
636  c1->Update();
637  // PadRow 17
638  fps->NewPage();
639  c1->cd(1);
640  sec49row1->Draw();
641  c1->cd(2);
642  sec50row1->Draw();
643  c1->cd(3);
644  sec51row1->Draw();
645  c1->cd(4);
646  sec52row1->Draw();
647  c1->cd(5);
648  sec53row1->Draw();
649  c1->cd(6);
650  sec54row1->Draw();
651  c1->Update();
652  // PadRow 18
653  fps->NewPage();
654  c1->cd(1);
655  sec49row2->Draw();
656  c1->cd(2);
657  sec50row2->Draw();
658  c1->cd(3);
659  sec51row2->Draw();
660  c1->cd(4);
661  sec52row2->Draw();
662  c1->cd(5);
663  sec53row2->Draw();
664  c1->cd(6);
665  sec54row2->Draw();
666  c1->Update();
667  // PadRow 19
668  fps->NewPage();
669  c1->cd(1);
670  sec55row1->Draw();
671  c1->cd(2);
672  sec56row1->Draw();
673  c1->cd(3);
674  sec57row1->Draw();
675  c1->cd(4);
676  sec58row1->Draw();
677  c1->cd(5);
678  sec59row1->Draw();
679  c1->cd(6);
680  sec60row1->Draw();
681  c1->Update();
682  // PadRow 20
683  fps->NewPage();
684  c1->cd(1);
685  sec55row2->Draw();
686  c1->cd(2);
687  sec56row2->Draw();
688  c1->cd(3);
689  sec57row2->Draw();
690  c1->cd(4);
691  sec58row2->Draw();
692  c1->cd(5);
693  sec59row2->Draw();
694  c1->cd(6);
695  sec60row2->Draw();
696  c1->Update();
697 
698  fps->Close();
699 }
700 
Generated by   doxygen 1.8.5