CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripDigiValid.cc
Go to the documentation of this file.
15 
16 using namespace std;
17 using namespace edm;
18 
20 
21  outputFile_ = ps.getUntrackedParameter<string>("outputFile", "stripdigihisto.root");
22  src_ = ps.getParameter<edm::InputTag>( "src" );
23 
25 
26 if ( dbe_ ) {
27  dbe_->setCurrentFolder("TrackerDigisV/TrackerDigis/Strip");
28 
29  for(int i = 0 ;i<3 ; i++) {
30  Char_t histo[200];
31  // Z Plus Side
32  sprintf(histo,"adc_tib_layer1_extmodule%d_zp",i+1);
33  meAdcTIBLayer1Extzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
34  sprintf(histo,"adc_tib_layer1_intmodule%d_zp",i+1);
35  meAdcTIBLayer1Intzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
36  sprintf(histo,"adc_tib_layer2_extmodule%d_zp",i+1);
37  meAdcTIBLayer2Extzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
38  sprintf(histo,"adc_tib_layer2_intmodule%d_zp",i+1);
39  meAdcTIBLayer2Intzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
40  sprintf(histo,"adc_tib_layer3_extmodule%d_zp",i+1);
41  meAdcTIBLayer3Extzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
42  sprintf(histo,"adc_tib_layer3_intmodule%d_zp",i+1);
43  meAdcTIBLayer3Intzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
44  sprintf(histo,"adc_tib_layer4_extmodule%d_zp",i+1);
45  meAdcTIBLayer4Extzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
46  sprintf(histo,"adc_tib_layer4_intmodule%d_zp",i+1);
47  meAdcTIBLayer4Intzp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
48  sprintf(histo,"strip_tib_layer1_extmodule%d_zp",i+1);
49  meStripTIBLayer1Extzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
50  sprintf(histo,"strip_tib_layer1_intmodule%d_zp",i+1);
51  meStripTIBLayer1Intzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
52  sprintf(histo,"strip_tib_layer2_extmodule%d_zp",i+1);
53  meStripTIBLayer2Extzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
54  sprintf(histo,"strip_tib_layer2_intmodule%d_zp",i+1);
55  meStripTIBLayer2Intzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
56  sprintf(histo,"strip_tib_layer3_extmodule%d_zp",i+1);
57  meStripTIBLayer3Extzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
58  sprintf(histo,"strip_tib_layer3_intmodule%d_zp",i+1);
59  meStripTIBLayer3Intzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
60  sprintf(histo,"strip_tib_layer4_extmodule%d_zp",i+1);
61  meStripTIBLayer4Extzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
62  sprintf(histo,"strip_tib_layer4_intmodule%d_zp",i+1);
63  meStripTIBLayer4Intzp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
64  // Z Minus Side
65  sprintf(histo,"adc_tib_layer1_extmodule%d_zm",i+1);
66  meAdcTIBLayer1Extzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
67  sprintf(histo,"adc_tib_layer1_intmodule%d_zm",i+1);
68  meAdcTIBLayer1Intzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
69  sprintf(histo,"adc_tib_layer2_extmodule%d_zm",i+1);
70  meAdcTIBLayer2Extzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
71  sprintf(histo,"adc_tib_layer2_intmodule%d_zm",i+1);
72  meAdcTIBLayer2Intzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
73  sprintf(histo,"adc_tib_layer3_extmodule%d_zm",i+1);
74  meAdcTIBLayer3Extzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
75  sprintf(histo,"adc_tib_layer3_intmodule%d_zm",i+1);
76  meAdcTIBLayer3Intzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
77  sprintf(histo,"adc_tib_layer4_extmodule%d_zm",i+1);
78  meAdcTIBLayer4Extzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
79  sprintf(histo,"adc_tib_layer4_intmodule%d_zm",i+1);
80  meAdcTIBLayer4Intzm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
81  sprintf(histo,"strip_tib_layer1_extmodule%d_zm",i+1);
82  meStripTIBLayer1Extzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
83  sprintf(histo,"strip_tib_layer1_intmodule%d_zm",i+1);
84  meStripTIBLayer1Intzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
85  sprintf(histo,"strip_tib_layer2_extmodule%d_zm",i+1);
86  meStripTIBLayer2Extzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
87  sprintf(histo,"strip_tib_layer2_intmodule%d_zm",i+1);
88  meStripTIBLayer2Intzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
89  sprintf(histo,"strip_tib_layer3_extmodule%d_zm",i+1);
90  meStripTIBLayer3Extzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
91  sprintf(histo,"strip_tib_layer3_intmodule%d_zm",i+1);
92  meStripTIBLayer3Intzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
93  sprintf(histo,"strip_tib_layer4_extmodule%d_zm",i+1);
94  meStripTIBLayer4Extzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
95  sprintf(histo,"strip_tib_layer4_intmodule%d_zm",i+1);
96  meStripTIBLayer4Intzm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
97 
98  }
99 
100  for(int i = 0 ;i<6 ; i++) {
101  Char_t histo[200];
102  // Z Plus Side
103  sprintf(histo,"adc_tob_layer1_module%d_zp",i+1);
104  meAdcTOBLayer1zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
105  sprintf(histo,"strip_tob_layer1_module%d_zp",i+1);
106  meStripTOBLayer1zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
107  sprintf(histo,"adc_tob_layer2_module%d_zp",i+1);
108  meAdcTOBLayer2zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
109  sprintf(histo,"strip_tob_layer2_module%d_zp",i+1);
110  meStripTOBLayer2zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
111  sprintf(histo,"adc_tob_layer3_module%d_zp",i+1);
112  meAdcTOBLayer3zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
113  sprintf(histo,"strip_tob_layer3_module%d_zp",i+1);
114  meStripTOBLayer3zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
115  sprintf(histo,"adc_tob_layer4_module%d_zp",i+1);
116  meAdcTOBLayer4zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
117  sprintf(histo,"strip_tob_layer4_module%d_zp",i+1);
118  meStripTOBLayer4zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
119  sprintf(histo,"adc_tob_layer5_module%d_zp",i+1);
120  meAdcTOBLayer5zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
121  sprintf(histo,"strip_tob_layer5_module%d_zp",i+1);
122  meStripTOBLayer5zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
123  sprintf(histo,"adc_tob_layer6_module%d_zp",i+1);
124  meAdcTOBLayer6zp_[i] = dbe_->book1D(histo,"Digis ADC",10,0.,300.);
125  sprintf(histo,"strip_tob_layer6_module%d_zp",i+1);
126  meStripTOBLayer6zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
127  // Z Minus Side
128  sprintf(histo,"adc_tob_layer1_module%d_zm",i+1);
129  meAdcTOBLayer1zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
130  sprintf(histo,"strip_tob_layer1_module%d_zm",i+1);
131  meStripTOBLayer1zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
132  sprintf(histo,"adc_tob_layer2_module%d_zm",i+1);
133  meAdcTOBLayer2zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
134  sprintf(histo,"strip_tob_layer2_module%d_zm",i+1);
135  meStripTOBLayer2zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
136  sprintf(histo,"adc_tob_layer3_module%d_zm",i+1);
137  meAdcTOBLayer3zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
138  sprintf(histo,"strip_tob_layer3_module%d_zm",i+1);
139  meStripTOBLayer3zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
140  sprintf(histo,"adc_tob_layer4_module%d_zm",i+1);
141  meAdcTOBLayer4zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
142  sprintf(histo,"strip_tob_layer4_module%d_zm",i+1);
143  meStripTOBLayer4zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
144  sprintf(histo,"adc_tob_layer5_module%d_zm",i+1);
145  meAdcTOBLayer5zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
146  sprintf(histo,"strip_tob_layer5_module%d_zm",i+1);
147  meStripTOBLayer5zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
148  sprintf(histo,"adc_tob_layer6_module%d_zm",i+1);
149  meAdcTOBLayer6zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
150  sprintf(histo,"strip_tob_layer6_module%d_zm",i+1);
151  meStripTOBLayer6zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
152  }
153 
154  for(int i = 0 ;i<3 ; i++) {
155  Char_t histo[200];
156  // Z Plus Side
157  sprintf(histo,"adc_tid_wheel1_ring%d_zp",i+1);
158  meAdcTIDWheel1zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
159  sprintf(histo,"adc_tid_wheel2_ring%d_zp",i+1);
160  meAdcTIDWheel2zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
161  sprintf(histo,"adc_tid_wheel3_ring%d_zp",i+1);
162  meAdcTIDWheel3zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
163  sprintf(histo,"strip_tid_wheel1_ring%d_zp",i+1);
164  meStripTIDWheel1zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
165  sprintf(histo,"strip_tid_wheel2_ring%d_zp",i+1);
166  meStripTIDWheel2zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
167  sprintf(histo,"strip_tid_wheel3_ring%d_zp",i+1);
168  meStripTIDWheel3zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
169  // Z minus Side
170  sprintf(histo,"adc_tid_wheel1_ring%d_zm",i+1);
171  meAdcTIDWheel1zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
172  sprintf(histo,"adc_tid_wheel2_ring%d_zm",i+1);
173  meAdcTIDWheel2zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
174  sprintf(histo,"adc_tid_wheel3_ring%d_zm",i+1);
175  meAdcTIDWheel3zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
176  sprintf(histo,"strip_tid_wheel1_ring%d_zm",i+1);
177  meStripTIDWheel1zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
178  sprintf(histo,"strip_tid_wheel2_ring%d_zm",i+1);
179  meStripTIDWheel2zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
180  sprintf(histo,"strip_tid_wheel3_ring%d_zm",i+1);
181  meStripTIDWheel3zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
182 
183  }
184 
185  for(int i = 0 ;i<7 ; i++) {
186  Char_t histo[200];
187  // Z Plus Side
188  sprintf(histo,"adc_tec_wheel1_ring%d_zp",i+1);
189  meAdcTECWheel1zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
190  sprintf(histo,"adc_tec_wheel2_ring%d_zp",i+1);
191  meAdcTECWheel2zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
192  sprintf(histo,"adc_tec_wheel3_ring%d_zp",i+1);
193  meAdcTECWheel3zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
194  sprintf(histo,"strip_tec_wheel1_ring%d_zp",i+1);
195  meStripTECWheel1zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
196  sprintf(histo,"strip_tec_wheel2_ring%d_zp",i+1);
197  meStripTECWheel2zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
198  sprintf(histo,"strip_tec_wheel3_ring%d_zp",i+1);
199  meStripTECWheel3zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
200 
201  // Z Minus Side
202  sprintf(histo,"adc_tec_wheel1_ring%d_zm",i+1);
203  meAdcTECWheel1zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
204  sprintf(histo,"adc_tec_wheel2_ring%d_zm",i+1);
205  meAdcTECWheel2zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
206  sprintf(histo,"adc_tec_wheel3_ring%d_zm",i+1);
207  meAdcTECWheel3zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
208  sprintf(histo,"strip_tec_wheel1_ring%d_zm",i+1);
209  meStripTECWheel1zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
210  sprintf(histo,"strip_tec_wheel2_ring%d_zm",i+1);
211  meStripTECWheel2zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
212  sprintf(histo,"strip_tec_wheel3_ring%d_zm",i+1);
213  meStripTECWheel3zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
214  }
215 
216  for(int i = 0 ;i<6 ; i++) {
217  Char_t histo[200];
218  // Z Plus Side
219  sprintf(histo,"adc_tec_wheel4_ring%d_zp",i+1);
220  meAdcTECWheel4zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
221  sprintf(histo,"adc_tec_wheel5_ring%d_zp",i+1);
222  meAdcTECWheel5zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
223  sprintf(histo,"adc_tec_wheel6_ring%d_zp",i+1);
224  meAdcTECWheel6zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
225  sprintf(histo,"strip_tec_wheel4_ring%d_zp",i+1);
226  meStripTECWheel4zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
227  sprintf(histo,"strip_tec_wheel5_ring%d_zp",i+1);
228  meStripTECWheel5zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
229  sprintf(histo,"strip_tec_wheel6_ring%d_zp",i+1);
230  meStripTECWheel6zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
231 
232  // Z Minus Side
233  sprintf(histo,"adc_tec_wheel4_ring%d_zm",i+1);
234  meAdcTECWheel4zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
235  sprintf(histo,"adc_tec_wheel5_ring%d_zm",i+1);
236  meAdcTECWheel5zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
237  sprintf(histo,"adc_tec_wheel6_ring%d_zm",i+1);
238  meAdcTECWheel6zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
239  sprintf(histo,"strip_tec_wheel4_ring%d_zm",i+1);
240  meStripTECWheel4zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
241  sprintf(histo,"strip_tec_wheel5_ring%d_zm",i+1);
242  meStripTECWheel5zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
243  sprintf(histo,"strip_tec_wheel6_ring%d_zm",i+1);
244  meStripTECWheel6zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
245  }
246 
247  for(int i = 0 ;i<5 ; i++) {
248  Char_t histo[200];
249  // Z Plus Side
250  sprintf(histo,"adc_tec_wheel7_ring%d_zp",i+1);
251  meAdcTECWheel7zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
252  sprintf(histo,"adc_tec_wheel8_ring%d_zp",i+1);
253  meAdcTECWheel8zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
254  sprintf(histo,"strip_tec_wheel7_ring%d_zp",i+1);
255  meStripTECWheel7zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
256  sprintf(histo,"strip_tec_wheel8_ring%d_zp",i+1);
257  meStripTECWheel8zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
258 
259  // Z Minus Side
260  sprintf(histo,"adc_tec_wheel7_ring%d_zm",i+1);
261  meAdcTECWheel7zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
262  sprintf(histo,"adc_tec_wheel8_ring%d_zm",i+1);
263  meAdcTECWheel8zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
264  sprintf(histo,"strip_tec_wheel7_ring%d_zm",i+1);
265  meStripTECWheel7zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
266  sprintf(histo,"strip_tec_wheel8_ring%d_zm",i+1);
267  meStripTECWheel8zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
268  }
269 
270  for(int i = 0 ;i<4 ; i++) {
271  Char_t histo[200];
272  // Z Plus Side
273  sprintf(histo,"adc_tec_wheel9_ring%d_zp",i+1);
274  meAdcTECWheel9zp_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
275  sprintf(histo,"strip_tec_wheel9_ring%d_zp",i+1);
276  meStripTECWheel9zp_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
277 
278  // Z Minus Side
279  sprintf(histo,"adc_tec_wheel9_ring%d_zm",i+1);
280  meAdcTECWheel9zm_[i] = dbe_->book1D(histo,"Digis ADC",50,0.,300.);
281  sprintf(histo,"strip_tec_wheel9_ring%d_zm",i+1);
282  meStripTECWheel9zm_[i] = dbe_->book1D(histo,"Digis Strip Num.",200,0.,800.);
283  }
284 
285 
286  for(int i = 0 ;i<4 ; i++) {
287  Char_t histo[200];
288  sprintf(histo,"ndigi_tib_layer_%d_zm",i+1);
289  meNDigiTIBLayerzm_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
290  sprintf(histo,"ndigi_tib_layer_%d_zp",i+1);
291  meNDigiTIBLayerzp_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
292  }
293 
294  for(int i = 0 ;i<6 ; i++) {
295  Char_t histo[200];
296  sprintf(histo,"ndigi_tob_layer_%d_zm",i+1);
297  meNDigiTOBLayerzm_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
298  sprintf(histo,"ndigi_tob_layer_%d_zp",i+1);
299  meNDigiTOBLayerzp_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
300 
301  }
302 
303  for(int i = 0 ;i<3 ; i++) {
304  Char_t histo[200];
305  sprintf(histo,"ndigi_tid_wheel_%d_zm",i+1);
306  meNDigiTIDWheelzm_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
307  sprintf(histo,"ndigi_tid_wheel_%d_zp",i+1);
308  meNDigiTIDWheelzp_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
309  }
310 
311  for(int i = 0 ;i<9 ; i++) {
312  Char_t histo[200];
313  sprintf(histo,"ndigi_tec_wheel_%d_zm",i+1);
314  meNDigiTECWheelzm_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
315  sprintf(histo,"ndigi_tec_wheel_%d_zp",i+1);
316  meNDigiTECWheelzp_[i] = dbe_->book1D(histo, "Digi Multiplicity",100,0.,500.);
317 
318  }
319  }
320 
321 }
322 
324 
325  //if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
326 }
327 
328 
330 
331 }
332 
334  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
335 
336 }
337 
338 
340  //Retrieve tracker topology from geometry
341  edm::ESHandle<TrackerTopology> tTopoHandle;
342  c.get<IdealGeometryRecord>().get(tTopoHandle);
343  const TrackerTopology* const tTopo = tTopoHandle.product();
344 
345 
346 
347  int ndigilayertibzp[4];
348  int ndigilayertibzm[4];
349 
350  for( int i = 0; i< 4; i++ ) {
351  ndigilayertibzp[i] = 0;
352  ndigilayertibzm[i] = 0;
353  }
354 
355  int ndigilayertobzp[6];
356  int ndigilayertobzm[6];
357 
358  for( int i = 0; i< 6; i++ ) {
359  ndigilayertobzp[i] = 0;
360  ndigilayertobzm[i] = 0;
361  }
362 
363  int ndigiwheeltidzp[3];
364  int ndigiwheeltidzm[3];
365 
366  for( int i = 0; i< 3; i++ ) {
367  ndigiwheeltidzp[i] = 0;
368  ndigiwheeltidzm[i] = 0;
369  }
370 
371  int ndigiwheelteczp[9];
372  int ndigiwheelteczm[9];
373 
374  for( int i = 0; i< 9; i++ ) {
375  ndigiwheelteczp[i] = 0;
376  ndigiwheelteczm[i] = 0;
377  }
378 
379 
380 
381  //LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
383  c.get<TrackerDigiGeometryRecord>().get( tracker );
384 
385  std::string digiProducer = "siStripDigis";
387  e.getByLabel(src_, stripDigis);
388  edm::DetSetVector<SiStripDigi>::const_iterator DSViter = stripDigis->begin();
389  for( ; DSViter != stripDigis->end(); DSViter++) {
390  unsigned int id = DSViter->id;
391  DetId detId(id);
395 
396  if(detId.subdetId()==StripSubdetector::TIB){
397 
398  for ( iter = begin ; iter != end; iter++ ) { // loop digis
399  if( tTopo->tibStringInfo(id)[0] == 1) {
400  ++ndigilayertibzm[tTopo->tibLayer(id)-1];
401  if( tTopo->tibLayer(id) == 1 ) {
402  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
403  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer1Intzm_[0] -> Fill((*iter).adc()); meStripTIBLayer1Intzm_[0] ->Fill((*iter).strip()); }
404  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer1Intzm_[1] -> Fill((*iter).adc()); meStripTIBLayer1Intzm_[1] ->Fill((*iter).strip()); }
405  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer1Intzm_[2] -> Fill((*iter).adc()); meStripTIBLayer1Intzm_[2] ->Fill((*iter).strip()); }
406  }else {
407  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer1Extzm_[0] -> Fill((*iter).adc()); meStripTIBLayer1Extzm_[0] ->Fill((*iter).strip()); }
408  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer1Extzm_[1] -> Fill((*iter).adc()); meStripTIBLayer1Extzm_[1] ->Fill((*iter).strip()); }
409  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer1Extzm_[2] -> Fill((*iter).adc()); meStripTIBLayer1Extzm_[2] ->Fill((*iter).strip()); }
410  }
411  }
412  if( tTopo->tibLayer(id) == 2 ) {
413  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
414  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer2Intzm_[0] -> Fill((*iter).adc()); meStripTIBLayer2Intzm_[0] ->Fill((*iter).strip()); }
415  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer2Intzm_[1] -> Fill((*iter).adc()); meStripTIBLayer2Intzm_[1] ->Fill((*iter).strip()); }
416  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer2Intzm_[2] -> Fill((*iter).adc()); meStripTIBLayer2Intzm_[2] ->Fill((*iter).strip()); }
417  }else {
418  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer2Extzm_[0] -> Fill((*iter).adc()); meStripTIBLayer2Extzm_[0] ->Fill((*iter).strip()); }
419  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer2Extzm_[1] -> Fill((*iter).adc()); meStripTIBLayer2Extzm_[1] ->Fill((*iter).strip()); }
420  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer2Extzm_[2] -> Fill((*iter).adc()); meStripTIBLayer2Extzm_[2] ->Fill((*iter).strip()); }
421  }
422  }
423  if( tTopo->tibLayer(id) == 3 ) {
424  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
425  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer3Intzm_[0] -> Fill((*iter).adc()); meStripTIBLayer3Intzm_[0] ->Fill((*iter).strip()); }
426  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer3Intzm_[1] -> Fill((*iter).adc()); meStripTIBLayer3Intzm_[1] ->Fill((*iter).strip()); }
427  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer3Intzm_[2] -> Fill((*iter).adc()); meStripTIBLayer3Intzm_[2] ->Fill((*iter).strip()); }
428  }else {
429  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer3Extzm_[0] -> Fill((*iter).adc()); meStripTIBLayer3Extzm_[0] ->Fill((*iter).strip()); }
430  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer3Extzm_[1] -> Fill((*iter).adc()); meStripTIBLayer3Extzm_[1] ->Fill((*iter).strip()); }
431  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer3Extzm_[2] -> Fill((*iter).adc()); meStripTIBLayer3Extzm_[2] ->Fill((*iter).strip()); }
432  }
433  }
434  if( tTopo->tibLayer(id) == 4 ) {
435  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
436  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer4Intzm_[0] -> Fill((*iter).adc()); meStripTIBLayer4Intzm_[0] ->Fill((*iter).strip()); }
437  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer4Intzm_[1] -> Fill((*iter).adc()); meStripTIBLayer4Intzm_[1] ->Fill((*iter).strip()); }
438  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer4Intzm_[2] -> Fill((*iter).adc()); meStripTIBLayer4Intzm_[2] ->Fill((*iter).strip()); }
439  }else {
440  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer4Extzm_[0] -> Fill((*iter).adc()); meStripTIBLayer4Extzm_[0] ->Fill((*iter).strip()); }
441  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer4Extzm_[1] -> Fill((*iter).adc()); meStripTIBLayer4Extzm_[1] ->Fill((*iter).strip()); }
442  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer4Extzm_[2] -> Fill((*iter).adc()); meStripTIBLayer4Extzm_[2] ->Fill((*iter).strip()); }
443  }
444  }
445  }else {
446  ++ndigilayertibzp[tTopo->tibLayer(id)-1];
447  if( tTopo->tibLayer(id) == 1 ) {
448  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
449  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer1Intzp_[0] -> Fill((*iter).adc()); meStripTIBLayer1Intzp_[0] ->Fill((*iter).strip()); }
450  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer1Intzp_[1] -> Fill((*iter).adc()); meStripTIBLayer1Intzp_[1] ->Fill((*iter).strip()); }
451  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer1Intzp_[2] -> Fill((*iter).adc()); meStripTIBLayer1Intzp_[2] ->Fill((*iter).strip()); }
452  }else {
453  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer1Extzp_[0] -> Fill((*iter).adc()); meStripTIBLayer1Extzp_[0] ->Fill((*iter).strip()); }
454  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer1Extzp_[1] -> Fill((*iter).adc()); meStripTIBLayer1Extzp_[1] ->Fill((*iter).strip()); }
455  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer1Extzp_[2] -> Fill((*iter).adc()); meStripTIBLayer1Extzp_[2] ->Fill((*iter).strip()); }
456  }
457  }
458  if( tTopo->tibLayer(id) == 2 ) {
459  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
460  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer2Intzp_[0] -> Fill((*iter).adc()); meStripTIBLayer2Intzp_[0] ->Fill((*iter).strip()); }
461  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer2Intzp_[1] -> Fill((*iter).adc()); meStripTIBLayer2Intzp_[1] ->Fill((*iter).strip()); }
462  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer2Intzp_[2] -> Fill((*iter).adc()); meStripTIBLayer2Intzp_[2] ->Fill((*iter).strip()); }
463  }else {
464  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer2Extzp_[0] -> Fill((*iter).adc()); meStripTIBLayer2Extzp_[0] ->Fill((*iter).strip()); }
465  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer2Extzp_[1] -> Fill((*iter).adc()); meStripTIBLayer2Extzp_[1] ->Fill((*iter).strip()); }
466  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer2Extzp_[2] -> Fill((*iter).adc()); meStripTIBLayer2Extzp_[2] ->Fill((*iter).strip()); }
467  }
468  }
469  if( tTopo->tibLayer(id) == 3 ) {
470  if ( tTopo->tibStringInfo(id)[1] == 1 ) {
471  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer3Intzp_[0] -> Fill((*iter).adc()); meStripTIBLayer3Intzp_[0] ->Fill((*iter).strip()); }
472  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer3Intzp_[1] -> Fill((*iter).adc()); meStripTIBLayer3Intzp_[1] ->Fill((*iter).strip()); }
473  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer3Intzp_[2] -> Fill((*iter).adc()); meStripTIBLayer3Intzp_[2] ->Fill((*iter).strip()); }
474  }else {
475  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer3Extzp_[0] -> Fill((*iter).adc()); meStripTIBLayer3Extzp_[0] ->Fill((*iter).strip()); }
476  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer3Extzp_[1] -> Fill((*iter).adc()); meStripTIBLayer3Extzp_[1] ->Fill((*iter).strip()); }
477  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer3Extzp_[2] -> Fill((*iter).adc()); meStripTIBLayer3Extzp_[2] ->Fill((*iter).strip()); }
478  }
479  }
480  if( tTopo->tibLayer(id) == 4 ) {
481  if ( tTopo->tibStringInfo(id)[1] == 1) {
482  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer4Intzp_[0] -> Fill((*iter).adc()); meStripTIBLayer4Intzp_[0] ->Fill((*iter).strip()); }
483  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer4Intzp_[1] -> Fill((*iter).adc()); meStripTIBLayer4Intzp_[1] ->Fill((*iter).strip()); }
484  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer4Intzp_[2] -> Fill((*iter).adc()); meStripTIBLayer4Intzp_[2] ->Fill((*iter).strip()); }
485  }else {
486  if( tTopo->tibModule(id) == 1 ) { meAdcTIBLayer4Extzp_[0] -> Fill((*iter).adc()); meStripTIBLayer4Extzp_[0] ->Fill((*iter).strip()); }
487  if( tTopo->tibModule(id) == 2 ) { meAdcTIBLayer4Extzp_[1] -> Fill((*iter).adc()); meStripTIBLayer4Extzp_[1] ->Fill((*iter).strip()); }
488  if( tTopo->tibModule(id) == 3 ) { meAdcTIBLayer4Extzp_[2] -> Fill((*iter).adc()); meStripTIBLayer4Extzp_[2] ->Fill((*iter).strip()); }
489  }
490  }
491 
492  }
493  }
494  }
495  if(detId.subdetId()==StripSubdetector::TOB){
496 
497  for ( iter = begin ; iter != end; iter++ ) { // loop digis
498  if( tTopo->tobRodInfo(id)[0] == 1) {
499  ++ndigilayertobzm[tTopo->tobLayer(id)-1];
500  if( tTopo->tobLayer(id) == 1 ) {
501  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer1zm_[0] -> Fill((*iter).adc()); meStripTOBLayer1zm_[0] ->Fill((*iter).strip()); }
502  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer1zm_[1] -> Fill((*iter).adc()); meStripTOBLayer1zm_[1] ->Fill((*iter).strip()); }
503  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer1zm_[2] -> Fill((*iter).adc()); meStripTOBLayer1zm_[2] ->Fill((*iter).strip()); }
504  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer1zm_[3] -> Fill((*iter).adc()); meStripTOBLayer1zm_[3] ->Fill((*iter).strip()); }
505  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer1zm_[4] -> Fill((*iter).adc()); meStripTOBLayer1zm_[4] ->Fill((*iter).strip()); }
506  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer1zm_[5] -> Fill((*iter).adc()); meStripTOBLayer1zm_[5] ->Fill((*iter).strip()); }
507  }
508  if( tTopo->tobLayer(id) == 2 ) {
509  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer2zm_[0] -> Fill((*iter).adc()); meStripTOBLayer2zm_[0] ->Fill((*iter).strip()); }
510  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer2zm_[1] -> Fill((*iter).adc()); meStripTOBLayer2zm_[1] ->Fill((*iter).strip()); }
511  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer2zm_[2] -> Fill((*iter).adc()); meStripTOBLayer2zm_[2] ->Fill((*iter).strip()); }
512  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer2zm_[3] -> Fill((*iter).adc()); meStripTOBLayer2zm_[3] ->Fill((*iter).strip()); }
513  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer2zm_[4] -> Fill((*iter).adc()); meStripTOBLayer2zm_[4] ->Fill((*iter).strip()); }
514  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer2zm_[5] -> Fill((*iter).adc()); meStripTOBLayer2zm_[5] ->Fill((*iter).strip()); }
515  }
516  if( tTopo->tobLayer(id) == 3 ) {
517  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer3zm_[0] -> Fill((*iter).adc()); meStripTOBLayer3zm_[0] ->Fill((*iter).strip()); }
518  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer3zm_[1] -> Fill((*iter).adc()); meStripTOBLayer3zm_[1] ->Fill((*iter).strip()); }
519  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer3zm_[2] -> Fill((*iter).adc()); meStripTOBLayer3zm_[2] ->Fill((*iter).strip()); }
520  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer3zm_[3] -> Fill((*iter).adc()); meStripTOBLayer3zm_[3] ->Fill((*iter).strip()); }
521  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer3zm_[4] -> Fill((*iter).adc()); meStripTOBLayer3zm_[4] ->Fill((*iter).strip()); }
522  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer3zm_[5] -> Fill((*iter).adc()); meStripTOBLayer3zm_[5] ->Fill((*iter).strip()); }
523  }
524  if( tTopo->tobLayer(id) == 4 ) {
525  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer4zm_[0] -> Fill((*iter).adc()); meStripTOBLayer4zm_[0] ->Fill((*iter).strip()); }
526  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer4zm_[1] -> Fill((*iter).adc()); meStripTOBLayer4zm_[1] ->Fill((*iter).strip()); }
527  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer4zm_[2] -> Fill((*iter).adc()); meStripTOBLayer4zm_[2] ->Fill((*iter).strip()); }
528  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer4zm_[3] -> Fill((*iter).adc()); meStripTOBLayer4zm_[3] ->Fill((*iter).strip()); }
529  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer4zm_[4] -> Fill((*iter).adc()); meStripTOBLayer4zm_[4] ->Fill((*iter).strip()); }
530  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer4zm_[5] -> Fill((*iter).adc()); meStripTOBLayer4zm_[5] ->Fill((*iter).strip()); }
531  }
532  if( tTopo->tobLayer(id) == 5 ) {
533  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer5zm_[0] -> Fill((*iter).adc()); meStripTOBLayer5zm_[0] ->Fill((*iter).strip()); }
534  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer5zm_[1] -> Fill((*iter).adc()); meStripTOBLayer5zm_[1] ->Fill((*iter).strip()); }
535  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer5zm_[2] -> Fill((*iter).adc()); meStripTOBLayer5zm_[2] ->Fill((*iter).strip()); }
536  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer5zm_[3] -> Fill((*iter).adc()); meStripTOBLayer5zm_[3] ->Fill((*iter).strip()); }
537  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer5zm_[4] -> Fill((*iter).adc()); meStripTOBLayer5zm_[4] ->Fill((*iter).strip()); }
538  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer5zm_[5] -> Fill((*iter).adc()); meStripTOBLayer5zm_[5] ->Fill((*iter).strip()); }
539  }
540  if( tTopo->tobLayer(id) == 6 ) {
541  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer6zm_[0] -> Fill((*iter).adc()); meStripTOBLayer6zm_[0] ->Fill((*iter).strip()); }
542  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer6zm_[1] -> Fill((*iter).adc()); meStripTOBLayer6zm_[1] ->Fill((*iter).strip()); }
543  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer6zm_[2] -> Fill((*iter).adc()); meStripTOBLayer6zm_[2] ->Fill((*iter).strip()); }
544  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer6zm_[3] -> Fill((*iter).adc()); meStripTOBLayer6zm_[3] ->Fill((*iter).strip()); }
545  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer6zm_[4] -> Fill((*iter).adc()); meStripTOBLayer6zm_[4] ->Fill((*iter).strip()); }
546  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer6zm_[5] -> Fill((*iter).adc()); meStripTOBLayer6zm_[5] ->Fill((*iter).strip()); }
547  }
548 
549  }else {
550  ++ndigilayertobzp[tTopo->tobLayer(id)-1];
551  if( tTopo->tobLayer(id) == 1 ) {
552  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer1zp_[0] -> Fill((*iter).adc()); meStripTOBLayer1zp_[0] ->Fill((*iter).strip()); }
553  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer1zp_[1] -> Fill((*iter).adc()); meStripTOBLayer1zp_[1] ->Fill((*iter).strip()); }
554  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer1zp_[2] -> Fill((*iter).adc()); meStripTOBLayer1zp_[2] ->Fill((*iter).strip()); }
555  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer1zp_[3] -> Fill((*iter).adc()); meStripTOBLayer1zp_[3] ->Fill((*iter).strip()); }
556  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer1zp_[4] -> Fill((*iter).adc()); meStripTOBLayer1zp_[4] ->Fill((*iter).strip()); }
557  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer1zp_[5] -> Fill((*iter).adc()); meStripTOBLayer1zp_[5] ->Fill((*iter).strip()); }
558  }
559  if( tTopo->tobLayer(id) == 2 ) {
560  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer2zp_[0] -> Fill((*iter).adc()); meStripTOBLayer2zp_[0] ->Fill((*iter).strip()); }
561  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer2zp_[1] -> Fill((*iter).adc()); meStripTOBLayer2zp_[1] ->Fill((*iter).strip()); }
562  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer2zp_[2] -> Fill((*iter).adc()); meStripTOBLayer2zp_[2] ->Fill((*iter).strip()); }
563  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer2zp_[3] -> Fill((*iter).adc()); meStripTOBLayer2zp_[3] ->Fill((*iter).strip()); }
564  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer2zp_[4] -> Fill((*iter).adc()); meStripTOBLayer2zp_[4] ->Fill((*iter).strip()); }
565  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer2zp_[5] -> Fill((*iter).adc()); meStripTOBLayer2zp_[5] ->Fill((*iter).strip()); }
566  }
567  if( tTopo->tobLayer(id) == 3 ) {
568  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer3zp_[0] -> Fill((*iter).adc()); meStripTOBLayer3zp_[0] ->Fill((*iter).strip()); }
569  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer3zp_[1] -> Fill((*iter).adc()); meStripTOBLayer3zp_[1] ->Fill((*iter).strip()); }
570  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer3zp_[2] -> Fill((*iter).adc()); meStripTOBLayer3zp_[2] ->Fill((*iter).strip()); }
571  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer3zp_[3] -> Fill((*iter).adc()); meStripTOBLayer3zp_[3] ->Fill((*iter).strip()); }
572  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer3zp_[4] -> Fill((*iter).adc()); meStripTOBLayer3zp_[4] ->Fill((*iter).strip()); }
573  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer3zp_[5] -> Fill((*iter).adc()); meStripTOBLayer3zp_[5] ->Fill((*iter).strip()); }
574  }
575  if( tTopo->tobLayer(id) == 4 ) {
576  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer4zp_[0] -> Fill((*iter).adc()); meStripTOBLayer4zp_[0] ->Fill((*iter).strip()); }
577  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer4zp_[1] -> Fill((*iter).adc()); meStripTOBLayer4zp_[1] ->Fill((*iter).strip()); }
578  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer4zp_[2] -> Fill((*iter).adc()); meStripTOBLayer4zp_[2] ->Fill((*iter).strip()); }
579  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer4zp_[3] -> Fill((*iter).adc()); meStripTOBLayer4zp_[3] ->Fill((*iter).strip()); }
580  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer4zp_[4] -> Fill((*iter).adc()); meStripTOBLayer4zp_[4] ->Fill((*iter).strip()); }
581  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer4zp_[5] -> Fill((*iter).adc()); meStripTOBLayer4zp_[5] ->Fill((*iter).strip()); }
582  }
583  if( tTopo->tobLayer(id) == 5 ) {
584  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer5zp_[0] -> Fill((*iter).adc()); meStripTOBLayer5zp_[0] ->Fill((*iter).strip()); }
585  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer5zp_[1] -> Fill((*iter).adc()); meStripTOBLayer5zp_[1] ->Fill((*iter).strip()); }
586  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer5zp_[2] -> Fill((*iter).adc()); meStripTOBLayer5zp_[2] ->Fill((*iter).strip()); }
587  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer5zp_[3] -> Fill((*iter).adc()); meStripTOBLayer5zp_[3] ->Fill((*iter).strip()); }
588  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer5zp_[4] -> Fill((*iter).adc()); meStripTOBLayer5zp_[4] ->Fill((*iter).strip()); }
589  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer5zp_[5] -> Fill((*iter).adc()); meStripTOBLayer5zp_[5] ->Fill((*iter).strip()); }
590  }
591  if( tTopo->tobLayer(id) == 6 ) {
592  if ( tTopo->tobModule(id) == 1 ) { meAdcTOBLayer6zp_[0] -> Fill((*iter).adc()); meStripTOBLayer6zp_[0] ->Fill((*iter).strip()); }
593  if ( tTopo->tobModule(id) == 2 ) { meAdcTOBLayer6zp_[1] -> Fill((*iter).adc()); meStripTOBLayer6zp_[1] ->Fill((*iter).strip()); }
594  if ( tTopo->tobModule(id) == 3 ) { meAdcTOBLayer6zp_[2] -> Fill((*iter).adc()); meStripTOBLayer6zp_[2] ->Fill((*iter).strip()); }
595  if ( tTopo->tobModule(id) == 4 ) { meAdcTOBLayer6zp_[3] -> Fill((*iter).adc()); meStripTOBLayer6zp_[3] ->Fill((*iter).strip()); }
596  if ( tTopo->tobModule(id) == 5 ) { meAdcTOBLayer6zp_[4] -> Fill((*iter).adc()); meStripTOBLayer6zp_[4] ->Fill((*iter).strip()); }
597  if ( tTopo->tobModule(id) == 6 ) { meAdcTOBLayer6zp_[5] -> Fill((*iter).adc()); meStripTOBLayer6zp_[5] ->Fill((*iter).strip()); }
598  }
599 
600 
601  }
602  }
603  }
604 
605  if (detId.subdetId()==StripSubdetector::TID) {
606 
607  for ( iter = begin ; iter != end; iter++ ) {
608  if( tTopo->tidSide(id) == 1){
609  ++ndigiwheeltidzm[tTopo->tidWheel(id)-1];
610  if( tTopo->tidWheel(id) == 1 ) {
611  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel1zm_[0] -> Fill((*iter).adc()); meStripTIDWheel1zm_[0] ->Fill((*iter).strip());}
612  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel1zm_[1] -> Fill((*iter).adc()); meStripTIDWheel1zm_[1] ->Fill((*iter).strip());}
613  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel1zm_[2] -> Fill((*iter).adc()); meStripTIDWheel1zm_[2] ->Fill((*iter).strip());}
614  }
615  if( tTopo->tidWheel(id) == 2 ) {
616  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel2zm_[0] -> Fill((*iter).adc()); meStripTIDWheel2zm_[0] ->Fill((*iter).strip());}
617  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel2zm_[1] -> Fill((*iter).adc()); meStripTIDWheel2zm_[1] ->Fill((*iter).strip());}
618  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel2zm_[2] -> Fill((*iter).adc()); meStripTIDWheel2zm_[2] ->Fill((*iter).strip());}
619  }
620  if( tTopo->tidWheel(id) == 3 ) {
621  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel3zm_[0] -> Fill((*iter).adc()); meStripTIDWheel3zm_[0] ->Fill((*iter).strip());}
622  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel3zm_[1] -> Fill((*iter).adc()); meStripTIDWheel3zm_[1] ->Fill((*iter).strip());}
623  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel3zm_[2] -> Fill((*iter).adc()); meStripTIDWheel3zm_[2] ->Fill((*iter).strip());}
624  }
625 
626  }else{
627  ++ndigiwheeltidzp[tTopo->tidWheel(id)-1];
628  if( tTopo->tidWheel(id) == 1 ) {
629  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel1zp_[0] -> Fill((*iter).adc()); meStripTIDWheel1zp_[0] ->Fill((*iter).strip());}
630  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel1zp_[1] -> Fill((*iter).adc()); meStripTIDWheel1zp_[1] ->Fill((*iter).strip());}
631  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel1zp_[2] -> Fill((*iter).adc()); meStripTIDWheel1zp_[2] ->Fill((*iter).strip());}
632  }
633  if( tTopo->tidWheel(id) == 2 ) {
634  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel2zp_[0] -> Fill((*iter).adc()); meStripTIDWheel2zp_[0] ->Fill((*iter).strip());}
635  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel2zp_[1] -> Fill((*iter).adc()); meStripTIDWheel2zp_[1] ->Fill((*iter).strip());}
636  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel2zp_[2] -> Fill((*iter).adc()); meStripTIDWheel2zp_[2] ->Fill((*iter).strip());}
637  }
638  if( tTopo->tidWheel(id) == 3 ) {
639  if(tTopo->tidRing(id)== 1) { meAdcTIDWheel3zp_[0] -> Fill((*iter).adc()); meStripTIDWheel3zp_[0] ->Fill((*iter).strip());}
640  if(tTopo->tidRing(id)== 2) { meAdcTIDWheel3zp_[1] -> Fill((*iter).adc()); meStripTIDWheel3zp_[1] ->Fill((*iter).strip());}
641  if(tTopo->tidRing(id)== 3) { meAdcTIDWheel3zp_[2] -> Fill((*iter).adc()); meStripTIDWheel3zp_[2] ->Fill((*iter).strip());}
642  }
643 
644  }
645  }
646  }
647  if (detId.subdetId()==StripSubdetector::TEC) {
648 
649  for ( iter = begin ; iter != end; iter++ ) {
650  if(tTopo->tecSide(id) == 1) {
651  ++ndigiwheelteczm[tTopo->tecWheel(id)-1];
652  if( tTopo->tecWheel(id) == 1 ) {
653  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel1zm_[0] -> Fill((*iter).adc()); meStripTECWheel1zm_[0] ->Fill((*iter).strip()); }
654  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel1zm_[1] -> Fill((*iter).adc()); meStripTECWheel1zm_[1] ->Fill((*iter).strip()); }
655  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel1zm_[2] -> Fill((*iter).adc()); meStripTECWheel1zm_[2] ->Fill((*iter).strip()); }
656  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel1zm_[3] -> Fill((*iter).adc()); meStripTECWheel1zm_[3] ->Fill((*iter).strip()); }
657  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel1zm_[4] -> Fill((*iter).adc()); meStripTECWheel1zm_[4] ->Fill((*iter).strip()); }
658  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel1zm_[5] -> Fill((*iter).adc()); meStripTECWheel1zm_[5] ->Fill((*iter).strip()); }
659  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel1zm_[6] -> Fill((*iter).adc()); meStripTECWheel1zm_[6] ->Fill((*iter).strip()); }
660  }
661  if( tTopo->tecWheel(id) == 2 ) {
662  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel2zm_[0] -> Fill((*iter).adc()); meStripTECWheel2zm_[0] ->Fill((*iter).strip()); }
663  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel2zm_[1] -> Fill((*iter).adc()); meStripTECWheel2zm_[1] ->Fill((*iter).strip()); }
664  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel2zm_[2] -> Fill((*iter).adc()); meStripTECWheel2zm_[2] ->Fill((*iter).strip()); }
665  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel2zm_[3] -> Fill((*iter).adc()); meStripTECWheel2zm_[3] ->Fill((*iter).strip()); }
666  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel2zm_[4] -> Fill((*iter).adc()); meStripTECWheel2zm_[4] ->Fill((*iter).strip()); }
667  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel2zm_[5] -> Fill((*iter).adc()); meStripTECWheel2zm_[5] ->Fill((*iter).strip()); }
668  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel2zm_[6] -> Fill((*iter).adc()); meStripTECWheel2zm_[6] ->Fill((*iter).strip()); }
669  }
670  if( tTopo->tecWheel(id) == 3 ) {
671  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel3zm_[0] -> Fill((*iter).adc()); meStripTECWheel3zm_[0] ->Fill((*iter).strip()); }
672  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel3zm_[1] -> Fill((*iter).adc()); meStripTECWheel3zm_[1] ->Fill((*iter).strip()); }
673  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel3zm_[2] -> Fill((*iter).adc()); meStripTECWheel3zm_[2] ->Fill((*iter).strip()); }
674  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel3zm_[3] -> Fill((*iter).adc()); meStripTECWheel3zm_[3] ->Fill((*iter).strip()); }
675  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel3zm_[4] -> Fill((*iter).adc()); meStripTECWheel3zm_[4] ->Fill((*iter).strip()); }
676  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel3zm_[5] -> Fill((*iter).adc()); meStripTECWheel3zm_[5] ->Fill((*iter).strip()); }
677  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel3zm_[6] -> Fill((*iter).adc()); meStripTECWheel3zm_[6] ->Fill((*iter).strip()); }
678  }
679  if( tTopo->tecWheel(id) == 4 ) {
680  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel4zm_[0] -> Fill((*iter).adc()); meStripTECWheel4zm_[0] ->Fill((*iter).strip()); }
681  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel4zm_[1] -> Fill((*iter).adc()); meStripTECWheel4zm_[1] ->Fill((*iter).strip()); }
682  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel4zm_[2] -> Fill((*iter).adc()); meStripTECWheel4zm_[2] ->Fill((*iter).strip()); }
683  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel4zm_[3] -> Fill((*iter).adc()); meStripTECWheel4zm_[3] ->Fill((*iter).strip()); }
684  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel4zm_[4] -> Fill((*iter).adc()); meStripTECWheel4zm_[4] ->Fill((*iter).strip()); }
685  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel4zm_[5] -> Fill((*iter).adc()); meStripTECWheel4zm_[5] ->Fill((*iter).strip()); }
686  }
687  if( tTopo->tecWheel(id) == 5 ) {
688  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel5zm_[0] -> Fill((*iter).adc()); meStripTECWheel5zm_[0] ->Fill((*iter).strip()); }
689  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel5zm_[1] -> Fill((*iter).adc()); meStripTECWheel5zm_[1] ->Fill((*iter).strip()); }
690  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel5zm_[2] -> Fill((*iter).adc()); meStripTECWheel5zm_[2] ->Fill((*iter).strip()); }
691  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel5zm_[3] -> Fill((*iter).adc()); meStripTECWheel5zm_[3] ->Fill((*iter).strip()); }
692  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel5zm_[4] -> Fill((*iter).adc()); meStripTECWheel5zm_[4] ->Fill((*iter).strip()); }
693  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel5zm_[5] -> Fill((*iter).adc()); meStripTECWheel5zm_[5] ->Fill((*iter).strip()); }
694  }
695  if( tTopo->tecWheel(id) == 6 ) {
696  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel6zm_[0] -> Fill((*iter).adc()); meStripTECWheel6zm_[0] ->Fill((*iter).strip()); }
697  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel6zm_[1] -> Fill((*iter).adc()); meStripTECWheel6zm_[1] ->Fill((*iter).strip()); }
698  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel6zm_[2] -> Fill((*iter).adc()); meStripTECWheel6zm_[2] ->Fill((*iter).strip()); }
699  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel6zm_[3] -> Fill((*iter).adc()); meStripTECWheel6zm_[3] ->Fill((*iter).strip()); }
700  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel6zm_[4] -> Fill((*iter).adc()); meStripTECWheel6zm_[4] ->Fill((*iter).strip()); }
701  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel6zm_[5] -> Fill((*iter).adc()); meStripTECWheel6zm_[5] ->Fill((*iter).strip()); }
702  }
703  if( tTopo->tecWheel(id) == 7 ) {
704  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel7zm_[0] -> Fill((*iter).adc()); meStripTECWheel7zm_[0] ->Fill((*iter).strip()); }
705  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel7zm_[1] -> Fill((*iter).adc()); meStripTECWheel7zm_[1] ->Fill((*iter).strip()); }
706  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel7zm_[2] -> Fill((*iter).adc()); meStripTECWheel7zm_[2] ->Fill((*iter).strip()); }
707  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel7zm_[3] -> Fill((*iter).adc()); meStripTECWheel7zm_[3] ->Fill((*iter).strip()); }
708  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel7zm_[4] -> Fill((*iter).adc()); meStripTECWheel7zm_[4] ->Fill((*iter).strip()); }
709  }
710  if( tTopo->tecWheel(id) == 8 ) {
711  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel8zm_[0] -> Fill((*iter).adc()); meStripTECWheel8zm_[0] ->Fill((*iter).strip()); }
712  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel8zm_[1] -> Fill((*iter).adc()); meStripTECWheel8zm_[1] ->Fill((*iter).strip()); }
713  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel8zm_[2] -> Fill((*iter).adc()); meStripTECWheel8zm_[2] ->Fill((*iter).strip()); }
714  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel8zm_[3] -> Fill((*iter).adc()); meStripTECWheel8zm_[3] ->Fill((*iter).strip()); }
715  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel8zm_[4] -> Fill((*iter).adc()); meStripTECWheel8zm_[4] ->Fill((*iter).strip()); }
716  }
717  if( tTopo->tecWheel(id) == 9 ) {
718  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel9zm_[0] -> Fill((*iter).adc()); meStripTECWheel9zm_[0] ->Fill((*iter).strip()); }
719  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel9zm_[1] -> Fill((*iter).adc()); meStripTECWheel9zm_[1] ->Fill((*iter).strip()); }
720  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel9zm_[2] -> Fill((*iter).adc()); meStripTECWheel9zm_[2] ->Fill((*iter).strip()); }
721  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel9zm_[3] -> Fill((*iter).adc()); meStripTECWheel9zm_[3] ->Fill((*iter).strip()); }
722  }
723  }else {
724  ++ndigiwheelteczp[tTopo->tecWheel(id)-1];
725  if( tTopo->tecWheel(id) == 1 ) {
726  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel1zp_[0] -> Fill((*iter).adc()); meStripTECWheel1zp_[0] ->Fill((*iter).strip()); }
727  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel1zp_[1] -> Fill((*iter).adc()); meStripTECWheel1zp_[1] ->Fill((*iter).strip()); }
728  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel1zp_[2] -> Fill((*iter).adc()); meStripTECWheel1zp_[2] ->Fill((*iter).strip()); }
729  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel1zp_[3] -> Fill((*iter).adc()); meStripTECWheel1zp_[3] ->Fill((*iter).strip()); }
730  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel1zp_[4] -> Fill((*iter).adc()); meStripTECWheel1zp_[4] ->Fill((*iter).strip()); }
731  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel1zp_[5] -> Fill((*iter).adc()); meStripTECWheel1zp_[5] ->Fill((*iter).strip()); }
732  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel1zp_[6] -> Fill((*iter).adc()); meStripTECWheel1zp_[6] ->Fill((*iter).strip()); }
733  }
734  if( tTopo->tecWheel(id) == 2 ) {
735  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel2zp_[0] -> Fill((*iter).adc()); meStripTECWheel2zp_[0] ->Fill((*iter).strip()); }
736  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel2zp_[1] -> Fill((*iter).adc()); meStripTECWheel2zp_[1] ->Fill((*iter).strip()); }
737  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel2zp_[2] -> Fill((*iter).adc()); meStripTECWheel2zp_[2] ->Fill((*iter).strip()); }
738  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel2zp_[3] -> Fill((*iter).adc()); meStripTECWheel2zp_[3] ->Fill((*iter).strip()); }
739  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel2zp_[4] -> Fill((*iter).adc()); meStripTECWheel2zp_[4] ->Fill((*iter).strip()); }
740  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel2zp_[5] -> Fill((*iter).adc()); meStripTECWheel2zp_[5] ->Fill((*iter).strip()); }
741  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel2zp_[6] -> Fill((*iter).adc()); meStripTECWheel2zp_[6] ->Fill((*iter).strip()); }
742  }
743  if( tTopo->tecWheel(id) == 3 ) {
744  if ( tTopo->tecRing(id) == 1 ) { meAdcTECWheel3zp_[0] -> Fill((*iter).adc()); meStripTECWheel3zp_[0] ->Fill((*iter).strip()); }
745  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel3zp_[1] -> Fill((*iter).adc()); meStripTECWheel3zp_[1] ->Fill((*iter).strip()); }
746  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel3zp_[2] -> Fill((*iter).adc()); meStripTECWheel3zp_[2] ->Fill((*iter).strip()); }
747  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel3zp_[3] -> Fill((*iter).adc()); meStripTECWheel3zp_[3] ->Fill((*iter).strip()); }
748  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel3zp_[4] -> Fill((*iter).adc()); meStripTECWheel3zp_[4] ->Fill((*iter).strip()); }
749  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel3zp_[5] -> Fill((*iter).adc()); meStripTECWheel3zp_[5] ->Fill((*iter).strip()); }
750  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel3zp_[6] -> Fill((*iter).adc()); meStripTECWheel3zp_[6] ->Fill((*iter).strip()); }
751  }
752  if( tTopo->tecWheel(id) == 4 ) {
753  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel4zp_[0] -> Fill((*iter).adc()); meStripTECWheel4zp_[0] ->Fill((*iter).strip()); }
754  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel4zp_[1] -> Fill((*iter).adc()); meStripTECWheel4zp_[1] ->Fill((*iter).strip()); }
755  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel4zp_[2] -> Fill((*iter).adc()); meStripTECWheel4zp_[2] ->Fill((*iter).strip()); }
756  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel4zp_[3] -> Fill((*iter).adc()); meStripTECWheel4zp_[3] ->Fill((*iter).strip()); }
757  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel4zp_[4] -> Fill((*iter).adc()); meStripTECWheel4zp_[4] ->Fill((*iter).strip()); }
758  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel4zp_[5] -> Fill((*iter).adc()); meStripTECWheel4zp_[5] ->Fill((*iter).strip()); }
759  }
760  if( tTopo->tecWheel(id) == 5 ) {
761  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel5zp_[0] -> Fill((*iter).adc()); meStripTECWheel5zp_[0] ->Fill((*iter).strip()); }
762  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel5zp_[1] -> Fill((*iter).adc()); meStripTECWheel5zp_[1] ->Fill((*iter).strip()); }
763  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel5zp_[2] -> Fill((*iter).adc()); meStripTECWheel5zp_[2] ->Fill((*iter).strip()); }
764  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel5zp_[3] -> Fill((*iter).adc()); meStripTECWheel5zp_[3] ->Fill((*iter).strip()); }
765  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel5zp_[4] -> Fill((*iter).adc()); meStripTECWheel5zp_[4] ->Fill((*iter).strip()); }
766  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel5zp_[5] -> Fill((*iter).adc()); meStripTECWheel5zp_[5] ->Fill((*iter).strip()); }
767  }
768  if( tTopo->tecWheel(id) == 6 ) {
769  if ( tTopo->tecRing(id) == 2 ) { meAdcTECWheel6zp_[0] -> Fill((*iter).adc()); meStripTECWheel6zp_[0] ->Fill((*iter).strip()); }
770  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel6zp_[1] -> Fill((*iter).adc()); meStripTECWheel6zp_[1] ->Fill((*iter).strip()); }
771  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel6zp_[2] -> Fill((*iter).adc()); meStripTECWheel6zp_[2] ->Fill((*iter).strip()); }
772  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel6zp_[3] -> Fill((*iter).adc()); meStripTECWheel6zp_[3] ->Fill((*iter).strip()); }
773  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel6zp_[4] -> Fill((*iter).adc()); meStripTECWheel6zp_[4] ->Fill((*iter).strip()); }
774  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel6zp_[5] -> Fill((*iter).adc()); meStripTECWheel6zp_[5] ->Fill((*iter).strip()); }
775  }
776  if( tTopo->tecWheel(id) == 7 ) {
777  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel7zp_[0] -> Fill((*iter).adc()); meStripTECWheel7zp_[0] ->Fill((*iter).strip()); }
778  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel7zp_[1] -> Fill((*iter).adc()); meStripTECWheel7zp_[1] ->Fill((*iter).strip()); }
779  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel7zp_[2] -> Fill((*iter).adc()); meStripTECWheel7zp_[2] ->Fill((*iter).strip()); }
780  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel7zp_[3] -> Fill((*iter).adc()); meStripTECWheel7zp_[3] ->Fill((*iter).strip()); }
781  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel7zp_[4] -> Fill((*iter).adc()); meStripTECWheel7zp_[4] ->Fill((*iter).strip()); }
782  }
783  if( tTopo->tecWheel(id) == 8 ) {
784  if ( tTopo->tecRing(id) == 3 ) { meAdcTECWheel8zp_[0] -> Fill((*iter).adc()); meStripTECWheel8zp_[0] ->Fill((*iter).strip()); }
785  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel8zp_[1] -> Fill((*iter).adc()); meStripTECWheel8zp_[1] ->Fill((*iter).strip()); }
786  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel8zp_[2] -> Fill((*iter).adc()); meStripTECWheel8zp_[2] ->Fill((*iter).strip()); }
787  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel8zp_[3] -> Fill((*iter).adc()); meStripTECWheel8zp_[3] ->Fill((*iter).strip()); }
788  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel8zp_[4] -> Fill((*iter).adc()); meStripTECWheel8zp_[4] ->Fill((*iter).strip()); }
789  }
790  if( tTopo->tecWheel(id) == 9 ) {
791  if ( tTopo->tecRing(id) == 4 ) { meAdcTECWheel9zp_[0] -> Fill((*iter).adc()); meStripTECWheel9zp_[0] ->Fill((*iter).strip()); }
792  if ( tTopo->tecRing(id) == 5 ) { meAdcTECWheel9zp_[1] -> Fill((*iter).adc()); meStripTECWheel9zp_[1] ->Fill((*iter).strip()); }
793  if ( tTopo->tecRing(id) == 6 ) { meAdcTECWheel9zp_[2] -> Fill((*iter).adc()); meStripTECWheel9zp_[2] ->Fill((*iter).strip()); }
794  if ( tTopo->tecRing(id) == 7 ) { meAdcTECWheel9zp_[3] -> Fill((*iter).adc()); meStripTECWheel9zp_[3] ->Fill((*iter).strip()); }
795  }
796  }
797  }
798  }
799 
800  }
801 
802  for ( int i =0; i< 4; i++ ) {
803  meNDigiTIBLayerzm_[i]->Fill(ndigilayertibzm[i]);
804  meNDigiTIBLayerzp_[i]->Fill(ndigilayertibzp[i]);
805  }
806 
807  for ( int i =0; i< 6; i++ ) {
808  meNDigiTOBLayerzm_[i]->Fill(ndigilayertobzm[i]);
809  meNDigiTOBLayerzp_[i]->Fill(ndigilayertobzp[i]);
810  }
811 
812 for ( int i =0; i< 3; i++ ) {
813  meNDigiTIDWheelzm_[i]->Fill(ndigiwheeltidzm[i]);
814  meNDigiTIDWheelzp_[i]->Fill(ndigiwheeltidzp[i]);
815  }
816 
817 for ( int i =0; i< 9; i++ ) {
818  meNDigiTECWheelzm_[i]->Fill(ndigiwheelteczm[i]);
819  meNDigiTECWheelzp_[i]->Fill(ndigiwheelteczp[i]);
820  }
821 
822 
823 
824 }
825 
MonitorElement * meAdcTOBLayer5zp_[6]
MonitorElement * meAdcTIBLayer4Extzm_[3]
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
int i
Definition: DBlmapReader.cc:9
MonitorElement * meStripTIBLayer4Intzm_[3]
MonitorElement * meAdcTIDWheel3zm_[3]
MonitorElement * meStripTECWheel1zm_[7]
MonitorElement * meAdcTECWheel6zp_[6]
MonitorElement * meStripTECWheel2zm_[7]
MonitorElement * meNDigiTIDWheelzm_[3]
MonitorElement * meAdcTIBLayer2Extzp_[3]
MonitorElement * meAdcTECWheel1zp_[7]
MonitorElement * meStripTECWheel8zp_[5]
MonitorElement * meStripTOBLayer2zm_[6]
MonitorElement * meStripTIBLayer2Extzm_[3]
unsigned int tibLayer(const DetId &id) const
MonitorElement * meStripTECWheel7zp_[5]
MonitorElement * meStripTIDWheel1zp_[3]
unsigned int tidRing(const DetId &id) const
std::string outputFile_
MonitorElement * meAdcTECWheel4zm_[6]
MonitorElement * meAdcTOBLayer4zp_[6]
MonitorElement * meStripTIBLayer4Intzp_[3]
MonitorElement * meAdcTIBLayer2Extzm_[3]
MonitorElement * meStripTOBLayer1zp_[6]
MonitorElement * meStripTECWheel9zp_[4]
MonitorElement * meStripTOBLayer3zp_[6]
MonitorElement * meNDigiTECWheelzm_[9]
MonitorElement * meAdcTECWheel9zp_[4]
MonitorElement * meAdcTOBLayer5zm_[6]
unsigned int tecRing(const DetId &id) const
ring id
MonitorElement * meStripTIBLayer2Intzp_[3]
MonitorElement * meAdcTECWheel3zm_[7]
MonitorElement * meStripTOBLayer4zm_[6]
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2118
MonitorElement * meStripTECWheel4zp_[6]
MonitorElement * meNDigiTIBLayerzp_[4]
MonitorElement * meNDigiTOBLayerzp_[6]
MonitorElement * meAdcTIBLayer1Intzp_[3]
MonitorElement * meStripTECWheel5zp_[6]
unsigned int tidWheel(const DetId &id) const
MonitorElement * meAdcTIDWheel2zm_[3]
MonitorElement * meStripTOBLayer4zp_[6]
MonitorElement * meAdcTIBLayer2Intzp_[3]
MonitorElement * meAdcTECWheel8zm_[5]
MonitorElement * meStripTIBLayer3Extzp_[3]
MonitorElement * meAdcTIBLayer2Intzm_[3]
MonitorElement * meAdcTECWheel5zm_[6]
std::vector< unsigned int > tibStringInfo(const DetId &id) const
MonitorElement * meNDigiTOBLayerzm_[6]
MonitorElement * meStripTOBLayer5zp_[6]
MonitorElement * meStripTECWheel3zp_[7]
void Fill(long long x)
MonitorElement * meStripTECWheel6zm_[6]
MonitorElement * meStripTIBLayer3Intzm_[3]
MonitorElement * meAdcTOBLayer2zp_[6]
MonitorElement * meAdcTIDWheel3zp_[3]
MonitorElement * meStripTIBLayer2Intzm_[3]
MonitorElement * meAdcTIBLayer1Extzp_[3]
MonitorElement * meStripTECWheel4zm_[6]
MonitorElement * meStripTIBLayer4Extzm_[3]
MonitorElement * meAdcTOBLayer1zp_[6]
unsigned int tidSide(const DetId &id) const
std::vector< unsigned int > tobRodInfo(const DetId &id) const
MonitorElement * meAdcTOBLayer4zm_[6]
void analyze(const edm::Event &e, const edm::EventSetup &c)
MonitorElement * meStripTOBLayer5zm_[6]
MonitorElement * meAdcTECWheel8zp_[5]
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
MonitorElement * meAdcTIDWheel1zp_[3]
MonitorElement * meAdcTIBLayer1Intzm_[3]
MonitorElement * meStripTECWheel9zm_[4]
MonitorElement * meAdcTECWheel7zp_[5]
MonitorElement * meAdcTECWheel1zm_[7]
MonitorElement * meStripTIBLayer1Intzp_[3]
#define end
Definition: vmac.h:38
MonitorElement * meStripTIDWheel2zm_[3]
MonitorElement * meStripTOBLayer3zm_[6]
MonitorElement * meStripTIDWheel1zm_[3]
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
MonitorElement * meAdcTIBLayer3Extzp_[3]
MonitorElement * meStripTIBLayer3Extzm_[3]
MonitorElement * meStripTECWheel7zm_[5]
MonitorElement * meStripTECWheel2zp_[7]
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:356
MonitorElement * meAdcTIBLayer3Intzm_[3]
MonitorElement * meNDigiTIBLayerzm_[4]
DQMStore * dbe_
MonitorElement * meAdcTECWheel3zp_[7]
unsigned int tibModule(const DetId &id) const
MonitorElement * meAdcTIDWheel1zm_[3]
MonitorElement * meAdcTECWheel7zm_[5]
MonitorElement * meAdcTECWheel5zp_[6]
MonitorElement * meAdcTOBLayer1zm_[6]
Definition: DetId.h:20
MonitorElement * meStripTECWheel8zm_[5]
MonitorElement * meStripTECWheel5zm_[6]
MonitorElement * meStripTIBLayer1Intzm_[3]
MonitorElement * meAdcTIBLayer4Intzp_[3]
MonitorElement * meStripTIBLayer1Extzp_[3]
MonitorElement * meStripTIBLayer1Extzm_[3]
MonitorElement * meNDigiTIDWheelzp_[3]
MonitorElement * meAdcTOBLayer3zp_[6]
const T & get() const
Definition: EventSetup.h:55
MonitorElement * meAdcTOBLayer6zm_[6]
SiStripDigiValid(const edm::ParameterSet &ps)
T const * product() const
Definition: ESHandle.h:62
MonitorElement * meStripTIDWheel3zp_[3]
MonitorElement * meStripTOBLayer6zp_[6]
MonitorElement * meAdcTOBLayer3zm_[6]
unsigned int tobModule(const DetId &id) const
MonitorElement * meAdcTIBLayer4Intzm_[3]
MonitorElement * meAdcTECWheel2zp_[7]
MonitorElement * meAdcTOBLayer6zp_[6]
MonitorElement * meAdcTIBLayer3Extzm_[3]
MonitorElement * meAdcTIDWheel2zp_[3]
#define begin
Definition: vmac.h:31
MonitorElement * meStripTOBLayer2zp_[6]
MonitorElement * meStripTIDWheel2zp_[3]
MonitorElement * meStripTIBLayer4Extzp_[3]
MonitorElement * meStripTOBLayer1zm_[6]
MonitorElement * meAdcTIBLayer1Extzm_[3]
MonitorElement * meStripTECWheel3zm_[7]
MonitorElement * meAdcTOBLayer2zm_[6]
MonitorElement * meAdcTIBLayer4Extzp_[3]
MonitorElement * meStripTECWheel6zp_[6]
iterator begin()
Return an iterator to the first DetSet.
Definition: DetSetVector.h:341
MonitorElement * meStripTOBLayer6zm_[6]
MonitorElement * meNDigiTECWheelzp_[9]
MonitorElement * meAdcTIBLayer3Intzp_[3]
collection_type::const_iterator const_iterator
Definition: DetSet.h:34
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:106
MonitorElement * meAdcTECWheel6zm_[6]
edm::InputTag src_
MonitorElement * meStripTECWheel1zp_[7]
unsigned int tecWheel(const DetId &id) const
MonitorElement * meStripTIBLayer3Intzp_[3]
MonitorElement * meAdcTECWheel4zp_[6]
MonitorElement * meAdcTECWheel9zm_[4]
MonitorElement * meStripTIBLayer2Extzp_[3]
unsigned int tobLayer(const DetId &id) const
unsigned int tecSide(const DetId &id) const
MonitorElement * meAdcTECWheel2zm_[7]
MonitorElement * meStripTIDWheel3zm_[3]