CMS 3D CMS Logo
HBHEIsolatedNoiseReflagger.cc
Go to the documentation of this file.
1 /*
2 Description: "Reflags" HB/HE hits based on their ECAL, HCAL, and tracking isolation.
3 
4 Original Author: John Paul Chou (Brown University)
5  Thursday, September 2, 2010
6 */
7 
8 #include "HBHEIsolatedNoiseReflagger.h"
9 
10 #include "FWCore/Framework/interface/ESHandle.h"
11 #include "FWCore/Framework/interface/EventSetup.h"
12 #include "Geometry/Records/interface/CaloGeometryRecord.h"
13 #include "DataFormats/JetReco/interface/TrackExtrapolation.h"
14 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
15 #include "RecoLocalCalo/HcalRecAlgos/interface/HcalSeverityLevelComputerRcd.h"
16 #include "DataFormats/METReco/interface/HcalCaloFlagLabels.h"
17 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
18 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
19 #include "CondFormats/HcalObjects/interface/HcalChannelQuality.h"
20 #include "CondFormats/DataRecord/interface/HcalChannelQualityRcd.h"
21 #include "CondFormats/DataRecord/interface/HcalFrontEndMapRcd.h"
22 #include "CondFormats/HcalObjects/interface/HcalFrontEndMap.h"
23 
24 HBHEIsolatedNoiseReflagger::HBHEIsolatedNoiseReflagger(const edm::ParameterSet& iConfig)
25  :
26 
27  LooseHcalIsol_(iConfig.getParameter<double>("LooseHcalIsol")),
28  LooseEcalIsol_(iConfig.getParameter<double>("LooseEcalIsol")),
29  LooseTrackIsol_(iConfig.getParameter<double>("LooseTrackIsol")),
30  TightHcalIsol_(iConfig.getParameter<double>("TightHcalIsol")),
31  TightEcalIsol_(iConfig.getParameter<double>("TightEcalIsol")),
32  TightTrackIsol_(iConfig.getParameter<double>("TightTrackIsol")),
33 
34  LooseRBXEne1_(iConfig.getParameter<double>("LooseRBXEne1")),
35  LooseRBXEne2_(iConfig.getParameter<double>("LooseRBXEne2")),
36  LooseRBXHits1_(iConfig.getParameter<int>("LooseRBXHits1")),
37  LooseRBXHits2_(iConfig.getParameter<int>("LooseRBXHits2")),
38  TightRBXEne1_(iConfig.getParameter<double>("TightRBXEne1")),
39  TightRBXEne2_(iConfig.getParameter<double>("TightRBXEne2")),
40  TightRBXHits1_(iConfig.getParameter<int>("TightRBXHits1")),
41  TightRBXHits2_(iConfig.getParameter<int>("TightRBXHits2")),
42 
43  LooseHPDEne1_(iConfig.getParameter<double>("LooseHPDEne1")),
44  LooseHPDEne2_(iConfig.getParameter<double>("LooseHPDEne2")),
45  LooseHPDHits1_(iConfig.getParameter<int>("LooseHPDHits1")),
46  LooseHPDHits2_(iConfig.getParameter<int>("LooseHPDHits2")),
47  TightHPDEne1_(iConfig.getParameter<double>("TightHPDEne1")),
48  TightHPDEne2_(iConfig.getParameter<double>("TightHPDEne2")),
49  TightHPDHits1_(iConfig.getParameter<int>("TightHPDHits1")),
50  TightHPDHits2_(iConfig.getParameter<int>("TightHPDHits2")),
51 
52  LooseDiHitEne_(iConfig.getParameter<double>("LooseDiHitEne")),
53  TightDiHitEne_(iConfig.getParameter<double>("TightDiHitEne")),
54  LooseMonoHitEne_(iConfig.getParameter<double>("LooseMonoHitEne")),
55  TightMonoHitEne_(iConfig.getParameter<double>("TightMonoHitEne")),
56 
57  RBXEneThreshold_(iConfig.getParameter<double>("RBXEneThreshold")),
58 
59  debug_(iConfig.getUntrackedParameter<bool>("debug", true)),
60  objvalidator_(iConfig) {
61  tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
62  tok_EB_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ebInput"));
63  tok_EE_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("eeInput"));
64  tok_trackExt_ =
65  consumes<std::vector<reco::TrackExtrapolation> >(iConfig.getParameter<edm::InputTag>("trackExtrapolationInput"));
66 
67  produces<HBHERecHitCollection>();
68 }
69 
70 HBHEIsolatedNoiseReflagger::~HBHEIsolatedNoiseReflagger() {}
71 
72 void HBHEIsolatedNoiseReflagger::produce(edm::Event& iEvent, const edm::EventSetup& evSetup) {
73  // get the ECAL channel status map
74  edm::ESHandle<EcalChannelStatus> ecalChStatus;
75  evSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
76  const EcalChannelStatus* dbEcalChStatus = ecalChStatus.product();
77 
78  // get the HCAL channel status map
79 
80  edm::ESHandle<HcalChannelQuality> hcalChStatus;
81  evSetup.get<HcalChannelQualityRcd>().get("withTopo", hcalChStatus);
82  const HcalChannelQuality* dbHcalChStatus = hcalChStatus.product();
83 
84  // get the severity level computers
85  edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl;
86  evSetup.get<HcalSeverityLevelComputerRcd>().get(hcalSevLvlComputerHndl);
87  const HcalSeverityLevelComputer* hcalSevLvlComputer = hcalSevLvlComputerHndl.product();
88 
89  edm::ESHandle<EcalSeverityLevelAlgo> ecalSevLvlAlgoHndl;
90  evSetup.get<EcalSeverityLevelAlgoRcd>().get(ecalSevLvlAlgoHndl);
91  const EcalSeverityLevelAlgo* ecalSevLvlAlgo = ecalSevLvlAlgoHndl.product();
92 
93  // get the calotower mappings
94  edm::ESHandle<CaloTowerConstituentsMap> ctcm;
95  evSetup.get<CaloGeometryRecord>().get(ctcm);
96 
97  // get hcal frontend map
98  edm::ESHandle<HcalFrontEndMap> hfemapHndl;
99  evSetup.get<HcalFrontEndMapRcd>().get(hfemapHndl);
100  hfemap = hfemapHndl.product();
101 
102  // get the HB/HE hits
103  edm::Handle<HBHERecHitCollection> hbhehits_h;
104  iEvent.getByToken(tok_hbhe_, hbhehits_h);
105 
106  // get the ECAL hits
107  edm::Handle<EcalRecHitCollection> ebhits_h;
108  iEvent.getByToken(tok_EB_, ebhits_h);
109  edm::Handle<EcalRecHitCollection> eehits_h;
110  iEvent.getByToken(tok_EE_, eehits_h);
111 
112  // get the tracks
113  edm::Handle<std::vector<reco::TrackExtrapolation> > trackextraps_h;
114  iEvent.getByToken(tok_trackExt_, trackextraps_h);
115 
116  // set the status maps and severity level computers for the hit validator
117  objvalidator_.setHcalChannelQuality(dbHcalChStatus);
118  objvalidator_.setEcalChannelStatus(dbEcalChStatus);
119  objvalidator_.setHcalSeverityLevelComputer(hcalSevLvlComputer);
120  objvalidator_.setEcalSeverityLevelAlgo(ecalSevLvlAlgo);
121  objvalidator_.setEBRecHitCollection(&(*ebhits_h));
122  objvalidator_.setEERecHitCollection(&(*eehits_h));
123 
124  // organizer the hits
125  PhysicsTowerOrganizer pto(
126  iEvent, evSetup, hbhehits_h, ebhits_h, eehits_h, trackextraps_h, objvalidator_, *(ctcm.product()));
127  HBHEHitMapOrganizer organizer(hbhehits_h, objvalidator_, pto, hfemap);
128 
129  // get the rbxs, hpds, dihits, and monohits
130  std::vector<HBHEHitMap> rbxs;
131  std::vector<HBHEHitMap> hpds;
132  std::vector<HBHEHitMap> dihits;
133  std::vector<HBHEHitMap> monohits;
134  organizer.getRBXs(rbxs, LooseRBXEne1_ < TightRBXEne1_ ? LooseRBXEne1_ : TightRBXEne1_);
135  organizer.getHPDs(hpds, LooseHPDEne1_ < TightHPDEne1_ ? LooseHPDEne1_ : TightHPDEne1_);
136  organizer.getDiHits(dihits, LooseDiHitEne_ < TightDiHitEne_ ? LooseDiHitEne_ : TightDiHitEne_);
137  organizer.getMonoHits(monohits, LooseMonoHitEne_ < TightMonoHitEne_ ? LooseMonoHitEne_ : TightMonoHitEne_);
138 
139  if (debug_ && (!rbxs.empty() || !hpds.empty() || !dihits.empty() || !monohits.empty())) {
140  edm::LogInfo("HBHEIsolatedNoiseReflagger") << "RBXs:" << std::endl;
141  DumpHBHEHitMap(rbxs);
142  edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nHPDs:" << std::endl;
143  DumpHBHEHitMap(hpds);
144  edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nDiHits:" << std::endl;
145  DumpHBHEHitMap(dihits);
146  edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nMonoHits:" << std::endl;
147  DumpHBHEHitMap(monohits);
148  }
149 
150  // bool result=true;
151 
152  // determine which hits are noisy
153  std::set<const HBHERecHit*> noisehits;
154  for (int i = 0; i < static_cast<int>(rbxs.size()); i++) {
155  int nhits = rbxs[i].nHits();
156  double ene = rbxs[i].hitEnergy();
157  double trkfide = rbxs[i].hitEnergyTrackFiducial();
158  double isolhcale = rbxs[i].hcalEnergySameTowers() + rbxs[i].hcalEnergyNeighborTowers();
159  double isolecale = rbxs[i].ecalEnergySameTowers();
160  double isoltrke = rbxs[i].trackEnergySameTowers() + rbxs[i].trackEnergyNeighborTowers();
161  //
162  // RBX mistag reduction
163  bool isLooseIso = false;
164  bool isTightIso = false;
165  if (ene > RBXEneThreshold_ && ene > 0) { // New absolute iso-cut for high energy RBX clusters
166  if (isolhcale < LooseHcalIsol_ * RBXEneThreshold_ && isolecale < LooseEcalIsol_ * RBXEneThreshold_ &&
167  isoltrke < LooseTrackIsol_ * RBXEneThreshold_)
168  isLooseIso = true;
169  if (isolhcale < TightHcalIsol_ * RBXEneThreshold_ && isolecale < TightEcalIsol_ * RBXEneThreshold_ &&
170  isoltrke < TightTrackIsol_ * RBXEneThreshold_)
171  isTightIso = true;
172  }
173  if (ene <= RBXEneThreshold_ && ene > 0) { // Old relative iso-cut for low energy RBX clusters
174  if (isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ && isoltrke / ene < LooseTrackIsol_)
175  isLooseIso = true;
176  if (isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ && isoltrke / ene < TightTrackIsol_)
177  isTightIso = true;
178  }
179  //
180  if ((isLooseIso && ((trkfide > LooseRBXEne1_ && nhits >= LooseRBXHits1_) ||
181  (trkfide > LooseRBXEne2_ && nhits >= LooseRBXHits2_))) ||
182  (isTightIso && ((trkfide > TightRBXEne1_ && nhits >= TightRBXHits1_) ||
183  (trkfide > TightRBXEne2_ && nhits >= TightRBXHits2_)))) {
184  for (HBHEHitMap::hitmap_const_iterator it = rbxs[i].beginHits(); it != rbxs[i].endHits(); ++it)
185  noisehits.insert(it->first);
186  // result=false;
187  }
188  }
189 
190  for (int i = 0; i < static_cast<int>(hpds.size()); i++) {
191  int nhits = hpds[i].nHits();
192  double ene = hpds[i].hitEnergy();
193  double trkfide = hpds[i].hitEnergyTrackFiducial();
194  double isolhcale = hpds[i].hcalEnergySameTowers() + hpds[i].hcalEnergyNeighborTowers();
195  double isolecale = hpds[i].ecalEnergySameTowers();
196  double isoltrke = hpds[i].trackEnergySameTowers() + hpds[i].trackEnergyNeighborTowers();
197  if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
198  isoltrke / ene < LooseTrackIsol_ &&
199  ((trkfide > LooseHPDEne1_ && nhits >= LooseHPDHits1_) ||
200  (trkfide > LooseHPDEne2_ && nhits >= LooseHPDHits2_))) ||
201  (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
202  isoltrke / ene < TightTrackIsol_ &&
203  ((trkfide > TightHPDEne1_ && nhits >= TightHPDHits1_) ||
204  (trkfide > TightHPDEne2_ && nhits >= TightHPDHits2_)))) {
205  for (HBHEHitMap::hitmap_const_iterator it = hpds[i].beginHits(); it != hpds[i].endHits(); ++it)
206  noisehits.insert(it->first);
207  // result=false;
208  }
209  }
210 
211  for (int i = 0; i < static_cast<int>(dihits.size()); i++) {
212  double ene = dihits[i].hitEnergy();
213  double trkfide = dihits[i].hitEnergyTrackFiducial();
214  double isolhcale = dihits[i].hcalEnergySameTowers() + dihits[i].hcalEnergyNeighborTowers();
215  double isolecale = dihits[i].ecalEnergySameTowers();
216  double isoltrke = dihits[i].trackEnergySameTowers() + dihits[i].trackEnergyNeighborTowers();
217  if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
218  isoltrke / ene < LooseTrackIsol_ && trkfide > 0.99 * ene && trkfide > LooseDiHitEne_) ||
219  (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
220  isoltrke / ene < TightTrackIsol_ && ene > TightDiHitEne_)) {
221  for (HBHEHitMap::hitmap_const_iterator it = dihits[i].beginHits(); it != dihits[i].endHits(); ++it)
222  noisehits.insert(it->first);
223  // result=false;
224  }
225  }
226 
227  for (int i = 0; i < static_cast<int>(monohits.size()); i++) {
228  double ene = monohits[i].hitEnergy();
229  double trkfide = monohits[i].hitEnergyTrackFiducial();
230  double isolhcale = monohits[i].hcalEnergySameTowers() + monohits[i].hcalEnergyNeighborTowers();
231  double isolecale = monohits[i].ecalEnergySameTowers();
232  double isoltrke = monohits[i].trackEnergySameTowers() + monohits[i].trackEnergyNeighborTowers();
233  if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
234  isoltrke / ene < LooseTrackIsol_ && trkfide > 0.99 * ene && trkfide > LooseMonoHitEne_) ||
235  (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
236  isoltrke / ene < TightTrackIsol_ && ene > TightMonoHitEne_)) {
237  for (HBHEHitMap::hitmap_const_iterator it = monohits[i].beginHits(); it != monohits[i].endHits(); ++it)
238  noisehits.insert(it->first);
239  // result=false;
240  }
241  }
242 
243  // prepare the output HBHE RecHit collection
244  auto pOut = std::make_unique<HBHERecHitCollection>();
245  // loop over rechits, and set the new bit you wish to use
246  for (HBHERecHitCollection::const_iterator it = hbhehits_h->begin(); it != hbhehits_h->end(); ++it) {
247  const HBHERecHit* hit = &(*it);
248  HBHERecHit newhit(*hit);
249  if (noisehits.end() != noisehits.find(hit)) {
250  newhit.setFlagField(1, HcalCaloFlagLabels::HBHEIsolatedNoise);
251  }
252  pOut->push_back(newhit);
253  }
254 
255  iEvent.put(std::move(pOut));
256 
257  return;
258 }
259 
260 void HBHEIsolatedNoiseReflagger::DumpHBHEHitMap(std::vector<HBHEHitMap>& i) const {
261  for (std::vector<HBHEHitMap>::const_iterator it = i.begin(); it != i.end(); ++it) {
262  edm::LogInfo("HBHEIsolatedNoiseReflagger")
263  << "hit energy=" << it->hitEnergy() << "; # hits=" << it->nHits()
264  << "; hcal energy same=" << it->hcalEnergySameTowers() << "; ecal energy same=" << it->ecalEnergySameTowers()
265  << "; track energy same=" << it->trackEnergySameTowers()
266  << "; neighbor hcal energy=" << it->hcalEnergyNeighborTowers() << std::endl;
267  edm::LogInfo("HBHEIsolatedNoiseReflagger") << "hits:" << std::endl;
268  for (HBHEHitMap::hitmap_const_iterator it2 = it->beginHits(); it2 != it->endHits(); ++it2) {
269  const HBHERecHit* hit = it2->first;
270  edm::LogInfo("HBHEIsolatedNoiseReflagger")
271  << "RBX #=" << hfemap->lookupRBX(hit->id()) << "; HPD #=" << hfemap->lookupRMIndex(hit->id()) << "; "
272  << (*hit) << std::endl;
273  }
274  }
275  return;
276 }
277 
278 //define this as a plug-in
279 DEFINE_FWK_MODULE(HBHEIsolatedNoiseReflagger);
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
EcalSeverityLevelAlgo
Definition: EcalSeverityLevelAlgo.h:33
electrons_cff.bool
bool
Definition: electrons_cff.py:372
mps_fire.i
i
Definition: mps_fire.py:355
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
hit::id
unsigned int id
Definition: SiStripHitEffFromCalibTree.cc:92
HBHEIsolatedNoiseReflagger::TightRBXHits1_
int TightRBXHits1_
Definition: HBHEIsolatedNoiseReflagger.h:45
HBHEIsolatedNoiseReflagger
Definition: HBHEIsolatedNoiseReflagger.h:18
ObjectValidator::setEcalChannelStatus
void setEcalChannelStatus(const EcalChannelStatus *q)
Definition: HBHEIsolatedNoiseAlgos.h:105
ESHandle.h
HBHEIsolatedNoiseReflagger::TightDiHitEne_
double TightDiHitEne_
Definition: HBHEIsolatedNoiseReflagger.h:51
HcalFrontEndMap::lookupRMIndex
const int lookupRMIndex(DetId fId) const
Definition: HcalFrontEndMap.cc:63
HcalChannelQualityRcd
Definition: HcalChannelQualityRcd.h:8
EcalSeverityLevelAlgoRcd.h
HBHEIsolatedNoiseReflagger::TightHPDHits1_
int TightHPDHits1_
Definition: HBHEIsolatedNoiseReflagger.h:49
HBHEIsolatedNoiseReflagger::LooseHPDEne1_
double LooseHPDEne1_
Definition: HBHEIsolatedNoiseReflagger.h:46
CaloGeometryRecord
Definition: CaloGeometryRecord.h:30
HBHEHitMapOrganizer::getHPDs
void getHPDs(std::vector< HBHEHitMap > &v, double energy) const
Definition: HBHEIsolatedNoiseAlgos.cc:816
HBHEIsolatedNoiseReflagger::RBXEneThreshold_
double RBXEneThreshold_
Definition: HBHEIsolatedNoiseReflagger.h:55
HBHERecHit
Definition: HBHERecHit.h:13
HBHEIsolatedNoiseReflagger::LooseRBXHits1_
int LooseRBXHits1_
Definition: HBHEIsolatedNoiseReflagger.h:43
edm::LogInfo
Definition: MessageLogger.h:254
HBHEIsolatedNoiseReflagger::TightRBXEne2_
double TightRBXEne2_
Definition: HBHEIsolatedNoiseReflagger.h:44
HBHEIsolatedNoiseReflagger::TightTrackIsol_
double TightTrackIsol_
Definition: HBHEIsolatedNoiseReflagger.h:40
HcalChannelQuality
Definition: HcalChannelQuality.h:17
HBHEIsolatedNoiseReflagger::TightMonoHitEne_
double TightMonoHitEne_
Definition: HBHEIsolatedNoiseReflagger.h:53
EcalCondObjectContainer< EcalChannelStatusCode >
edm::Handle
Definition: AssociativeIterator.h:50
HBHEIsolatedNoiseReflagger::hfemap
const HcalFrontEndMap * hfemap
Definition: HBHEIsolatedNoiseReflagger.h:33
HBHEIsolatedNoiseReflagger::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HBHEIsolatedNoiseReflagger.h:31
TrackExtrapolation.h
HBHEHitMapOrganizer::getRBXs
void getRBXs(std::vector< HBHEHitMap > &v, double energy) const
Definition: HBHEIsolatedNoiseAlgos.cc:807
HBHEIsolatedNoiseReflagger::LooseEcalIsol_
double LooseEcalIsol_
Definition: HBHEIsolatedNoiseReflagger.h:36
HBHEHitMapOrganizer
Definition: HBHEIsolatedNoiseAlgos.h:330
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
ObjectValidator::setEERecHitCollection
void setEERecHitCollection(const EcalRecHitCollection *q)
Definition: HBHEIsolatedNoiseAlgos.h:109
HBHEIsolatedNoiseReflagger::TightHcalIsol_
double TightHcalIsol_
Definition: HBHEIsolatedNoiseReflagger.h:38
HBHEIsolatedNoiseReflagger::LooseRBXHits2_
int LooseRBXHits2_
Definition: HBHEIsolatedNoiseReflagger.h:43
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
HcalSeverityLevelComputer
Definition: HcalSeverityLevelComputer.h:24
edm::ESHandle
Definition: DTSurvey.h:22
EcalSeverityLevelAlgoRcd
Definition: EcalSeverityLevelAlgoRcd.h:12
HBHEHitMapOrganizer::getDiHits
void getDiHits(std::vector< HBHEHitMap > &v, double energy) const
Definition: HBHEIsolatedNoiseAlgos.cc:825
PhysicsTowerOrganizer
Definition: HBHEIsolatedNoiseAlgos.h:170
HBHEIsolatedNoiseReflagger::~HBHEIsolatedNoiseReflagger
~HBHEIsolatedNoiseReflagger() override
Definition: HBHEIsolatedNoiseReflagger.cc:70
ObjectValidator::setHcalSeverityLevelComputer
void setHcalSeverityLevelComputer(const HcalSeverityLevelComputer *q)
Definition: HBHEIsolatedNoiseAlgos.h:106
HBHEIsolatedNoiseReflagger::TightHPDEne2_
double TightHPDEne2_
Definition: HBHEIsolatedNoiseReflagger.h:48
HcalCaloFlagLabels::HBHEIsolatedNoise
Definition: HcalCaloFlagLabels.h:21
HBHEIsolatedNoiseReflagger::produce
void produce(edm::Event &, const edm::EventSetup &) override
Definition: HBHEIsolatedNoiseReflagger.cc:72
nhits
Definition: HIMultiTrackSelector.h:42
EcalSeverityLevelAlgo.h
HBHEIsolatedNoiseReflagger::LooseHPDHits1_
int LooseHPDHits1_
Definition: HBHEIsolatedNoiseReflagger.h:47
CaloGeometryRecord.h
HcalSeverityLevelComputerRcd.h
HcalFrontEndMap::lookupRBX
const std::string lookupRBX(DetId fId) const
brief lookup the RBX associated with the given logical id
Definition: HcalFrontEndMap.cc:71
HcalFrontEndMap.h
funct::true
true
Definition: Factorize.h:173
HBHEIsolatedNoiseReflagger::TightRBXHits2_
int TightRBXHits2_
Definition: HBHEIsolatedNoiseReflagger.h:45
edm::ParameterSet
Definition: ParameterSet.h:36
HBHEIsolatedNoiseReflagger::LooseHcalIsol_
double LooseHcalIsol_
Definition: HBHEIsolatedNoiseReflagger.h:35
HBHEIsolatedNoiseReflagger.h
ObjectValidator::setEcalSeverityLevelAlgo
void setEcalSeverityLevelAlgo(const EcalSeverityLevelAlgo *q)
Definition: HBHEIsolatedNoiseAlgos.h:107
HcalCaloFlagLabels.h
HBHEIsolatedNoiseReflagger::HBHEIsolatedNoiseReflagger
HBHEIsolatedNoiseReflagger(const edm::ParameterSet &)
Definition: HBHEIsolatedNoiseReflagger.cc:24
EcalChannelStatusRcd
Definition: EcalChannelStatusRcd.h:5
HBHEIsolatedNoiseReflagger::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HBHEIsolatedNoiseReflagger.h:30
HBHEIsolatedNoiseReflagger::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HBHEIsolatedNoiseReflagger.h:29
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:267
createfilelist.int
int
Definition: createfilelist.py:10
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalFrontEndMapRcd
Definition: HcalFrontEndMapRcd.h:7
HBHEIsolatedNoiseReflagger::TightRBXEne1_
double TightRBXEne1_
Definition: HBHEIsolatedNoiseReflagger.h:44
edm::EventSetup
Definition: EventSetup.h:57
ObjectValidator::setEBRecHitCollection
void setEBRecHitCollection(const EcalRecHitCollection *q)
Definition: HBHEIsolatedNoiseAlgos.h:108
get
#define get
HcalChannelQualityRcd.h
HcalChannelQuality.h
HBHEIsolatedNoiseReflagger::LooseHPDEne2_
double LooseHPDEne2_
Definition: HBHEIsolatedNoiseReflagger.h:46
HBHEIsolatedNoiseReflagger::TightHPDHits2_
int TightHPDHits2_
Definition: HBHEIsolatedNoiseReflagger.h:49
HBHEIsolatedNoiseReflagger::DumpHBHEHitMap
void DumpHBHEHitMap(std::vector< HBHEHitMap > &i) const
Definition: HBHEIsolatedNoiseReflagger.cc:260
HBHEIsolatedNoiseReflagger::LooseTrackIsol_
double LooseTrackIsol_
Definition: HBHEIsolatedNoiseReflagger.h:37
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
HBHEIsolatedNoiseReflagger::TightEcalIsol_
double TightEcalIsol_
Definition: HBHEIsolatedNoiseReflagger.h:39
HBHEIsolatedNoiseReflagger::tok_trackExt_
edm::EDGetTokenT< std::vector< reco::TrackExtrapolation > > tok_trackExt_
Definition: HBHEIsolatedNoiseReflagger.h:32
eostools.move
def move(src, dest)
Definition: eostools.py:511
HBHEIsolatedNoiseReflagger::LooseMonoHitEne_
double LooseMonoHitEne_
Definition: HBHEIsolatedNoiseReflagger.h:52
EventSetup.h
HBHEIsolatedNoiseReflagger::LooseRBXEne1_
double LooseRBXEne1_
Definition: HBHEIsolatedNoiseReflagger.h:42
HBHEIsolatedNoiseReflagger::LooseDiHitEne_
double LooseDiHitEne_
Definition: HBHEIsolatedNoiseReflagger.h:50
ObjectValidator::setHcalChannelQuality
void setHcalChannelQuality(const HcalChannelQuality *q)
Definition: HBHEIsolatedNoiseAlgos.h:104
HBHEIsolatedNoiseReflagger::debug_
bool debug_
Definition: HBHEIsolatedNoiseReflagger.h:57
CaloRecHit::setFlagField
constexpr void setFlagField(uint32_t value, int base, int width=1)
Definition: CaloRecHit.h:36
HBHEHitMap::hitmap_const_iterator
std::map< const HBHERecHit *, const PhysicsTower * >::const_iterator hitmap_const_iterator
Definition: HBHEIsolatedNoiseAlgos.h:224
HBHEIsolatedNoiseReflagger::TightHPDEne1_
double TightHPDEne1_
Definition: HBHEIsolatedNoiseReflagger.h:48
HBHEIsolatedNoiseReflagger::objvalidator_
ObjectValidator objvalidator_
Definition: HBHEIsolatedNoiseReflagger.h:60
edm::Event
Definition: Event.h:73
HBHEHitMapOrganizer::getMonoHits
void getMonoHits(std::vector< HBHEHitMap > &v, double energy) const
Definition: HBHEIsolatedNoiseAlgos.cc:833
HBHEIsolatedNoiseReflagger::LooseRBXEne2_
double LooseRBXEne2_
Definition: HBHEIsolatedNoiseReflagger.h:42
HcalFrontEndMapRcd.h
edm::InputTag
Definition: InputTag.h:15
HcalSeverityLevelComputerRcd
Definition: HcalSeverityLevelComputerRcd.h:23
HBHEIsolatedNoiseReflagger::LooseHPDHits2_
int LooseHPDHits2_
Definition: HBHEIsolatedNoiseReflagger.h:47
EcalChannelStatusRcd.h
hit
Definition: SiStripHitEffFromCalibTree.cc:88
Generated for CMSSW Reference Manual by doxygen 1.8.16