CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
SiStripDigiValid.cc
Go to the documentation of this file.
21 #include "SiStripDigiValid.h"
22 
24  : m_topoToken(esConsumes()),
25  dbe_(nullptr),
26  runStandalone(ps.getParameter<bool>("runStandalone")),
27  outputFile_(ps.getUntrackedParameter<std::string>("outputFile", "stripdigihisto.root")),
28  edmDetSetVector_SiStripDigi_Token_(
29  consumes<edm::DetSetVector<SiStripDigi>>(ps.getParameter<edm::InputTag>("src"))) {}
30 
32 
35 
36  if (dbe_) {
37  ibooker.setCurrentFolder("TrackerDigisV/TrackerDigis/Strip");
38 
39  for (int i = 0; i < 3; i++) {
40  Char_t histo[200];
41  // Z Plus Side
42  sprintf(histo, "adc_tib_layer1_extmodule%d_zp", i + 1);
43  meAdcTIBLayer1Extzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
44  sprintf(histo, "adc_tib_layer1_intmodule%d_zp", i + 1);
45  meAdcTIBLayer1Intzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
46  sprintf(histo, "adc_tib_layer2_extmodule%d_zp", i + 1);
47  meAdcTIBLayer2Extzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
48  sprintf(histo, "adc_tib_layer2_intmodule%d_zp", i + 1);
49  meAdcTIBLayer2Intzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
50  sprintf(histo, "adc_tib_layer3_extmodule%d_zp", i + 1);
51  meAdcTIBLayer3Extzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
52  sprintf(histo, "adc_tib_layer3_intmodule%d_zp", i + 1);
53  meAdcTIBLayer3Intzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
54  sprintf(histo, "adc_tib_layer4_extmodule%d_zp", i + 1);
55  meAdcTIBLayer4Extzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
56  sprintf(histo, "adc_tib_layer4_intmodule%d_zp", i + 1);
57  meAdcTIBLayer4Intzp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
58  sprintf(histo, "strip_tib_layer1_extmodule%d_zp", i + 1);
59  meStripTIBLayer1Extzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
60  sprintf(histo, "strip_tib_layer1_intmodule%d_zp", i + 1);
61  meStripTIBLayer1Intzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
62  sprintf(histo, "strip_tib_layer2_extmodule%d_zp", i + 1);
63  meStripTIBLayer2Extzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
64  sprintf(histo, "strip_tib_layer2_intmodule%d_zp", i + 1);
65  meStripTIBLayer2Intzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
66  sprintf(histo, "strip_tib_layer3_extmodule%d_zp", i + 1);
67  meStripTIBLayer3Extzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
68  sprintf(histo, "strip_tib_layer3_intmodule%d_zp", i + 1);
69  meStripTIBLayer3Intzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
70  sprintf(histo, "strip_tib_layer4_extmodule%d_zp", i + 1);
71  meStripTIBLayer4Extzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
72  sprintf(histo, "strip_tib_layer4_intmodule%d_zp", i + 1);
73  meStripTIBLayer4Intzp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
74  // Z Minus Side
75  sprintf(histo, "adc_tib_layer1_extmodule%d_zm", i + 1);
76  meAdcTIBLayer1Extzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
77  sprintf(histo, "adc_tib_layer1_intmodule%d_zm", i + 1);
78  meAdcTIBLayer1Intzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
79  sprintf(histo, "adc_tib_layer2_extmodule%d_zm", i + 1);
80  meAdcTIBLayer2Extzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
81  sprintf(histo, "adc_tib_layer2_intmodule%d_zm", i + 1);
82  meAdcTIBLayer2Intzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
83  sprintf(histo, "adc_tib_layer3_extmodule%d_zm", i + 1);
84  meAdcTIBLayer3Extzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
85  sprintf(histo, "adc_tib_layer3_intmodule%d_zm", i + 1);
86  meAdcTIBLayer3Intzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
87  sprintf(histo, "adc_tib_layer4_extmodule%d_zm", i + 1);
88  meAdcTIBLayer4Extzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
89  sprintf(histo, "adc_tib_layer4_intmodule%d_zm", i + 1);
90  meAdcTIBLayer4Intzm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
91  sprintf(histo, "strip_tib_layer1_extmodule%d_zm", i + 1);
92  meStripTIBLayer1Extzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
93  sprintf(histo, "strip_tib_layer1_intmodule%d_zm", i + 1);
94  meStripTIBLayer1Intzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
95  sprintf(histo, "strip_tib_layer2_extmodule%d_zm", i + 1);
96  meStripTIBLayer2Extzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
97  sprintf(histo, "strip_tib_layer2_intmodule%d_zm", i + 1);
98  meStripTIBLayer2Intzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
99  sprintf(histo, "strip_tib_layer3_extmodule%d_zm", i + 1);
100  meStripTIBLayer3Extzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
101  sprintf(histo, "strip_tib_layer3_intmodule%d_zm", i + 1);
102  meStripTIBLayer3Intzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
103  sprintf(histo, "strip_tib_layer4_extmodule%d_zm", i + 1);
104  meStripTIBLayer4Extzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
105  sprintf(histo, "strip_tib_layer4_intmodule%d_zm", i + 1);
106  meStripTIBLayer4Intzm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
107  }
108 
109  for (int i = 0; i < 6; i++) {
110  Char_t histo[200];
111  // Z Plus Side
112  sprintf(histo, "adc_tob_layer1_module%d_zp", i + 1);
113  meAdcTOBLayer1zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
114  sprintf(histo, "strip_tob_layer1_module%d_zp", i + 1);
115  meStripTOBLayer1zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
116  sprintf(histo, "adc_tob_layer2_module%d_zp", i + 1);
117  meAdcTOBLayer2zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
118  sprintf(histo, "strip_tob_layer2_module%d_zp", i + 1);
119  meStripTOBLayer2zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
120  sprintf(histo, "adc_tob_layer3_module%d_zp", i + 1);
121  meAdcTOBLayer3zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
122  sprintf(histo, "strip_tob_layer3_module%d_zp", i + 1);
123  meStripTOBLayer3zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
124  sprintf(histo, "adc_tob_layer4_module%d_zp", i + 1);
125  meAdcTOBLayer4zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
126  sprintf(histo, "strip_tob_layer4_module%d_zp", i + 1);
127  meStripTOBLayer4zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
128  sprintf(histo, "adc_tob_layer5_module%d_zp", i + 1);
129  meAdcTOBLayer5zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
130  sprintf(histo, "strip_tob_layer5_module%d_zp", i + 1);
131  meStripTOBLayer5zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
132  sprintf(histo, "adc_tob_layer6_module%d_zp", i + 1);
133  meAdcTOBLayer6zp_[i] = ibooker.book1D(histo, "Digis ADC", 10, 0., 300.);
134  sprintf(histo, "strip_tob_layer6_module%d_zp", i + 1);
135  meStripTOBLayer6zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
136  // Z Minus Side
137  sprintf(histo, "adc_tob_layer1_module%d_zm", i + 1);
138  meAdcTOBLayer1zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
139  sprintf(histo, "strip_tob_layer1_module%d_zm", i + 1);
140  meStripTOBLayer1zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
141  sprintf(histo, "adc_tob_layer2_module%d_zm", i + 1);
142  meAdcTOBLayer2zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
143  sprintf(histo, "strip_tob_layer2_module%d_zm", i + 1);
144  meStripTOBLayer2zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
145  sprintf(histo, "adc_tob_layer3_module%d_zm", i + 1);
146  meAdcTOBLayer3zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
147  sprintf(histo, "strip_tob_layer3_module%d_zm", i + 1);
148  meStripTOBLayer3zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
149  sprintf(histo, "adc_tob_layer4_module%d_zm", i + 1);
150  meAdcTOBLayer4zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
151  sprintf(histo, "strip_tob_layer4_module%d_zm", i + 1);
152  meStripTOBLayer4zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
153  sprintf(histo, "adc_tob_layer5_module%d_zm", i + 1);
154  meAdcTOBLayer5zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
155  sprintf(histo, "strip_tob_layer5_module%d_zm", i + 1);
156  meStripTOBLayer5zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
157  sprintf(histo, "adc_tob_layer6_module%d_zm", i + 1);
158  meAdcTOBLayer6zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
159  sprintf(histo, "strip_tob_layer6_module%d_zm", i + 1);
160  meStripTOBLayer6zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
161  }
162 
163  for (int i = 0; i < 3; i++) {
164  Char_t histo[200];
165  // Z Plus Side
166  sprintf(histo, "adc_tid_wheel1_ring%d_zp", i + 1);
167  meAdcTIDWheel1zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
168  sprintf(histo, "adc_tid_wheel2_ring%d_zp", i + 1);
169  meAdcTIDWheel2zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
170  sprintf(histo, "adc_tid_wheel3_ring%d_zp", i + 1);
171  meAdcTIDWheel3zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
172  sprintf(histo, "strip_tid_wheel1_ring%d_zp", i + 1);
173  meStripTIDWheel1zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
174  sprintf(histo, "strip_tid_wheel2_ring%d_zp", i + 1);
175  meStripTIDWheel2zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
176  sprintf(histo, "strip_tid_wheel3_ring%d_zp", i + 1);
177  meStripTIDWheel3zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
178  // Z minus Side
179  sprintf(histo, "adc_tid_wheel1_ring%d_zm", i + 1);
180  meAdcTIDWheel1zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
181  sprintf(histo, "adc_tid_wheel2_ring%d_zm", i + 1);
182  meAdcTIDWheel2zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
183  sprintf(histo, "adc_tid_wheel3_ring%d_zm", i + 1);
184  meAdcTIDWheel3zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
185  sprintf(histo, "strip_tid_wheel1_ring%d_zm", i + 1);
186  meStripTIDWheel1zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
187  sprintf(histo, "strip_tid_wheel2_ring%d_zm", i + 1);
188  meStripTIDWheel2zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
189  sprintf(histo, "strip_tid_wheel3_ring%d_zm", i + 1);
190  meStripTIDWheel3zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
191  }
192 
193  for (int i = 0; i < 7; i++) {
194  Char_t histo[200];
195  // Z Plus Side
196  sprintf(histo, "adc_tec_wheel1_ring%d_zp", i + 1);
197  meAdcTECWheel1zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
198  sprintf(histo, "adc_tec_wheel2_ring%d_zp", i + 1);
199  meAdcTECWheel2zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
200  sprintf(histo, "adc_tec_wheel3_ring%d_zp", i + 1);
201  meAdcTECWheel3zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
202  sprintf(histo, "strip_tec_wheel1_ring%d_zp", i + 1);
203  meStripTECWheel1zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
204  sprintf(histo, "strip_tec_wheel2_ring%d_zp", i + 1);
205  meStripTECWheel2zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
206  sprintf(histo, "strip_tec_wheel3_ring%d_zp", i + 1);
207  meStripTECWheel3zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
208 
209  // Z Minus Side
210  sprintf(histo, "adc_tec_wheel1_ring%d_zm", i + 1);
211  meAdcTECWheel1zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
212  sprintf(histo, "adc_tec_wheel2_ring%d_zm", i + 1);
213  meAdcTECWheel2zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
214  sprintf(histo, "adc_tec_wheel3_ring%d_zm", i + 1);
215  meAdcTECWheel3zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
216  sprintf(histo, "strip_tec_wheel1_ring%d_zm", i + 1);
217  meStripTECWheel1zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
218  sprintf(histo, "strip_tec_wheel2_ring%d_zm", i + 1);
219  meStripTECWheel2zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
220  sprintf(histo, "strip_tec_wheel3_ring%d_zm", i + 1);
221  meStripTECWheel3zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
222  }
223 
224  for (int i = 0; i < 6; i++) {
225  Char_t histo[200];
226  // Z Plus Side
227  sprintf(histo, "adc_tec_wheel4_ring%d_zp", i + 1);
228  meAdcTECWheel4zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
229  sprintf(histo, "adc_tec_wheel5_ring%d_zp", i + 1);
230  meAdcTECWheel5zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
231  sprintf(histo, "adc_tec_wheel6_ring%d_zp", i + 1);
232  meAdcTECWheel6zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
233  sprintf(histo, "strip_tec_wheel4_ring%d_zp", i + 1);
234  meStripTECWheel4zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
235  sprintf(histo, "strip_tec_wheel5_ring%d_zp", i + 1);
236  meStripTECWheel5zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
237  sprintf(histo, "strip_tec_wheel6_ring%d_zp", i + 1);
238  meStripTECWheel6zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
239 
240  // Z Minus Side
241  sprintf(histo, "adc_tec_wheel4_ring%d_zm", i + 1);
242  meAdcTECWheel4zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
243  sprintf(histo, "adc_tec_wheel5_ring%d_zm", i + 1);
244  meAdcTECWheel5zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
245  sprintf(histo, "adc_tec_wheel6_ring%d_zm", i + 1);
246  meAdcTECWheel6zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
247  sprintf(histo, "strip_tec_wheel4_ring%d_zm", i + 1);
248  meStripTECWheel4zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
249  sprintf(histo, "strip_tec_wheel5_ring%d_zm", i + 1);
250  meStripTECWheel5zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
251  sprintf(histo, "strip_tec_wheel6_ring%d_zm", i + 1);
252  meStripTECWheel6zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
253  }
254 
255  for (int i = 0; i < 5; i++) {
256  Char_t histo[200];
257  // Z Plus Side
258  sprintf(histo, "adc_tec_wheel7_ring%d_zp", i + 1);
259  meAdcTECWheel7zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
260  sprintf(histo, "adc_tec_wheel8_ring%d_zp", i + 1);
261  meAdcTECWheel8zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
262  sprintf(histo, "strip_tec_wheel7_ring%d_zp", i + 1);
263  meStripTECWheel7zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
264  sprintf(histo, "strip_tec_wheel8_ring%d_zp", i + 1);
265  meStripTECWheel8zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
266 
267  // Z Minus Side
268  sprintf(histo, "adc_tec_wheel7_ring%d_zm", i + 1);
269  meAdcTECWheel7zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
270  sprintf(histo, "adc_tec_wheel8_ring%d_zm", i + 1);
271  meAdcTECWheel8zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
272  sprintf(histo, "strip_tec_wheel7_ring%d_zm", i + 1);
273  meStripTECWheel7zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
274  sprintf(histo, "strip_tec_wheel8_ring%d_zm", i + 1);
275  meStripTECWheel8zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
276  }
277 
278  for (int i = 0; i < 4; i++) {
279  Char_t histo[200];
280  // Z Plus Side
281  sprintf(histo, "adc_tec_wheel9_ring%d_zp", i + 1);
282  meAdcTECWheel9zp_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
283  sprintf(histo, "strip_tec_wheel9_ring%d_zp", i + 1);
284  meStripTECWheel9zp_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
285 
286  // Z Minus Side
287  sprintf(histo, "adc_tec_wheel9_ring%d_zm", i + 1);
288  meAdcTECWheel9zm_[i] = ibooker.book1D(histo, "Digis ADC", 50, 0., 300.);
289  sprintf(histo, "strip_tec_wheel9_ring%d_zm", i + 1);
290  meStripTECWheel9zm_[i] = ibooker.book1D(histo, "Digis Strip Num.", 200, 0., 800.);
291  }
292 
293  for (int i = 0; i < 4; i++) {
294  Char_t histo[200];
295  sprintf(histo, "ndigi_tib_layer_%d_zm", i + 1);
296  meNDigiTIBLayerzm_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
297  sprintf(histo, "ndigi_tib_layer_%d_zp", i + 1);
298  meNDigiTIBLayerzp_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
299  }
300 
301  for (int i = 0; i < 6; i++) {
302  Char_t histo[200];
303  sprintf(histo, "ndigi_tob_layer_%d_zm", i + 1);
304  meNDigiTOBLayerzm_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
305  sprintf(histo, "ndigi_tob_layer_%d_zp", i + 1);
306  meNDigiTOBLayerzp_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
307  }
308 
309  for (int i = 0; i < 3; i++) {
310  Char_t histo[200];
311  sprintf(histo, "ndigi_tid_wheel_%d_zm", i + 1);
312  meNDigiTIDWheelzm_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
313  sprintf(histo, "ndigi_tid_wheel_%d_zp", i + 1);
314  meNDigiTIDWheelzp_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
315  }
316 
317  for (int i = 0; i < 9; i++) {
318  Char_t histo[200];
319  sprintf(histo, "ndigi_tec_wheel_%d_zm", i + 1);
320  meNDigiTECWheelzm_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
321  sprintf(histo, "ndigi_tec_wheel_%d_zp", i + 1);
322  meNDigiTECWheelzp_[i] = ibooker.book1D(histo, "Digi Multiplicity", 100, 0., 500.);
323  }
324  }
325 }
326 
328  // Retrieve tracker topology from geometry
329  const TrackerTopology *const tTopo = &c.getData(m_topoToken);
330 
331  int ndigilayertibzp[4];
332  int ndigilayertibzm[4];
333 
334  for (int i = 0; i < 4; i++) {
335  ndigilayertibzp[i] = 0;
336  ndigilayertibzm[i] = 0;
337  }
338 
339  int ndigilayertobzp[6];
340  int ndigilayertobzm[6];
341 
342  for (int i = 0; i < 6; i++) {
343  ndigilayertobzp[i] = 0;
344  ndigilayertobzm[i] = 0;
345  }
346 
347  int ndigiwheeltidzp[3];
348  int ndigiwheeltidzm[3];
349 
350  for (int i = 0; i < 3; i++) {
351  ndigiwheeltidzp[i] = 0;
352  ndigiwheeltidzm[i] = 0;
353  }
354 
355  int ndigiwheelteczp[9];
356  int ndigiwheelteczm[9];
357 
358  for (int i = 0; i < 9; i++) {
359  ndigiwheelteczp[i] = 0;
360  ndigiwheelteczm[i] = 0;
361  }
362 
363  // LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " <<
364  // e.id().event();
365 
366  std::string digiProducer = "siStripDigis";
369  edm::DetSetVector<SiStripDigi>::const_iterator DSViter = stripDigis->begin();
370  for (; DSViter != stripDigis->end(); DSViter++) {
371  unsigned int id = DSViter->id;
372  DetId detId(id);
376 
377  if (detId.subdetId() == StripSubdetector::TIB) {
378  for (iter = begin; iter != end; iter++) { // loop digis
379  if (tTopo->tibStringInfo(id)[0] == 1) {
380  ++ndigilayertibzm[tTopo->tibLayer(id) - 1];
381  if (tTopo->tibLayer(id) == 1) {
382  if (tTopo->tibStringInfo(id)[1] == 1) {
383  if (tTopo->tibModule(id) == 1) {
384  meAdcTIBLayer1Intzm_[0]->Fill((*iter).adc());
385  meStripTIBLayer1Intzm_[0]->Fill((*iter).strip());
386  }
387  if (tTopo->tibModule(id) == 2) {
388  meAdcTIBLayer1Intzm_[1]->Fill((*iter).adc());
389  meStripTIBLayer1Intzm_[1]->Fill((*iter).strip());
390  }
391  if (tTopo->tibModule(id) == 3) {
392  meAdcTIBLayer1Intzm_[2]->Fill((*iter).adc());
393  meStripTIBLayer1Intzm_[2]->Fill((*iter).strip());
394  }
395  } else {
396  if (tTopo->tibModule(id) == 1) {
397  meAdcTIBLayer1Extzm_[0]->Fill((*iter).adc());
398  meStripTIBLayer1Extzm_[0]->Fill((*iter).strip());
399  }
400  if (tTopo->tibModule(id) == 2) {
401  meAdcTIBLayer1Extzm_[1]->Fill((*iter).adc());
402  meStripTIBLayer1Extzm_[1]->Fill((*iter).strip());
403  }
404  if (tTopo->tibModule(id) == 3) {
405  meAdcTIBLayer1Extzm_[2]->Fill((*iter).adc());
406  meStripTIBLayer1Extzm_[2]->Fill((*iter).strip());
407  }
408  }
409  }
410  if (tTopo->tibLayer(id) == 2) {
411  if (tTopo->tibStringInfo(id)[1] == 1) {
412  if (tTopo->tibModule(id) == 1) {
413  meAdcTIBLayer2Intzm_[0]->Fill((*iter).adc());
414  meStripTIBLayer2Intzm_[0]->Fill((*iter).strip());
415  }
416  if (tTopo->tibModule(id) == 2) {
417  meAdcTIBLayer2Intzm_[1]->Fill((*iter).adc());
418  meStripTIBLayer2Intzm_[1]->Fill((*iter).strip());
419  }
420  if (tTopo->tibModule(id) == 3) {
421  meAdcTIBLayer2Intzm_[2]->Fill((*iter).adc());
422  meStripTIBLayer2Intzm_[2]->Fill((*iter).strip());
423  }
424  } else {
425  if (tTopo->tibModule(id) == 1) {
426  meAdcTIBLayer2Extzm_[0]->Fill((*iter).adc());
427  meStripTIBLayer2Extzm_[0]->Fill((*iter).strip());
428  }
429  if (tTopo->tibModule(id) == 2) {
430  meAdcTIBLayer2Extzm_[1]->Fill((*iter).adc());
431  meStripTIBLayer2Extzm_[1]->Fill((*iter).strip());
432  }
433  if (tTopo->tibModule(id) == 3) {
434  meAdcTIBLayer2Extzm_[2]->Fill((*iter).adc());
435  meStripTIBLayer2Extzm_[2]->Fill((*iter).strip());
436  }
437  }
438  }
439  if (tTopo->tibLayer(id) == 3) {
440  if (tTopo->tibStringInfo(id)[1] == 1) {
441  if (tTopo->tibModule(id) == 1) {
442  meAdcTIBLayer3Intzm_[0]->Fill((*iter).adc());
443  meStripTIBLayer3Intzm_[0]->Fill((*iter).strip());
444  }
445  if (tTopo->tibModule(id) == 2) {
446  meAdcTIBLayer3Intzm_[1]->Fill((*iter).adc());
447  meStripTIBLayer3Intzm_[1]->Fill((*iter).strip());
448  }
449  if (tTopo->tibModule(id) == 3) {
450  meAdcTIBLayer3Intzm_[2]->Fill((*iter).adc());
451  meStripTIBLayer3Intzm_[2]->Fill((*iter).strip());
452  }
453  } else {
454  if (tTopo->tibModule(id) == 1) {
455  meAdcTIBLayer3Extzm_[0]->Fill((*iter).adc());
456  meStripTIBLayer3Extzm_[0]->Fill((*iter).strip());
457  }
458  if (tTopo->tibModule(id) == 2) {
459  meAdcTIBLayer3Extzm_[1]->Fill((*iter).adc());
460  meStripTIBLayer3Extzm_[1]->Fill((*iter).strip());
461  }
462  if (tTopo->tibModule(id) == 3) {
463  meAdcTIBLayer3Extzm_[2]->Fill((*iter).adc());
464  meStripTIBLayer3Extzm_[2]->Fill((*iter).strip());
465  }
466  }
467  }
468  if (tTopo->tibLayer(id) == 4) {
469  if (tTopo->tibStringInfo(id)[1] == 1) {
470  if (tTopo->tibModule(id) == 1) {
471  meAdcTIBLayer4Intzm_[0]->Fill((*iter).adc());
472  meStripTIBLayer4Intzm_[0]->Fill((*iter).strip());
473  }
474  if (tTopo->tibModule(id) == 2) {
475  meAdcTIBLayer4Intzm_[1]->Fill((*iter).adc());
476  meStripTIBLayer4Intzm_[1]->Fill((*iter).strip());
477  }
478  if (tTopo->tibModule(id) == 3) {
479  meAdcTIBLayer4Intzm_[2]->Fill((*iter).adc());
480  meStripTIBLayer4Intzm_[2]->Fill((*iter).strip());
481  }
482  } else {
483  if (tTopo->tibModule(id) == 1) {
484  meAdcTIBLayer4Extzm_[0]->Fill((*iter).adc());
485  meStripTIBLayer4Extzm_[0]->Fill((*iter).strip());
486  }
487  if (tTopo->tibModule(id) == 2) {
488  meAdcTIBLayer4Extzm_[1]->Fill((*iter).adc());
489  meStripTIBLayer4Extzm_[1]->Fill((*iter).strip());
490  }
491  if (tTopo->tibModule(id) == 3) {
492  meAdcTIBLayer4Extzm_[2]->Fill((*iter).adc());
493  meStripTIBLayer4Extzm_[2]->Fill((*iter).strip());
494  }
495  }
496  }
497  } else {
498  ++ndigilayertibzp[tTopo->tibLayer(id) - 1];
499  if (tTopo->tibLayer(id) == 1) {
500  if (tTopo->tibStringInfo(id)[1] == 1) {
501  if (tTopo->tibModule(id) == 1) {
502  meAdcTIBLayer1Intzp_[0]->Fill((*iter).adc());
503  meStripTIBLayer1Intzp_[0]->Fill((*iter).strip());
504  }
505  if (tTopo->tibModule(id) == 2) {
506  meAdcTIBLayer1Intzp_[1]->Fill((*iter).adc());
507  meStripTIBLayer1Intzp_[1]->Fill((*iter).strip());
508  }
509  if (tTopo->tibModule(id) == 3) {
510  meAdcTIBLayer1Intzp_[2]->Fill((*iter).adc());
511  meStripTIBLayer1Intzp_[2]->Fill((*iter).strip());
512  }
513  } else {
514  if (tTopo->tibModule(id) == 1) {
515  meAdcTIBLayer1Extzp_[0]->Fill((*iter).adc());
516  meStripTIBLayer1Extzp_[0]->Fill((*iter).strip());
517  }
518  if (tTopo->tibModule(id) == 2) {
519  meAdcTIBLayer1Extzp_[1]->Fill((*iter).adc());
520  meStripTIBLayer1Extzp_[1]->Fill((*iter).strip());
521  }
522  if (tTopo->tibModule(id) == 3) {
523  meAdcTIBLayer1Extzp_[2]->Fill((*iter).adc());
524  meStripTIBLayer1Extzp_[2]->Fill((*iter).strip());
525  }
526  }
527  }
528  if (tTopo->tibLayer(id) == 2) {
529  if (tTopo->tibStringInfo(id)[1] == 1) {
530  if (tTopo->tibModule(id) == 1) {
531  meAdcTIBLayer2Intzp_[0]->Fill((*iter).adc());
532  meStripTIBLayer2Intzp_[0]->Fill((*iter).strip());
533  }
534  if (tTopo->tibModule(id) == 2) {
535  meAdcTIBLayer2Intzp_[1]->Fill((*iter).adc());
536  meStripTIBLayer2Intzp_[1]->Fill((*iter).strip());
537  }
538  if (tTopo->tibModule(id) == 3) {
539  meAdcTIBLayer2Intzp_[2]->Fill((*iter).adc());
540  meStripTIBLayer2Intzp_[2]->Fill((*iter).strip());
541  }
542  } else {
543  if (tTopo->tibModule(id) == 1) {
544  meAdcTIBLayer2Extzp_[0]->Fill((*iter).adc());
545  meStripTIBLayer2Extzp_[0]->Fill((*iter).strip());
546  }
547  if (tTopo->tibModule(id) == 2) {
548  meAdcTIBLayer2Extzp_[1]->Fill((*iter).adc());
549  meStripTIBLayer2Extzp_[1]->Fill((*iter).strip());
550  }
551  if (tTopo->tibModule(id) == 3) {
552  meAdcTIBLayer2Extzp_[2]->Fill((*iter).adc());
553  meStripTIBLayer2Extzp_[2]->Fill((*iter).strip());
554  }
555  }
556  }
557  if (tTopo->tibLayer(id) == 3) {
558  if (tTopo->tibStringInfo(id)[1] == 1) {
559  if (tTopo->tibModule(id) == 1) {
560  meAdcTIBLayer3Intzp_[0]->Fill((*iter).adc());
561  meStripTIBLayer3Intzp_[0]->Fill((*iter).strip());
562  }
563  if (tTopo->tibModule(id) == 2) {
564  meAdcTIBLayer3Intzp_[1]->Fill((*iter).adc());
565  meStripTIBLayer3Intzp_[1]->Fill((*iter).strip());
566  }
567  if (tTopo->tibModule(id) == 3) {
568  meAdcTIBLayer3Intzp_[2]->Fill((*iter).adc());
569  meStripTIBLayer3Intzp_[2]->Fill((*iter).strip());
570  }
571  } else {
572  if (tTopo->tibModule(id) == 1) {
573  meAdcTIBLayer3Extzp_[0]->Fill((*iter).adc());
574  meStripTIBLayer3Extzp_[0]->Fill((*iter).strip());
575  }
576  if (tTopo->tibModule(id) == 2) {
577  meAdcTIBLayer3Extzp_[1]->Fill((*iter).adc());
578  meStripTIBLayer3Extzp_[1]->Fill((*iter).strip());
579  }
580  if (tTopo->tibModule(id) == 3) {
581  meAdcTIBLayer3Extzp_[2]->Fill((*iter).adc());
582  meStripTIBLayer3Extzp_[2]->Fill((*iter).strip());
583  }
584  }
585  }
586  if (tTopo->tibLayer(id) == 4) {
587  if (tTopo->tibStringInfo(id)[1] == 1) {
588  if (tTopo->tibModule(id) == 1) {
589  meAdcTIBLayer4Intzp_[0]->Fill((*iter).adc());
590  meStripTIBLayer4Intzp_[0]->Fill((*iter).strip());
591  }
592  if (tTopo->tibModule(id) == 2) {
593  meAdcTIBLayer4Intzp_[1]->Fill((*iter).adc());
594  meStripTIBLayer4Intzp_[1]->Fill((*iter).strip());
595  }
596  if (tTopo->tibModule(id) == 3) {
597  meAdcTIBLayer4Intzp_[2]->Fill((*iter).adc());
598  meStripTIBLayer4Intzp_[2]->Fill((*iter).strip());
599  }
600  } else {
601  if (tTopo->tibModule(id) == 1) {
602  meAdcTIBLayer4Extzp_[0]->Fill((*iter).adc());
603  meStripTIBLayer4Extzp_[0]->Fill((*iter).strip());
604  }
605  if (tTopo->tibModule(id) == 2) {
606  meAdcTIBLayer4Extzp_[1]->Fill((*iter).adc());
607  meStripTIBLayer4Extzp_[1]->Fill((*iter).strip());
608  }
609  if (tTopo->tibModule(id) == 3) {
610  meAdcTIBLayer4Extzp_[2]->Fill((*iter).adc());
611  meStripTIBLayer4Extzp_[2]->Fill((*iter).strip());
612  }
613  }
614  }
615  }
616  }
617  }
618  if (detId.subdetId() == StripSubdetector::TOB) {
619  for (iter = begin; iter != end; iter++) { // loop digis
620  if (tTopo->tobRodInfo(id)[0] == 1) {
621  ++ndigilayertobzm[tTopo->tobLayer(id) - 1];
622  if (tTopo->tobLayer(id) == 1) {
623  if (tTopo->tobModule(id) == 1) {
624  meAdcTOBLayer1zm_[0]->Fill((*iter).adc());
625  meStripTOBLayer1zm_[0]->Fill((*iter).strip());
626  }
627  if (tTopo->tobModule(id) == 2) {
628  meAdcTOBLayer1zm_[1]->Fill((*iter).adc());
629  meStripTOBLayer1zm_[1]->Fill((*iter).strip());
630  }
631  if (tTopo->tobModule(id) == 3) {
632  meAdcTOBLayer1zm_[2]->Fill((*iter).adc());
633  meStripTOBLayer1zm_[2]->Fill((*iter).strip());
634  }
635  if (tTopo->tobModule(id) == 4) {
636  meAdcTOBLayer1zm_[3]->Fill((*iter).adc());
637  meStripTOBLayer1zm_[3]->Fill((*iter).strip());
638  }
639  if (tTopo->tobModule(id) == 5) {
640  meAdcTOBLayer1zm_[4]->Fill((*iter).adc());
641  meStripTOBLayer1zm_[4]->Fill((*iter).strip());
642  }
643  if (tTopo->tobModule(id) == 6) {
644  meAdcTOBLayer1zm_[5]->Fill((*iter).adc());
645  meStripTOBLayer1zm_[5]->Fill((*iter).strip());
646  }
647  }
648  if (tTopo->tobLayer(id) == 2) {
649  if (tTopo->tobModule(id) == 1) {
650  meAdcTOBLayer2zm_[0]->Fill((*iter).adc());
651  meStripTOBLayer2zm_[0]->Fill((*iter).strip());
652  }
653  if (tTopo->tobModule(id) == 2) {
654  meAdcTOBLayer2zm_[1]->Fill((*iter).adc());
655  meStripTOBLayer2zm_[1]->Fill((*iter).strip());
656  }
657  if (tTopo->tobModule(id) == 3) {
658  meAdcTOBLayer2zm_[2]->Fill((*iter).adc());
659  meStripTOBLayer2zm_[2]->Fill((*iter).strip());
660  }
661  if (tTopo->tobModule(id) == 4) {
662  meAdcTOBLayer2zm_[3]->Fill((*iter).adc());
663  meStripTOBLayer2zm_[3]->Fill((*iter).strip());
664  }
665  if (tTopo->tobModule(id) == 5) {
666  meAdcTOBLayer2zm_[4]->Fill((*iter).adc());
667  meStripTOBLayer2zm_[4]->Fill((*iter).strip());
668  }
669  if (tTopo->tobModule(id) == 6) {
670  meAdcTOBLayer2zm_[5]->Fill((*iter).adc());
671  meStripTOBLayer2zm_[5]->Fill((*iter).strip());
672  }
673  }
674  if (tTopo->tobLayer(id) == 3) {
675  if (tTopo->tobModule(id) == 1) {
676  meAdcTOBLayer3zm_[0]->Fill((*iter).adc());
677  meStripTOBLayer3zm_[0]->Fill((*iter).strip());
678  }
679  if (tTopo->tobModule(id) == 2) {
680  meAdcTOBLayer3zm_[1]->Fill((*iter).adc());
681  meStripTOBLayer3zm_[1]->Fill((*iter).strip());
682  }
683  if (tTopo->tobModule(id) == 3) {
684  meAdcTOBLayer3zm_[2]->Fill((*iter).adc());
685  meStripTOBLayer3zm_[2]->Fill((*iter).strip());
686  }
687  if (tTopo->tobModule(id) == 4) {
688  meAdcTOBLayer3zm_[3]->Fill((*iter).adc());
689  meStripTOBLayer3zm_[3]->Fill((*iter).strip());
690  }
691  if (tTopo->tobModule(id) == 5) {
692  meAdcTOBLayer3zm_[4]->Fill((*iter).adc());
693  meStripTOBLayer3zm_[4]->Fill((*iter).strip());
694  }
695  if (tTopo->tobModule(id) == 6) {
696  meAdcTOBLayer3zm_[5]->Fill((*iter).adc());
697  meStripTOBLayer3zm_[5]->Fill((*iter).strip());
698  }
699  }
700  if (tTopo->tobLayer(id) == 4) {
701  if (tTopo->tobModule(id) == 1) {
702  meAdcTOBLayer4zm_[0]->Fill((*iter).adc());
703  meStripTOBLayer4zm_[0]->Fill((*iter).strip());
704  }
705  if (tTopo->tobModule(id) == 2) {
706  meAdcTOBLayer4zm_[1]->Fill((*iter).adc());
707  meStripTOBLayer4zm_[1]->Fill((*iter).strip());
708  }
709  if (tTopo->tobModule(id) == 3) {
710  meAdcTOBLayer4zm_[2]->Fill((*iter).adc());
711  meStripTOBLayer4zm_[2]->Fill((*iter).strip());
712  }
713  if (tTopo->tobModule(id) == 4) {
714  meAdcTOBLayer4zm_[3]->Fill((*iter).adc());
715  meStripTOBLayer4zm_[3]->Fill((*iter).strip());
716  }
717  if (tTopo->tobModule(id) == 5) {
718  meAdcTOBLayer4zm_[4]->Fill((*iter).adc());
719  meStripTOBLayer4zm_[4]->Fill((*iter).strip());
720  }
721  if (tTopo->tobModule(id) == 6) {
722  meAdcTOBLayer4zm_[5]->Fill((*iter).adc());
723  meStripTOBLayer4zm_[5]->Fill((*iter).strip());
724  }
725  }
726  if (tTopo->tobLayer(id) == 5) {
727  if (tTopo->tobModule(id) == 1) {
728  meAdcTOBLayer5zm_[0]->Fill((*iter).adc());
729  meStripTOBLayer5zm_[0]->Fill((*iter).strip());
730  }
731  if (tTopo->tobModule(id) == 2) {
732  meAdcTOBLayer5zm_[1]->Fill((*iter).adc());
733  meStripTOBLayer5zm_[1]->Fill((*iter).strip());
734  }
735  if (tTopo->tobModule(id) == 3) {
736  meAdcTOBLayer5zm_[2]->Fill((*iter).adc());
737  meStripTOBLayer5zm_[2]->Fill((*iter).strip());
738  }
739  if (tTopo->tobModule(id) == 4) {
740  meAdcTOBLayer5zm_[3]->Fill((*iter).adc());
741  meStripTOBLayer5zm_[3]->Fill((*iter).strip());
742  }
743  if (tTopo->tobModule(id) == 5) {
744  meAdcTOBLayer5zm_[4]->Fill((*iter).adc());
745  meStripTOBLayer5zm_[4]->Fill((*iter).strip());
746  }
747  if (tTopo->tobModule(id) == 6) {
748  meAdcTOBLayer5zm_[5]->Fill((*iter).adc());
749  meStripTOBLayer5zm_[5]->Fill((*iter).strip());
750  }
751  }
752  if (tTopo->tobLayer(id) == 6) {
753  if (tTopo->tobModule(id) == 1) {
754  meAdcTOBLayer6zm_[0]->Fill((*iter).adc());
755  meStripTOBLayer6zm_[0]->Fill((*iter).strip());
756  }
757  if (tTopo->tobModule(id) == 2) {
758  meAdcTOBLayer6zm_[1]->Fill((*iter).adc());
759  meStripTOBLayer6zm_[1]->Fill((*iter).strip());
760  }
761  if (tTopo->tobModule(id) == 3) {
762  meAdcTOBLayer6zm_[2]->Fill((*iter).adc());
763  meStripTOBLayer6zm_[2]->Fill((*iter).strip());
764  }
765  if (tTopo->tobModule(id) == 4) {
766  meAdcTOBLayer6zm_[3]->Fill((*iter).adc());
767  meStripTOBLayer6zm_[3]->Fill((*iter).strip());
768  }
769  if (tTopo->tobModule(id) == 5) {
770  meAdcTOBLayer6zm_[4]->Fill((*iter).adc());
771  meStripTOBLayer6zm_[4]->Fill((*iter).strip());
772  }
773  if (tTopo->tobModule(id) == 6) {
774  meAdcTOBLayer6zm_[5]->Fill((*iter).adc());
775  meStripTOBLayer6zm_[5]->Fill((*iter).strip());
776  }
777  }
778 
779  } else {
780  ++ndigilayertobzp[tTopo->tobLayer(id) - 1];
781  if (tTopo->tobLayer(id) == 1) {
782  if (tTopo->tobModule(id) == 1) {
783  meAdcTOBLayer1zp_[0]->Fill((*iter).adc());
784  meStripTOBLayer1zp_[0]->Fill((*iter).strip());
785  }
786  if (tTopo->tobModule(id) == 2) {
787  meAdcTOBLayer1zp_[1]->Fill((*iter).adc());
788  meStripTOBLayer1zp_[1]->Fill((*iter).strip());
789  }
790  if (tTopo->tobModule(id) == 3) {
791  meAdcTOBLayer1zp_[2]->Fill((*iter).adc());
792  meStripTOBLayer1zp_[2]->Fill((*iter).strip());
793  }
794  if (tTopo->tobModule(id) == 4) {
795  meAdcTOBLayer1zp_[3]->Fill((*iter).adc());
796  meStripTOBLayer1zp_[3]->Fill((*iter).strip());
797  }
798  if (tTopo->tobModule(id) == 5) {
799  meAdcTOBLayer1zp_[4]->Fill((*iter).adc());
800  meStripTOBLayer1zp_[4]->Fill((*iter).strip());
801  }
802  if (tTopo->tobModule(id) == 6) {
803  meAdcTOBLayer1zp_[5]->Fill((*iter).adc());
804  meStripTOBLayer1zp_[5]->Fill((*iter).strip());
805  }
806  }
807  if (tTopo->tobLayer(id) == 2) {
808  if (tTopo->tobModule(id) == 1) {
809  meAdcTOBLayer2zp_[0]->Fill((*iter).adc());
810  meStripTOBLayer2zp_[0]->Fill((*iter).strip());
811  }
812  if (tTopo->tobModule(id) == 2) {
813  meAdcTOBLayer2zp_[1]->Fill((*iter).adc());
814  meStripTOBLayer2zp_[1]->Fill((*iter).strip());
815  }
816  if (tTopo->tobModule(id) == 3) {
817  meAdcTOBLayer2zp_[2]->Fill((*iter).adc());
818  meStripTOBLayer2zp_[2]->Fill((*iter).strip());
819  }
820  if (tTopo->tobModule(id) == 4) {
821  meAdcTOBLayer2zp_[3]->Fill((*iter).adc());
822  meStripTOBLayer2zp_[3]->Fill((*iter).strip());
823  }
824  if (tTopo->tobModule(id) == 5) {
825  meAdcTOBLayer2zp_[4]->Fill((*iter).adc());
826  meStripTOBLayer2zp_[4]->Fill((*iter).strip());
827  }
828  if (tTopo->tobModule(id) == 6) {
829  meAdcTOBLayer2zp_[5]->Fill((*iter).adc());
830  meStripTOBLayer2zp_[5]->Fill((*iter).strip());
831  }
832  }
833  if (tTopo->tobLayer(id) == 3) {
834  if (tTopo->tobModule(id) == 1) {
835  meAdcTOBLayer3zp_[0]->Fill((*iter).adc());
836  meStripTOBLayer3zp_[0]->Fill((*iter).strip());
837  }
838  if (tTopo->tobModule(id) == 2) {
839  meAdcTOBLayer3zp_[1]->Fill((*iter).adc());
840  meStripTOBLayer3zp_[1]->Fill((*iter).strip());
841  }
842  if (tTopo->tobModule(id) == 3) {
843  meAdcTOBLayer3zp_[2]->Fill((*iter).adc());
844  meStripTOBLayer3zp_[2]->Fill((*iter).strip());
845  }
846  if (tTopo->tobModule(id) == 4) {
847  meAdcTOBLayer3zp_[3]->Fill((*iter).adc());
848  meStripTOBLayer3zp_[3]->Fill((*iter).strip());
849  }
850  if (tTopo->tobModule(id) == 5) {
851  meAdcTOBLayer3zp_[4]->Fill((*iter).adc());
852  meStripTOBLayer3zp_[4]->Fill((*iter).strip());
853  }
854  if (tTopo->tobModule(id) == 6) {
855  meAdcTOBLayer3zp_[5]->Fill((*iter).adc());
856  meStripTOBLayer3zp_[5]->Fill((*iter).strip());
857  }
858  }
859  if (tTopo->tobLayer(id) == 4) {
860  if (tTopo->tobModule(id) == 1) {
861  meAdcTOBLayer4zp_[0]->Fill((*iter).adc());
862  meStripTOBLayer4zp_[0]->Fill((*iter).strip());
863  }
864  if (tTopo->tobModule(id) == 2) {
865  meAdcTOBLayer4zp_[1]->Fill((*iter).adc());
866  meStripTOBLayer4zp_[1]->Fill((*iter).strip());
867  }
868  if (tTopo->tobModule(id) == 3) {
869  meAdcTOBLayer4zp_[2]->Fill((*iter).adc());
870  meStripTOBLayer4zp_[2]->Fill((*iter).strip());
871  }
872  if (tTopo->tobModule(id) == 4) {
873  meAdcTOBLayer4zp_[3]->Fill((*iter).adc());
874  meStripTOBLayer4zp_[3]->Fill((*iter).strip());
875  }
876  if (tTopo->tobModule(id) == 5) {
877  meAdcTOBLayer4zp_[4]->Fill((*iter).adc());
878  meStripTOBLayer4zp_[4]->Fill((*iter).strip());
879  }
880  if (tTopo->tobModule(id) == 6) {
881  meAdcTOBLayer4zp_[5]->Fill((*iter).adc());
882  meStripTOBLayer4zp_[5]->Fill((*iter).strip());
883  }
884  }
885  if (tTopo->tobLayer(id) == 5) {
886  if (tTopo->tobModule(id) == 1) {
887  meAdcTOBLayer5zp_[0]->Fill((*iter).adc());
888  meStripTOBLayer5zp_[0]->Fill((*iter).strip());
889  }
890  if (tTopo->tobModule(id) == 2) {
891  meAdcTOBLayer5zp_[1]->Fill((*iter).adc());
892  meStripTOBLayer5zp_[1]->Fill((*iter).strip());
893  }
894  if (tTopo->tobModule(id) == 3) {
895  meAdcTOBLayer5zp_[2]->Fill((*iter).adc());
896  meStripTOBLayer5zp_[2]->Fill((*iter).strip());
897  }
898  if (tTopo->tobModule(id) == 4) {
899  meAdcTOBLayer5zp_[3]->Fill((*iter).adc());
900  meStripTOBLayer5zp_[3]->Fill((*iter).strip());
901  }
902  if (tTopo->tobModule(id) == 5) {
903  meAdcTOBLayer5zp_[4]->Fill((*iter).adc());
904  meStripTOBLayer5zp_[4]->Fill((*iter).strip());
905  }
906  if (tTopo->tobModule(id) == 6) {
907  meAdcTOBLayer5zp_[5]->Fill((*iter).adc());
908  meStripTOBLayer5zp_[5]->Fill((*iter).strip());
909  }
910  }
911  if (tTopo->tobLayer(id) == 6) {
912  if (tTopo->tobModule(id) == 1) {
913  meAdcTOBLayer6zp_[0]->Fill((*iter).adc());
914  meStripTOBLayer6zp_[0]->Fill((*iter).strip());
915  }
916  if (tTopo->tobModule(id) == 2) {
917  meAdcTOBLayer6zp_[1]->Fill((*iter).adc());
918  meStripTOBLayer6zp_[1]->Fill((*iter).strip());
919  }
920  if (tTopo->tobModule(id) == 3) {
921  meAdcTOBLayer6zp_[2]->Fill((*iter).adc());
922  meStripTOBLayer6zp_[2]->Fill((*iter).strip());
923  }
924  if (tTopo->tobModule(id) == 4) {
925  meAdcTOBLayer6zp_[3]->Fill((*iter).adc());
926  meStripTOBLayer6zp_[3]->Fill((*iter).strip());
927  }
928  if (tTopo->tobModule(id) == 5) {
929  meAdcTOBLayer6zp_[4]->Fill((*iter).adc());
930  meStripTOBLayer6zp_[4]->Fill((*iter).strip());
931  }
932  if (tTopo->tobModule(id) == 6) {
933  meAdcTOBLayer6zp_[5]->Fill((*iter).adc());
934  meStripTOBLayer6zp_[5]->Fill((*iter).strip());
935  }
936  }
937  }
938  }
939  }
940 
941  if (detId.subdetId() == StripSubdetector::TID) {
942  for (iter = begin; iter != end; iter++) {
943  if (tTopo->tidSide(id) == 1) {
944  ++ndigiwheeltidzm[tTopo->tidWheel(id) - 1];
945  if (tTopo->tidWheel(id) == 1) {
946  if (tTopo->tidRing(id) == 1) {
947  meAdcTIDWheel1zm_[0]->Fill((*iter).adc());
948  meStripTIDWheel1zm_[0]->Fill((*iter).strip());
949  }
950  if (tTopo->tidRing(id) == 2) {
951  meAdcTIDWheel1zm_[1]->Fill((*iter).adc());
952  meStripTIDWheel1zm_[1]->Fill((*iter).strip());
953  }
954  if (tTopo->tidRing(id) == 3) {
955  meAdcTIDWheel1zm_[2]->Fill((*iter).adc());
956  meStripTIDWheel1zm_[2]->Fill((*iter).strip());
957  }
958  }
959  if (tTopo->tidWheel(id) == 2) {
960  if (tTopo->tidRing(id) == 1) {
961  meAdcTIDWheel2zm_[0]->Fill((*iter).adc());
962  meStripTIDWheel2zm_[0]->Fill((*iter).strip());
963  }
964  if (tTopo->tidRing(id) == 2) {
965  meAdcTIDWheel2zm_[1]->Fill((*iter).adc());
966  meStripTIDWheel2zm_[1]->Fill((*iter).strip());
967  }
968  if (tTopo->tidRing(id) == 3) {
969  meAdcTIDWheel2zm_[2]->Fill((*iter).adc());
970  meStripTIDWheel2zm_[2]->Fill((*iter).strip());
971  }
972  }
973  if (tTopo->tidWheel(id) == 3) {
974  if (tTopo->tidRing(id) == 1) {
975  meAdcTIDWheel3zm_[0]->Fill((*iter).adc());
976  meStripTIDWheel3zm_[0]->Fill((*iter).strip());
977  }
978  if (tTopo->tidRing(id) == 2) {
979  meAdcTIDWheel3zm_[1]->Fill((*iter).adc());
980  meStripTIDWheel3zm_[1]->Fill((*iter).strip());
981  }
982  if (tTopo->tidRing(id) == 3) {
983  meAdcTIDWheel3zm_[2]->Fill((*iter).adc());
984  meStripTIDWheel3zm_[2]->Fill((*iter).strip());
985  }
986  }
987 
988  } else {
989  ++ndigiwheeltidzp[tTopo->tidWheel(id) - 1];
990  if (tTopo->tidWheel(id) == 1) {
991  if (tTopo->tidRing(id) == 1) {
992  meAdcTIDWheel1zp_[0]->Fill((*iter).adc());
993  meStripTIDWheel1zp_[0]->Fill((*iter).strip());
994  }
995  if (tTopo->tidRing(id) == 2) {
996  meAdcTIDWheel1zp_[1]->Fill((*iter).adc());
997  meStripTIDWheel1zp_[1]->Fill((*iter).strip());
998  }
999  if (tTopo->tidRing(id) == 3) {
1000  meAdcTIDWheel1zp_[2]->Fill((*iter).adc());
1001  meStripTIDWheel1zp_[2]->Fill((*iter).strip());
1002  }
1003  }
1004  if (tTopo->tidWheel(id) == 2) {
1005  if (tTopo->tidRing(id) == 1) {
1006  meAdcTIDWheel2zp_[0]->Fill((*iter).adc());
1007  meStripTIDWheel2zp_[0]->Fill((*iter).strip());
1008  }
1009  if (tTopo->tidRing(id) == 2) {
1010  meAdcTIDWheel2zp_[1]->Fill((*iter).adc());
1011  meStripTIDWheel2zp_[1]->Fill((*iter).strip());
1012  }
1013  if (tTopo->tidRing(id) == 3) {
1014  meAdcTIDWheel2zp_[2]->Fill((*iter).adc());
1015  meStripTIDWheel2zp_[2]->Fill((*iter).strip());
1016  }
1017  }
1018  if (tTopo->tidWheel(id) == 3) {
1019  if (tTopo->tidRing(id) == 1) {
1020  meAdcTIDWheel3zp_[0]->Fill((*iter).adc());
1021  meStripTIDWheel3zp_[0]->Fill((*iter).strip());
1022  }
1023  if (tTopo->tidRing(id) == 2) {
1024  meAdcTIDWheel3zp_[1]->Fill((*iter).adc());
1025  meStripTIDWheel3zp_[1]->Fill((*iter).strip());
1026  }
1027  if (tTopo->tidRing(id) == 3) {
1028  meAdcTIDWheel3zp_[2]->Fill((*iter).adc());
1029  meStripTIDWheel3zp_[2]->Fill((*iter).strip());
1030  }
1031  }
1032  }
1033  }
1034  }
1035  if (detId.subdetId() == StripSubdetector::TEC) {
1036  for (iter = begin; iter != end; iter++) {
1037  if (tTopo->tecSide(id) == 1) {
1038  ++ndigiwheelteczm[tTopo->tecWheel(id) - 1];
1039  if (tTopo->tecWheel(id) == 1) {
1040  if (tTopo->tecRing(id) == 1) {
1041  meAdcTECWheel1zm_[0]->Fill((*iter).adc());
1042  meStripTECWheel1zm_[0]->Fill((*iter).strip());
1043  }
1044  if (tTopo->tecRing(id) == 2) {
1045  meAdcTECWheel1zm_[1]->Fill((*iter).adc());
1046  meStripTECWheel1zm_[1]->Fill((*iter).strip());
1047  }
1048  if (tTopo->tecRing(id) == 3) {
1049  meAdcTECWheel1zm_[2]->Fill((*iter).adc());
1050  meStripTECWheel1zm_[2]->Fill((*iter).strip());
1051  }
1052  if (tTopo->tecRing(id) == 4) {
1053  meAdcTECWheel1zm_[3]->Fill((*iter).adc());
1054  meStripTECWheel1zm_[3]->Fill((*iter).strip());
1055  }
1056  if (tTopo->tecRing(id) == 5) {
1057  meAdcTECWheel1zm_[4]->Fill((*iter).adc());
1058  meStripTECWheel1zm_[4]->Fill((*iter).strip());
1059  }
1060  if (tTopo->tecRing(id) == 6) {
1061  meAdcTECWheel1zm_[5]->Fill((*iter).adc());
1062  meStripTECWheel1zm_[5]->Fill((*iter).strip());
1063  }
1064  if (tTopo->tecRing(id) == 7) {
1065  meAdcTECWheel1zm_[6]->Fill((*iter).adc());
1066  meStripTECWheel1zm_[6]->Fill((*iter).strip());
1067  }
1068  }
1069  if (tTopo->tecWheel(id) == 2) {
1070  if (tTopo->tecRing(id) == 1) {
1071  meAdcTECWheel2zm_[0]->Fill((*iter).adc());
1072  meStripTECWheel2zm_[0]->Fill((*iter).strip());
1073  }
1074  if (tTopo->tecRing(id) == 2) {
1075  meAdcTECWheel2zm_[1]->Fill((*iter).adc());
1076  meStripTECWheel2zm_[1]->Fill((*iter).strip());
1077  }
1078  if (tTopo->tecRing(id) == 3) {
1079  meAdcTECWheel2zm_[2]->Fill((*iter).adc());
1080  meStripTECWheel2zm_[2]->Fill((*iter).strip());
1081  }
1082  if (tTopo->tecRing(id) == 4) {
1083  meAdcTECWheel2zm_[3]->Fill((*iter).adc());
1084  meStripTECWheel2zm_[3]->Fill((*iter).strip());
1085  }
1086  if (tTopo->tecRing(id) == 5) {
1087  meAdcTECWheel2zm_[4]->Fill((*iter).adc());
1088  meStripTECWheel2zm_[4]->Fill((*iter).strip());
1089  }
1090  if (tTopo->tecRing(id) == 6) {
1091  meAdcTECWheel2zm_[5]->Fill((*iter).adc());
1092  meStripTECWheel2zm_[5]->Fill((*iter).strip());
1093  }
1094  if (tTopo->tecRing(id) == 7) {
1095  meAdcTECWheel2zm_[6]->Fill((*iter).adc());
1096  meStripTECWheel2zm_[6]->Fill((*iter).strip());
1097  }
1098  }
1099  if (tTopo->tecWheel(id) == 3) {
1100  if (tTopo->tecRing(id) == 1) {
1101  meAdcTECWheel3zm_[0]->Fill((*iter).adc());
1102  meStripTECWheel3zm_[0]->Fill((*iter).strip());
1103  }
1104  if (tTopo->tecRing(id) == 2) {
1105  meAdcTECWheel3zm_[1]->Fill((*iter).adc());
1106  meStripTECWheel3zm_[1]->Fill((*iter).strip());
1107  }
1108  if (tTopo->tecRing(id) == 3) {
1109  meAdcTECWheel3zm_[2]->Fill((*iter).adc());
1110  meStripTECWheel3zm_[2]->Fill((*iter).strip());
1111  }
1112  if (tTopo->tecRing(id) == 4) {
1113  meAdcTECWheel3zm_[3]->Fill((*iter).adc());
1114  meStripTECWheel3zm_[3]->Fill((*iter).strip());
1115  }
1116  if (tTopo->tecRing(id) == 5) {
1117  meAdcTECWheel3zm_[4]->Fill((*iter).adc());
1118  meStripTECWheel3zm_[4]->Fill((*iter).strip());
1119  }
1120  if (tTopo->tecRing(id) == 6) {
1121  meAdcTECWheel3zm_[5]->Fill((*iter).adc());
1122  meStripTECWheel3zm_[5]->Fill((*iter).strip());
1123  }
1124  if (tTopo->tecRing(id) == 7) {
1125  meAdcTECWheel3zm_[6]->Fill((*iter).adc());
1126  meStripTECWheel3zm_[6]->Fill((*iter).strip());
1127  }
1128  }
1129  if (tTopo->tecWheel(id) == 4) {
1130  if (tTopo->tecRing(id) == 2) {
1131  meAdcTECWheel4zm_[0]->Fill((*iter).adc());
1132  meStripTECWheel4zm_[0]->Fill((*iter).strip());
1133  }
1134  if (tTopo->tecRing(id) == 3) {
1135  meAdcTECWheel4zm_[1]->Fill((*iter).adc());
1136  meStripTECWheel4zm_[1]->Fill((*iter).strip());
1137  }
1138  if (tTopo->tecRing(id) == 4) {
1139  meAdcTECWheel4zm_[2]->Fill((*iter).adc());
1140  meStripTECWheel4zm_[2]->Fill((*iter).strip());
1141  }
1142  if (tTopo->tecRing(id) == 5) {
1143  meAdcTECWheel4zm_[3]->Fill((*iter).adc());
1144  meStripTECWheel4zm_[3]->Fill((*iter).strip());
1145  }
1146  if (tTopo->tecRing(id) == 6) {
1147  meAdcTECWheel4zm_[4]->Fill((*iter).adc());
1148  meStripTECWheel4zm_[4]->Fill((*iter).strip());
1149  }
1150  if (tTopo->tecRing(id) == 7) {
1151  meAdcTECWheel4zm_[5]->Fill((*iter).adc());
1152  meStripTECWheel4zm_[5]->Fill((*iter).strip());
1153  }
1154  }
1155  if (tTopo->tecWheel(id) == 5) {
1156  if (tTopo->tecRing(id) == 2) {
1157  meAdcTECWheel5zm_[0]->Fill((*iter).adc());
1158  meStripTECWheel5zm_[0]->Fill((*iter).strip());
1159  }
1160  if (tTopo->tecRing(id) == 3) {
1161  meAdcTECWheel5zm_[1]->Fill((*iter).adc());
1162  meStripTECWheel5zm_[1]->Fill((*iter).strip());
1163  }
1164  if (tTopo->tecRing(id) == 4) {
1165  meAdcTECWheel5zm_[2]->Fill((*iter).adc());
1166  meStripTECWheel5zm_[2]->Fill((*iter).strip());
1167  }
1168  if (tTopo->tecRing(id) == 5) {
1169  meAdcTECWheel5zm_[3]->Fill((*iter).adc());
1170  meStripTECWheel5zm_[3]->Fill((*iter).strip());
1171  }
1172  if (tTopo->tecRing(id) == 6) {
1173  meAdcTECWheel5zm_[4]->Fill((*iter).adc());
1174  meStripTECWheel5zm_[4]->Fill((*iter).strip());
1175  }
1176  if (tTopo->tecRing(id) == 7) {
1177  meAdcTECWheel5zm_[5]->Fill((*iter).adc());
1178  meStripTECWheel5zm_[5]->Fill((*iter).strip());
1179  }
1180  }
1181  if (tTopo->tecWheel(id) == 6) {
1182  if (tTopo->tecRing(id) == 2) {
1183  meAdcTECWheel6zm_[0]->Fill((*iter).adc());
1184  meStripTECWheel6zm_[0]->Fill((*iter).strip());
1185  }
1186  if (tTopo->tecRing(id) == 3) {
1187  meAdcTECWheel6zm_[1]->Fill((*iter).adc());
1188  meStripTECWheel6zm_[1]->Fill((*iter).strip());
1189  }
1190  if (tTopo->tecRing(id) == 4) {
1191  meAdcTECWheel6zm_[2]->Fill((*iter).adc());
1192  meStripTECWheel6zm_[2]->Fill((*iter).strip());
1193  }
1194  if (tTopo->tecRing(id) == 5) {
1195  meAdcTECWheel6zm_[3]->Fill((*iter).adc());
1196  meStripTECWheel6zm_[3]->Fill((*iter).strip());
1197  }
1198  if (tTopo->tecRing(id) == 6) {
1199  meAdcTECWheel6zm_[4]->Fill((*iter).adc());
1200  meStripTECWheel6zm_[4]->Fill((*iter).strip());
1201  }
1202  if (tTopo->tecRing(id) == 7) {
1203  meAdcTECWheel6zm_[5]->Fill((*iter).adc());
1204  meStripTECWheel6zm_[5]->Fill((*iter).strip());
1205  }
1206  }
1207  if (tTopo->tecWheel(id) == 7) {
1208  if (tTopo->tecRing(id) == 3) {
1209  meAdcTECWheel7zm_[0]->Fill((*iter).adc());
1210  meStripTECWheel7zm_[0]->Fill((*iter).strip());
1211  }
1212  if (tTopo->tecRing(id) == 4) {
1213  meAdcTECWheel7zm_[1]->Fill((*iter).adc());
1214  meStripTECWheel7zm_[1]->Fill((*iter).strip());
1215  }
1216  if (tTopo->tecRing(id) == 5) {
1217  meAdcTECWheel7zm_[2]->Fill((*iter).adc());
1218  meStripTECWheel7zm_[2]->Fill((*iter).strip());
1219  }
1220  if (tTopo->tecRing(id) == 6) {
1221  meAdcTECWheel7zm_[3]->Fill((*iter).adc());
1222  meStripTECWheel7zm_[3]->Fill((*iter).strip());
1223  }
1224  if (tTopo->tecRing(id) == 7) {
1225  meAdcTECWheel7zm_[4]->Fill((*iter).adc());
1226  meStripTECWheel7zm_[4]->Fill((*iter).strip());
1227  }
1228  }
1229  if (tTopo->tecWheel(id) == 8) {
1230  if (tTopo->tecRing(id) == 3) {
1231  meAdcTECWheel8zm_[0]->Fill((*iter).adc());
1232  meStripTECWheel8zm_[0]->Fill((*iter).strip());
1233  }
1234  if (tTopo->tecRing(id) == 4) {
1235  meAdcTECWheel8zm_[1]->Fill((*iter).adc());
1236  meStripTECWheel8zm_[1]->Fill((*iter).strip());
1237  }
1238  if (tTopo->tecRing(id) == 5) {
1239  meAdcTECWheel8zm_[2]->Fill((*iter).adc());
1240  meStripTECWheel8zm_[2]->Fill((*iter).strip());
1241  }
1242  if (tTopo->tecRing(id) == 6) {
1243  meAdcTECWheel8zm_[3]->Fill((*iter).adc());
1244  meStripTECWheel8zm_[3]->Fill((*iter).strip());
1245  }
1246  if (tTopo->tecRing(id) == 7) {
1247  meAdcTECWheel8zm_[4]->Fill((*iter).adc());
1248  meStripTECWheel8zm_[4]->Fill((*iter).strip());
1249  }
1250  }
1251  if (tTopo->tecWheel(id) == 9) {
1252  if (tTopo->tecRing(id) == 4) {
1253  meAdcTECWheel9zm_[0]->Fill((*iter).adc());
1254  meStripTECWheel9zm_[0]->Fill((*iter).strip());
1255  }
1256  if (tTopo->tecRing(id) == 5) {
1257  meAdcTECWheel9zm_[1]->Fill((*iter).adc());
1258  meStripTECWheel9zm_[1]->Fill((*iter).strip());
1259  }
1260  if (tTopo->tecRing(id) == 6) {
1261  meAdcTECWheel9zm_[2]->Fill((*iter).adc());
1262  meStripTECWheel9zm_[2]->Fill((*iter).strip());
1263  }
1264  if (tTopo->tecRing(id) == 7) {
1265  meAdcTECWheel9zm_[3]->Fill((*iter).adc());
1266  meStripTECWheel9zm_[3]->Fill((*iter).strip());
1267  }
1268  }
1269  } else {
1270  ++ndigiwheelteczp[tTopo->tecWheel(id) - 1];
1271  if (tTopo->tecWheel(id) == 1) {
1272  if (tTopo->tecRing(id) == 1) {
1273  meAdcTECWheel1zp_[0]->Fill((*iter).adc());
1274  meStripTECWheel1zp_[0]->Fill((*iter).strip());
1275  }
1276  if (tTopo->tecRing(id) == 2) {
1277  meAdcTECWheel1zp_[1]->Fill((*iter).adc());
1278  meStripTECWheel1zp_[1]->Fill((*iter).strip());
1279  }
1280  if (tTopo->tecRing(id) == 3) {
1281  meAdcTECWheel1zp_[2]->Fill((*iter).adc());
1282  meStripTECWheel1zp_[2]->Fill((*iter).strip());
1283  }
1284  if (tTopo->tecRing(id) == 4) {
1285  meAdcTECWheel1zp_[3]->Fill((*iter).adc());
1286  meStripTECWheel1zp_[3]->Fill((*iter).strip());
1287  }
1288  if (tTopo->tecRing(id) == 5) {
1289  meAdcTECWheel1zp_[4]->Fill((*iter).adc());
1290  meStripTECWheel1zp_[4]->Fill((*iter).strip());
1291  }
1292  if (tTopo->tecRing(id) == 6) {
1293  meAdcTECWheel1zp_[5]->Fill((*iter).adc());
1294  meStripTECWheel1zp_[5]->Fill((*iter).strip());
1295  }
1296  if (tTopo->tecRing(id) == 7) {
1297  meAdcTECWheel1zp_[6]->Fill((*iter).adc());
1298  meStripTECWheel1zp_[6]->Fill((*iter).strip());
1299  }
1300  }
1301  if (tTopo->tecWheel(id) == 2) {
1302  if (tTopo->tecRing(id) == 1) {
1303  meAdcTECWheel2zp_[0]->Fill((*iter).adc());
1304  meStripTECWheel2zp_[0]->Fill((*iter).strip());
1305  }
1306  if (tTopo->tecRing(id) == 2) {
1307  meAdcTECWheel2zp_[1]->Fill((*iter).adc());
1308  meStripTECWheel2zp_[1]->Fill((*iter).strip());
1309  }
1310  if (tTopo->tecRing(id) == 3) {
1311  meAdcTECWheel2zp_[2]->Fill((*iter).adc());
1312  meStripTECWheel2zp_[2]->Fill((*iter).strip());
1313  }
1314  if (tTopo->tecRing(id) == 4) {
1315  meAdcTECWheel2zp_[3]->Fill((*iter).adc());
1316  meStripTECWheel2zp_[3]->Fill((*iter).strip());
1317  }
1318  if (tTopo->tecRing(id) == 5) {
1319  meAdcTECWheel2zp_[4]->Fill((*iter).adc());
1320  meStripTECWheel2zp_[4]->Fill((*iter).strip());
1321  }
1322  if (tTopo->tecRing(id) == 6) {
1323  meAdcTECWheel2zp_[5]->Fill((*iter).adc());
1324  meStripTECWheel2zp_[5]->Fill((*iter).strip());
1325  }
1326  if (tTopo->tecRing(id) == 7) {
1327  meAdcTECWheel2zp_[6]->Fill((*iter).adc());
1328  meStripTECWheel2zp_[6]->Fill((*iter).strip());
1329  }
1330  }
1331  if (tTopo->tecWheel(id) == 3) {
1332  if (tTopo->tecRing(id) == 1) {
1333  meAdcTECWheel3zp_[0]->Fill((*iter).adc());
1334  meStripTECWheel3zp_[0]->Fill((*iter).strip());
1335  }
1336  if (tTopo->tecRing(id) == 2) {
1337  meAdcTECWheel3zp_[1]->Fill((*iter).adc());
1338  meStripTECWheel3zp_[1]->Fill((*iter).strip());
1339  }
1340  if (tTopo->tecRing(id) == 3) {
1341  meAdcTECWheel3zp_[2]->Fill((*iter).adc());
1342  meStripTECWheel3zp_[2]->Fill((*iter).strip());
1343  }
1344  if (tTopo->tecRing(id) == 4) {
1345  meAdcTECWheel3zp_[3]->Fill((*iter).adc());
1346  meStripTECWheel3zp_[3]->Fill((*iter).strip());
1347  }
1348  if (tTopo->tecRing(id) == 5) {
1349  meAdcTECWheel3zp_[4]->Fill((*iter).adc());
1350  meStripTECWheel3zp_[4]->Fill((*iter).strip());
1351  }
1352  if (tTopo->tecRing(id) == 6) {
1353  meAdcTECWheel3zp_[5]->Fill((*iter).adc());
1354  meStripTECWheel3zp_[5]->Fill((*iter).strip());
1355  }
1356  if (tTopo->tecRing(id) == 7) {
1357  meAdcTECWheel3zp_[6]->Fill((*iter).adc());
1358  meStripTECWheel3zp_[6]->Fill((*iter).strip());
1359  }
1360  }
1361  if (tTopo->tecWheel(id) == 4) {
1362  if (tTopo->tecRing(id) == 2) {
1363  meAdcTECWheel4zp_[0]->Fill((*iter).adc());
1364  meStripTECWheel4zp_[0]->Fill((*iter).strip());
1365  }
1366  if (tTopo->tecRing(id) == 3) {
1367  meAdcTECWheel4zp_[1]->Fill((*iter).adc());
1368  meStripTECWheel4zp_[1]->Fill((*iter).strip());
1369  }
1370  if (tTopo->tecRing(id) == 4) {
1371  meAdcTECWheel4zp_[2]->Fill((*iter).adc());
1372  meStripTECWheel4zp_[2]->Fill((*iter).strip());
1373  }
1374  if (tTopo->tecRing(id) == 5) {
1375  meAdcTECWheel4zp_[3]->Fill((*iter).adc());
1376  meStripTECWheel4zp_[3]->Fill((*iter).strip());
1377  }
1378  if (tTopo->tecRing(id) == 6) {
1379  meAdcTECWheel4zp_[4]->Fill((*iter).adc());
1380  meStripTECWheel4zp_[4]->Fill((*iter).strip());
1381  }
1382  if (tTopo->tecRing(id) == 7) {
1383  meAdcTECWheel4zp_[5]->Fill((*iter).adc());
1384  meStripTECWheel4zp_[5]->Fill((*iter).strip());
1385  }
1386  }
1387  if (tTopo->tecWheel(id) == 5) {
1388  if (tTopo->tecRing(id) == 2) {
1389  meAdcTECWheel5zp_[0]->Fill((*iter).adc());
1390  meStripTECWheel5zp_[0]->Fill((*iter).strip());
1391  }
1392  if (tTopo->tecRing(id) == 3) {
1393  meAdcTECWheel5zp_[1]->Fill((*iter).adc());
1394  meStripTECWheel5zp_[1]->Fill((*iter).strip());
1395  }
1396  if (tTopo->tecRing(id) == 4) {
1397  meAdcTECWheel5zp_[2]->Fill((*iter).adc());
1398  meStripTECWheel5zp_[2]->Fill((*iter).strip());
1399  }
1400  if (tTopo->tecRing(id) == 5) {
1401  meAdcTECWheel5zp_[3]->Fill((*iter).adc());
1402  meStripTECWheel5zp_[3]->Fill((*iter).strip());
1403  }
1404  if (tTopo->tecRing(id) == 6) {
1405  meAdcTECWheel5zp_[4]->Fill((*iter).adc());
1406  meStripTECWheel5zp_[4]->Fill((*iter).strip());
1407  }
1408  if (tTopo->tecRing(id) == 7) {
1409  meAdcTECWheel5zp_[5]->Fill((*iter).adc());
1410  meStripTECWheel5zp_[5]->Fill((*iter).strip());
1411  }
1412  }
1413  if (tTopo->tecWheel(id) == 6) {
1414  if (tTopo->tecRing(id) == 2) {
1415  meAdcTECWheel6zp_[0]->Fill((*iter).adc());
1416  meStripTECWheel6zp_[0]->Fill((*iter).strip());
1417  }
1418  if (tTopo->tecRing(id) == 3) {
1419  meAdcTECWheel6zp_[1]->Fill((*iter).adc());
1420  meStripTECWheel6zp_[1]->Fill((*iter).strip());
1421  }
1422  if (tTopo->tecRing(id) == 4) {
1423  meAdcTECWheel6zp_[2]->Fill((*iter).adc());
1424  meStripTECWheel6zp_[2]->Fill((*iter).strip());
1425  }
1426  if (tTopo->tecRing(id) == 5) {
1427  meAdcTECWheel6zp_[3]->Fill((*iter).adc());
1428  meStripTECWheel6zp_[3]->Fill((*iter).strip());
1429  }
1430  if (tTopo->tecRing(id) == 6) {
1431  meAdcTECWheel6zp_[4]->Fill((*iter).adc());
1432  meStripTECWheel6zp_[4]->Fill((*iter).strip());
1433  }
1434  if (tTopo->tecRing(id) == 7) {
1435  meAdcTECWheel6zp_[5]->Fill((*iter).adc());
1436  meStripTECWheel6zp_[5]->Fill((*iter).strip());
1437  }
1438  }
1439  if (tTopo->tecWheel(id) == 7) {
1440  if (tTopo->tecRing(id) == 3) {
1441  meAdcTECWheel7zp_[0]->Fill((*iter).adc());
1442  meStripTECWheel7zp_[0]->Fill((*iter).strip());
1443  }
1444  if (tTopo->tecRing(id) == 4) {
1445  meAdcTECWheel7zp_[1]->Fill((*iter).adc());
1446  meStripTECWheel7zp_[1]->Fill((*iter).strip());
1447  }
1448  if (tTopo->tecRing(id) == 5) {
1449  meAdcTECWheel7zp_[2]->Fill((*iter).adc());
1450  meStripTECWheel7zp_[2]->Fill((*iter).strip());
1451  }
1452  if (tTopo->tecRing(id) == 6) {
1453  meAdcTECWheel7zp_[3]->Fill((*iter).adc());
1454  meStripTECWheel7zp_[3]->Fill((*iter).strip());
1455  }
1456  if (tTopo->tecRing(id) == 7) {
1457  meAdcTECWheel7zp_[4]->Fill((*iter).adc());
1458  meStripTECWheel7zp_[4]->Fill((*iter).strip());
1459  }
1460  }
1461  if (tTopo->tecWheel(id) == 8) {
1462  if (tTopo->tecRing(id) == 3) {
1463  meAdcTECWheel8zp_[0]->Fill((*iter).adc());
1464  meStripTECWheel8zp_[0]->Fill((*iter).strip());
1465  }
1466  if (tTopo->tecRing(id) == 4) {
1467  meAdcTECWheel8zp_[1]->Fill((*iter).adc());
1468  meStripTECWheel8zp_[1]->Fill((*iter).strip());
1469  }
1470  if (tTopo->tecRing(id) == 5) {
1471  meAdcTECWheel8zp_[2]->Fill((*iter).adc());
1472  meStripTECWheel8zp_[2]->Fill((*iter).strip());
1473  }
1474  if (tTopo->tecRing(id) == 6) {
1475  meAdcTECWheel8zp_[3]->Fill((*iter).adc());
1476  meStripTECWheel8zp_[3]->Fill((*iter).strip());
1477  }
1478  if (tTopo->tecRing(id) == 7) {
1479  meAdcTECWheel8zp_[4]->Fill((*iter).adc());
1480  meStripTECWheel8zp_[4]->Fill((*iter).strip());
1481  }
1482  }
1483  if (tTopo->tecWheel(id) == 9) {
1484  if (tTopo->tecRing(id) == 4) {
1485  meAdcTECWheel9zp_[0]->Fill((*iter).adc());
1486  meStripTECWheel9zp_[0]->Fill((*iter).strip());
1487  }
1488  if (tTopo->tecRing(id) == 5) {
1489  meAdcTECWheel9zp_[1]->Fill((*iter).adc());
1490  meStripTECWheel9zp_[1]->Fill((*iter).strip());
1491  }
1492  if (tTopo->tecRing(id) == 6) {
1493  meAdcTECWheel9zp_[2]->Fill((*iter).adc());
1494  meStripTECWheel9zp_[2]->Fill((*iter).strip());
1495  }
1496  if (tTopo->tecRing(id) == 7) {
1497  meAdcTECWheel9zp_[3]->Fill((*iter).adc());
1498  meStripTECWheel9zp_[3]->Fill((*iter).strip());
1499  }
1500  }
1501  }
1502  }
1503  }
1504  }
1505 
1506  for (int i = 0; i < 4; i++) {
1507  meNDigiTIBLayerzm_[i]->Fill(ndigilayertibzm[i]);
1508  meNDigiTIBLayerzp_[i]->Fill(ndigilayertibzp[i]);
1509  }
1510 
1511  for (int i = 0; i < 6; i++) {
1512  meNDigiTOBLayerzm_[i]->Fill(ndigilayertobzm[i]);
1513  meNDigiTOBLayerzp_[i]->Fill(ndigilayertobzp[i]);
1514  }
1515 
1516  for (int i = 0; i < 3; i++) {
1517  meNDigiTIDWheelzm_[i]->Fill(ndigiwheeltidzm[i]);
1518  meNDigiTIDWheelzp_[i]->Fill(ndigiwheeltidzp[i]);
1519  }
1520 
1521  for (int i = 0; i < 9; i++) {
1522  meNDigiTECWheelzm_[i]->Fill(ndigiwheelteczm[i]);
1523  meNDigiTECWheelzp_[i]->Fill(ndigiwheelteczp[i]);
1524  }
1525 }
MonitorElement * meAdcTOBLayer5zp_[6]
MonitorElement * meAdcTIBLayer4Extzm_[3]
static constexpr auto TEC
MonitorElement * meStripTIBLayer4Intzm_[3]
MonitorElement * meAdcTIDWheel3zm_[3]
MonitorElement * meStripTECWheel1zm_[7]
const edm::EventSetup & c
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
MonitorElement * meAdcTECWheel4zm_[6]
MonitorElement * meAdcTOBLayer4zp_[6]
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
MonitorElement * meStripTIBLayer4Intzp_[3]
edm::EDGetTokenT< edm::DetSetVector< SiStripDigi > > edmDetSetVector_SiStripDigi_Token_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
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]
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]
MonitorElement * meStripTECWheel6zm_[6]
void Fill(long long x)
bool getData(T &iHolder) const
Definition: EventSetup.h:128
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]
~SiStripDigiValid() override
MonitorElement * meStripTOBLayer5zm_[6]
MonitorElement * meAdcTECWheel8zp_[5]
MonitorElement * meAdcTIDWheel1zp_[3]
MonitorElement * meAdcTIBLayer1Intzm_[3]
void bookHistograms(DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) override
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
MonitorElement * meStripTECWheel9zm_[4]
MonitorElement * meAdcTECWheel7zp_[5]
MonitorElement * meAdcTECWheel1zm_[7]
MonitorElement * meStripTIBLayer1Intzp_[3]
MonitorElement * meStripTIDWheel2zm_[3]
A Digi for the silicon strip detector, containing both strip and adc information, and suitable for st...
Definition: SiStripDigi.h:12
static constexpr auto TOB
MonitorElement * meStripTOBLayer3zm_[6]
MonitorElement * meStripTIDWheel1zm_[3]
MonitorElement * meAdcTIBLayer3Extzp_[3]
MonitorElement * meStripTIBLayer3Extzm_[3]
MonitorElement * meStripTECWheel7zm_[5]
MonitorElement * meStripTECWheel2zp_[7]
iterator end()
Return the off-the-end iterator.
Definition: DetSetVector.h:325
MonitorElement * meAdcTIBLayer3Intzm_[3]
MonitorElement * meNDigiTIBLayerzm_[4]
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:17
MonitorElement * meStripTECWheel8zm_[5]
MonitorElement * meStripTECWheel5zm_[6]
MonitorElement * meStripTIBLayer1Intzm_[3]
static constexpr auto TIB
MonitorElement * meAdcTIBLayer4Intzp_[3]
MonitorElement * meStripTIBLayer1Extzp_[3]
MonitorElement * meStripTIBLayer1Extzm_[3]
MonitorElement * meNDigiTIDWheelzp_[3]
MonitorElement * meAdcTOBLayer3zp_[6]
MonitorElement * meAdcTOBLayer6zm_[6]
SiStripDigiValid(const edm::ParameterSet &ps)
void analyze(const edm::Event &e, const edm::EventSetup &c) override
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]
MonitorElement * meStripTOBLayer2zp_[6]
MonitorElement * meStripTIDWheel2zp_[3]
string end
Definition: dataset.py:937
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:314
MonitorElement * meStripTOBLayer6zm_[6]
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > m_topoToken
MonitorElement * meNDigiTECWheelzp_[9]
MonitorElement * meAdcTIBLayer3Intzp_[3]
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
collection_type::const_iterator const_iterator
Definition: DetSet.h:31
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:102
MonitorElement * meAdcTECWheel6zm_[6]
MonitorElement * meStripTECWheel1zp_[7]
unsigned int tecWheel(const DetId &id) const
static constexpr auto TID
MonitorElement * meStripTIBLayer3Intzp_[3]
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
MonitorElement * meAdcTECWheel4zp_[6]
MonitorElement * meAdcTECWheel9zm_[4]
MonitorElement * meStripTIBLayer2Extzp_[3]
unsigned int tobLayer(const DetId &id) const
Definition: Run.h:45
unsigned int tecSide(const DetId &id) const
MonitorElement * meAdcTECWheel2zm_[7]
MonitorElement * meStripTIDWheel3zm_[3]