CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EBClusterTask.cc
Go to the documentation of this file.
1 /*
2  * \file EBClusterTask.cc
3  *
4  * $Date: 2012/04/27 13:46:01 $
5  * $Revision: 1.95 $
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 "CondFormats/EcalObjects/interface/EcalADCToGeVConstant.h"
34 #include "CondFormats/DataRecord/interface/EcalADCToGeVConstantRcd.h"
35 
36 #include "DQM/EcalCommon/interface/Numbers.h"
37 
38 #include "DQM/EcalBarrelMonitorTasks/interface/EBClusterTask.h"
39 
40 #include "DataFormats/EcalDetId/interface/EBDetId.h"
41 
42 #include "TLorentzVector.h"
43 
44 EBClusterTask::EBClusterTask(const edm::ParameterSet& ps){
45 
46  init_ = false;
47 
48  dqmStore_ = edm::Service<DQMStore>().operator->();
49 
50  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
51 
52  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
53 
54  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
55 
56  // parameters...
57  EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
58  BasicClusterCollection_ = ps.getParameter<edm::InputTag>("BasicClusterCollection");
59  SuperClusterCollection_ = ps.getParameter<edm::InputTag>("SuperClusterCollection");
60  EcalRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalRecHitCollection");
61 
62  // histograms...
63  meBCEne_ = 0;
64  meBCNum_ = 0;
65  meBCSiz_ = 0;
66 
67  meBCEneMap_ = 0;
68  meBCNumMap_ = 0;
69  meBCETMap_ = 0;
70  meBCSizMap_ = 0;
71 
72  meBCEneMapProjEta_ = 0;
73  meBCNumMapProjEta_ = 0;
74  meBCETMapProjEta_ = 0;
75  meBCSizMapProjEta_ = 0;
76 
77  meBCEneMapProjPhi_ = 0;
78  meBCNumMapProjPhi_ = 0;
79  meBCETMapProjPhi_ = 0;
80  meBCSizMapProjPhi_ = 0;
81 
82  meSCEne_ = 0;
83  meSCNum_ = 0;
84  meSCSiz_ = 0;
85 
86  meSCCrystalSiz_ = 0;
87  meSCSeedEne_ = 0;
88  meSCEne2_ = 0;
89  meSCEneVsEMax_ = 0;
90  meSCEneLowScale_ = 0;
91  meSCSeedMapOcc_ = 0;
92  meSCMapSingleCrystal_ = 0;
93 
94  mes1s9_ = 0;
95  mes1s9thr_ = 0;
96  mes9s25_ = 0;
97 
98  meInvMassPi0_ = 0;
99  meInvMassJPsi_ = 0;
100  meInvMassZ0_ = 0;
101  meInvMassHigh_ = 0;
102 
103  meInvMassPi0Sel_ = 0;
104  meInvMassJPsiSel_ = 0;
105  meInvMassZ0Sel_ = 0;
106  meInvMassHighSel_ = 0;
107 
108  thrS4S9_ = 0.85;
109  thrClusEt_ = 0.200;
110  thrCandEt_ = 0.650;
111 
112 }
113 
114 EBClusterTask::~EBClusterTask(){
115 
116 }
117 
118 void EBClusterTask::beginJob(void){
119 
120  ievt_ = 0;
121 
122  if ( dqmStore_ ) {
123  dqmStore_->setCurrentFolder(prefixME_ + "/EBClusterTask");
124  dqmStore_->rmdir(prefixME_ + "/EBClusterTask");
125  }
126 
127 }
128 
129 void EBClusterTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
130 
131  Numbers::initGeometry(c, false);
132 
133  if ( ! mergeRuns_ ) this->reset();
134 
135 }
136 
137 void EBClusterTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
138 
139 }
140 
141 void EBClusterTask::reset(void) {
142 
143  if ( meBCEne_ ) meBCEne_->Reset();
144 
145  if ( meBCNum_ ) meBCNum_->Reset();
146 
147  if ( meBCSiz_ ) meBCSiz_->Reset();
148 
149  if ( meBCEneMap_ ) meBCEneMap_->Reset();
150 
151  if ( meBCNumMap_ ) meBCNumMap_->Reset();
152 
153  if ( meBCETMap_ ) meBCETMap_->Reset();
154 
155  if ( meBCSizMap_ ) meBCSizMap_->Reset();
156 
157  if ( meBCEneMapProjEta_ ) meBCEneMapProjEta_->Reset();
158 
159  if ( meBCEneMapProjPhi_ ) meBCEneMapProjPhi_->Reset();
160 
161  if ( meBCNumMapProjEta_ ) meBCNumMapProjEta_->Reset();
162 
163  if ( meBCNumMapProjPhi_ ) meBCNumMapProjPhi_->Reset();
164 
165  if ( meBCETMapProjEta_ ) meBCETMapProjEta_->Reset();
166 
167  if ( meBCETMapProjPhi_ ) meBCETMapProjPhi_->Reset();
168 
169  if ( meBCSizMapProjEta_ ) meBCSizMapProjEta_->Reset();
170 
171  if ( meBCSizMapProjPhi_ ) meBCSizMapProjPhi_->Reset();
172 
173  if ( meSCEne_ ) meSCEne_->Reset();
174 
175  if ( meSCNum_ ) meSCNum_->Reset();
176 
177  if ( meSCSiz_ ) meSCSiz_->Reset();
178 
179  if ( meSCCrystalSiz_ ) meSCCrystalSiz_->Reset();
180 
181  if ( meSCSeedEne_ ) meSCSeedEne_->Reset();
182 
183  if ( meSCEne2_ ) meSCEne2_->Reset();
184 
185  if ( meSCEneVsEMax_ ) meSCEneVsEMax_->Reset();
186 
187  if ( meSCEneLowScale_ ) meSCEneLowScale_->Reset();
188 
189  if ( meSCSeedMapOcc_ ) meSCSeedMapOcc_->Reset();
190 
191  if ( meSCMapSingleCrystal_ ) meSCMapSingleCrystal_->Reset();
192 
193  if ( mes1s9_ ) mes1s9_->Reset();
194 
195  if ( mes1s9thr_ ) mes1s9thr_->Reset();
196 
197  if ( mes9s25_ ) mes9s25_->Reset();
198 
199  if ( meInvMassPi0_ ) meInvMassPi0_->Reset();
200 
201  if ( meInvMassJPsi_ ) meInvMassJPsi_->Reset();
202 
203  if ( meInvMassZ0_ ) meInvMassZ0_->Reset();
204 
205  if ( meInvMassHigh_ ) meInvMassHigh_->Reset();
206 
207  if ( meInvMassPi0Sel_ ) meInvMassPi0Sel_->Reset();
208 
209  if ( meInvMassJPsiSel_ ) meInvMassJPsiSel_->Reset();
210 
211  if ( meInvMassZ0Sel_ ) meInvMassZ0Sel_->Reset();
212 
213  if ( meInvMassHighSel_ ) meInvMassHighSel_->Reset();
214 
215 }
216 
217 void EBClusterTask::setup(void){
218 
219  init_ = true;
220 
221  std::string name;
222 
223  if ( dqmStore_ ) {
224  dqmStore_->setCurrentFolder(prefixME_ + "/EBClusterTask");
225 
226  name = "EBCLT BC energy";
227  meBCEne_ = dqmStore_->book1D(name, name, 100, 0., 150.);
228  meBCEne_->setAxisTitle("energy (GeV)", 1);
229 
230  name = "EBCLT BC number";
231  meBCNum_ = dqmStore_->book1D(name, name, 100, 0., 100.);
232  meBCNum_->setAxisTitle("number of clusters", 1);
233 
234  name = "EBCLT BC size";
235  meBCSiz_ = dqmStore_->book1D(name, name, 100, 0., 100.);
236  meBCSiz_->setAxisTitle("cluster size", 1);
237 
238  name = "EBCLT BC energy map";
239  meBCEneMap_ = dqmStore_->bookProfile2D(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 500., "s");
240  meBCEneMap_->setAxisTitle("phi", 1);
241  meBCEneMap_->setAxisTitle("eta", 2);
242 
243  name = "EBCLT BC number map";
244  meBCNumMap_ = dqmStore_->book2D(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479);
245  meBCNumMap_->setAxisTitle("phi", 1);
246  meBCNumMap_->setAxisTitle("eta", 2);
247 
248  name = "EBCLT BC ET map";
249  meBCETMap_ = dqmStore_->bookProfile2D(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 500., "s");
250  meBCETMap_->setAxisTitle("phi", 1);
251  meBCETMap_->setAxisTitle("eta", 2);
252 
253  name = "EBCLT BC size map";
254  meBCSizMap_ = dqmStore_->bookProfile2D(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 100., "s");
255  meBCSizMap_->setAxisTitle("phi", 1);
256  meBCSizMap_->setAxisTitle("eta", 2);
257 
258  name = "EBCLT BC energy projection eta";
259  meBCEneMapProjEta_ = dqmStore_->bookProfile(name, name, 34, -1.479, 1.479, 100, 0., 500., "s");
260  meBCEneMapProjEta_->setAxisTitle("eta", 1);
261  meBCEneMapProjEta_->setAxisTitle("energy (GeV)", 2);
262 
263  name = "EBCLT BC energy projection phi";
264  meBCEneMapProjPhi_ = dqmStore_->bookProfile(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 500., "s");
265  meBCEneMapProjPhi_->setAxisTitle("phi", 1);
266  meBCEneMapProjPhi_->setAxisTitle("energy (GeV)", 2);
267 
268  name = "EBCLT BC number projection eta";
269  meBCNumMapProjEta_ = dqmStore_->book1D(name, name, 34, -1.479, 1.479);
270  meBCNumMapProjEta_->setAxisTitle("eta", 1);
271  meBCNumMapProjEta_->setAxisTitle("number of clusters", 2);
272 
273  name = "EBCLT BC number projection phi";
274  meBCNumMapProjPhi_ = dqmStore_->book1D(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9);
275  meBCNumMapProjPhi_->setAxisTitle("phi", 1);
276  meBCNumMapProjPhi_->setAxisTitle("number of clusters", 2);
277 
278  name = "EBCLT BC ET projection eta";
279  meBCETMapProjEta_ = dqmStore_->bookProfile(name, name, 34, -1.479, 1.479, 100, 0., 500., "s");
280  meBCETMapProjEta_->setAxisTitle("eta", 1);
281  meBCETMapProjEta_->setAxisTitle("transverse energy (GeV)", 2);
282 
283  name = "EBCLT BC ET projection phi";
284  meBCETMapProjPhi_ = dqmStore_->bookProfile(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 500., "s");
285  meBCETMapProjPhi_->setAxisTitle("phi", 1);
286  meBCETMapProjPhi_->setAxisTitle("transverse energy (GeV)", 2);
287 
288  name = "EBCLT BC size projection eta";
289  meBCSizMapProjEta_ = dqmStore_->bookProfile(name, name, 34, -1.479, 1.479, 100, 0., 100., "s");
290  meBCSizMapProjEta_->setAxisTitle("eta", 1);
291  meBCSizMapProjEta_->setAxisTitle("cluster size", 2);
292 
293  name = "EBCLT BC size projection phi";
294  meBCSizMapProjPhi_ = dqmStore_->bookProfile(name, name, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 100., "s");
295  meBCSizMapProjPhi_->setAxisTitle("phi", 1);
296  meBCSizMapProjPhi_->setAxisTitle("cluster size", 2);
297 
298  name = "EBCLT SC energy";
299  meSCEne_ = dqmStore_->book1D(name, name, 100, 0., 150.);
300  meSCEne_->setAxisTitle("energy (GeV)", 1);
301 
302  name = "EBCLT SC number";
303  meSCNum_ = dqmStore_->book1D(name, name, 50, 0., 50.);
304  meSCNum_->setAxisTitle("number of clusters", 1);
305 
306  name = "EBCLT SC size";
307  meSCSiz_ = dqmStore_->book1D(name, name, 50, 0., 50.);
308  meSCSiz_->setAxisTitle("cluster size", 1);
309 
310  name = "EBCLT SC size (crystal)";
311  meSCCrystalSiz_ = dqmStore_->book1D(name, name, 150, 0, 150);
312  meSCCrystalSiz_->setAxisTitle("cluster size in crystals", 1);
313 
314  name = "EBCLT SC seed crystal energy";
315  meSCSeedEne_ = dqmStore_->book1D(name, name, 100, 0., 10.);
316  meSCSeedEne_->setAxisTitle("seed crystal energy (GeV)", 1);
317 
318  name = "EBCLT SC e2";
319  meSCEne2_ = dqmStore_->book1D(name, name, 100, 0., 10.);
320  meSCEne2_->setAxisTitle("seed + highest neighbor crystal energy (GeV)", 1);
321 
322  name = "EBCLT SC energy vs seed crystal energy";
323  meSCEneVsEMax_ = dqmStore_->book2D(name, name, 50, 0., 10., 50, 0., 10.);
324  meSCEneVsEMax_->setAxisTitle("seed crystal energy (GeV)", 1);
325  meSCEneVsEMax_->setAxisTitle("cluster energy (GeV)", 2);
326 
327  name = "EBCLT SC energy (low scale)";
328  meSCEneLowScale_ = dqmStore_->book1D(name, name, 200, 0., 10.);
329  meSCEneLowScale_->setAxisTitle("cluster energy (GeV)", 1);
330 
331  name = "EBCLT SC seed occupancy map";
332  meSCSeedMapOcc_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85, 85);
333  meSCSeedMapOcc_->setAxisTitle("jphi", 1);
334  meSCSeedMapOcc_->setAxisTitle("jeta", 2);
335 
336  name = "EBCLT SC single crystal cluster seed occupancy map";
337  meSCMapSingleCrystal_ = dqmStore_->book2D(name, name, 72, 0., 360., 34, -85, 85);
338  meSCMapSingleCrystal_->setAxisTitle("jphi", 1);
339  meSCMapSingleCrystal_->setAxisTitle("jeta", 2);
340 
341  name = "EBCLT s1s9";
342  mes1s9_ = dqmStore_->book1D(name, name, 50, 0., 1.5);
343  mes1s9_->setAxisTitle("s1/s9", 1);
344 
345  name = "EBCLT s1s9 thr";
346  mes1s9thr_ = dqmStore_->book1D(name, name, 50, 0., 1.5);
347  mes1s9thr_->setAxisTitle("s1/s9", 1);
348 
349  name = "EBCLT s9s25";
350  mes9s25_ = dqmStore_->book1D(name, name, 75, 0., 1.5);
351  mes9s25_->setAxisTitle("s9/s25", 1);
352 
353  name = "EBCLT dicluster invariant mass Pi0";
354  meInvMassPi0_ = dqmStore_->book1D(name, name, 50, 0.0, 0.500);
355  meInvMassPi0_->setAxisTitle("mass (GeV)", 1);
356 
357  name = "EBCLT dicluster invariant mass JPsi";
358  meInvMassJPsi_ = dqmStore_->book1D(name, name, 50, 2.9, 3.3);
359  meInvMassJPsi_->setAxisTitle("mass (GeV)", 1);
360 
361  name = "EBCLT dicluster invariant mass Z0";
362  meInvMassZ0_ = dqmStore_->book1D(name, name, 50, 40, 110);
363  meInvMassZ0_->setAxisTitle("mass (GeV)", 1);
364 
365  name = "EBCLT dicluster invariant mass high";
366  meInvMassHigh_ = dqmStore_->book1D(name, name, 500, 110, 3000);
367  meInvMassHigh_->setAxisTitle("mass (GeV)", 1);
368 
369  name = "EBCLT dicluster invariant mass Pi0 sel";
370  meInvMassPi0Sel_ = dqmStore_->book1D(name, name, 50, 0.00, 0.500);
371  meInvMassPi0Sel_->setAxisTitle("mass (GeV)", 1);
372 
373  name = "EBCLT dicluster invariant mass JPsi sel";
374  meInvMassJPsiSel_ = dqmStore_->book1D(name, name, 50, 2.9, 3.3);
375  meInvMassJPsiSel_->setAxisTitle("mass (GeV)", 1);
376 
377  name = "EBCLT dicluster invariant mass Z0 sel";
378  meInvMassZ0Sel_ = dqmStore_->book1D(name, name, 50, 40, 110);
379  meInvMassZ0Sel_->setAxisTitle("mass (GeV)", 1);
380 
381  name = "EBCLT dicluster invariant mass high sel";
382  meInvMassHighSel_ = dqmStore_->book1D(name, name, 500, 110, 3000);
383  meInvMassHighSel_->setAxisTitle("mass (GeV)", 1);
384 
385  }
386 
387 }
388 
389 void EBClusterTask::cleanup(void){
390 
391  if ( ! init_ ) return;
392 
393  if ( dqmStore_ ) {
394  dqmStore_->setCurrentFolder(prefixME_ + "/EBClusterTask");
395 
396  if ( meBCEne_ ) dqmStore_->removeElement( meBCEne_->getName() );
397  meBCEne_ = 0;
398 
399  if ( meBCNum_ ) dqmStore_->removeElement( meBCNum_->getName() );
400  meBCNum_ = 0;
401 
402  if ( meBCSiz_ ) dqmStore_->removeElement( meBCSiz_->getName() );
403  meBCSiz_ = 0;
404 
405  if ( meBCEneMap_ ) dqmStore_->removeElement( meBCEneMap_->getName() );
406  meBCEneMap_ = 0;
407 
408  if ( meBCNumMap_ ) dqmStore_->removeElement( meBCNumMap_->getName() );
409  meBCNumMap_ = 0;
410 
411  if ( meBCETMap_ ) dqmStore_->removeElement( meBCETMap_->getName() );
412  meBCETMap_ = 0;
413 
414  if ( meBCSizMap_ ) dqmStore_->removeElement( meBCSizMap_->getName() );
415  meBCSizMap_ = 0;
416 
417  if ( meBCEneMapProjEta_ ) dqmStore_->removeElement( meBCEneMapProjEta_->getName() );
418  meBCEneMapProjEta_ = 0;
419 
420  if ( meBCEneMapProjPhi_ ) dqmStore_->removeElement( meBCEneMapProjPhi_->getName() );
421  meBCEneMapProjPhi_ = 0;
422 
423  if ( meBCNumMapProjEta_ ) dqmStore_->removeElement( meBCNumMapProjEta_->getName() );
424  meBCNumMapProjEta_ = 0;
425 
426  if ( meBCNumMapProjPhi_ ) dqmStore_->removeElement( meBCNumMapProjPhi_->getName() );
427  meBCNumMapProjPhi_ = 0;
428 
429  if ( meBCETMapProjEta_ ) dqmStore_->removeElement( meBCETMapProjEta_->getName() );
430  meBCETMapProjEta_ = 0;
431 
432  if ( meBCETMapProjPhi_ ) dqmStore_->removeElement( meBCETMapProjPhi_->getName() );
433  meBCETMapProjPhi_ = 0;
434 
435  if ( meBCSizMapProjEta_ ) dqmStore_->removeElement( meBCSizMapProjEta_->getName() );
436  meBCSizMapProjEta_ = 0;
437 
438  if ( meBCSizMapProjPhi_ ) dqmStore_->removeElement( meBCSizMapProjPhi_->getName() );
439  meBCSizMapProjPhi_ = 0;
440 
441  if ( meSCEne_ ) dqmStore_->removeElement( meSCEne_->getName() );
442  meSCEne_ = 0;
443 
444  if ( meSCNum_ ) dqmStore_->removeElement( meSCNum_->getName() );
445  meSCNum_ = 0;
446 
447  if ( meSCSiz_ ) dqmStore_->removeElement( meSCSiz_->getName() );
448  meSCSiz_ = 0;
449 
450  if ( meSCCrystalSiz_ ) dqmStore_->removeElement( meSCCrystalSiz_->getName() );
451  meSCCrystalSiz_ = 0;
452 
453  if ( meSCSeedEne_ ) dqmStore_->removeElement( meSCSeedEne_->getName() );
454  meSCSeedEne_ = 0;
455 
456  if ( meSCEne2_ ) dqmStore_->removeElement( meSCEne2_->getName() );
457  meSCEne2_ = 0;
458 
459  if ( meSCEneVsEMax_ ) dqmStore_->removeElement( meSCEneVsEMax_->getName() );
460  meSCEneVsEMax_ = 0;
461 
462  if ( meSCEneLowScale_ ) dqmStore_->removeElement( meSCEneLowScale_->getName() );
463  meSCEneLowScale_ = 0;
464 
465  if ( meSCSeedMapOcc_ ) dqmStore_->removeElement( meSCSeedMapOcc_->getName() );
466  meSCSeedMapOcc_ = 0;
467 
468  if ( meSCMapSingleCrystal_ ) dqmStore_->removeElement( meSCMapSingleCrystal_->getName() );
469  meSCMapSingleCrystal_ = 0;
470 
471  if ( mes1s9_ ) dqmStore_->removeElement( mes1s9_->getName() );
472  mes1s9_ = 0;
473 
474  if ( mes1s9thr_ ) dqmStore_->removeElement( mes1s9thr_->getName() );
475  mes1s9thr_ = 0;
476 
477  if ( mes9s25_ ) dqmStore_->removeElement( mes9s25_->getName() );
478  mes9s25_ = 0;
479 
480  if ( meInvMassPi0_ ) dqmStore_->removeElement( meInvMassPi0_->getName() );
481  meInvMassPi0_ = 0;
482 
483  if ( meInvMassJPsi_ ) dqmStore_->removeElement( meInvMassJPsi_->getName() );
484  meInvMassJPsi_ = 0;
485 
486  if ( meInvMassZ0_ ) dqmStore_->removeElement( meInvMassZ0_->getName() );
487  meInvMassZ0_ = 0;
488 
489  if ( meInvMassHigh_ ) dqmStore_->removeElement( meInvMassHigh_->getName() );
490  meInvMassHigh_ = 0;
491 
492  if ( meInvMassPi0Sel_ ) dqmStore_->removeElement( meInvMassPi0Sel_->getName() );
493  meInvMassPi0Sel_ = 0;
494 
495  if ( meInvMassJPsiSel_ ) dqmStore_->removeElement( meInvMassJPsiSel_->getName() );
496  meInvMassJPsiSel_ = 0;
497 
498  if ( meInvMassZ0Sel_ ) dqmStore_->removeElement( meInvMassZ0Sel_->getName() );
499  meInvMassZ0Sel_ = 0;
500 
501  if ( meInvMassHighSel_ ) dqmStore_->removeElement( meInvMassHighSel_->getName() );
502  meInvMassHighSel_ = 0;
503 
504  }
505 
506  init_ = false;
507 
508 }
509 
510 void EBClusterTask::endJob(void){
511 
512  edm::LogInfo("EBClusterTask") << "analyzed " << ievt_ << " events";
513 
514  if ( enableCleanup_ ) this->cleanup();
515 
516 }
517 
518 void EBClusterTask::analyze(const edm::Event& e, const edm::EventSetup& c){
519 
520  bool enable = false;
521 
522  edm::Handle<EcalRawDataCollection> dcchs;
523 
524  if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
525 
526  for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
527 
528  if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue;
529 
530  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH4 ||
531  dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH2 ) enable = true;
532 
533  if ( dcchItr->getRunType() == EcalDCCHeaderBlock::COSMIC ||
534  dcchItr->getRunType() == EcalDCCHeaderBlock::MTCC ||
535  dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
536  dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
537  dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_LOCAL ||
538  dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_LOCAL ) enable = true;
539 
540  break;
541 
542  }
543 
544  } else {
545 
546  enable = true;
547  edm::LogWarning("EBClusterTask") << EcalRawDataCollection_ << " not available";
548 
549  }
550 
551  if ( ! enable ) return;
552 
553  if ( ! init_ ) this->setup();
554 
555  ievt_++;
556 
557  // ECAL topology
558  edm::ESHandle<CaloTopology> pTopology;
559  c.get<CaloTopologyRecord>().get(pTopology);
560  if ( !pTopology.isValid() ) {
561  edm::LogWarning("EBClusterTask") << "Topology not valid";
562  return;
563  }
564  const CaloTopology* topology = pTopology.product();
565 
566  // recHits
567  edm::Handle< EcalRecHitCollection > pEBRecHits;
568  e.getByLabel( EcalRecHitCollection_, pEBRecHits );
569  if ( !pEBRecHits.isValid() ) {
570  edm::LogWarning("EBClusterTask") << "RecHit collection " << EcalRecHitCollection_ << " not available.";
571  return;
572  }
573  const EcalRecHitCollection* ebRecHits = pEBRecHits.product();
574 
575  reco::BasicClusterCollection bcSel;
576 
577  // --- Barrel Basic Clusters ---
578  edm::Handle<reco::BasicClusterCollection> pBasicClusters;
579  if ( e.getByLabel(BasicClusterCollection_, pBasicClusters) ) {
580 
581  int nbcc = pBasicClusters->size();
582  if ( nbcc > 0 ) meBCNum_->Fill(float(nbcc));
583 
584  for ( reco::BasicClusterCollection::const_iterator bCluster = pBasicClusters->begin(); bCluster != pBasicClusters->end(); ++bCluster ) {
585 
586  meBCEne_->Fill(bCluster->energy());
587  meBCSiz_->Fill(float(bCluster->size()));
588 
589  float xphi = bCluster->phi();
590  if ( xphi > M_PI*(9-1.5)/9 ) xphi = xphi - M_PI*2;
591 
592  meBCEneMap_->Fill(xphi, bCluster->eta(), bCluster->energy());
593  meBCEneMapProjEta_->Fill(bCluster->eta(), bCluster->energy());
594  meBCEneMapProjPhi_->Fill(xphi, bCluster->energy());
595 
596  meBCNumMap_->Fill(xphi, bCluster->eta());
597  meBCNumMapProjEta_->Fill(bCluster->eta());
598  meBCNumMapProjPhi_->Fill(xphi);
599 
600  meBCSizMap_->Fill(xphi, bCluster->eta(), float(bCluster->size()));
601  meBCSizMapProjEta_->Fill(bCluster->eta(), float(bCluster->size()));
602  meBCSizMapProjPhi_->Fill(xphi, float(bCluster->size()));
603 
604  meBCETMap_->Fill(xphi, bCluster->eta(), float(bCluster->energy()) * sin(bCluster->position().theta()));
605  meBCETMapProjEta_->Fill(bCluster->eta(), float(bCluster->energy()) * sin(bCluster->position().theta()));
606  meBCETMapProjPhi_->Fill(xphi, float(bCluster->energy()) * sin(bCluster->position().theta()));
607 
608  float e2x2 = EcalClusterTools::e2x2( *bCluster, ebRecHits, topology );
609  float e3x3 = EcalClusterTools::e3x3( *bCluster, ebRecHits, topology );
610 
611  // fill the selected cluster collection
612  float pt = std::abs( bCluster->energy()*sin(bCluster->position().theta()) );
613  if ( pt > thrClusEt_ && e2x2/e3x3 > thrS4S9_ ) bcSel.push_back(*bCluster);
614  }
615 
616  } else {
617  edm::LogWarning("EBClusterTask") << BasicClusterCollection_ << " not available";
618  }
619 
620  for ( reco::BasicClusterCollection::const_iterator bc1 = bcSel.begin(); bc1 != bcSel.end(); ++bc1 ) {
621  TLorentzVector bc1P;
622  bc1P.SetPtEtaPhiE(std::abs(bc1->energy()*sin(bc1->position().theta())),
623  bc1->eta(), bc1->phi(), bc1->energy());
624  for ( reco::BasicClusterCollection::const_iterator bc2 = bc1+1; bc2 != bcSel.end(); ++bc2 ) {
625  TLorentzVector bc2P;
626  bc2P.SetPtEtaPhiE(std::abs(bc2->energy()*sin(bc2->position().theta())),
627  bc2->eta(), bc2->phi(), bc2->energy());
628 
629  TLorentzVector candP = bc1P + bc2P;
630 
631  if ( candP.Pt() > thrCandEt_ ) {
632  float mass = candP.M();
633  if ( mass < 0.500 ) {
634  meInvMassPi0Sel_->Fill( mass );
635  } else if ( mass > 2.9 && mass < 3.3 ) {
636  meInvMassJPsiSel_->Fill( mass );
637  } else if ( mass > 40 && mass < 110 ) {
638  meInvMassZ0Sel_->Fill( mass );
639  } else if ( mass > 110 ) {
640  meInvMassHighSel_->Fill( mass );
641  }
642 
643  }
644 
645  }
646  }
647 
648  // --- Barrel Super Clusters ---
649  edm::Handle<reco::SuperClusterCollection> pSuperClusters;
650  if ( e.getByLabel(SuperClusterCollection_, pSuperClusters) ) {
651 
652  int nscc = pSuperClusters->size();
653  if ( nscc > 0 ) meSCNum_->Fill(float(nscc));
654 
655  TLorentzVector sc1_p(0,0,0,0);
656  TLorentzVector sc2_p(0,0,0,0);
657 
658  for ( reco::SuperClusterCollection::const_iterator sCluster = pSuperClusters->begin(); sCluster != pSuperClusters->end(); ++sCluster ) {
659 
660  // energy, size
661  meSCEne_->Fill( sCluster->energy() );
662  meSCSiz_->Fill( float(sCluster->clustersSize()) );
663 
664  reco::CaloClusterPtr theSeed = sCluster->seed();
665 
666  // Find the seed rec hit
667  std::vector< std::pair<DetId,float> > sIds = sCluster->hitsAndFractions();
668 
669  float eMax, e2nd;
670  EcalRecHitCollection::const_iterator seedItr = ebRecHits->begin();
671  EcalRecHitCollection::const_iterator secondItr = ebRecHits->begin();
672 
673  for(std::vector< std::pair<DetId,float> >::const_iterator idItr = sIds.begin(); idItr != sIds.end(); ++idItr) {
674  DetId id = idItr->first;
675  if(id.det() != DetId::Ecal) { continue; }
676  EcalRecHitCollection::const_iterator hitItr = ebRecHits->find(id);
677  if(hitItr == ebRecHits->end()) { continue; }
678  if(hitItr->energy() > secondItr->energy()) { secondItr = hitItr; }
679  if(hitItr->energy() > seedItr->energy()) { std::swap(seedItr,secondItr); }
680  }
681 
682  eMax = seedItr->energy();
683  e2nd = secondItr->energy();
684  EBDetId seedId = (EBDetId) seedItr->id();
685 
686  float e3x3 = EcalClusterTools::e3x3( *theSeed, ebRecHits, topology );
687  float e5x5 = EcalClusterTools::e5x5( *theSeed, ebRecHits, topology );
688 
689  meSCCrystalSiz_->Fill(sIds.size());
690  meSCSeedEne_->Fill(eMax);
691  meSCEne2_->Fill(eMax+e2nd);
692  meSCEneVsEMax_->Fill(eMax,sCluster->energy());
693  meSCEneLowScale_->Fill(sCluster->energy());
694 
695  // Prepare to fill maps
696  int ebeta = seedId.ieta();
697  int ebphi = seedId.iphi();
698  float xebeta = ebeta - 0.5 * seedId.zside();
699  float xebphi = ebphi - 0.5;
700 
701  meSCSeedMapOcc_->Fill(xebphi,xebeta);
702 
703  if(sIds.size() == 1) meSCMapSingleCrystal_->Fill(xebphi,xebeta);
704 
705  mes1s9_->Fill( eMax/e3x3 );
706  if ( eMax > 3.0 ) mes1s9thr_->Fill( eMax/e3x3 );
707  mes9s25_->Fill( e3x3/e5x5 );
708 
709  if ( nscc >= 2 ) {
710  // look for the two most energetic super clusters
711  if ( sCluster->energy() > sc1_p.Energy() ) {
712  sc2_p=sc1_p;
713  sc1_p.SetPtEtaPhiE(std::abs(sCluster->energy()*sin(sCluster->position().theta())),
714  sCluster->eta(), sCluster->phi(), sCluster->energy());
715  } else if ( sCluster->energy() > sc2_p.Energy() ) {
716  sc2_p.SetPtEtaPhiE(std::abs(sCluster->energy()*sin(sCluster->position().theta())),
717  sCluster->eta(), sCluster->phi(), sCluster->energy());
718  }
719  }
720 
721  }
722  // Get the invariant mass of the two most energetic super clusters
723  if ( nscc >= 2 ) {
724  TLorentzVector sum = sc1_p+sc2_p;
725  float mass = sum.M();
726  if ( mass < 0.500 ) {
727  meInvMassPi0_->Fill( mass );
728  } else if ( mass > 2.9 && mass < 3.3 ) {
729  meInvMassJPsi_->Fill( mass );
730  } else if ( mass > 40 && mass < 110 ) {
731  meInvMassZ0_->Fill( mass );
732  } else if ( mass > 110 ) {
733  meInvMassHigh_->Fill( mass );
734  }
735  }
736 
737  } else {
738 
739  edm::LogWarning("EBClusterTask") << SuperClusterCollection_ << " not available";
740 
741  }
742 
743 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EBClusterTask::meBCSizMap_
MonitorElement * meBCSizMap_
Definition: EBClusterTask.h:80
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
EBClusterTask::EcalRawDataCollection_
edm::InputTag EcalRawDataCollection_
Definition: EBClusterTask.h:68
EBClusterTask::meBCEneMapProjPhi_
MonitorElement * meBCEneMapProjPhi_
Definition: EBClusterTask.h:77
EBClusterTask::meBCNumMapProjPhi_
MonitorElement * meBCNumMapProjPhi_
Definition: EBClusterTask.h:78
EBClusterTask::meSCSiz_
MonitorElement * meSCSiz_
Definition: EBClusterTask.h:84
EBClusterTask::meSCSeedEne_
MonitorElement * meSCSeedEne_
Definition: EBClusterTask.h:87
EBClusterTask::setup
void setup(void)
Setup.
Definition: EBClusterTask.cc:217
EBClusterTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EBClusterTask.cc:129
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
EBClusterTask::meBCSiz_
MonitorElement * meBCSiz_
Definition: EBClusterTask.h:75
EBClusterTask::meBCSizMapProjEta_
MonitorElement * meBCSizMapProjEta_
Definition: EBClusterTask.h:80
EBClusterTask::meInvMassPi0Sel_
MonitorElement * meInvMassPi0Sel_
Definition: EBClusterTask.h:101
Numbers.h
Some &quot;id&quot; conversions.
EBClusterTask::mes9s25_
MonitorElement * mes9s25_
Definition: EBClusterTask.h:96
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
EBClusterTask::mes1s9thr_
MonitorElement * mes1s9thr_
Definition: EBClusterTask.h:95
EBClusterTask::meSCEne_
MonitorElement * meSCEne_
Definition: EBClusterTask.h:82
EBClusterTask::meBCEneMap_
MonitorElement * meBCEneMap_
Definition: EBClusterTask.h:77
abs
#define abs(x)
Definition: mlp_lapack.h:159
EBClusterTask::meSCCrystalSiz_
MonitorElement * meSCCrystalSiz_
Definition: EBClusterTask.h:86
EBClusterTask::meBCEneMapProjEta_
MonitorElement * meBCEneMapProjEta_
Definition: EBClusterTask.h:77
EBClusterTask::meInvMassJPsiSel_
MonitorElement * meInvMassJPsiSel_
Definition: EBClusterTask.h:102
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
EBClusterTask::meInvMassZ0_
MonitorElement * meInvMassZ0_
Definition: EBClusterTask.h:99
EBDetId::iphi
int iphi() const
get the crystal iphi
Definition: EBDetId.h:46
EBClusterTask::meSCNum_
MonitorElement * meSCNum_
Definition: EBClusterTask.h:83
EBClusterTask::meInvMassPi0_
MonitorElement * meInvMassPi0_
Definition: EBClusterTask.h:97
EBClusterTask::meBCSizMapProjPhi_
MonitorElement * meBCSizMapProjPhi_
Definition: EBClusterTask.h:80
EBClusterTask::EBClusterTask
EBClusterTask(const edm::ParameterSet &ps)
Constructor.
Definition: EBClusterTask.cc:44
EBClusterTask::BasicClusterCollection_
edm::InputTag BasicClusterCollection_
Definition: EBClusterTask.h:69
EBClusterTask::meBCEne_
MonitorElement * meBCEne_
Definition: EBClusterTask.h:73
std::swap
void swap(edm::DataFrameContainer &lhs, edm::DataFrameContainer &rhs)
Definition: DataFrameContainer.h:223
EBClusterTask::meSCSeedMapOcc_
MonitorElement * meSCSeedMapOcc_
Definition: EBClusterTask.h:91
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2572
EcalClusterTools::e3x3
static float e3x3(const reco::BasicCluster &cluster, const EcalRecHitCollection *recHits, const CaloTopology *topology)
Definition: EcalClusterTools.cc:240
EBClusterTask::meBCETMap_
MonitorElement * meBCETMap_
Definition: EBClusterTask.h:79
EBClusterTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EBClusterTask.cc:518
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
EBClusterTask::dqmStore_
DQMStore * dqmStore_
Definition: EBClusterTask.h:60
EBDetId::ieta
int ieta() const
get the crystal ieta
Definition: EBDetId.h:44
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
EBClusterTask::meSCEneVsEMax_
MonitorElement * meSCEneVsEMax_
Definition: EBClusterTask.h:89
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
EBClusterTask::reset
void reset(void)
Reset.
Definition: EBClusterTask.cc:141
jptDQMConfig_cff.eMax
int eMax
Definition: jptDQMConfig_cff.py:21
EBClusterTask::meSCEneLowScale_
MonitorElement * meSCEneLowScale_
Definition: EBClusterTask.h:90
EBClusterTask::meSCMapSingleCrystal_
MonitorElement * meSCMapSingleCrystal_
Definition: EBClusterTask.h:92
EBClusterTask::mes1s9_
MonitorElement * mes1s9_
Definition: EBClusterTask.h:94
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:301
EBClusterTask::meInvMassHigh_
MonitorElement * meInvMassHigh_
Definition: EBClusterTask.h:100
DetId
Definition: DetId.h:20
EBClusterTask::meInvMassHighSel_
MonitorElement * meInvMassHighSel_
Definition: EBClusterTask.h:104
EBClusterTask::meBCNumMap_
MonitorElement * meBCNumMap_
Definition: EBClusterTask.h:78
EBClusterTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EBClusterTask.cc:137
M_PI
#define M_PI
Definition: BFit3D.cc:3
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:50
EBClusterTask::beginJob
void beginJob(void)
BeginJob.
Definition: EBClusterTask.cc:118
EBClusterTask::meBCNumMapProjEta_
MonitorElement * meBCNumMapProjEta_
Definition: EBClusterTask.h:78
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
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
EBClusterTask::EcalRecHitCollection_
edm::InputTag EcalRecHitCollection_
Definition: EBClusterTask.h:71
EBClusterTask::meSCEne2_
MonitorElement * meSCEne2_
Definition: EBClusterTask.h:88
EBClusterTask::meBCNum_
MonitorElement * meBCNum_
Definition: EBClusterTask.h:74
edm::SortedCollection::find
iterator find(key_type k)
Definition: SortedCollection.h:262
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
EBClusterTask::meBCETMapProjEta_
MonitorElement * meBCETMapProjEta_
Definition: EBClusterTask.h:79
EBClusterTask::SuperClusterCollection_
edm::InputTag SuperClusterCollection_
Definition: EBClusterTask.h:70
EBClusterTask::prefixME_
std::string prefixME_
Definition: EBClusterTask.h:62
Numbers::subDet
static EcalSubdetector subDet(const EBDetId &id)
Definition: Numbers.cc:145
EBClusterTask::endJob
void endJob(void)
EndJob.
Definition: EBClusterTask.cc:510
EBClusterTask::meInvMassJPsi_
MonitorElement * meInvMassJPsi_
Definition: EBClusterTask.h:98
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
EBClusterTask::cleanup
void cleanup(void)
Cleanup.
Definition: EBClusterTask.cc:389
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
EBClusterTask::meInvMassZ0Sel_
MonitorElement * meInvMassZ0Sel_
Definition: EBClusterTask.h:103
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:294
edm::Run
Definition: Run.h:33
EBClusterTask::meBCETMapProjPhi_
MonitorElement * meBCETMapProjPhi_
Definition: EBClusterTask.h:79
EBDetId::zside
int zside() const
get the z-side of the crystal (1/-1)
Definition: EBDetId.h:40
EBClusterTask::~EBClusterTask
virtual ~EBClusterTask()
Destructor.
Definition: EBClusterTask.cc:114
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
Generated for CMSSW Reference Manual by doxygen 1.8.5