CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEClusterTask.cc
Go to the documentation of this file.
1 /*
2  * \file EEClusterTask.cc
3  *
4  * $Date: 2012/04/27 13:46:14 $
5  * $Revision: 1.86 $
6  * \author G. Della Ricca
7  * \author E. Di Marco
8  *
9 */
10 
11 #include <iostream>
12 #include <fstream>
13 #include <vector>
14 #include <math.h>
15 
16 #include "FWCore/ServiceRegistry/interface/Service.h"
17 #include "FWCore/MessageLogger/interface/MessageLogger.h"
18 
19 #include "DQMServices/Core/interface/MonitorElement.h"
20 
21 #include "DQMServices/Core/interface/DQMStore.h"
22 
23 #include "DataFormats/EcalRawData/interface/EcalRawDataCollections.h"
24 #include "DataFormats/EgammaReco/interface/BasicCluster.h"
25 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
26 #include "DataFormats/EgammaReco/interface/SuperCluster.h"
27 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
28 #include "FWCore/Framework/interface/ESHandle.h"
29 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"
30 #include "Geometry/CaloTopology/interface/CaloTopology.h"
31 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
32 #include "DataFormats/Math/interface/Point3D.h"
33 #include "DataFormats/EcalDetId/interface/EEDetId.h"
34 #include "CondFormats/EcalObjects/interface/EcalADCToGeVConstant.h"
35 #include "CondFormats/DataRecord/interface/EcalADCToGeVConstantRcd.h"
36 
37 #include "DQM/EcalCommon/interface/Numbers.h"
38 
39 #include "DQM/EcalEndcapMonitorTasks/interface/EEClusterTask.h"
40 
41 #include "TLorentzVector.h"
42 
43 EEClusterTask::EEClusterTask(const edm::ParameterSet& ps){
44 
45  init_ = false;
46 
47  dqmStore_ = edm::Service<DQMStore>().operator->();
48 
49  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
50 
51  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
52 
53  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
54 
55  // parameters...
56  EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
57  BasicClusterCollection_ = ps.getParameter<edm::InputTag>("BasicClusterCollection");
58  SuperClusterCollection_ = ps.getParameter<edm::InputTag>("SuperClusterCollection");
59  EcalRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalRecHitCollection");
60 
61  // histograms...
62  meBCEne_ = 0;
63  meBCNum_ = 0;
64  meBCSiz_ = 0;
65 
66  meBCEneFwdMap_ = 0;
67  meBCNumFwdMap_ = 0;
68  meBCETFwdMap_ = 0;
69  meBCSizFwdMap_ = 0;
70 
71  meBCEneFwdMapProjEta_ = 0;
72  meBCEneFwdMapProjPhi_ = 0;
73  meBCNumFwdMapProjEta_ = 0;
74  meBCNumFwdMapProjPhi_ = 0;
75  meBCETFwdMapProjEta_ = 0;
76  meBCETFwdMapProjPhi_ = 0;
77  meBCSizFwdMapProjEta_ = 0;
78  meBCSizFwdMapProjPhi_ = 0;
79 
80  meBCEneBwdMap_ = 0;
81  meBCNumBwdMap_ = 0;
82  meBCETBwdMap_ = 0;
83  meBCSizBwdMap_ = 0;
84 
85  meBCEneBwdMapProjEta_ = 0;
86  meBCEneBwdMapProjPhi_ = 0;
87  meBCNumBwdMapProjEta_ = 0;
88  meBCNumBwdMapProjPhi_ = 0;
89  meBCETBwdMapProjEta_ = 0;
90  meBCETBwdMapProjPhi_ = 0;
91  meBCSizBwdMapProjEta_ = 0;
92  meBCSizBwdMapProjPhi_ = 0;
93 
94  meSCEne_ = 0;
95  meSCNum_ = 0;
96  meSCSiz_ = 0;
97 
98  meSCCrystalSiz_ = 0;
99  meSCSeedEne_ = 0;
100  meSCEne2_ = 0;
101  meSCEneVsEMax_ = 0;
102  meSCEneLowScale_ = 0;
103  meSCSeedMapOcc_[0] = 0;
104  meSCSeedMapOcc_[1] = 0;
105  meSCMapSingleCrystal_[0] = 0;
106  meSCMapSingleCrystal_[1] = 0;
107 
108  mes1s9_ = 0;
109  mes1s9thr_ = 0;
110  mes9s25_ = 0;
111 
112  meInvMassPi0_ = 0;
113  meInvMassJPsi_ = 0;
114  meInvMassZ0_ = 0;
115  meInvMassHigh_ = 0;
116 
117  meInvMassPi0Sel_ = 0;
118  meInvMassJPsiSel_ = 0;
119  meInvMassZ0Sel_ = 0;
120  meInvMassHighSel_ = 0;
121 
122  thrS4S9_ = 0.85;
123  thrClusEt_ = 0.250;
124  thrCandEt_ = 0.800;
125 
126 }
127 
128 EEClusterTask::~EEClusterTask(){
129 
130 }
131 
132 void EEClusterTask::beginJob(void){
133 
134  ievt_ = 0;
135 
136  if ( dqmStore_ ) {
137  dqmStore_->setCurrentFolder(prefixME_ + "/EEClusterTask");
138  dqmStore_->rmdir(prefixME_ + "/EEClusterTask");
139  }
140 
141 }
142 
143 void EEClusterTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
144 
145  Numbers::initGeometry(c, false);
146 
147  if ( ! mergeRuns_ ) this->reset();
148 
149 }
150 
151 void EEClusterTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
152 
153 }
154 
155 void EEClusterTask::reset(void) {
156 
157  if ( meBCEne_ ) meBCEne_->Reset();
158 
159  if ( meBCNum_ ) meBCNum_->Reset();
160 
161  if ( meBCSiz_ ) meBCSiz_->Reset();
162 
163  if ( meBCEneFwdMap_ ) meBCEneFwdMap_->Reset();
164 
165  if ( meBCNumFwdMap_ ) meBCNumFwdMap_->Reset();
166 
167  if ( meBCETFwdMap_ ) meBCETFwdMap_->Reset();
168 
169  if ( meBCSizFwdMap_ ) meBCSizFwdMap_->Reset();
170 
171  if ( meBCEneFwdMapProjEta_ ) meBCEneFwdMapProjEta_->Reset();
172 
173  if ( meBCEneFwdMapProjPhi_ ) meBCEneFwdMapProjPhi_->Reset();
174 
175  if ( meBCNumFwdMapProjEta_ ) meBCNumFwdMapProjEta_->Reset();
176 
177  if ( meBCNumFwdMapProjPhi_ ) meBCNumFwdMapProjPhi_->Reset();
178 
179  if ( meBCETFwdMapProjEta_ ) meBCETFwdMapProjEta_->Reset();
180 
181  if ( meBCETFwdMapProjPhi_ ) meBCETFwdMapProjPhi_->Reset();
182 
183  if ( meBCSizFwdMapProjEta_ ) meBCSizFwdMapProjEta_->Reset();
184 
185  if ( meBCSizFwdMapProjPhi_ ) meBCSizFwdMapProjPhi_->Reset();
186 
187  if ( meBCEneBwdMap_ ) meBCEneBwdMap_->Reset();
188 
189  if ( meBCNumBwdMap_ ) meBCNumBwdMap_->Reset();
190 
191  if ( meBCETBwdMap_ ) meBCETBwdMap_->Reset();
192 
193  if ( meBCSizBwdMap_ ) meBCSizBwdMap_->Reset();
194 
195  if ( meBCEneBwdMapProjEta_ ) meBCEneBwdMapProjEta_->Reset();
196 
197  if ( meBCEneBwdMapProjPhi_ ) meBCEneBwdMapProjPhi_->Reset();
198 
199  if ( meBCNumBwdMapProjEta_ ) meBCNumBwdMapProjEta_->Reset();
200 
201  if ( meBCNumBwdMapProjPhi_ ) meBCNumBwdMapProjPhi_->Reset();
202 
203  if ( meBCETBwdMapProjEta_ ) meBCETBwdMapProjEta_->Reset();
204 
205  if ( meBCETBwdMapProjPhi_ ) meBCETBwdMapProjPhi_->Reset();
206 
207  if ( meBCSizBwdMapProjEta_ ) meBCSizBwdMapProjEta_->Reset();
208 
209  if ( meBCSizBwdMapProjPhi_ ) meBCSizBwdMapProjPhi_->Reset();
210 
211  if ( meSCEne_ ) meSCEne_->Reset();
212 
213  if ( meSCNum_ ) meSCNum_->Reset();
214 
215  if ( meSCSiz_ ) meSCSiz_->Reset();
216 
217  if ( meSCCrystalSiz_ ) meSCCrystalSiz_->Reset();
218 
219  if ( meSCSeedEne_ ) meSCSeedEne_->Reset();
220 
221  if ( meSCEne2_ ) meSCEne2_->Reset();
222 
223  if ( meSCEneVsEMax_ ) meSCEneVsEMax_->Reset();
224 
225  if ( meSCEneLowScale_ ) meSCEneLowScale_->Reset();
226 
227  if ( meSCSeedMapOcc_[0] ) meSCSeedMapOcc_[0]->Reset();
228 
229  if ( meSCSeedMapOcc_[1] ) meSCSeedMapOcc_[1]->Reset();
230 
231  if ( meSCMapSingleCrystal_[0] ) meSCMapSingleCrystal_[0]->Reset();
232 
233  if ( meSCMapSingleCrystal_[1] ) meSCMapSingleCrystal_[1]->Reset();
234 
235  if ( mes1s9_ ) mes1s9_->Reset();
236 
237  if ( mes1s9thr_ ) mes1s9thr_->Reset();
238 
239  if ( mes9s25_ ) mes9s25_->Reset();
240 
241  if ( meInvMassPi0_ ) meInvMassPi0_->Reset();
242 
243  if ( meInvMassJPsi_ ) meInvMassJPsi_->Reset();
244 
245  if ( meInvMassZ0_ ) meInvMassZ0_->Reset();
246 
247  if ( meInvMassHigh_ ) meInvMassHigh_->Reset();
248 
249  if ( meInvMassPi0Sel_ ) meInvMassPi0Sel_->Reset();
250 
251  if ( meInvMassJPsiSel_ ) meInvMassJPsiSel_->Reset();
252 
253  if ( meInvMassZ0Sel_ ) meInvMassZ0Sel_->Reset();
254 
255  if ( meInvMassHighSel_ ) meInvMassHighSel_->Reset();
256 
257 }
258 
259 void EEClusterTask::setup(void){
260 
261  init_ = true;
262 
263  std::string name;
264 
265  if ( dqmStore_ ) {
266  dqmStore_->setCurrentFolder(prefixME_ + "/EEClusterTask");
267 
268  name = "EECLT BC energy";
269  meBCEne_ = dqmStore_->book1D(name, name, 100, 0., 150.);
270  meBCEne_->setAxisTitle("energy (GeV)", 1);
271 
272  name = "EECLT BC number";
273  meBCNum_ = dqmStore_->book1D(name, name, 100, 0., 150.);
274  meBCNum_->setAxisTitle("number of clusters", 1);
275 
276  name = "EECLT BC size";
277  meBCSiz_ = dqmStore_->book1D(name, name, 100, 0., 150.);
278  meBCSiz_->setAxisTitle("cluster size", 1);
279 
280  name = "EECLT BC energy map EE +";
281  meBCEneFwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 500., "s");
282  meBCEneFwdMap_->setAxisTitle("x", 1);
283  meBCEneFwdMap_->setAxisTitle("y", 2);
284 
285  name = "EECLT BC number map EE +";
286  meBCNumFwdMap_ = dqmStore_->book2D(name, name, 20, -150., 150., 20, -150., 150.);
287  meBCNumFwdMap_->setAxisTitle("x", 1);
288  meBCNumFwdMap_->setAxisTitle("y", 2);
289 
290  name = "EECLT BC ET map EE +";
291  meBCETFwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 500., "s");
292  meBCETFwdMap_->setAxisTitle("x", 1);
293  meBCETFwdMap_->setAxisTitle("y", 2);
294 
295  name = "EECLT BC size map EE +";
296  meBCSizFwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 100., "s");
297  meBCSizFwdMap_->setAxisTitle("x", 1);
298  meBCSizFwdMap_->setAxisTitle("y", 2);
299 
300  name = "EECLT BC energy projection eta EE +";
301  meBCEneFwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, 1.479, 3.0, 100, 0., 500., "s");
302  meBCEneFwdMapProjEta_->setAxisTitle("eta", 1);
303  meBCEneFwdMapProjEta_->setAxisTitle("energy (GeV)", 2);
304 
305  name = "EECLT BC energy projection phi EE +";
306  meBCEneFwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 500., "s");
307  meBCEneFwdMapProjPhi_->setAxisTitle("phi", 1);
308  meBCEneFwdMapProjPhi_->setAxisTitle("energy (GeV)", 2);
309 
310  name = "EECLT BC number projection eta EE +";
311  meBCNumFwdMapProjEta_ = dqmStore_->book1D(name, name, 20, 1.479, 3.0);
312  meBCNumFwdMapProjEta_->setAxisTitle("eta", 1);
313  meBCNumFwdMapProjEta_->setAxisTitle("number of clusters", 2);
314 
315  name = "EECLT BC number projection phi EE +";
316  meBCNumFwdMapProjPhi_ = dqmStore_->book1D(name, name, 50, -M_PI, M_PI);
317  meBCNumFwdMapProjPhi_->setAxisTitle("phi", 1);
318  meBCNumFwdMapProjPhi_->setAxisTitle("number of clusters", 2);
319 
320  name = "EECLT BC ET projection eta EE +";
321  meBCETFwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, 1.479, 3.0, 100, 0., 500., "s");
322  meBCETFwdMapProjEta_->setAxisTitle("eta", 1);
323  meBCETFwdMapProjEta_->setAxisTitle("transverse energy (GeV)", 2);
324 
325  name = "EECLT BC ET projection phi EE +";
326  meBCETFwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 500., "s");
327  meBCETFwdMapProjPhi_->setAxisTitle("phi", 1);
328  meBCETFwdMapProjPhi_->setAxisTitle("transverse energy (GeV)", 2);
329 
330  name = "EECLT BC size projection eta EE +";
331  meBCSizFwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, 1.479, 3.0, 100, 0., 100., "s");
332  meBCSizFwdMapProjEta_->setAxisTitle("eta", 1);
333  meBCSizFwdMapProjEta_->setAxisTitle("cluster size", 2);
334 
335  name = "EECLT BC size projection phi EE +";
336  meBCSizFwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 100., "s");
337  meBCSizFwdMapProjPhi_->setAxisTitle("phi", 1);
338  meBCSizFwdMapProjPhi_->setAxisTitle("cluster size", 2);
339 
340  name = "EECLT BC energy map EE -";
341  meBCEneBwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 500., "s");
342  meBCEneBwdMap_->setAxisTitle("x", 1);
343  meBCEneBwdMap_->setAxisTitle("y", 2);
344 
345  name = "EECLT BC number map EE -";
346  meBCNumBwdMap_ = dqmStore_->book2D(name, name, 20, -150., 150., 20, -150., 150.);
347  meBCNumBwdMap_->setAxisTitle("x", 1);
348  meBCNumBwdMap_->setAxisTitle("y", 2);
349 
350  name = "EECLT BC ET map EE -";
351  meBCETBwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 500., "s");
352  meBCETBwdMap_->setAxisTitle("x", 1);
353  meBCETBwdMap_->setAxisTitle("y", 2);
354 
355  name = "EECLT BC size map EE -";
356  meBCSizBwdMap_ = dqmStore_->bookProfile2D(name, name, 20, -150., 150., 20, -150., 150., 100, 0., 100., "s");
357  meBCSizBwdMap_->setAxisTitle("x", 1);
358  meBCSizBwdMap_->setAxisTitle("y", 2);
359 
360  name = "EECLT BC energy projection eta EE -";
361  meBCEneBwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, -3.0, -1.479, 100, 0., 500., "s");
362  meBCEneBwdMapProjEta_->setAxisTitle("eta", 1);
363  meBCEneBwdMapProjEta_->setAxisTitle("energy (GeV)", 2);
364 
365  name = "EECLT BC energy projection phi EE -";
366  meBCEneBwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 500., "s");
367  meBCEneBwdMapProjPhi_->setAxisTitle("phi", 1);
368  meBCEneBwdMapProjPhi_->setAxisTitle("energy (GeV)", 2);
369 
370  name = "EECLT BC number projection eta EE -";
371  meBCNumBwdMapProjEta_ = dqmStore_->book1D(name, name, 20, -3.0, -1.479);
372  meBCNumBwdMapProjEta_->setAxisTitle("eta", 1);
373  meBCNumBwdMapProjEta_->setAxisTitle("number of clusters", 2);
374 
375  name = "EECLT BC number projection phi EE -";
376  meBCNumBwdMapProjPhi_ = dqmStore_->book1D(name, name, 50, -M_PI, M_PI);
377  meBCNumBwdMapProjPhi_->setAxisTitle("phi", 1);
378  meBCNumBwdMapProjPhi_->setAxisTitle("number of clusters", 2);
379 
380  name = "EECLT BC ET projection eta EE -";
381  meBCETBwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, -3.0, -1.479, 100, 0., 500., "s");
382  meBCETBwdMapProjEta_->setAxisTitle("eta", 1);
383  meBCETBwdMapProjEta_->setAxisTitle("transverse energy (GeV)", 2);
384 
385  name = "EECLT BC ET projection phi EE -";
386  meBCETBwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 500., "s");
387  meBCETBwdMapProjPhi_->setAxisTitle("phi", 1);
388  meBCETBwdMapProjPhi_->setAxisTitle("transverse energy (GeV)", 2);
389 
390  name = "EECLT BC size projection eta EE -";
391  meBCSizBwdMapProjEta_ = dqmStore_->bookProfile(name, name, 20, -3.0, -1.479, 100, 0., 100., "s");
392  meBCSizBwdMapProjEta_->setAxisTitle("eta", 1);
393  meBCSizBwdMapProjEta_->setAxisTitle("cluster size", 2);
394 
395  name = "EECLT BC size projection phi EE -";
396  meBCSizBwdMapProjPhi_ = dqmStore_->bookProfile(name, name, 50, -M_PI, M_PI, 100, 0., 100., "s");
397  meBCSizBwdMapProjPhi_->setAxisTitle("phi", 1);
398  meBCSizBwdMapProjPhi_->setAxisTitle("cluster size", 2);
399 
400  name = "EECLT SC energy";
401  meSCEne_ = dqmStore_->book1D(name, name, 100, 0., 150.);
402  meSCEne_->setAxisTitle("energy (GeV)", 1);
403 
404  name = "EECLT SC number";
405  meSCNum_ = dqmStore_->book1D(name, name, 50, 0., 50.);
406  meSCNum_->setAxisTitle("number of clusters", 1);
407 
408  name = "EECLT SC size";
409  meSCSiz_ = dqmStore_->book1D(name, name, 50, 0., 50.);
410  meSCSiz_->setAxisTitle("cluster size", 1);
411 
412  name = "EECLT SC size (crystal)";
413  meSCCrystalSiz_ = dqmStore_->book1D(name, name, 150, 0, 150);
414  meSCCrystalSiz_->setAxisTitle("cluster size in crystals", 1);
415 
416  name = "EECLT SC seed crystal energy";
417  meSCSeedEne_ = dqmStore_->book1D(name, name, 100, 0., 10.);
418  meSCSeedEne_->setAxisTitle("seed crystal energy (GeV)", 1);
419 
420  name = "EECLT SC e2";
421  meSCEne2_ = dqmStore_->book1D(name, name, 100, 0., 10.);
422  meSCEne2_->setAxisTitle("seed + highest neighbor crystal energy (GeV)", 1);
423 
424  name = "EECLT SC energy vs seed crystal energy";
425  meSCEneVsEMax_ = dqmStore_->book2D(name, name, 50, 0., 10., 50, 0., 10.);
426  meSCEneVsEMax_->setAxisTitle("seed crystal energy (GeV)", 1);
427  meSCEneVsEMax_->setAxisTitle("cluster energy (GeV)", 2);
428 
429  name = "EECLT SC energy (low scale)";
430  meSCEneLowScale_ = dqmStore_->book1D(name, name, 100, 0., 10.);
431  meSCEneLowScale_->setAxisTitle("cluster energy (GeV)", 1);
432 
433  name = "EECLT SC seed occupancy map EE -";
434  meSCSeedMapOcc_[0] = dqmStore_->book2D(name, name, 20, 0., 100., 20, 0., 100.);
435  meSCSeedMapOcc_[0]->setAxisTitle("jx'", 1);
436  meSCSeedMapOcc_[0]->setAxisTitle("jy'", 2);
437 
438  name = "EECLT SC seed occupancy map EE +";
439  meSCSeedMapOcc_[1] = dqmStore_->book2D(name, name, 20, 0., 100., 20, 0., 100.);
440  meSCSeedMapOcc_[1]->setAxisTitle("jx'", 1);
441  meSCSeedMapOcc_[1]->setAxisTitle("jy'", 2);
442 
443  name = "EECLT SC single crystal cluster seed occupancy map EE -";
444  meSCMapSingleCrystal_[0] = dqmStore_->book2D(name, name, 20, 0., 100., 20, 0., 100.);
445  meSCMapSingleCrystal_[0]->setAxisTitle("jx'", 1);
446  meSCMapSingleCrystal_[0]->setAxisTitle("jy'", 2);
447 
448  name = "EECLT SC single crystal cluster seed occupancy map EE +";
449  meSCMapSingleCrystal_[1] = dqmStore_->book2D(name, name, 20, 0., 100., 20, 0., 100.);
450  meSCMapSingleCrystal_[1]->setAxisTitle("jx'", 1);
451  meSCMapSingleCrystal_[1]->setAxisTitle("jy'", 2);
452 
453  name = "EECLT s1s9";
454  mes1s9_ = dqmStore_->book1D(name, name, 50, 0., 1.5);
455  mes1s9_->setAxisTitle("s1/s9", 1);
456 
457  name = "EECLT s1s9 thr";
458  mes1s9thr_ = dqmStore_->book1D(name, name, 50, 0., 1.5);
459  mes1s9thr_->setAxisTitle("s1/s9", 1);
460 
461  name = "EECLT s9s25";
462  mes9s25_ = dqmStore_->book1D(name, name, 75, 0., 1.5);
463  mes9s25_->setAxisTitle("s9/s25", 1);
464 
465  name = "EECLT dicluster invariant mass Pi0";
466  meInvMassPi0_ = dqmStore_->book1D(name, name, 50, 0.0, 0.500);
467  meInvMassPi0_->setAxisTitle("mass (GeV)", 1);
468 
469  name = "EECLT dicluster invariant mass JPsi";
470  meInvMassJPsi_ = dqmStore_->book1D(name, name, 50, 2.9, 3.3);
471  meInvMassJPsi_->setAxisTitle("mass (GeV)", 1);
472 
473  name = "EECLT dicluster invariant mass Z0";
474  meInvMassZ0_ = dqmStore_->book1D(name, name, 50, 40, 110);
475  meInvMassZ0_->setAxisTitle("mass (GeV)", 1);
476 
477  name = "EECLT dicluster invariant mass high";
478  meInvMassHigh_ = dqmStore_->book1D(name, name, 500, 110, 3000);
479  meInvMassHigh_->setAxisTitle("mass (GeV)", 1);
480 
481  name = "EECLT dicluster invariant mass Pi0 sel";
482  meInvMassPi0Sel_ = dqmStore_->book1D(name, name, 50, 0.00, 0.500);
483  meInvMassPi0Sel_->setAxisTitle("mass (GeV)", 1);
484 
485  name = "EECLT dicluster invariant mass JPsi sel";
486  meInvMassJPsiSel_ = dqmStore_->book1D(name, name, 50, 2.9, 3.3);
487  meInvMassJPsiSel_->setAxisTitle("mass (GeV)", 1);
488 
489  name = "EECLT dicluster invariant mass Z0 sel";
490  meInvMassZ0Sel_ = dqmStore_->book1D(name, name, 50, 40, 110);
491  meInvMassZ0Sel_->setAxisTitle("mass (GeV)", 1);
492 
493  name = "EECLT dicluster invariant mass high sel";
494  meInvMassHighSel_ = dqmStore_->book1D(name, name, 500, 110, 3000);
495  meInvMassHighSel_->setAxisTitle("mass (GeV)", 1);
496 
497  }
498 
499 }
500 
501 void EEClusterTask::cleanup(void){
502 
503  if ( ! init_ ) return;
504 
505  if ( dqmStore_ ) {
506  dqmStore_->setCurrentFolder(prefixME_ + "/EEClusterTask");
507 
508  if ( meBCEne_ ) dqmStore_->removeElement( meBCEne_->getName() );
509  meBCEne_ = 0;
510 
511  if ( meBCNum_ ) dqmStore_->removeElement( meBCNum_->getName() );
512  meBCNum_ = 0;
513 
514  if ( meBCSiz_ ) dqmStore_->removeElement( meBCSiz_->getName() );
515  meBCSiz_ = 0;
516 
517  if ( meBCEneFwdMap_ ) dqmStore_->removeElement( meBCEneFwdMap_->getName() );
518  meBCEneFwdMap_ = 0;
519 
520  if ( meBCNumFwdMap_ ) dqmStore_->removeElement( meBCNumFwdMap_->getName() );
521  meBCNumFwdMap_ = 0;
522 
523  if ( meBCETFwdMap_ ) dqmStore_->removeElement( meBCETFwdMap_->getName() );
524  meBCETFwdMap_ = 0;
525 
526  if ( meBCSizFwdMap_ ) dqmStore_->removeElement( meBCSizFwdMap_->getName() );
527  meBCSizFwdMap_ = 0;
528 
529  if ( meBCEneFwdMapProjEta_ ) dqmStore_->removeElement( meBCEneFwdMapProjEta_->getName() );
530  meBCEneFwdMapProjEta_ = 0;
531 
532  if ( meBCEneFwdMapProjPhi_ ) dqmStore_->removeElement( meBCEneFwdMapProjPhi_->getName() );
533  meBCEneFwdMapProjPhi_ = 0;
534 
535  if ( meBCNumFwdMapProjEta_ ) dqmStore_->removeElement( meBCNumFwdMapProjEta_->getName() );
536  meBCNumFwdMapProjEta_ = 0;
537 
538  if ( meBCNumFwdMapProjPhi_ ) dqmStore_->removeElement( meBCNumFwdMapProjPhi_->getName() );
539  meBCNumFwdMapProjPhi_ = 0;
540 
541  if ( meBCETFwdMapProjEta_ ) dqmStore_->removeElement( meBCETFwdMapProjEta_->getName() );
542  meBCETFwdMapProjEta_ = 0;
543 
544  if ( meBCETFwdMapProjPhi_ ) dqmStore_->removeElement( meBCETFwdMapProjPhi_->getName() );
545  meBCETFwdMapProjPhi_ = 0;
546 
547  if ( meBCSizFwdMapProjEta_ ) dqmStore_->removeElement( meBCSizFwdMapProjEta_->getName() );
548  meBCSizFwdMapProjEta_ = 0;
549 
550  if ( meBCSizFwdMapProjPhi_ ) dqmStore_->removeElement( meBCSizFwdMapProjPhi_->getName() );
551  meBCSizFwdMapProjPhi_ = 0;
552 
553  if ( meBCEneBwdMap_ ) dqmStore_->removeElement( meBCEneBwdMap_->getName() );
554  meBCEneBwdMap_ = 0;
555 
556  if ( meBCNumBwdMap_ ) dqmStore_->removeElement( meBCNumBwdMap_->getName() );
557  meBCNumBwdMap_ = 0;
558 
559  if ( meBCETBwdMap_ ) dqmStore_->removeElement( meBCETBwdMap_->getName() );
560  meBCETBwdMap_ = 0;
561 
562  if ( meBCSizBwdMap_ ) dqmStore_->removeElement( meBCSizBwdMap_->getName() );
563  meBCSizBwdMap_ = 0;
564 
565  if ( meBCEneBwdMapProjEta_ ) dqmStore_->removeElement( meBCEneBwdMapProjEta_->getName() );
566  meBCEneBwdMapProjEta_ = 0;
567 
568  if ( meBCEneBwdMapProjPhi_ ) dqmStore_->removeElement( meBCEneBwdMapProjPhi_->getName() );
569  meBCEneBwdMapProjPhi_ = 0;
570 
571  if ( meBCNumBwdMapProjEta_ ) dqmStore_->removeElement( meBCNumBwdMapProjEta_->getName() );
572  meBCNumBwdMapProjEta_ = 0;
573 
574  if ( meBCNumBwdMapProjPhi_ ) dqmStore_->removeElement( meBCNumBwdMapProjPhi_->getName() );
575  meBCNumBwdMapProjPhi_ = 0;
576 
577  if ( meBCETBwdMapProjEta_ ) dqmStore_->removeElement( meBCETBwdMapProjEta_->getName() );
578  meBCETBwdMapProjEta_ = 0;
579 
580  if ( meBCETBwdMapProjPhi_ ) dqmStore_->removeElement( meBCETBwdMapProjPhi_->getName() );
581  meBCETBwdMapProjPhi_ = 0;
582 
583  if ( meBCSizBwdMapProjEta_ ) dqmStore_->removeElement( meBCSizBwdMapProjEta_->getName() );
584  meBCSizBwdMapProjEta_ = 0;
585 
586  if ( meBCSizBwdMapProjPhi_ ) dqmStore_->removeElement( meBCSizBwdMapProjPhi_->getName() );
587  meBCSizBwdMapProjPhi_ = 0;
588 
589  if ( meSCEne_ ) dqmStore_->removeElement( meSCEne_->getName() );
590  meSCEne_ = 0;
591 
592  if ( meSCNum_ ) dqmStore_->removeElement( meSCNum_->getName() );
593  meSCNum_ = 0;
594 
595  if ( meSCSiz_ ) dqmStore_->removeElement( meSCSiz_->getName() );
596  meSCSiz_ = 0;
597 
598  if ( meSCCrystalSiz_ ) dqmStore_->removeElement( meSCCrystalSiz_->getName() );
599  meSCCrystalSiz_ = 0;
600 
601  if ( meSCSeedEne_ ) dqmStore_->removeElement( meSCSeedEne_->getName() );
602  meSCSeedEne_ = 0;
603 
604  if ( meSCEne2_ ) dqmStore_->removeElement( meSCEne2_->getName() );
605  meSCEne2_ = 0;
606 
607  if ( meSCEneVsEMax_ ) dqmStore_->removeElement( meSCEneVsEMax_->getName() );
608  meSCEneVsEMax_ = 0;
609 
610  if ( meSCEneLowScale_ ) dqmStore_->removeElement( meSCEneLowScale_->getName() );
611  meSCEneLowScale_ = 0;
612 
613  if ( meSCSeedMapOcc_[0] ) dqmStore_->removeElement( meSCSeedMapOcc_[0]->getName() );
614  meSCSeedMapOcc_[0] = 0;
615 
616  if ( meSCSeedMapOcc_[1] ) dqmStore_->removeElement( meSCSeedMapOcc_[1]->getName() );
617  meSCSeedMapOcc_[1] = 0;
618 
619  if ( meSCMapSingleCrystal_[0] ) dqmStore_->removeElement( meSCMapSingleCrystal_[0]->getName() );
620  meSCMapSingleCrystal_[0] = 0;
621 
622  if ( meSCMapSingleCrystal_[1] ) dqmStore_->removeElement( meSCMapSingleCrystal_[1]->getName() );
623  meSCMapSingleCrystal_[1] = 0;
624 
625  if ( mes1s9_ ) dqmStore_->removeElement( mes1s9_->getName() );
626  mes1s9_ = 0;
627 
628  if ( mes1s9thr_ ) dqmStore_->removeElement( mes1s9thr_->getName() );
629  mes1s9thr_ = 0;
630 
631  if ( mes9s25_ ) dqmStore_->removeElement( mes9s25_->getName() );
632  mes9s25_ = 0;
633 
634  if ( meInvMassPi0_ ) dqmStore_->removeElement( meInvMassPi0_->getName() );
635  meInvMassPi0_ = 0;
636 
637  if ( meInvMassJPsi_ ) dqmStore_->removeElement( meInvMassJPsi_->getName() );
638  meInvMassJPsi_ = 0;
639 
640  if ( meInvMassZ0_ ) dqmStore_->removeElement( meInvMassZ0_->getName() );
641  meInvMassZ0_ = 0;
642 
643  if ( meInvMassHigh_ ) dqmStore_->removeElement( meInvMassHigh_->getName() );
644  meInvMassHigh_ = 0;
645 
646  if ( meInvMassPi0Sel_ ) dqmStore_->removeElement( meInvMassPi0Sel_->getName() );
647  meInvMassPi0Sel_ = 0;
648 
649  if ( meInvMassJPsiSel_ ) dqmStore_->removeElement( meInvMassJPsiSel_->getName() );
650  meInvMassJPsiSel_ = 0;
651 
652  if ( meInvMassZ0Sel_ ) dqmStore_->removeElement( meInvMassZ0Sel_->getName() );
653  meInvMassZ0Sel_ = 0;
654 
655  if ( meInvMassHighSel_ ) dqmStore_->removeElement( meInvMassHighSel_->getName() );
656  meInvMassHighSel_ = 0;
657 
658  }
659 
660  init_ = false;
661 
662 }
663 
664 void EEClusterTask::endJob(void){
665 
666  edm::LogInfo("EEClusterTask") << "analyzed " << ievt_ << " events";
667 
668  if ( enableCleanup_ ) this->cleanup();
669 
670 }
671 
672 void EEClusterTask::analyze(const edm::Event& e, const edm::EventSetup& c){
673 
674  bool enable = false;
675 
676  edm::Handle<EcalRawDataCollection> dcchs;
677 
678  if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
679 
680  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
681 
682  if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
683 
684  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH4 ||
685  dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH2 ) enable = true;
686 
687  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::COSMIC ||
688  dcchItr->getRunType() == EcalDCCHeaderBlock::MTCC ||
689  dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
690  dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
691  dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_LOCAL ||
692  dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_LOCAL ) enable = true;
693 
694  break;
695 
696  }
697 
698  } else {
699 
700  enable = true;
701  edm::LogWarning("EEClusterTask") << EcalRawDataCollection_ << " not available";
702 
703  }
704 
705  if ( ! enable ) return;
706 
707  if ( ! init_ ) this->setup();
708 
709  ievt_++;
710 
711  // ECAL topology
712  edm::ESHandle<CaloTopology> pTopology;
713  c.get<CaloTopologyRecord>().get(pTopology);
714  if ( !pTopology.isValid() ) {
715  edm::LogWarning("EEClusterTask") << "Topology not valid";
716  return;
717  }
718  const CaloTopology* topology = pTopology.product();
719 
720  // recHits
721  edm::Handle< EcalRecHitCollection > pEERecHits;
722  e.getByLabel( EcalRecHitCollection_, pEERecHits );
723  if ( !pEERecHits.isValid() ) {
724  edm::LogWarning("EEClusterTask") << "RecHit collection " << EcalRecHitCollection_ << " not available.";
725  return;
726  }
727  const EcalRecHitCollection* eeRecHits = pEERecHits.product();
728 
729  reco::BasicClusterCollection bcSel;
730 
731  // --- Endcap Basic Clusters ---
732  edm::Handle<reco::BasicClusterCollection> pBasicClusters;
733  if ( e.getByLabel(BasicClusterCollection_, pBasicClusters) ) {
734 
735  int nbcc = pBasicClusters->size();
736  if (nbcc>0) meBCNum_->Fill(float(nbcc));
737 
738  for ( reco::BasicClusterCollection::const_iterator bCluster = pBasicClusters->begin(); bCluster != pBasicClusters->end(); ++bCluster ) {
739 
740  meBCEne_->Fill(bCluster->energy());
741  meBCSiz_->Fill(float(bCluster->size()));
742 
743  if ( bCluster->eta() > 0 ) {
744  meBCEneFwdMap_->Fill(bCluster->x(), bCluster->y(), bCluster->energy());
745  meBCEneFwdMapProjEta_->Fill( bCluster->eta(), bCluster->energy() );
746  meBCEneFwdMapProjPhi_->Fill( bCluster->phi(), bCluster->energy() );
747 
748  meBCNumFwdMap_->Fill(bCluster->x(), bCluster->y());
749  meBCNumFwdMapProjEta_->Fill( bCluster->eta() );
750  meBCNumFwdMapProjPhi_->Fill( bCluster->phi() );
751 
752  meBCETFwdMap_->Fill(bCluster->x(), bCluster->y(), bCluster->energy() * sin(bCluster->position().theta()) );
753  meBCETFwdMapProjEta_->Fill( bCluster->eta(), bCluster->energy() * sin(bCluster->position().theta()) );
754  meBCETFwdMapProjPhi_->Fill( bCluster->phi(), bCluster->energy() * sin(bCluster->position().theta()) );
755 
756  meBCSizFwdMap_->Fill(bCluster->x(), bCluster->y(), float(bCluster->size()) );
757  meBCSizFwdMapProjEta_->Fill( bCluster->eta(), float(bCluster->size()) );
758  meBCSizFwdMapProjPhi_->Fill( bCluster->phi(), float(bCluster->size()) );
759  } else {
760  meBCEneBwdMap_->Fill(bCluster->x(), bCluster->y(), bCluster->energy());
761  meBCEneBwdMapProjEta_->Fill( bCluster->eta(), bCluster->energy() );
762  meBCEneBwdMapProjPhi_->Fill( bCluster->phi(), bCluster->energy() );
763 
764  meBCNumBwdMap_->Fill(bCluster->x(), bCluster->y());
765  meBCNumBwdMapProjEta_->Fill( bCluster->eta() );
766  meBCNumBwdMapProjPhi_->Fill( bCluster->phi() );
767 
768  meBCETBwdMap_->Fill(bCluster->x(), bCluster->y(), bCluster->energy() * sin(bCluster->position().theta()) );
769  meBCETBwdMapProjEta_->Fill( bCluster->eta(), bCluster->energy() * sin(bCluster->position().theta()) );
770  meBCETBwdMapProjPhi_->Fill( bCluster->phi(), bCluster->energy() * sin(bCluster->position().theta()) );
771 
772  meBCSizBwdMap_->Fill(bCluster->x(), bCluster->y(), float(bCluster->size()) );
773  meBCSizBwdMapProjEta_->Fill( bCluster->eta(), float(bCluster->size()) );
774  meBCSizBwdMapProjPhi_->Fill( bCluster->phi(), float(bCluster->size()) );
775 
776  float e2x2 = EcalClusterTools::e2x2( *bCluster, eeRecHits, topology );
777  float e3x3 = EcalClusterTools::e3x3( *bCluster, eeRecHits, topology );
778 
779  // fill the selected cluster collection
780  float pt = std::abs( bCluster->energy()*sin(bCluster->position().theta()) );
781  if ( pt > thrClusEt_ && e2x2/e3x3 > thrS4S9_ ) bcSel.push_back(*bCluster);
782  }
783 
784  }
785 
786  } else {
787 
788  edm::LogWarning("EEClusterTask") << BasicClusterCollection_ << " not available";
789 
790  }
791 
792  for ( reco::BasicClusterCollection::const_iterator bc1 = bcSel.begin(); bc1 != bcSel.end(); ++bc1 ) {
793  TLorentzVector bc1P;
794  bc1P.SetPtEtaPhiE(std::abs(bc1->energy()*sin(bc1->position().theta())),
795  bc1->eta(), bc1->phi(), bc1->energy());
796  for ( reco::BasicClusterCollection::const_iterator bc2 = bc1+1; bc2 != bcSel.end(); ++bc2 ) {
797  TLorentzVector bc2P;
798  bc2P.SetPtEtaPhiE(std::abs(bc2->energy()*sin(bc2->position().theta())),
799  bc2->eta(), bc2->phi(), bc2->energy());
800 
801  TLorentzVector candP = bc1P + bc2P;
802 
803  if ( candP.Pt() > thrCandEt_ ) {
804  float mass = candP.M();
805  if ( mass < 0.500 ) {
806  meInvMassPi0Sel_->Fill( mass );
807  } else if ( mass > 2.9 && mass < 3.3 ) {
808  meInvMassJPsiSel_->Fill( mass );
809  } else if ( mass > 40 && mass < 110 ) {
810  meInvMassZ0Sel_->Fill( mass );
811  } else if ( mass > 110 ) {
812  meInvMassHighSel_->Fill( mass );
813  }
814 
815  }
816 
817  }
818  }
819 
820  // --- Endcap Super Clusters ----
821  edm::Handle<reco::SuperClusterCollection> pSuperClusters;
822  if ( e.getByLabel(SuperClusterCollection_, pSuperClusters) ) {
823 
824  int nscc = pSuperClusters->size();
825  if ( nscc > 0 ) meSCNum_->Fill(float(nscc));
826 
827  TLorentzVector sc1_p(0,0,0,0);
828  TLorentzVector sc2_p(0,0,0,0);
829 
830  reco::SuperClusterCollection scSel;
831 
832  for ( reco::SuperClusterCollection::const_iterator sCluster = pSuperClusters->begin(); sCluster != pSuperClusters->end(); sCluster++ ) {
833 
834  // energy, size
835  meSCEne_->Fill(sCluster->energy());
836  meSCSiz_->Fill(float(sCluster->clustersSize()));
837 
838  reco::CaloClusterPtr theSeed = sCluster->seed();
839 
840  // Find the seed rec hit
841  std::vector< std::pair<DetId,float> > sIds = sCluster->hitsAndFractions();
842 
843  float eMax, e2nd;
844  EcalRecHitCollection::const_iterator seedItr = eeRecHits->begin();
845  EcalRecHitCollection::const_iterator secondItr = eeRecHits->begin();
846 
847  for(std::vector< std::pair<DetId,float> >::const_iterator idItr = sIds.begin(); idItr != sIds.end(); ++idItr) {
848  DetId id = idItr->first;
849  if(id.det() != DetId::Ecal) { continue; }
850  EcalRecHitCollection::const_iterator hitItr = eeRecHits->find(id);
851  if(hitItr == eeRecHits->end()) { continue; }
852  if(hitItr->energy() > secondItr->energy()) { secondItr = hitItr; }
853  if(hitItr->energy() > seedItr->energy()) { std::swap(seedItr,secondItr); }
854  }
855 
856  eMax = seedItr->energy();
857  e2nd = secondItr->energy();
858  EEDetId seedId = (EEDetId) seedItr->id();
859 
860  float e3x3 = EcalClusterTools::e3x3( *theSeed, eeRecHits, topology );
861  float e5x5 = EcalClusterTools::e5x5( *theSeed, eeRecHits, topology );
862 
863  meSCCrystalSiz_->Fill(sIds.size());
864  meSCSeedEne_->Fill(eMax);
865  meSCEne2_->Fill(eMax+e2nd);
866  meSCEneVsEMax_->Fill(eMax,sCluster->energy());
867  meSCEneLowScale_->Fill(sCluster->energy());
868 
869  // Prepare to fill maps
870  int ism = Numbers::iSM(seedId);
871  int eeSide;
872  if( ism >= 1 && ism <= 9)
873  eeSide = 0;
874  else
875  eeSide = 1;
876  int eex = seedId.ix();
877  int eey = seedId.iy();
878  float xeex = eex - 0.5;
879  float xeey = eey - 0.5;
880 
881  meSCSeedMapOcc_[eeSide]->Fill(xeex, xeey);
882 
883  if(sIds.size() == 1) meSCMapSingleCrystal_[eeSide]->Fill(xeex, xeey);
884 
885  mes1s9_->Fill( eMax/e3x3 );
886  if ( eMax > 3.0 ) mes1s9thr_->Fill( eMax/e3x3 );
887  mes9s25_->Fill( e3x3/e5x5 );
888 
889  // look for the two most energetic super clusters
890  if ( sCluster->energy() > sc1_p.Energy() ) {
891  sc2_p=sc1_p;
892  sc1_p.SetPtEtaPhiE(sCluster->energy()*sin(sCluster->position().theta()),
893  sCluster->eta(), sCluster->phi(), sCluster->energy());
894  } else if ( sCluster->energy() > sc2_p.Energy() ) {
895  sc2_p.SetPtEtaPhiE(sCluster->energy()*sin(sCluster->position().theta()),
896  sCluster->eta(), sCluster->phi(), sCluster->energy());
897  }
898 
899  }
900  // Get the invariant mass of the two most energetic super clusters
901  if ( nscc >= 2) {
902  TLorentzVector sum = sc1_p+sc2_p;
903  float mass = sum.M();
904  if ( mass < 0.500 ) {
905  meInvMassPi0_->Fill( mass );
906  } else if ( mass > 2.9 && mass < 3.3 ) {
907  meInvMassJPsi_->Fill( mass );
908  } else if ( mass > 40 && mass < 110 ) {
909  meInvMassZ0_->Fill( mass );
910  } else if ( mass > 110 ) {
911  meInvMassHigh_->Fill( mass );
912  }
913  }
914 
915  } else {
916 
917  edm::LogWarning("EEClusterTask") << SuperClusterCollection_ << " not available";
918 
919  }
920 
921 }
EEClusterTask::meBCETFwdMap_
MonitorElement * meBCETFwdMap_
Definition: EEClusterTask.h:77
EEClusterTask::mes9s25_
MonitorElement * mes9s25_
Definition: EEClusterTask.h:99
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:99
EEClusterTask::meBCETBwdMapProjEta_
MonitorElement * meBCETBwdMapProjEta_
Definition: EEClusterTask.h:82
EEClusterTask::reset
void reset(void)
Reset.
Definition: EEClusterTask.cc:155
EEClusterTask::endJob
void endJob(void)
EndJob.
Definition: EEClusterTask.cc:664
EEDetId::ix
int ix() const
Definition: EEDetId.h:71
EEClusterTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EEClusterTask.cc:672
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:717
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:2530
EEClusterTask::meBCNum_
MonitorElement * meBCNum_
Definition: EEClusterTask.h:74
Numbers.h
Some &quot;id&quot; conversions.
EEClusterTask::meBCNumFwdMap_
MonitorElement * meBCNumFwdMap_
Definition: EEClusterTask.h:77
EEClusterTask::meInvMassPi0_
MonitorElement * meInvMassPi0_
Definition: EEClusterTask.h:100
EEClusterTask::meInvMassHigh_
MonitorElement * meInvMassHigh_
Definition: EEClusterTask.h:103
EEClusterTask::beginJob
void beginJob(void)
BeginJob.
Definition: EEClusterTask.cc:132
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
abs
#define abs(x)
Definition: mlp_lapack.h:159
EEClusterTask::meBCEneBwdMap_
MonitorElement * meBCEneBwdMap_
Definition: EEClusterTask.h:81
EEClusterTask::meBCSizFwdMapProjEta_
MonitorElement * meBCSizFwdMapProjEta_
Definition: EEClusterTask.h:78
EEClusterTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EEClusterTask.cc:143
EEClusterTask::cleanup
void cleanup(void)
Cleanup.
Definition: EEClusterTask.cc:501
EEClusterTask::EEClusterTask
EEClusterTask(const edm::ParameterSet &ps)
Constructor.
Definition: EEClusterTask.cc:43
EEClusterTask::meSCNum_
MonitorElement * meSCNum_
Definition: EEClusterTask.h:86
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
EEClusterTask::meBCNumBwdMapProjPhi_
MonitorElement * meBCNumBwdMapProjPhi_
Definition: EEClusterTask.h:83
EEClusterTask::meBCETBwdMapProjPhi_
MonitorElement * meBCETBwdMapProjPhi_
Definition: EEClusterTask.h:83
EEClusterTask::meBCEneFwdMapProjPhi_
MonitorElement * meBCEneFwdMapProjPhi_
Definition: EEClusterTask.h:79
EEClusterTask::meBCEneFwdMapProjEta_
MonitorElement * meBCEneFwdMapProjEta_
Definition: EEClusterTask.h:78
EEClusterTask::SuperClusterCollection_
edm::InputTag SuperClusterCollection_
Definition: EEClusterTask.h:70
reco::SuperClusterCollection
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
Definition: SuperClusterFwd.h:9
EEClusterTask::meBCNumFwdMapProjPhi_
MonitorElement * meBCNumFwdMapProjPhi_
Definition: EEClusterTask.h:79
EEClusterTask::meBCNumBwdMapProjEta_
MonitorElement * meBCNumBwdMapProjEta_
Definition: EEClusterTask.h:82
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
std::swap
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
Definition: DataFrameContainer.h:223
EEClusterTask::meInvMassZ0Sel_
MonitorElement * meInvMassZ0Sel_
Definition: EEClusterTask.h:106
EEClusterTask::meBCEneBwdMapProjEta_
MonitorElement * meBCEneBwdMapProjEta_
Definition: EEClusterTask.h:82
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2572
EEClusterTask::meBCSizBwdMap_
MonitorElement * meBCSizBwdMap_
Definition: EEClusterTask.h:81
EcalClusterTools::e3x3
static float e3x3(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.cc:240
EEClusterTask::meBCEneBwdMapProjPhi_
MonitorElement * meBCEneBwdMapProjPhi_
Definition: EEClusterTask.h:83
EEClusterTask::meInvMassJPsiSel_
MonitorElement * meInvMassJPsiSel_
Definition: EEClusterTask.h:105
EEClusterTask::meBCSizFwdMapProjPhi_
MonitorElement * meBCSizFwdMapProjPhi_
Definition: EEClusterTask.h:79
EEClusterTask::meBCETFwdMapProjPhi_
MonitorElement * meBCETFwdMapProjPhi_
Definition: EEClusterTask.h:79
EEClusterTask::meBCSiz_
MonitorElement * meBCSiz_
Definition: EEClusterTask.h:75
EEDetId::iy
int iy() const
Definition: EEDetId.h:77
EEClusterTask::meBCNumFwdMapProjEta_
MonitorElement * meBCNumFwdMapProjEta_
Definition: EEClusterTask.h:78
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1031
edm::getName
std::string getName(Reflex::Type &cc)
Definition: ClassFiller.cc:18
EEClusterTask::meSCEneVsEMax_
MonitorElement * meSCEneVsEMax_
Definition: EEClusterTask.h:92
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
EEClusterTask::meSCSiz_
MonitorElement * meSCSiz_
Definition: EEClusterTask.h:87
EEClusterTask::meBCSizFwdMap_
MonitorElement * meBCSizFwdMap_
Definition: EEClusterTask.h:77
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
EEClusterTask::meSCEne_
MonitorElement * meSCEne_
Definition: EEClusterTask.h:85
EEClusterTask::meSCMapSingleCrystal_
MonitorElement * meSCMapSingleCrystal_[2]
Definition: EEClusterTask.h:95
EEClusterTask::meSCSeedEne_
MonitorElement * meSCSeedEne_
Definition: EEClusterTask.h:90
jptDQMConfig_cff.eMax
int eMax
Definition: jptDQMConfig_cff.py:21
EEClusterTask::meInvMassJPsi_
MonitorElement * meInvMassJPsi_
Definition: EEClusterTask.h:101
EEClusterTask::meBCETFwdMapProjEta_
MonitorElement * meBCETFwdMapProjEta_
Definition: EEClusterTask.h:78
EEClusterTask::meSCEneLowScale_
MonitorElement * meSCEneLowScale_
Definition: EEClusterTask.h:93
EEClusterTask::setup
void setup(void)
Setup.
Definition: EEClusterTask.cc:259
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:301
EEClusterTask::EcalRecHitCollection_
edm::InputTag EcalRecHitCollection_
Definition: EEClusterTask.h:71
DetId
Definition: DetId.h:20
EEClusterTask::meSCCrystalSiz_
MonitorElement * meSCCrystalSiz_
Definition: EEClusterTask.h:89
M_PI
#define M_PI
Definition: BFit3D.cc:3
EEClusterTask::meBCNumBwdMap_
MonitorElement * meBCNumBwdMap_
Definition: EEClusterTask.h:81
EEClusterTask::meBCSizBwdMapProjPhi_
MonitorElement * meBCSizBwdMapProjPhi_
Definition: EEClusterTask.h:83
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:50
EEClusterTask::meInvMassPi0Sel_
MonitorElement * meInvMassPi0Sel_
Definition: EEClusterTask.h:104
EEClusterTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EEClusterTask.cc:151
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
EEClusterTask::meBCEneFwdMap_
MonitorElement * meBCEneFwdMap_
Definition: EEClusterTask.h:77
reco::BasicClusterCollection
std::vector< BasicCluster > BasicClusterCollection
collection of BasicCluster objects
Definition: BasicClusterFwd.h:27
edm::Handle::product
T const * product() const
Definition: Handle.h:74
EEClusterTask::meSCSeedMapOcc_
MonitorElement * meSCSeedMapOcc_[2]
Definition: EEClusterTask.h:94
edm::SortedCollection::find
iterator find(key_type k)
Definition: SortedCollection.h:262
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:246
scaleCards.mass
tuple mass
Definition: scaleCards.py:27
EcalClusterTools::e5x5
static float e5x5(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.cc:278
EEClusterTask::BasicClusterCollection_
edm::InputTag BasicClusterCollection_
Definition: EEClusterTask.h:69
EEClusterTask::meBCEne_
MonitorElement * meBCEne_
Definition: EEClusterTask.h:73
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:145
EEClusterTask::meBCETBwdMap_
MonitorElement * meBCETBwdMap_
Definition: EEClusterTask.h:81
EEClusterTask::meSCEne2_
MonitorElement * meSCEne2_
Definition: EEClusterTask.h:91
EEClusterTask::meInvMassHighSel_
MonitorElement * meInvMassHighSel_
Definition: EEClusterTask.h:107
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:37
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:845
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:916
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:422
EEClusterTask::EcalRawDataCollection_
edm::InputTag EcalRawDataCollection_
Definition: EEClusterTask.h:68
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:294
EEClusterTask::prefixME_
std::string prefixME_
Definition: EEClusterTask.h:62
edm::Run
Definition: Run.h:33
EEClusterTask::dqmStore_
DQMStore * dqmStore_
Definition: EEClusterTask.h:60
EEClusterTask::~EEClusterTask
virtual ~EEClusterTask()
Destructor.
Definition: EEClusterTask.cc:128
EEClusterTask::mes1s9_
MonitorElement * mes1s9_
Definition: EEClusterTask.h:97
EEClusterTask::mes1s9thr_
MonitorElement * mes1s9thr_
Definition: EEClusterTask.h:98
EEClusterTask::meInvMassZ0_
MonitorElement * meInvMassZ0_
Definition: EEClusterTask.h:102
DQMStore::bookProfile2D
MonitorElement * bookProfile2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, int nchZ, double lowZ, double highZ, const char *option="s")
Definition: DQMStore.cc:1175
EcalClusterTools::e2x2
static float e2x2(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.cc:188
EEClusterTask::meBCSizBwdMapProjEta_
MonitorElement * meBCSizBwdMapProjEta_
Definition: EEClusterTask.h:82
Generated for CMSSW Reference Manual by doxygen 1.8.5