CMS 3D CMS Logo

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