74 #ifndef EECLUSTERTASKEXTRAS_DQMOFFLINE
82 meSCEneSingleCrystal_ = 0;
84 for(
int i=0;
i!=2;++) {
85 meSCSeedMapOccSingleCrystal_[
i] = 0;
86 meSCSeedMapOccSC_[
i] = 0;
87 meSCSeedMapOccHighEne_[
i] = 0;
88 meSCSeedMapTimeSC_[
i] = 0;
89 for(
int j=0;
j!=5;++
j) {
90 meSCSeedMapOccTrg_[
i][
j] = 0;
91 meSCSeedMapOccTrgExcl_[
i][
j] = 0;
96 meSCSeedTimeVsAmp_ = 0;
99 for(
int i=0;
i!=18;++
i)
100 meSCSeedTimePerFed_[
i] = 0;
103 meSCSizCrystalVsEne_ = 0;
105 for(
int i=0;
i!=2;++
i) {
106 meSCSeedMapOcc_[
i] = 0;
107 meSCSeedMapOccHighEneSC_[
i] = 0;
108 for(
int j=0;
j!=5;++
j) {
109 meSCSeedMapTimeTrgSC_[
i][
j] = 0;
110 meSCSeedMapOccTrgSC_[
i][
j] = 0;
111 meSCSeedMapOccTrgExclSC_[
i][
j] = 0;
126 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras");
127 dqmStore_->rmdir(prefixME_ +
"/EEClusterTaskExtras");
136 if ( ! mergeRuns_ ) this->
reset();
145 #ifndef EECLUSTERTASKEXTRAS_DQMOFFLINE
146 if ( meSCSizCrystal_ ) meSCSizCrystal_->Reset();
147 if ( meSCSizBC_ ) meSCSizBC_->Reset();
149 if ( meSCSeedEne_ ) meSCSeedEne_->Reset();
150 if ( meSCEne2_ ) meSCEne2_->Reset();
151 if ( meSCEneLow_ ) meSCEneLow_->Reset();
152 if ( meSCEneHigh_ ) meSCEneHigh_->Reset();
153 if ( meSCEneSingleCrystal_ ) meSCEneSingleCrystal_->Reset();
155 for(
int i=0;
i!=2;++
i) {
156 if ( meSCSeedMapOccSingleCrystal_[
i] ) meSCSeedMapOccSingleCrystal_[
i]->Reset();
157 if ( meSCSeedMapOccSC_[i] ) meSCSeedMapOccSC_[
i]->Reset();
158 if ( meSCSeedMapOccHighEne_[i] ) meSCSeedMapOccHighEne_[
i]->Reset();
159 if ( meSCSeedMapTimeSC_[i] ) meSCSeedMapTimeSC_[
i]->Reset();
160 for(
int j=0;
j!=5;++
j) {
161 if ( meSCSeedMapOccTrg_[i][
j] ) meSCSeedMapOccTrg_[
i][
j]->Reset();
162 if ( meSCSeedMapOccTrgExcl_[i][j] ) meSCSeedMapOccTrgExcl_[
i][
j]->Reset();
166 if ( meSCSeedTime_ ) meSCSeedTime_->Reset();
167 if ( meSCSeedTimeVsAmp_ ) meSCSeedTimeVsAmp_->Reset();
168 if ( meSCSeedTimeEEM_ ) meSCSeedTimeEEM_->Reset();
169 if ( meSCSeedTimeEEP_ ) meSCSeedTimeEEP_->Reset();
171 for(
int i=0;
i!=18;++
i)
172 if ( meSCSeedTimePerFed_[
i] ) meSCSeedTimePerFed_[
i]->Reset();
175 if ( meSCSizCrystalVsEne_ ) meSCSizCrystalVsEne_->Reset();
177 for(
int i=0;i!=2;++
i) {
178 if ( meSCSeedMapOcc_[i] ) meSCSeedMapOcc_[
i]->Reset();
179 if ( meSCSeedMapOccHighEneSC_[i] ) meSCSeedMapOccHighEneSC_[
i]->Reset();
180 for(
int j=0;
j!=5; ++
j) {
181 if ( meSCSeedMapOccTrgSC_[i][
j] ) meSCSeedMapOccTrgSC_[
i][
j]->Reset();
182 if ( meSCSeedMapOccTrgExclSC_[i][j] ) meSCSeedMapOccTrgExclSC_[
i][
j]->Reset();
183 if ( meSCSeedMapTimeTrgSC_[i][j] ) meSCSeedMapTimeTrgSC_[
i][
j]->Reset();
196 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras");
198 #ifndef EECLUSTERTASKEXTRAS_DQMOFFLINE
200 histo =
"EECLTE SC size (crystal)";
201 meSCSizCrystal_ = dqmStore_->book1D(histo,histo,150,0,150);
202 meSCSizCrystal_->setAxisTitle(
"super cluster size (crystal)", 1);
204 histo =
"EECLTE SC size (basic clusters)";
205 meSCSizBC_ = dqmStore_->book1D(histo,histo,20,0,20);
206 meSCSizBC_->setAxisTitle(
"super cluster size (basic clusters)", 1);
208 histo =
"EECLTE SC energy";
209 meSCSeedEne_ = dqmStore_->book1D(histo,histo,200,0,1.8);
210 meSCSeedEne_->setAxisTitle(
"energy (GeV)", 1);
212 histo =
"EECLTE SC + highest neighbor energy";
213 meSCEne2_ = dqmStore_->book1D(histo,histo,200,0,1.8);
214 meSCEne2_->setAxisTitle(
"+ highest neighbor energy (GeV)", 1);
216 histo =
"EECLTE SC energy low scale";
217 meSCEneLow_ = dqmStore_->book1D(histo,histo,200,0,1.8);
218 meSCEneLow_->setAxisTitle(
"energy (GeV)", 1);
220 histo =
"EECLTE SC energy high scale";
221 meSCEneHigh_ = dqmStore_->book1D(histo,histo,200,0,200);
222 meSCEneHigh_->setAxisTitle(
"energy (GeV)", 1);
224 histo =
"EECLTE SC single crystal cluster energy (GeV)";
225 meSCEneSingleCrystal_ = dqmStore_->book1D(histo,histo,200,0,200);
226 meSCEneSingleCrystal_->setAxisTitle(
"energy (GeV)", 1);
228 histo =
"EECLTE SC seed occupancy map super crystal binned EE -";
229 meSCSeedMapOccSC_[0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
230 meSCSeedMapOccSC_[0]->setAxisTitle(
"jx", 1);
231 meSCSeedMapOccSC_[0]->setAxisTitle(
"jy", 2);
233 histo =
"EECLTE SC seed occupancy map EE - (high energy clusters)";
234 meSCSeedMapOccHighEne_[0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
235 meSCSeedMapOccHighEne_[0]->setAxisTitle(
"jx", 1);
236 meSCSeedMapOccHighEne_[0]->setAxisTitle(
"jy", 2);
238 histo =
"EECLTE SC single crystal cluster seed occupancy map EE -";
239 meSCSeedMapOccSingleCrystal_[0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
240 meSCSeedMapOccSingleCrystal_[0]->setAxisTitle(
"jx", 1);
241 meSCSeedMapOccSingleCrystal_[0]->setAxisTitle(
"jy", 2);
243 histo =
"EECLTE SC seed occupancy map EE - (CSC triggered)";
244 meSCSeedMapOccTrg_[0][0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
245 meSCSeedMapOccTrg_[0][0]->setAxisTitle(
"jx", 1);
246 meSCSeedMapOccTrg_[0][0]->setAxisTitle(
"jy", 2);
248 histo =
"EECLTE SC seed occupancy map EE - (DT triggered)";
249 meSCSeedMapOccTrg_[0][1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
250 meSCSeedMapOccTrg_[0][1]->setAxisTitle(
"jx", 1);
251 meSCSeedMapOccTrg_[0][1]->setAxisTitle(
"jy", 2);
253 histo =
"EECLTE SC seed occupancy map EE - (ECAL triggered)";
254 meSCSeedMapOccTrg_[0][2] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
255 meSCSeedMapOccTrg_[0][2]->setAxisTitle(
"jx", 1);
256 meSCSeedMapOccTrg_[0][2]->setAxisTitle(
"jy", 2);
258 histo =
"EECLTE SC seed occupancy map EE - (HCAL triggered)";
259 meSCSeedMapOccTrg_[0][3] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
260 meSCSeedMapOccTrg_[0][3]->setAxisTitle(
"jx", 1);
261 meSCSeedMapOccTrg_[0][3]->setAxisTitle(
"jy", 2);
263 histo =
"EECLTE SC seed occupancy map EE - (RPC triggered)";
264 meSCSeedMapOccTrg_[0][4] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
265 meSCSeedMapOccTrg_[0][4]->setAxisTitle(
"jx", 1);
266 meSCSeedMapOccTrg_[0][4]->setAxisTitle(
"jy", 2);
268 histo =
"EECLTE SC seed occupancy map EE - (CSC exclusive triggered)";
269 meSCSeedMapOccTrgExcl_[0][0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
270 meSCSeedMapOccTrgExcl_[0][0]->setAxisTitle(
"jx", 1);
271 meSCSeedMapOccTrgExcl_[0][0]->setAxisTitle(
"jy", 2);
273 histo =
"EECLTE SC seed occupancy map EE - (DT exclusive triggered)";
274 meSCSeedMapOccTrgExcl_[0][1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
275 meSCSeedMapOccTrgExcl_[0][1]->setAxisTitle(
"jx", 1);
276 meSCSeedMapOccTrgExcl_[0][1]->setAxisTitle(
"jy", 2);
278 histo =
"EECLTE SC seed occupancy map EE - (ECAL exclusive triggered)";
279 meSCSeedMapOccTrgExcl_[0][2] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
280 meSCSeedMapOccTrgExcl_[0][2]->setAxisTitle(
"jx", 1);
281 meSCSeedMapOccTrgExcl_[0][2]->setAxisTitle(
"jy", 2);
283 histo =
"EECLTE SC seed occupancy map EE - (HCAL exclusive triggered)";
284 meSCSeedMapOccTrgExcl_[0][3] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
285 meSCSeedMapOccTrgExcl_[0][3]->setAxisTitle(
"jx", 1);
286 meSCSeedMapOccTrgExcl_[0][3]->setAxisTitle(
"jy", 2);
288 histo =
"EECLTE SC seed occupancy map EE - (RPC exclusive triggered)";
289 meSCSeedMapOccTrgExcl_[0][4] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
290 meSCSeedMapOccTrgExcl_[0][4]->setAxisTitle(
"jx", 1);
291 meSCSeedMapOccTrgExcl_[0][4]->setAxisTitle(
"jy", 2);
293 histo =
"EECLTE SC seed occupancy map super crystal binned EE +";
294 meSCSeedMapOccSC_[1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
295 meSCSeedMapOccSC_[1]->setAxisTitle(
"jx", 1);
296 meSCSeedMapOccSC_[1]->setAxisTitle(
"jy", 2);
298 histo =
"EECLTE SC seed occupancy map EE + (high energy clusters)";
299 meSCSeedMapOccHighEne_[1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
300 meSCSeedMapOccHighEne_[1]->setAxisTitle(
"jx", 1);
301 meSCSeedMapOccHighEne_[1]->setAxisTitle(
"jy", 2);
303 histo =
"EECLTE SC single crystal cluster seed occupancy map EE +";
304 meSCSeedMapOccSingleCrystal_[1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
305 meSCSeedMapOccSingleCrystal_[1]->setAxisTitle(
"jx", 1);
306 meSCSeedMapOccSingleCrystal_[1]->setAxisTitle(
"jy", 2);
308 histo =
"EECLTE SC seed occupancy map EE + (CSC triggered)";
309 meSCSeedMapOccTrg_[1][0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
310 meSCSeedMapOccTrg_[1][0]->setAxisTitle(
"jx", 1);
311 meSCSeedMapOccTrg_[1][0]->setAxisTitle(
"jy", 2);
313 histo =
"EECLTE SC seed occupancy map EE + (DT triggered)";
314 meSCSeedMapOccTrg_[1][1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
315 meSCSeedMapOccTrg_[1][1]->setAxisTitle(
"jx", 1);
316 meSCSeedMapOccTrg_[1][1]->setAxisTitle(
"jy", 2);
318 histo =
"EECLTE SC seed occupancy map EE + (ECAL triggered)";
319 meSCSeedMapOccTrg_[1][2] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
320 meSCSeedMapOccTrg_[1][2]->setAxisTitle(
"jx", 1);
321 meSCSeedMapOccTrg_[1][2]->setAxisTitle(
"jy", 2);
323 histo =
"EECLTE SC seed occupancy map EE + (HCAL triggered)";
324 meSCSeedMapOccTrg_[1][3] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
325 meSCSeedMapOccTrg_[1][3]->setAxisTitle(
"jx", 1);
326 meSCSeedMapOccTrg_[1][3]->setAxisTitle(
"jy", 2);
328 histo =
"EECLTE SC seed occupancy map EE + (RPC triggered)";
329 meSCSeedMapOccTrg_[1][4] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
330 meSCSeedMapOccTrg_[1][4]->setAxisTitle(
"jx", 1);
331 meSCSeedMapOccTrg_[1][4]->setAxisTitle(
"jy", 2);
333 histo =
"EECLTE SC seed occupancy map EE + (CSC exclusive triggered)";
334 meSCSeedMapOccTrgExcl_[1][0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
335 meSCSeedMapOccTrgExcl_[1][0]->setAxisTitle(
"jx", 1);
336 meSCSeedMapOccTrgExcl_[1][0]->setAxisTitle(
"jy", 2);
338 histo =
"EECLTE SC seed occupancy map EE + (DT exclusive triggered)";
339 meSCSeedMapOccTrgExcl_[1][1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
340 meSCSeedMapOccTrgExcl_[1][1]->setAxisTitle(
"jx", 1);
341 meSCSeedMapOccTrgExcl_[1][1]->setAxisTitle(
"jy", 2);
343 histo =
"EECLTE SC seed occupancy map EE + (ECAL exclusive triggered)";
344 meSCSeedMapOccTrgExcl_[1][2] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
345 meSCSeedMapOccTrgExcl_[1][2]->setAxisTitle(
"jx", 1);
346 meSCSeedMapOccTrgExcl_[1][2]->setAxisTitle(
"jy", 2);
348 histo =
"EECLTE SC seed occupancy map EE + (HCAL exclusive triggered)";
349 meSCSeedMapOccTrgExcl_[1][3] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
350 meSCSeedMapOccTrgExcl_[1][3]->setAxisTitle(
"jx", 1);
351 meSCSeedMapOccTrgExcl_[1][3]->setAxisTitle(
"jy", 2);
353 histo =
"EECLTE SC seed occupancy map EE + (RPC exclusive triggered)";
354 meSCSeedMapOccTrgExcl_[1][4] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
355 meSCSeedMapOccTrgExcl_[1][4]->setAxisTitle(
"jx", 1);
356 meSCSeedMapOccTrgExcl_[1][4]->setAxisTitle(
"jy", 2);
358 histo =
"EECLTE SC seed crystal timing map EE + super crystal binned";
359 meSCSeedMapTimeSC_[1] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
360 meSCSeedMapTimeSC_[1]->setAxisTitle(
"jx", 1);
361 meSCSeedMapTimeSC_[1]->setAxisTitle(
"jy", 2);
363 histo =
"EECLTE SC seed crystal timing map EE + super crystal binned";
364 meSCSeedMapTimeSC_[1] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
365 meSCSeedMapTimeSC_[1]->setAxisTitle(
"jx", 1);
366 meSCSeedMapTimeSC_[1]->setAxisTitle(
"jy", 2);
368 histo =
"EECLTE SC relative timing";
369 meSCSeedTime_ = dqmStore_->book1D(histo,histo,78,0.,10.);
370 meSCSeedTime_->setAxisTitle(
"seed crystal timing");
372 histo =
"EECLTE SC relative timing vs amplitude";
373 meSCSeedTimeVsAmp_ = dqmStore_->bookProfile(histo, histo, 78, -7, 7, 200, 0, 1.8);
374 meSCSeedTimeVsAmp_->setAxisTitle(
"seed crystal timing", 1);
375 meSCSeedTimeVsAmp_->setAxisTitle(
"energy (GeV)", 2);
377 histo =
"EECLTE SC relative timing EE -";
378 meSCSeedTimeEEM_ = dqmStore_->book1D(histo,histo,78,0.,10.);
379 meSCSeedTimeEEM_->setAxisTitle(
"seed crystal timing");
381 histo =
"EECLTE SC relative timing EE +";
382 meSCSeedTimeEEP_ = dqmStore_->book1D(histo,histo,78,0.,10.);
383 meSCSeedTimeEEP_->setAxisTitle(
"seed crystal timing");
385 std::stringstream ss;
386 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras/EECLTE timing per super module");
387 for(
int i=0;
i!=18;++
i) {
390 ss <<
"EECLTE SC timing EE - " <<
i+1;
394 ss <<
"EECLTE SC timing EE + " <<
i+1-9;
397 meSCSeedTimePerFed_[
i] = dqmStore_->book1D(histo,histo,78,0.,10.);
398 meSCSeedTimePerFed_[
i]->setAxisTitle(
"seed crystal timing", 1);
401 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras");
405 histo =
"EECLTE SC size (crystal) vs energy (GeV)";
406 meSCSizCrystalVsEne_ = dqmStore_->bookProfile(histo,histo,200,0.,10.,150,0,150);
407 meSCSizCrystalVsEne_->setAxisTitle(
"energy (GeV)", 1);
408 meSCSizCrystalVsEne_->setAxisTitle(
"super cluster size (crystal)", 2);
410 histo =
"EECLTE SC seed occupancy map EE -";
411 meSCSeedMapOcc_[0] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
412 meSCSeedMapOcc_[0]->setAxisTitle(
"jx", 1);
413 meSCSeedMapOcc_[0]->setAxisTitle(
"jy", 2);
415 histo =
"EECLTE SC seed occupancy map EE - (high energy clusters) super crystal binned";
416 meSCSeedMapOccHighEneSC_[0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
417 meSCSeedMapOccHighEneSC_[0]->setAxisTitle(
"jx", 1);
418 meSCSeedMapOccHighEneSC_[0]->setAxisTitle(
"jy", 2);
420 histo =
"EECLTE SC seed occupancy map EE - (CSC triggered) super crystal binned";
421 meSCSeedMapOccTrgSC_[0][0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
422 meSCSeedMapOccTrgSC_[0][0]->setAxisTitle(
"jx", 1);
423 meSCSeedMapOccTrgSC_[0][0]->setAxisTitle(
"jy", 2);
425 histo =
"EECLTE SC seed occupancy map EE - (DT triggered) super crystal binned";
426 meSCSeedMapOccTrgSC_[0][1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
427 meSCSeedMapOccTrgSC_[0][1]->setAxisTitle(
"jx", 1);
428 meSCSeedMapOccTrgSC_[0][1]->setAxisTitle(
"jy", 2);
430 histo =
"EECLTE SC seed occupancy map EE - (ECAL triggered) super crystal binned";
431 meSCSeedMapOccTrgSC_[0][2] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
432 meSCSeedMapOccTrgSC_[0][2]->setAxisTitle(
"jx", 1);
433 meSCSeedMapOccTrgSC_[0][2]->setAxisTitle(
"jy", 2);
435 histo =
"EECLTE SC seed occupancy map EE - (HCAL triggered) super crystal binned";
436 meSCSeedMapOccTrgSC_[0][3] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
437 meSCSeedMapOccTrgSC_[0][3]->setAxisTitle(
"jx", 1);
438 meSCSeedMapOccTrgSC_[0][3]->setAxisTitle(
"jy", 2);
440 histo =
"EECLTE SC seed occupancy map EE - (RPC triggered) super crystal binned";
441 meSCSeedMapOccTrgSC_[0][4] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
442 meSCSeedMapOccTrgSC_[0][4]->setAxisTitle(
"jx", 1);
443 meSCSeedMapOccTrgSC_[0][4]->setAxisTitle(
"jy", 2);
445 histo =
"EECLTE SC seed occupancy map EE - (CSC exclusive triggered) super crystal binned";
446 meSCSeedMapOccTrgExclSC_[0][0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
447 meSCSeedMapOccTrgExclSC_[0][0]->setAxisTitle(
"jx", 1);
448 meSCSeedMapOccTrgExclSC_[0][0]->setAxisTitle(
"jy", 2);
450 histo =
"EECLTE SC seed occupancy map EE - (DT exclusive triggered) super crystal binned";
451 meSCSeedMapOccTrgExclSC_[0][1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
452 meSCSeedMapOccTrgExclSC_[0][1]->setAxisTitle(
"jx", 1);
453 meSCSeedMapOccTrgExclSC_[0][1]->setAxisTitle(
"jy", 2);
455 histo =
"EECLTE SC seed occupancy map EE - (ECAL exclusive triggered) super crystal binned";
456 meSCSeedMapOccTrgExclSC_[0][2] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
457 meSCSeedMapOccTrgExclSC_[0][2]->setAxisTitle(
"jx", 1);
458 meSCSeedMapOccTrgExclSC_[0][2]->setAxisTitle(
"jy", 2);
460 histo =
"EECLTE SC seed occupancy map EE - (HCAL exclusive triggered) super crystal binned";
461 meSCSeedMapOccTrgExclSC_[0][3] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
462 meSCSeedMapOccTrgExclSC_[0][3]->setAxisTitle(
"jx", 1);
463 meSCSeedMapOccTrgExclSC_[0][3]->setAxisTitle(
"jy", 2);
465 histo =
"EECLTE SC seed occupancy map EE - (RPC exclusive triggered) super crystal binned";
466 meSCSeedMapOccTrgExclSC_[0][4] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
467 meSCSeedMapOccTrgExclSC_[0][4]->setAxisTitle(
"jx", 1);
468 meSCSeedMapOccTrgExclSC_[0][4]->setAxisTitle(
"jy", 2);
470 histo =
"EECLTE SC seed occupancy map EE +";
471 meSCSeedMapOcc_[1] = dqmStore_->book2D(histo,histo,100,0,100,100,0,100);
472 meSCSeedMapOcc_[1]->setAxisTitle(
"jx", 1);
473 meSCSeedMapOcc_[1]->setAxisTitle(
"jy", 2);
475 histo =
"EECLTE SC seed occupancy map EE + (high energy clusters) super crystal binned";
476 meSCSeedMapOccHighEneSC_[1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
477 meSCSeedMapOccHighEneSC_[1]->setAxisTitle(
"jx", 1);
478 meSCSeedMapOccHighEneSC_[1]->setAxisTitle(
"jy", 2);
480 histo =
"EECLTE SC seed occupancy map EE + (CSC triggered) super crystal binned";
481 meSCSeedMapOccTrgSC_[1][0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
482 meSCSeedMapOccTrgSC_[1][0]->setAxisTitle(
"jx", 1);
483 meSCSeedMapOccTrgSC_[1][0]->setAxisTitle(
"jy", 2);
485 histo =
"EECLTE SC seed occupancy map EE + (DT triggered) super crystal binned";
486 meSCSeedMapOccTrgSC_[1][1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
487 meSCSeedMapOccTrgSC_[1][1]->setAxisTitle(
"jx", 1);
488 meSCSeedMapOccTrgSC_[1][1]->setAxisTitle(
"jy", 2);
490 histo =
"EECLTE SC seed occupancy map EE + (ECAL triggered) super crystal binned";
491 meSCSeedMapOccTrgSC_[1][2] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
492 meSCSeedMapOccTrgSC_[1][2]->setAxisTitle(
"jx", 1);
493 meSCSeedMapOccTrgSC_[1][2]->setAxisTitle(
"jy", 2);
495 histo =
"EECLTE SC seed occupancy map EE + (HCAL triggered) super crystal binned";
496 meSCSeedMapOccTrgSC_[1][3] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
497 meSCSeedMapOccTrgSC_[1][3]->setAxisTitle(
"jx", 1);
498 meSCSeedMapOccTrgSC_[1][3]->setAxisTitle(
"jy", 2);
500 histo =
"EECLTE SC seed occupancy map EE + (RPC triggered) super crystal binned";
501 meSCSeedMapOccTrgSC_[1][4] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
502 meSCSeedMapOccTrgSC_[1][4]->setAxisTitle(
"jx", 1);
503 meSCSeedMapOccTrgSC_[1][4]->setAxisTitle(
"jy", 2);
505 histo =
"EECLTE SC seed occupancy map EE + (CSC exclusive triggered) super crystal binned";
506 meSCSeedMapOccTrgExclSC_[1][0] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
507 meSCSeedMapOccTrgExclSC_[1][0]->setAxisTitle(
"jx", 1);
508 meSCSeedMapOccTrgExclSC_[1][0]->setAxisTitle(
"jy", 2);
510 histo =
"EECLTE SC seed occupancy map EE + (DT exclusive triggered) super crystal binned";
511 meSCSeedMapOccTrgExclSC_[1][1] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
512 meSCSeedMapOccTrgExclSC_[1][1]->setAxisTitle(
"jx", 1);
513 meSCSeedMapOccTrgExclSC_[1][1]->setAxisTitle(
"jy", 2);
515 histo =
"EECLTE SC seed occupancy map EE + (ECAL exclusive triggered) super crystal binned";
516 meSCSeedMapOccTrgExclSC_[1][2] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
517 meSCSeedMapOccTrgExclSC_[1][2]->setAxisTitle(
"jx", 1);
518 meSCSeedMapOccTrgExclSC_[1][2]->setAxisTitle(
"jy", 2);
520 histo =
"EECLTE SC seed occupancy map EE + (HCAL exclusive triggered) super crystal binned";
521 meSCSeedMapOccTrgExclSC_[1][3] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
522 meSCSeedMapOccTrgExclSC_[1][3]->setAxisTitle(
"jx", 1);
523 meSCSeedMapOccTrgExclSC_[1][3]->setAxisTitle(
"jy", 2);
525 histo =
"EECLTE SC seed occupancy map EE + (RPC exclusive triggered) super crystal binned";
526 meSCSeedMapOccTrgExclSC_[1][4] = dqmStore_->book2D(histo,histo,20,0,100,20,0,100);
527 meSCSeedMapOccTrgExclSC_[1][4]->setAxisTitle(
"jx", 1);
528 meSCSeedMapOccTrgExclSC_[1][4]->setAxisTitle(
"jy", 2);
530 histo =
"EECLTE SC seed crystal timing map EE - (CSC exclusive triggered) super crystal binned";
531 meSCSeedMapTimeTrgSC_[0][0] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
532 meSCSeedMapTimeTrgSC_[0][0]->setAxisTitle(
"jx", 1);
533 meSCSeedMapTimeTrgSC_[0][0]->setAxisTitle(
"jy", 2);
535 histo =
"EECLTE SC seed crystal timing map EE - (DT exclusive triggered) super crystal binned";
536 meSCSeedMapTimeTrgSC_[0][1] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
537 meSCSeedMapTimeTrgSC_[0][1]->setAxisTitle(
"jx", 1);
538 meSCSeedMapTimeTrgSC_[0][1]->setAxisTitle(
"jy", 2);
540 histo =
"EECLTE SC seed crystal timing map EE - (ECAL exclusive triggered) super crystal binned";
541 meSCSeedMapTimeTrgSC_[0][2] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
542 meSCSeedMapTimeTrgSC_[0][2]->setAxisTitle(
"jx", 1);
543 meSCSeedMapTimeTrgSC_[0][2]->setAxisTitle(
"jy", 2);
545 histo =
"EECLTE SC seed crystal timing map EE - (HCAL exclusive triggered) super crystal binned";
546 meSCSeedMapTimeTrgSC_[0][3] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
547 meSCSeedMapTimeTrgSC_[0][3]->setAxisTitle(
"jx", 1);
548 meSCSeedMapTimeTrgSC_[0][3]->setAxisTitle(
"jy", 2);
550 histo =
"EECLTE SC seed crystal timing map EE - (RPC exclusive triggered) super crystal binned";
551 meSCSeedMapTimeTrgSC_[0][4] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
552 meSCSeedMapTimeTrgSC_[0][4]->setAxisTitle(
"jx", 1);
553 meSCSeedMapTimeTrgSC_[0][4]->setAxisTitle(
"jy", 2);
555 histo =
"EECLTE SC seed crystal timing map EE + (CSC exclusive triggered) super crystal binned";
556 meSCSeedMapTimeTrgSC_[1][0] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
557 meSCSeedMapTimeTrgSC_[1][0]->setAxisTitle(
"jx", 1);
558 meSCSeedMapTimeTrgSC_[1][0]->setAxisTitle(
"jy", 2);
560 histo =
"EECLTE SC seed crystal timing map EE + (DT exclusive triggered) super crystal binned";
561 meSCSeedMapTimeTrgSC_[1][1] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
562 meSCSeedMapTimeTrgSC_[1][1]->setAxisTitle(
"jx", 1);
563 meSCSeedMapTimeTrgSC_[1][1]->setAxisTitle(
"jy", 2);
565 histo =
"EECLTE SC seed crystal timing map EE + (ECAL exclusive triggered) super crystal binned";
566 meSCSeedMapTimeTrgSC_[1][2] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
567 meSCSeedMapTimeTrgSC_[1][2]->setAxisTitle(
"jx", 1);
568 meSCSeedMapTimeTrgSC_[1][2]->setAxisTitle(
"jy", 2);
570 histo =
"EECLTE SC seed crystal timing map EE + (HCAL exclusive triggered) super crystal binned";
571 meSCSeedMapTimeTrgSC_[1][3] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
572 meSCSeedMapTimeTrgSC_[1][3]->setAxisTitle(
"jx", 1);
573 meSCSeedMapTimeTrgSC_[1][3]->setAxisTitle(
"jy", 2);
575 histo =
"EECLTE SC seed crystal timing map EE + (RPC exclusive triggered) super crystal binned";
576 meSCSeedMapTimeTrgSC_[1][4] = dqmStore_->bookProfile2D(histo,histo,20,0,100,20,0,100,78,0.,10.,
"s");
577 meSCSeedMapTimeTrgSC_[1][4]->setAxisTitle(
"jx", 1);
578 meSCSeedMapTimeTrgSC_[1][4]->setAxisTitle(
"jy", 2);
584 if ( ! init_ )
return;
587 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras");
589 #ifndef EECLUSTERTASKEXTRAS_DQMOFFLINE
590 if ( meSCSizCrystal_ ) dqmStore_->removeElement( meSCSizCrystal_->getName() );
592 if ( meSCSizBC_ ) dqmStore_->removeElement( meSCSizBC_->getName() );
595 if ( meSCSeedEne_ ) dqmStore_->removeElement( meSCSeedEne_->getName() );
597 if ( meSCEne2_ ) dqmStore_->removeElement( meSCEne2_->getName() );
599 if ( meSCEneLow_ ) dqmStore_->removeElement( meSCEneLow_->getName() );
601 if ( meSCEneHigh_ ) dqmStore_->removeElement( meSCEneHigh_->getName() );
603 if ( meSCEneSingleCrystal_ ) dqmStore_->removeElement( meSCEneSingleCrystal_->getName() );
604 meSCEneSingleCrystal_ = 0;
606 for(
int i=0;
i!=2;++
i) {
607 if ( meSCSeedMapOccSC_[
i] ) dqmStore_->removeElement( meSCSeedMapOccSC_[i]->getName() );
608 meSCSeedMapOccSC_[
i] = 0;
609 if ( meSCSeedMapOccHighEne_[i] ) dqmStore_->removeElement( meSCSeedMapOccHighEne_[i]->getName() );
610 meSCSeedMapOccHighEne_[
i] = 0;
611 if ( meSCSeedMapOccSingleCrystal_[i] ) dqmStore_->removeElement( meSCSeedMapOccSingleCrystal_[i]->getName() );
612 meSCSeedMapOccSingleCrystal_[
i] = 0;
613 if ( meSCSeedMapTimeSC_[i] ) dqmStore_->removeElement( meSCSeedMapTimeSC_[i]->getName() );
614 meSCSeedMapTimeSC_[
i] = 0;
615 for(
int j=0;
j!=5;++
j) {
616 if ( meSCSeedMapOccTrg_[i][
j] ) dqmStore_->removeElement( meSCSeedMapOccTrg_[i][j]->getName() );
617 meSCSeedMapOccTrg_[
i][
j] = 0;
618 if ( meSCSeedMapOccTrgExcl_[i][j] ) dqmStore_->removeElement( meSCSeedMapOccTrgExcl_[i][j]->getName() );
619 meSCSeedMapOccTrgExcl_[
i][
j] = 0;
623 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras/EECLTE timing per super module");
624 for(
int i=0;
i!=18; ++
i) {
625 if( meSCSeedTimePerFed_[
i] ) dqmStore_->removeElement( meSCSeedTimePerFed_[i]->getName() );
626 meSCSeedTimePerFed_[
i] = 0;
628 dqmStore_->setCurrentFolder(prefixME_ +
"/EEClusterTaskExtras");
630 if ( meSCSeedTime_ ) dqmStore_->removeElement( meSCSeedTime_->getName() );
632 if ( meSCSeedTimeVsAmp_ ) dqmStore_->removeElement( meSCSeedTimeVsAmp_->getName() );
633 meSCSeedTimeVsAmp_ = 0;
634 if ( meSCSeedTimeEEM_ ) dqmStore_->removeElement( meSCSeedTimeEEM_->getName() );
635 meSCSeedTimeEEM_ = 0;
636 if ( meSCSeedTimeEEP_ ) dqmStore_->removeElement( meSCSeedTimeEEP_->getName() );
637 meSCSeedTimeEEP_ = 0;
641 if (meSCSizCrystalVsEne_ ) dqmStore_->removeElement(meSCSizCrystalVsEne_->getName() );
642 meSCSizCrystalVsEne_ = 0;
644 for(
int i=0;
i!=2;++
i) {
645 if ( meSCSeedMapOcc_[
i] ) dqmStore_->removeElement( meSCSeedMapOcc_[i]->getName() );
646 meSCSeedMapOcc_[
i] = 0;
647 if ( meSCSeedMapOccHighEneSC_[i] ) dqmStore_->removeElement( meSCSeedMapOccHighEneSC_[i]->getName() );
648 meSCSeedMapOccHighEneSC_[
i] = 0;
650 for(
int j=0;
j!=5; ++
j) {
651 if ( meSCSeedMapOccTrgSC_[i][
j] ) dqmStore_->removeElement( meSCSeedMapOccTrgSC_[i][j]->getName() );
652 meSCSeedMapOccTrgSC_[
i][
j] = 0;
653 if ( meSCSeedMapOccTrgExclSC_[i][j] ) dqmStore_->removeElement( meSCSeedMapOccTrgExclSC_[i][j]->getName() );
654 meSCSeedMapOccTrgExclSC_[
i][
j] = 0;
656 if ( meSCSeedMapTimeTrgSC_[i][j] ) dqmStore_->removeElement( meSCSeedMapTimeTrgSC_[i][j]->getName() );
657 meSCSeedMapTimeTrgSC_[
i][
j] = 0;
669 LogInfo(
"EEClusterTaskExtras") <<
"analyzed " << ievt_ <<
" events";
671 if ( enableCleanup_ ) this->
cleanup();
679 if ( ! init_ ) this->
setup();
687 if ( e.
getByLabel(SuperClusterCollection_, pSuperClusters) ) {
694 for ( SuperClusterCollection::const_iterator sCluster = pSuperClusters->begin(); sCluster != pSuperClusters->end(); ++sCluster ) {
698 e.
getByLabel( EcalRecHitCollection_, pEERecHits );
706 std::vector< std::pair<DetId,float> > sIds = sCluster->hitsAndFractions();
716 for(std::vector< std::pair<DetId,float> >::const_iterator idItr = sIds.begin(); idItr != sIds.end(); ++idItr) {
717 DetId id = idItr->first;
720 if(hitItr == eeRecHits->
end()) {
continue; }
721 if(hitItr->energy() > secondItr->energy()) { secondItr = hitItr; }
722 if(hitItr->energy() > seedItr->energy()) {
std::swap(seedItr,secondItr); }
729 int eey = seedId.
iy();
730 int eex = seedId.
ix();
731 float xeey = eey - 0.5;
732 float xeex = eex - 0.5;
734 int side = (ism >=1 && ism <= 9) ? 0 : 1;
739 vector<bool>
triggers = determineTriggers(e,c);
741 #ifndef EECLUSTERTASKEXTRAS_DQMOFFLINE
743 eMax = seedItr->energy();
744 e2nd = secondItr->energy();
746 if(meSCEneLow_) meSCEneLow_->Fill( sCluster->energy() );
747 if(meSCEneHigh_) meSCEneHigh_->Fill( sCluster->energy() );
748 if(meSCSizBC_) meSCSizBC_->Fill(
float(sCluster->clustersSize()) );
750 if(meSCSizCrystal_) meSCSizCrystal_->Fill(sIds.size());
751 if(meSCSeedEne_) meSCSeedEne_->Fill(eMax);
752 if(meSCEne2_) meSCEne2_->Fill(eMax+e2nd);
755 if(meSCSeedMapOccSC_[side]) meSCSeedMapOccSC_[side]->Fill(xeex, xeey);
757 if(sCluster->energy() > 2) {
758 if(meSCSeedMapOccHighEne_[side]) meSCSeedMapOccHighEne_[side]->Fill(xeex, xeey);
759 if(meSCSeedMapOccHighEneSC_[side]) meSCSeedMapOccHighEneSC_[side]->Fill(xeex, xeey);
761 if(sIds.size() == 1) {
762 if(meSCEneSingleCrystal_) meSCEneSingleCrystal_->Fill(sCluster->energy());
763 if(meSCSeedMapOccSingleCrystal_[side]) meSCSeedMapOccSingleCrystal_[side]->Fill(xeex, xeey);
766 if(meSCSeedMapOcc_[side]) meSCSeedMapOcc_[side]->Fill(xeex, xeey);
770 if(seedItr->energy() / agc->
getEBValue() > 12) {
771 if(meSCSeedTime_) meSCSeedTime_->Fill(seedItr->time());
772 if(meSCSeedTimeVsAmp_) meSCSeedTimeVsAmp_->Fill(seedItr->time(),sCluster->energy());
774 if(meSCSeedTimeEEM_) meSCSeedTimeEEM_->Fill(seedItr->time());
776 if(meSCSeedTimeEEP_) meSCSeedTimeEEP_->Fill(seedItr->time());
777 if(meSCSeedTimePerFed_[ism-1]) meSCSeedTimePerFed_[ism-1]->Fill(seedItr->time());
778 if(meSCSeedMapTimeSC_[side]) meSCSeedMapTimeSC_[side]->Fill(xeex,xeey,seedItr->time());
783 LogWarning(
"EBClusterTaskExtras") <<
"EcalADCToGeVConstant not valid";
785 for(
int i=0;
i!=5;++
i) {
787 if(meSCSeedMapOccTrg_[side][i]) meSCSeedMapOccTrg_[side][
i]->Fill(xeex, xeey);
788 bool isExclusive =
true;
789 for(
int j=0;
j!=5;++
j) {
790 if(
j != i && triggers[
j])
794 if(meSCSeedMapOccTrgExcl_[side][i]) meSCSeedMapOccTrgExcl_[side][
i]->Fill(xeex, xeey);
799 if(meSCSizCrystalVsEne_) meSCSizCrystalVsEne_->Fill(sCluster->energy(),sIds.size());
801 for(
int i=0;
i!=5;++
i) {
803 if(meSCSeedMapOccTrgSC_[side][i]) meSCSeedMapOccTrgSC_[side][
i]->Fill(xeex, xeey);
807 if(seedItr->energy() / agc->
getEBValue() > 12) {
808 if(meSCSeedMapTimeTrgSC_[side][i]) meSCSeedMapTimeTrgSC_[side][
i]->Fill(xeex, xeey, seedItr->time());
812 LogWarning(
"EBClusterTaskExtras") <<
"EcalADCToGeVConstant not valid";
815 bool isExclusive =
true;
816 for(
int j=0;
j!=5;++
j) {
817 if(
j != i && triggers[
j])
821 if(meSCSeedMapOccTrgExclSC_[side][i]) meSCSeedMapOccTrgExclSC_[side][
i]->Fill(xeex, xeey);
826 LogWarning(
"EEClusterTaskExtras") << pEERecHits <<
" not available";
832 LogWarning(
"EEClusterTaskExtras") << SuperClusterCollection_ <<
" not available";
842 std::vector<bool> l1Triggers;
845 l1Triggers.push_back(
false);
849 iEvent.
getByLabel(l1GMTReadoutRecTag_,gmtrc_handle);
851 if (!(gmtrc_handle.isValid()))
853 LogWarning(
"EcalCosmicsHists") <<
"l1MuGMTReadoutCollection" <<
" not available";
869 bool l1SingleEG2 =
menu->gtAlgorithmResult(
"L1_SingleEG2", dWord);
870 bool l1SingleEG5 =
menu->gtAlgorithmResult(
"L1_SingleEG5", dWord);
871 bool l1SingleEG8 =
menu->gtAlgorithmResult(
"L1_SingleEG8", dWord);
872 bool l1SingleEG10 =
menu->gtAlgorithmResult(
"L1_SingleEG10", dWord);
873 bool l1SingleEG12 =
menu->gtAlgorithmResult(
"L1_SingleEG12", dWord);
874 bool l1SingleEG15 =
menu->gtAlgorithmResult(
"L1_SingleEG15", dWord);
875 bool l1SingleEG20 =
menu->gtAlgorithmResult(
"L1_SingleEG20", dWord);
876 bool l1SingleEG25 =
menu->gtAlgorithmResult(
"L1_SingleEG25", dWord);
877 bool l1DoubleNoIsoEGBTBtight =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEG_BTB_tight", dWord);
878 bool l1DoubleNoIsoEGBTBloose =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEG_BTB_loose ", dWord);
879 bool l1DoubleNoIsoEGTopBottom =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEGTopBottom", dWord);
880 bool l1DoubleNoIsoEGTopBottomCen =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEGTopBottomCen", dWord);
881 bool l1DoubleNoIsoEGTopBottomCen2 =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEGTopBottomCen2", dWord);
882 bool l1DoubleNoIsoEGTopBottomCenVert =
menu->gtAlgorithmResult(
"L1_DoubleNoIsoEGTopBottomCenVert", dWord);
884 l1Triggers[ECAL_TRIGGER] = l1SingleEG2 || l1SingleEG5 || l1SingleEG8 || l1SingleEG10 || l1SingleEG12 || l1SingleEG15
885 || l1SingleEG20 || l1SingleEG25 || l1DoubleNoIsoEGBTBtight || l1DoubleNoIsoEGBTBloose
886 || l1DoubleNoIsoEGTopBottom || l1DoubleNoIsoEGTopBottomCen || l1DoubleNoIsoEGTopBottomCen2
887 || l1DoubleNoIsoEGTopBottomCenVert;
889 std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
890 std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
891 for(igmtrr=gmt_records.begin(); igmtrr!=gmt_records.end(); igmtrr++) {
892 std::vector<L1MuRegionalCand>::const_iterator iter1;
893 std::vector<L1MuRegionalCand> rmc;
897 rmc = igmtrr->getDTBXCands();
898 for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
899 if ( !(*iter1).empty() ) {
905 if(igmtrr->getBxInEvent()==0 && idt>0) l1Triggers[DT_TRIGGER] =
true;
909 rmc = igmtrr->getBrlRPCCands();
910 for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
911 if ( !(*iter1).empty() ) {
917 if(igmtrr->getBxInEvent()==0 && irpcb>0) l1Triggers[RPC_TRIGGER] =
true;
921 rmc = igmtrr->getCSCCands();
922 for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
923 if ( !(*iter1).empty() ) {
929 if(igmtrr->getBxInEvent()==0 && icsc>0) l1Triggers[CSC_TRIGGER] =
true;
934 for(
int ibx=-1; ibx<=1; ibx++) {
935 bool hcal_top =
false;
936 bool hcal_bot =
false;
937 const L1GtPsbWord psb = gtrr->gtPsbWord(0xbb0d,ibx);
938 std::vector<int> valid_x;
939 if((psb.
aData(4)&0x3f) >= 1) {valid_x.push_back( (psb.
aData(4)>>10)&0x1f ); }
940 if((psb.
bData(4)&0x3f) >= 1) {valid_x.push_back( (psb.
bData(4)>>10)&0x1f ); }
941 if((psb.
aData(5)&0x3f) >= 1) {valid_x.push_back( (psb.
aData(5)>>10)&0x1f ); }
942 if((psb.
bData(5)&0x3f) >= 1) {valid_x.push_back( (psb.
bData(5)>>10)&0x1f ); }
943 std::vector<int>::const_iterator ix;
944 for(ix=valid_x.begin(); ix!=valid_x.end(); ix++) {
946 if(*ix<9) hcal_top=
true;
947 if(*ix>8) hcal_bot=
true;
949 if(ibx==0 && hcal_top && hcal_bot) l1Triggers[HCAL_TRIGGER]=
true;
T getParameter(std::string const &) const
const cms_uint16_t bData(int iB) const
get/set B_DATA_CH_IB
T getUntrackedParameter(std::string const &, T const &) const
Some "id" conversions.
std::vector< EcalRecHit >::const_iterator const_iterator
static void cleanup(const Factory::MakerMap::value_type &v)
const cms_uint16_t aData(int iA) const
get/set A_DATA_CH_IA
std::vector< bool > DecisionWord
typedefs
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
const_iterator end() const
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
T const * product() const
T const * product() const
iterator find(key_type k)
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
int ism(int ieta, int iphi)
void reset(double vett[256])
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
const_iterator begin() const