|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
|
// -*- C++ -*-
//
// Package: ESRecoSummary
// Class: ESRecoSummary
// Original Author: Martina Malberti
//
// system include files
#include <cmath>
#include <iostream>
#include <memory>
// user include files
#include "DQMOffline/Ecal/interface/ESRecoSummary.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "Geometry/Records/interface/CaloTopologyRecord.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "DataFormats/EcalDetId/interface/ESDetId.h"
#include "DataFormats/EgammaReco/interface/BasicCluster.h"
#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
#include "DataFormats/EgammaReco/interface/PreshowerCluster.h"
#include "DataFormats/EgammaReco/interface/SuperCluster.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalTools.h"
#include "RecoLocalCalo/EcalRecAlgos/interface/EcalCleaningAlgo.h"
#include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
//
// constructors and destructor
//
ESRecoSummary::ESRecoSummary(const edm::ParameterSet &ps) {
prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
// now do what ever initialization is needed
esRecHitCollection_ = consumes<ESRecHitCollection>(ps.getParameter<edm::InputTag>("recHitCollection_ES"));
esClusterCollectionX_ =
consumes<reco::PreshowerClusterCollection>(ps.getParameter<edm::InputTag>("ClusterCollectionX_ES"));
esClusterCollectionY_ =
consumes<reco::PreshowerClusterCollection>(ps.getParameter<edm::InputTag>("ClusterCollectionY_ES"));
superClusterCollection_EE_ =
consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("superClusterCollection_EE"));
}
void ESRecoSummary::bookHistograms(DQMStore::IBooker &iBooker, edm::Run const &, edm::EventSetup const &) {
// Monitor Elements (ex THXD)
iBooker.setCurrentFolder(prefixME_ + "/ESRecoSummary"); // to organise the histos in folders
// Preshower ----------------------------------------------
h_recHits_ES_energyMax = iBooker.book1D("recHits_ES_energyMax", "recHits_ES_energyMax", 200, 0., 0.01);
h_recHits_ES_time = iBooker.book1D("recHits_ES_time", "recHits_ES_time", 200, -100., 100.);
h_esClusters_energy_plane1 = iBooker.book1D("esClusters_energy_plane1", "esClusters_energy_plane1", 200, 0., 0.01);
h_esClusters_energy_plane2 = iBooker.book1D("esClusters_energy_plane2", "esClusters_energy_plane2", 200, 0., 0.01);
h_esClusters_energy_ratio = iBooker.book1D("esClusters_energy_ratio", "esClusters_energy_ratio", 200, 0., 20.);
}
//
// member functions
//
// ------------ method called to for each event ------------
void ESRecoSummary::analyze(const edm::Event &ev, const edm::EventSetup &) {
// Preshower RecHits
edm::Handle<ESRecHitCollection> recHitsES;
ev.getByToken(esRecHitCollection_, recHitsES);
const ESRecHitCollection *thePreShowerRecHits = recHitsES.product();
if (!recHitsES.isValid()) {
std::cerr << "ESRecoSummary::analyze --> recHitsES not found" << std::endl;
}
float maxRecHitEnergyES = -999.;
for (ESRecHitCollection::const_iterator esItr = thePreShowerRecHits->begin(); esItr != thePreShowerRecHits->end();
++esItr) {
h_recHits_ES_time->Fill(esItr->time());
if (esItr->energy() > maxRecHitEnergyES)
maxRecHitEnergyES = esItr->energy();
} // end loop over ES rec Hits
h_recHits_ES_energyMax->Fill(maxRecHitEnergyES);
// ES clusters in X plane
edm::Handle<reco::PreshowerClusterCollection> esClustersX;
ev.getByToken(esClusterCollectionX_, esClustersX);
const reco::PreshowerClusterCollection *ESclustersX = esClustersX.product();
// ES clusters in Y plane
edm::Handle<reco::PreshowerClusterCollection> esClustersY;
ev.getByToken(esClusterCollectionY_, esClustersY);
const reco::PreshowerClusterCollection *ESclustersY = esClustersY.product();
// ... endcap
edm::Handle<reco::SuperClusterCollection> superClusters_EE_h;
ev.getByToken(superClusterCollection_EE_, superClusters_EE_h);
const reco::SuperClusterCollection *theEndcapSuperClusters = superClusters_EE_h.product();
if (!superClusters_EE_h.isValid()) {
std::cerr << "EcalRecHitSummary::analyze --> superClusters_EE_h not found" << std::endl;
}
// loop over all super clusters
for (reco::SuperClusterCollection::const_iterator itSC = theEndcapSuperClusters->begin();
itSC != theEndcapSuperClusters->end();
++itSC) {
if (fabs(itSC->eta()) < 1.65 || fabs(itSC->eta()) > 2.6)
continue;
float ESenergyPlane1 = 0.;
float ESenergyPlane2 = 0.;
// Loop over all ECAL Basic clusters in the supercluster
for (reco::CaloCluster_iterator ecalBasicCluster = itSC->clustersBegin(); ecalBasicCluster != itSC->clustersEnd();
ecalBasicCluster++) {
const reco::CaloClusterPtr ecalBasicClusterPtr = *(ecalBasicCluster);
for (reco::PreshowerClusterCollection::const_iterator iESClus = ESclustersX->begin();
iESClus != ESclustersX->end();
++iESClus) {
const reco::CaloClusterPtr preshBasicCluster = iESClus->basicCluster();
const reco::PreshowerCluster *esCluster = &*iESClus;
if (preshBasicCluster == ecalBasicClusterPtr) {
ESenergyPlane1 += esCluster->energy();
}
} // end of x loop
for (reco::PreshowerClusterCollection::const_iterator iESClus = ESclustersY->begin();
iESClus != ESclustersY->end();
++iESClus) {
const reco::CaloClusterPtr preshBasicCluster = iESClus->basicCluster();
const reco::PreshowerCluster *esCluster = &*iESClus;
if (preshBasicCluster == ecalBasicClusterPtr) {
ESenergyPlane2 += esCluster->energy();
}
} // end of y loop
} // end loop over all basic clusters in the supercluster
// cout<<"DQM : "<<ESenergyPlane1<<" "<<ESenergyPlane2<<endl;
h_esClusters_energy_plane1->Fill(ESenergyPlane1);
h_esClusters_energy_plane2->Fill(ESenergyPlane2);
if (ESenergyPlane1 > 0 && ESenergyPlane2 > 0)
h_esClusters_energy_ratio->Fill(ESenergyPlane1 / ESenergyPlane2);
} // end loop over superclusters
}
// define this as a plug-in
DEFINE_FWK_MODULE(ESRecoSummary);
|