CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
reco::HcalNoiseInfoProducer Class Reference

#include <HcalNoiseInfoProducer.h>

Inheritance diagram for reco::HcalNoiseInfoProducer:
edm::stream::EDProducer<>

Public Member Functions

 HcalNoiseInfoProducer (const edm::ParameterSet &)
 
 ~HcalNoiseInfoProducer () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProducer & operator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Member Functions

void fillcalotwrs (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const
 
void filldigis (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &)
 
void filljetinfo (edm::Event &, const edm::EventSetup &, HcalNoiseSummary &) const
 
void fillOtherSummaryVariables (HcalNoiseSummary &summary, const CommonHcalNoiseRBXData &data) const
 
void fillrechits (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const
 
void filltracks (edm::Event &, const edm::EventSetup &, HcalNoiseSummary &) const
 
void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

std::vector< float > adc2fC
 
std::vector< float > adc2fCHF
 
HcalNoiseAlgo algo_
 
double calibdigiHBHEthreshold_
 
std::vector< int > calibdigiHBHEtimeslices_
 
double calibdigiHFthreshold_
 
std::vector< int > calibdigiHFtimeslices_
 
edm::ESGetToken< CaloGeometry,
CaloGeometryRecord > 		calogeometry_token_
 
edm::EDGetTokenT
< CaloTowerCollection > 		calotower_token_
 
std::string caloTowerCollName_
 
std::string digiCollName_
 
bool fillCaloTowers_
 
bool fillDigis_
 
bool fillLaserMonitor_
 
bool fillRecHits_
 
bool fillTracks_
 
edm::EDGetTokenT
< HBHEDigiCollection > 		hbhedigi_token_
 
edm::EDGetTokenT
< HBHERecHitCollection > 		hbherechit_token_
 
uint32_t HcalAcceptSeverityLevel_
 
edm::EDGetTokenT
< HcalCalibDigiCollection > 		hcalcalibdigi_token_
 
std::vector< int > HcalRecHitFlagsToBeExcluded_
 
edm::ESGetToken< HcalTopology,
HcalRecNumberingRecord > 		hcaltopo_token_
 
edm::EDGetTokenT
< reco::PFJetCollection > 		jet_token_
 
std::string jetCollName_
 
std::vector< int > laserMonCBoxList_
 
edm::EDGetTokenT
< QIE10DigiCollection > 		lasermondigi_token_
 
std::vector< int > laserMonIEtaList_
 
std::vector< int > laserMonIPhiList_
 
unsigned laserMonitorSamples_
 
int laserMonitorTSEnd_
 
int laserMonitorTSStart_
 
int maxCaloTowerIEta_
 
int maxjetindex_
 
double maxNHF_
 
int maxProblemRBXs_
 
double maxTrackEta_
 
double minHighHitE_
 
double minLowHitE_
 
double minR45HitE_
 
double minRecHitE_
 
double minTrackPt_
 
edm::ESGetToken
< HcalChannelQuality,
HcalChannelQualityRcd > 		quality_token_
 
std::string recHitCollName_
 
edm::ESGetToken< HcalDbService,
HcalDbRecord > 		service_token_
 
edm::ESGetToken
< HcalSeverityLevelComputer,
HcalSeverityLevelComputerRcd > 		severitycomputer_token_
 
const HcalTopology * theHcalTopology_
 
double totalCalibCharge
 
double totalLasmonCharge
 
edm::EDGetTokenT
< reco::TrackCollection > 		track_token_
 
std::string trackCollName_
 
double TS4TS5EnergyThreshold_
 
std::vector< std::pair< double,
double > > 		TS4TS5LowerCut_
 
std::vector< std::pair< double,
double > > 		TS4TS5UpperCut_
 
bool useCalibDigi_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
using CacheTypes = CacheContexts< T...>
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T...>
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 

Detailed Description

Definition at line 56 of file HcalNoiseInfoProducer.h.

Constructor & Destructor Documentation

HcalNoiseInfoProducer::HcalNoiseInfoProducer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 26 of file HcalNoiseInfoProducer.cc.

References adc2fC, adc2fCHF, calibdigiHBHEthreshold_, calibdigiHBHEtimeslices_, calibdigiHFthreshold_, calibdigiHFtimeslices_, calogeometry_token_, calotower_token_, caloTowerCollName_, digiCollName_, fillCaloTowers_, fillDigis_, fillLaserMonitor_, fillRecHits_, fillTracks_, edm::ParameterSet::getParameter(), hbhedigi_token_, hbherechit_token_, HcalAcceptSeverityLevel_, hcalcalibdigi_token_, HcalRecHitFlagsToBeExcluded_, hcaltopo_token_, mps_fire::i, HLT_FULL_cff::InputTag, jet_token_, jetCollName_, laserMonCBoxList_, lasermondigi_token_, laserMonIEtaList_, laserMonIPhiList_, laserMonitorSamples_, laserMonitorTSEnd_, laserMonitorTSStart_, maxCaloTowerIEta_, maxjetindex_, maxNHF_, maxProblemRBXs_, maxTrackEta_, minHighHitE_, minLowHitE_, minR45HitE_, minRecHitE_, minTrackPt_, quality_token_, recHitCollName_, service_token_, severitycomputer_token_, AlCaHLTBitMon_QueryRunRegistry::string, track_token_, trackCollName_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

26  : algo_(iConfig) {
27  // set the parameters
28  fillDigis_ = iConfig.getParameter<bool>("fillDigis");
29  fillRecHits_ = iConfig.getParameter<bool>("fillRecHits");
30  fillCaloTowers_ = iConfig.getParameter<bool>("fillCaloTowers");
31  fillTracks_ = iConfig.getParameter<bool>("fillTracks");
32  fillLaserMonitor_ = iConfig.getParameter<bool>("fillLaserMonitor");
33 
34  maxProblemRBXs_ = iConfig.getParameter<int>("maxProblemRBXs");
35 
36  maxCaloTowerIEta_ = iConfig.getParameter<int>("maxCaloTowerIEta");
37  maxTrackEta_ = iConfig.getParameter<double>("maxTrackEta");
38  minTrackPt_ = iConfig.getParameter<double>("minTrackPt");
39 
40  digiCollName_ = iConfig.getParameter<std::string>("digiCollName");
41  recHitCollName_ = iConfig.getParameter<std::string>("recHitCollName");
42  caloTowerCollName_ = iConfig.getParameter<std::string>("caloTowerCollName");
43  trackCollName_ = iConfig.getParameter<std::string>("trackCollName");
44 
45  jetCollName_ = iConfig.getParameter<std::string>("jetCollName");
46  maxNHF_ = iConfig.getParameter<double>("maxNHF");
47  maxjetindex_ = iConfig.getParameter<int>("maxjetindex");
48  jet_token_ = consumes<reco::PFJetCollection>(edm::InputTag(jetCollName_));
49 
50  minRecHitE_ = iConfig.getParameter<double>("minRecHitE");
51  minLowHitE_ = iConfig.getParameter<double>("minLowHitE");
52  minHighHitE_ = iConfig.getParameter<double>("minHighHitE");
53 
54  minR45HitE_ = iConfig.getParameter<double>("minR45HitE");
55 
56  HcalAcceptSeverityLevel_ = iConfig.getParameter<uint32_t>("HcalAcceptSeverityLevel");
57  HcalRecHitFlagsToBeExcluded_ = iConfig.getParameter<std::vector<int>>("HcalRecHitFlagsToBeExcluded");
58 
59  // Digi threshold and time slices to use for HBHE and HF calibration digis
60  calibdigiHBHEthreshold_ = 0;
61  calibdigiHBHEtimeslices_ = std::vector<int>();
62  calibdigiHFthreshold_ = 0;
63  calibdigiHFtimeslices_ = std::vector<int>();
64 
65  calibdigiHBHEthreshold_ = iConfig.getParameter<double>("calibdigiHBHEthreshold");
66  calibdigiHBHEtimeslices_ = iConfig.getParameter<std::vector<int>>("calibdigiHBHEtimeslices");
67  calibdigiHFthreshold_ = iConfig.getParameter<double>("calibdigiHFthreshold");
68  calibdigiHFtimeslices_ = iConfig.getParameter<std::vector<int>>("calibdigiHFtimeslices");
69 
70  TS4TS5EnergyThreshold_ = iConfig.getParameter<double>("TS4TS5EnergyThreshold");
71 
72  std::vector<double> TS4TS5UpperThresholdTemp = iConfig.getParameter<std::vector<double>>("TS4TS5UpperThreshold");
73  std::vector<double> TS4TS5UpperCutTemp = iConfig.getParameter<std::vector<double>>("TS4TS5UpperCut");
74  std::vector<double> TS4TS5LowerThresholdTemp = iConfig.getParameter<std::vector<double>>("TS4TS5LowerThreshold");
75  std::vector<double> TS4TS5LowerCutTemp = iConfig.getParameter<std::vector<double>>("TS4TS5LowerCut");
76 
77  for (int i = 0; i < (int)TS4TS5UpperThresholdTemp.size() && i < (int)TS4TS5UpperCutTemp.size(); i++)
78  TS4TS5UpperCut_.push_back(std::pair<double, double>(TS4TS5UpperThresholdTemp[i], TS4TS5UpperCutTemp[i]));
79  sort(TS4TS5UpperCut_.begin(), TS4TS5UpperCut_.end());
80 
81  for (int i = 0; i < (int)TS4TS5LowerThresholdTemp.size() && i < (int)TS4TS5LowerCutTemp.size(); i++)
82  TS4TS5LowerCut_.push_back(std::pair<double, double>(TS4TS5LowerThresholdTemp[i], TS4TS5LowerCutTemp[i]));
83  sort(TS4TS5LowerCut_.begin(), TS4TS5LowerCut_.end());
84 
85  // if digis are filled, then rechits must also be filled
86  if (fillDigis_ && !fillRecHits_) {
87  fillRecHits_ = true;
88  edm::LogWarning("HCalNoiseInfoProducer") << " forcing fillRecHits to be true if fillDigis is true.\n";
89  }
90 
91  // get the fiber configuration vectors
92  laserMonCBoxList_ = iConfig.getParameter<std::vector<int>>("laserMonCBoxList");
93  laserMonIPhiList_ = iConfig.getParameter<std::vector<int>>("laserMonIPhiList");
94  laserMonIEtaList_ = iConfig.getParameter<std::vector<int>>("laserMonIEtaList");
95 
96  // check that the vectors have the same size, if not
97  // disable the laser monitor
98  if (!((laserMonCBoxList_.size() == laserMonIEtaList_.size()) &&
99  (laserMonCBoxList_.size() == laserMonIPhiList_.size()))) {
100  edm::LogWarning("MisConfiguration") << "Must provide equally sized lists for laserMonCBoxList, laserMonIEtaList, "
101  "and laserMonIPhiList. Will not fill LaserMon\n";
102  fillLaserMonitor_ = false;
103  }
104 
105  // get the integration region with defaults
106  laserMonitorTSStart_ = iConfig.getParameter<int>("laserMonTSStart");
107  laserMonitorTSEnd_ = iConfig.getParameter<int>("laserMonTSEnd");
108  laserMonitorSamples_ = iConfig.getParameter<unsigned>("laserMonSamples");
109 
110  adc2fC = std::vector<float>{
111  -0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5,
112  12.5, 13.5, 15., 17., 19., 21., 23., 25., 27., 29.5, 32.5, 35.5, 38.5,
113  42., 46., 50., 54.5, 59.5, 64.5, 59.5, 64.5, 69.5, 74.5, 79.5, 84.5, 89.5,
114  94.5, 99.5, 104.5, 109.5, 114.5, 119.5, 124.5, 129.5, 137., 147., 157., 167., 177.,
115  187., 197., 209.5, 224.5, 239.5, 254.5, 272., 292., 312., 334.5, 359.5, 384.5, 359.5,
116  384.5, 409.5, 434.5, 459.5, 484.5, 509.5, 534.5, 559.5, 584.5, 609.5, 634.5, 659.5, 684.5,
117  709.5, 747., 797., 847., 897., 947., 997., 1047., 1109.5, 1184.5, 1259.5, 1334.5, 1422.,
118  1522., 1622., 1734.5, 1859.5, 1984.5, 1859.5, 1984.5, 2109.5, 2234.5, 2359.5, 2484.5, 2609.5, 2734.5,
119  2859.5, 2984.5, 3109.5, 3234.5, 3359.5, 3484.5, 3609.5, 3797., 4047., 4297., 4547., 4797., 5047.,
120  5297., 5609.5, 5984.5, 6359.5, 6734.5, 7172., 7672., 8172., 8734.5, 9359.5, 9984.5};
121 
122  // adc -> fC for qie10, for laser monitor
123  // Taken from page 3 in
124  // https://cms-docdb.cern.ch/cgi-bin/DocDB/RetrieveFile?docid=12570&filename=QIE10_final.pdf&version=5
125  adc2fCHF = std::vector<float>{// - - - - - - - range 0 - - - - - - - -
126  //subrange0
127  1.58,
128  4.73,
129  7.88,
130  11.0,
131  14.2,
132  17.3,
133  20.5,
134  23.6,
135  26.8,
136  29.9,
137  33.1,
138  36.2,
139  39.4,
140  42.5,
141  45.7,
142  48.8,
143  //subrange1
144  53.6,
145  60.1,
146  66.6,
147  73.0,
148  79.5,
149  86.0,
150  92.5,
151  98.9,
152  105,
153  112,
154  118,
155  125,
156  131,
157  138,
158  144,
159  151,
160  //subrange2
161  157,
162  164,
163  170,
164  177,
165  186,
166  199,
167  212,
168  225,
169  238,
170  251,
171  264,
172  277,
173  289,
174  302,
175  315,
176  328,
177  //subrange3
178  341,
179  354,
180  367,
181  380,
182  393,
183  406,
184  418,
185  431,
186  444,
187  464,
188  490,
189  516,
190  542,
191  568,
192  594,
193  620,
194 
195  // - - - - - - - range 1 - - - - - - - -
196  //subrange0
197  569,
198  594,
199  619,
200  645,
201  670,
202  695,
203  720,
204  745,
205  771,
206  796,
207  821,
208  846,
209  871,
210  897,
211  922,
212  947,
213  //subrange1
214  960,
215  1010,
216  1060,
217  1120,
218  1170,
219  1220,
220  1270,
221  1320,
222  1370,
223  1430,
224  1480,
225  1530,
226  1580,
227  1630,
228  1690,
229  1740,
230  //subrange2
231  1790,
232  1840,
233  1890,
234  1940,
235  2020,
236  2120,
237  2230,
238  2330,
239  2430,
240  2540,
241  2640,
242  2740,
243  2850,
244  2950,
245  3050,
246  3150,
247  //subrange3
248  3260,
249  3360,
250  3460,
251  3570,
252  3670,
253  3770,
254  3880,
255  3980,
256  4080,
257  4240,
258  4450,
259  4650,
260  4860,
261  5070,
262  5280,
263  5490,
264 
265  // - - - - - - - range 2 - - - - - - - -
266  //subrange0
267  5080,
268  5280,
269  5480,
270  5680,
271  5880,
272  6080,
273  6280,
274  6480,
275  6680,
276  6890,
277  7090,
278  7290,
279  7490,
280  7690,
281  7890,
282  8090,
283  //subrange1
284  8400,
285  8810,
286  9220,
287  9630,
288  10000,
289  10400,
290  10900,
291  11300,
292  11700,
293  12100,
294  12500,
295  12900,
296  13300,
297  13700,
298  14100,
299  14500,
300  //subrange2
301  15000,
302  15400,
303  15800,
304  16200,
305  16800,
306  17600,
307  18400,
308  19300,
309  20100,
310  20900,
311  21700,
312  22500,
313  23400,
314  24200,
315  25000,
316  25800,
317  //subrange3
318  26600,
319  27500,
320  28300,
321  29100,
322  29900,
323  30700,
324  31600,
325  32400,
326  33200,
327  34400,
328  36100,
329  37700,
330  39400,
331  41000,
332  42700,
333  44300,
334 
335  // - - - - - - - range 3 - - - - - - - - -
336  //subrange0
337  41100,
338  42700,
339  44300,
340  45900,
341  47600,
342  49200,
343  50800,
344  52500,
345  54100,
346  55700,
347  57400,
348  59000,
349  60600,
350  62200,
351  63900,
352  65500,
353  //subrange1
354  68000,
355  71300,
356  74700,
357  78000,
358  81400,
359  84700,
360  88000,
361  91400,
362  94700,
363  98100,
364  101000,
365  105000,
366  108000,
367  111000,
368  115000,
369  118000,
370  //subrange2
371  121000,
372  125000,
373  128000,
374  131000,
375  137000,
376  145000,
377  152000,
378  160000,
379  168000,
380  176000,
381  183000,
382  191000,
383  199000,
384  206000,
385  214000,
386  222000,
387  //subrange3
388  230000,
389  237000,
390  245000,
391  253000,
392  261000,
393  268000,
394  276000,
395  284000,
396  291000,
397  302000,
398  316000,
399  329000,
400  343000,
401  356000,
402  370000,
403  384000};
404 
405  hbhedigi_token_ = consumes<HBHEDigiCollection>(edm::InputTag(digiCollName_));
406  hcalcalibdigi_token_ = consumes<HcalCalibDigiCollection>(edm::InputTag("hcalDigis"));
407  lasermondigi_token_ = consumes<QIE10DigiCollection>(iConfig.getParameter<edm::InputTag>("lasermonDigis"));
408  hbherechit_token_ = consumes<HBHERecHitCollection>(edm::InputTag(recHitCollName_));
409  calotower_token_ = consumes<CaloTowerCollection>(edm::InputTag(caloTowerCollName_));
410  track_token_ = consumes<reco::TrackCollection>(edm::InputTag(trackCollName_));
411  hcaltopo_token_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
412  service_token_ = esConsumes<HcalDbService, HcalDbRecord>();
413  quality_token_ = esConsumes<HcalChannelQuality, HcalChannelQualityRcd>(edm::ESInputTag("", "withTopo"));
414  severitycomputer_token_ = esConsumes<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd>();
415  calogeometry_token_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
416 
417  // we produce a vector of HcalNoiseRBXs
418  produces<HcalNoiseRBXCollection>();
419  // we also produce a noise summary
420  produces<HcalNoiseSummary>();
421 }
reco::HcalNoiseInfoProducer::lasermondigi_token_
edm::EDGetTokenT< QIE10DigiCollection > lasermondigi_token_
Definition: HcalNoiseInfoProducer.h:116
reco::HcalNoiseInfoProducer::calibdigiHFtimeslices_
std::vector< int > calibdigiHFtimeslices_
Definition: HcalNoiseInfoProducer.h:140
reco::HcalNoiseInfoProducer::HcalRecHitFlagsToBeExcluded_
std::vector< int > HcalRecHitFlagsToBeExcluded_
Definition: HcalNoiseInfoProducer.h:147
mps_fire.i
i
Definition: mps_fire.py:428
reco::HcalNoiseInfoProducer::TS4TS5LowerCut_
std::vector< std::pair< double, double > > TS4TS5LowerCut_
Definition: HcalNoiseInfoProducer.h:144
reco::HcalNoiseInfoProducer::severitycomputer_token_
edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd > severitycomputer_token_
Definition: HcalNoiseInfoProducer.h:124
reco::HcalNoiseInfoProducer::hcaltopo_token_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > hcaltopo_token_
Definition: HcalNoiseInfoProducer.h:121
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::HcalNoiseInfoProducer::jet_token_
edm::EDGetTokenT< reco::PFJetCollection > jet_token_
Definition: HcalNoiseInfoProducer.h:120
reco::HcalNoiseInfoProducer::quality_token_
edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd > quality_token_
Definition: HcalNoiseInfoProducer.h:123
reco::HcalNoiseInfoProducer::hcalcalibdigi_token_
edm::EDGetTokenT< HcalCalibDigiCollection > hcalcalibdigi_token_
Definition: HcalNoiseInfoProducer.h:115
reco::HcalNoiseInfoProducer::hbhedigi_token_
edm::EDGetTokenT< HBHEDigiCollection > hbhedigi_token_
Definition: HcalNoiseInfoProducer.h:114
reco::HcalNoiseInfoProducer::calogeometry_token_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > calogeometry_token_
Definition: HcalNoiseInfoProducer.h:125
reco::HcalNoiseInfoProducer::TS4TS5UpperCut_
std::vector< std::pair< double, double > > TS4TS5UpperCut_
Definition: HcalNoiseInfoProducer.h:143
reco::HcalNoiseInfoProducer::track_token_
edm::EDGetTokenT< reco::TrackCollection > track_token_
Definition: HcalNoiseInfoProducer.h:119
reco::HcalNoiseInfoProducer::hbherechit_token_
edm::EDGetTokenT< HBHERecHitCollection > hbherechit_token_
Definition: HcalNoiseInfoProducer.h:117
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
reco::HcalNoiseInfoProducer::service_token_
edm::ESGetToken< HcalDbService, HcalDbRecord > service_token_
Definition: HcalNoiseInfoProducer.h:122
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:71427
reco::HcalNoiseInfoProducer::calibdigiHBHEtimeslices_
std::vector< int > calibdigiHBHEtimeslices_
Definition: HcalNoiseInfoProducer.h:137
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
reco::HcalNoiseInfoProducer::calotower_token_
edm::EDGetTokenT< CaloTowerCollection > calotower_token_
Definition: HcalNoiseInfoProducer.h:118
HcalNoiseInfoProducer::~HcalNoiseInfoProducer ( )
override

Definition at line 423 of file HcalNoiseInfoProducer.cc.

423 {}

Member Function Documentation

void HcalNoiseInfoProducer::fillcalotwrs ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
HcalNoiseRBXArray &  array,
HcalNoiseSummary &  summary 
) const
private

Definition at line 1283 of file HcalNoiseInfoProducer.cc.

References calotower_token_, caloTowerCollName_, CaloTower::emEnergy(), HcalNoiseSummary::emenergy_, Exception, reco::HcalNoiseRBXArray::findHPD(), edm::Event::getByToken(), CaloTower::hadEnergy(), HcalNoiseSummary::hadenergy_, patZpeak::handle, CaloTower::ietaAbs(), edm::HandleBase::isValid(), maxCaloTowerIEta_, edm::errors::ProductNotFound, and trackerHitRTTI::vector.

Referenced by produce().

1286  {
1287  // get the calotowers
1288  edm::Handle<CaloTowerCollection> handle;
1289  // iEvent.getByLabel(caloTowerCollName_, handle);
1290  iEvent.getByToken(calotower_token_, handle);
1291 
1292  if (!handle.isValid()) {
1293  throw edm::Exception(edm::errors::ProductNotFound)
1294  << " could not find CaloTowerCollection named " << caloTowerCollName_ << "\n.";
1295  return;
1296  }
1297 
1298  summary.emenergy_ = summary.hadenergy_ = 0.0;
1299 
1300  // loop over all of the calotower information
1301  for (CaloTowerCollection::const_iterator it = handle->begin(); it != handle->end(); ++it) {
1302  const CaloTower& twr = (*it);
1303 
1304  // create a persistent reference to the tower
1305  edm::Ref<CaloTowerCollection> myRef(handle, it - handle->begin());
1306 
1307  // get all of the hpd's that are pointed to by the calotower
1308  std::vector<std::vector<HcalNoiseHPD>::iterator> hpditervec;
1309  array.findHPD(twr, hpditervec);
1310 
1311  // loop over the hpd's and add the reference to the RefVectors
1312  for (std::vector<std::vector<HcalNoiseHPD>::iterator>::iterator it = hpditervec.begin(); it != hpditervec.end();
1313  ++it)
1314  (*it)->calotowers_.push_back(myRef);
1315 
1316  // skip over anything with |ieta|>maxCaloTowerIEta
1317  if (twr.ietaAbs() > maxCaloTowerIEta_) {
1318  summary.emenergy_ += twr.emEnergy();
1319  summary.hadenergy_ += twr.hadEnergy();
1320  }
1321  }
1322 
1323  return;
1324 }
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:80
CaloTower::emEnergy
double emEnergy() const
Definition: CaloTower.h:130
patZpeak.handle
tuple handle
Definition: patZpeak.py:25
reco::HcalNoiseRBXArray::findHPD
std::vector< HcalNoiseHPD >::iterator findHPD(int hpdindex)
Definition: HcalNoiseRBXArray.cc:46
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
CaloTower::hadEnergy
double hadEnergy() const
Definition: CaloTower.h:131
CaloTower::ietaAbs
int ietaAbs() const
Definition: CaloTower.h:201
reco::HcalNoiseInfoProducer::calotower_token_
edm::EDGetTokenT< CaloTowerCollection > calotower_token_
Definition: HcalNoiseInfoProducer.h:118
void HcalNoiseInfoProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 425 of file HcalNoiseInfoProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

425  {
426  edm::ParameterSetDescription desc;
427  // define hit energy thesholds
428  desc.add<double>("minRecHitE", 1.5);
429  desc.add<double>("minLowHitE", 10.0);
430  desc.add<double>("minHighHitE", 25.0);
431  desc.add<double>("minR45HitE", 5.0);
432 
433  // define energy threshold for "problematic" cuts
434  desc.add<double>("pMinERatio", 25.0);
435  desc.add<double>("pMinEZeros", 5.0);
436  desc.add<double>("pMinEEMF", 10.0);
437 
438  // define energy threshold for loose/tight/high level cuts
439  desc.add<double>("minERatio", 50.0);
440  desc.add<double>("minEZeros", 10.0);
441  desc.add<double>("minEEMF", 50.0);
442 
443  // define problematic RBX
444  desc.add<double>("pMinE", 40.0);
445  desc.add<double>("pMinRatio", 0.75);
446  desc.add<double>("pMaxRatio", 0.85);
447  desc.add<int>("pMinHPDHits", 10);
448  desc.add<int>("pMinRBXHits", 20);
449  desc.add<int>("pMinHPDNoOtherHits", 7);
450  desc.add<int>("pMinZeros", 4);
451  desc.add<double>("pMinLowEHitTime", -6.0);
452  desc.add<double>("pMaxLowEHitTime", 6.0);
453  desc.add<double>("pMinHighEHitTime", -4.0);
454  desc.add<double>("pMaxHighEHitTime", 5.0);
455  desc.add<double>("pMaxHPDEMF", -0.02);
456  desc.add<double>("pMaxRBXEMF", 0.02);
457  desc.add<int>("pMinRBXRechitR45Count", 1);
458  desc.add<double>("pMinRBXRechitR45Fraction", 0.1);
459  desc.add<double>("pMinRBXRechitR45EnergyFraction", 0.1);
460 
461  // define loose noise cuts
462  desc.add<double>("lMinRatio", -999.0);
463  desc.add<double>("lMaxRatio", 999.0);
464  desc.add<int>("lMinHPDHits", 17);
465  desc.add<int>("lMinRBXHits", 999);
466  desc.add<int>("lMinHPDNoOtherHits", 10);
467  desc.add<int>("lMinZeros", 10);
468  desc.add<double>("lMinLowEHitTime", -9999.0);
469  desc.add<double>("lMaxLowEHitTime", 9999.0);
470  desc.add<double>("lMinHighEHitTime", -9999.0);
471  desc.add<double>("lMaxHighEHitTime", 9999.0);
472 
473  // define tight noise cuts
474  desc.add<double>("tMinRatio", -999.0);
475  desc.add<double>("tMaxRatio", 999.0);
476  desc.add<int>("tMinHPDHits", 16);
477  desc.add<int>("tMinRBXHits", 50);
478  desc.add<int>("tMinHPDNoOtherHits", 9);
479  desc.add<int>("tMinZeros", 8);
480  desc.add<double>("tMinLowEHitTime", -9999.0);
481  desc.add<double>("tMaxLowEHitTime", 9999.0);
482  desc.add<double>("tMinHighEHitTime", -7.0);
483  desc.add<double>("tMaxHighEHitTime", 6.0);
484 
485  // define high level noise cuts
486  desc.add<double>("hlMaxHPDEMF", -9999.0);
487  desc.add<double>("hlMaxRBXEMF", 0.01);
488 
489  // Calibration digi noise variables (used for finding laser noise events)
490  desc.add<double>("calibdigiHBHEthreshold", 15)
491  ->setComment(
492  "minimum threshold in fC of any HBHE \
493  calib digi to be counted in summary");
494  desc.add<std::vector<int>>("calibdigiHBHEtimeslices",
495  {
496  3,
497  4,
498  5,
499  6,
500  })
501  ->setComment("time slices to use when determining charge of HBHE calib digis");
502  desc.add<double>("calibdigiHFthreshold", -999)
503  ->setComment("minimum threshold in fC of any HF calib digi to be counted in summary");
504  desc.add<std::vector<int>>("calibdigiHFtimeslices",
505  {
506  0,
507  1,
508  2,
509  3,
510  4,
511  5,
512  6,
513  7,
514  8,
515  9,
516  })
517  ->setComment("time slices to use when determining charge of HF calib digis");
518 
519  // RBX-wide TS4TS5 variable
520  desc.add<double>("TS4TS5EnergyThreshold", 50);
521  desc.add<std::vector<double>>("TS4TS5UpperThreshold",
522  {
523  70,
524  90,
525  100,
526  400,
527  4000,
528  });
529  desc.add<std::vector<double>>("TS4TS5UpperCut",
530  {
531  1,
532  0.8,
533  0.75,
534  0.72,
535  0.72,
536  });
537  desc.add<std::vector<double>>("TS4TS5LowerThreshold",
538  {
539  100,
540  120,
541  150,
542  200,
543  300,
544  400,
545  500,
546  });
547  desc.add<std::vector<double>>("TS4TS5LowerCut",
548  {
549  -1,
550  -0.7,
551  -0.4,
552  -0.2,
553  -0.08,
554  0,
555  0.1,
556  });
557 
558  // rechit R45 population filter variables
559  // this comes in groups of four: (a_Count, a_Fraction, a_EnergyFraction, const)
560  // flag as noise if (count * a_count + fraction * a_fraction + energyfraction * a_energyfraction + const) > 0
561  desc.add<std::vector<double>>("lRBXRecHitR45Cuts",
562  {
563  0.0,
564  1.0,
565  0.0,
566  -0.5,
567  0.0,
568  0.0,
569  1.0,
570  -0.5,
571  })
572  ->setComment(
573  "first 4 entries : equivalent to 'fraction > 0.5' \
574  last 4 entries : equivalent to 'energy fraction > 0.5'");
575  desc.add<std::vector<double>>("tRBXRecHitR45Cuts",
576  {
577  0.0,
578  1.0,
579  0.0,
580  -0.2,
581  0.0,
582  0.0,
583  1.0,
584  -0.2,
585  })
586  ->setComment(
587  "first 4 entries : equivalent to 'fraction > 0.2' \
588  last 4 entries : equivalent to 'energy fraction > 0.2'");
589 
590  // define the channels used for laser monitoring
591  // note that the order here indicates the time order
592  // of the channels
593  desc.add<std::vector<int>>("laserMonCBoxList",
594  {
595  5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5,
596  })
597  ->setComment("time ordered list of the cBox values of laser monitor channels");
598  desc.add<std::vector<int>>("laserMonIPhiList",
599  {23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0})
600  ->setComment("time ordered list of the iPhi values of laser monitor channels");
601  desc.add<std::vector<int>>("laserMonIEtaList",
602  {
603  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
604  })
605  ->setComment("time ordered list of the iEta values of laser monitor channels");
606 
607  // boundaries for total charge integration
608  desc.add<int>("laserMonTSStart", 0)->setComment("lower bound of laser monitor charge integration window");
609  desc.add<int>("laserMonTSEnd", -1)
610  ->setComment("upper bound of laser monitor charge integration window (-1 = no bound)");
611  desc.add<unsigned>("laserMonSamples", 4)->setComment("Number of laser monitor samples to take per channel");
612 
613  // what to fill
614  desc.add<bool>("fillDigis", true);
615  desc.add<bool>("fillRecHits", true);
616  desc.add<bool>("fillCaloTowers", true);
617  desc.add<bool>("fillTracks", true);
618  desc.add<bool>("fillLaserMonitor", true);
619 
620  // maximum number of RBXs to fill
621  // if you want to record all RBXs above some energy threshold,
622  // change maxProblemRBXs to 999 and pMinE (above) to the threshold you want
623  desc.add<int>("maxProblemRBXs", 72)
624  ->setComment(
625  "maximum number of RBXs to fill. if you want to record \
626  all RBXs above some energy threshold,change maxProblemRBXs to \
627  999 and pMinE (above) to the threshold you want");
628  ;
629 
630  // parameters for calculating summary variables
631  desc.add<int>("maxCaloTowerIEta", 20);
632  desc.add<double>("maxTrackEta", 2.0);
633  desc.add<double>("minTrackPt", 1.0);
634  desc.add<double>("maxNHF", 0.9);
635  desc.add<int>("maxjetindex", 0);
636 
637  // collection names
638  desc.add<std::string>("digiCollName", "hcalDigis");
639  desc.add<std::string>("recHitCollName", "hbhereco");
640  desc.add<std::string>("caloTowerCollName", "towerMaker");
641  desc.add<std::string>("trackCollName", "generalTracks");
642  desc.add<std::string>("jetCollName", "ak4PFJets");
643  desc.add<edm::InputTag>("lasermonDigis", edm::InputTag("hcalDigis", "LASERMON"));
644 
645  // severity level
646  desc.add<unsigned int>("HcalAcceptSeverityLevel", 9);
647 
648  // which hcal calo flags to mask
649  // (HBHEIsolatedNoise=11, HBHEFlatNoise=12, HBHESpikeNoise=13,
650  // HBHETriangleNoise=14, HBHETS4TS5Noise=15, HBHENegativeNoise=27)
651  desc.add<std::vector<int>>("HcalRecHitFlagsToBeExcluded",
652  {
653  11,
654  12,
655  13,
656  14,
657  15,
658  27,
659  })
660  ->setComment(
661  "which hcal calo flags to mask (HBHEIsolatedNoise=11, \
662  HBHEFlatNoise=12, HBHESpikeNoise=13, \
663  HBHETriangleNoise=14, HBHETS4TS5Noise=15, HBHENegativeNoise=27)");
664  ;
665 
666  descriptions.add("hcalnoise", desc);
667 }
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:71427
void HcalNoiseInfoProducer::filldigis ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
HcalNoiseRBXArray &  array,
HcalNoiseSummary &  summary 
)
private

Definition at line 860 of file HcalNoiseInfoProducer.cc.

References gpuClustering::adc, HcalCoderDb::adc2fC(), adc2fC, adc2fCHF, reco::HcalNoiseRBX::allCharge_, edm::RefVector< C, T, F >::begin(), reco::HcalNoiseHPD::big5Charge_, reco::HcalNoiseHPD::bigCharge_, HcalNoiseSummary::calibCharge_, HcalNoiseSummary::calibChargegt15TS45_, HcalNoiseSummary::calibChargeTS45_, HcalNoiseSummary::calibCountgt15TS45_, HcalNoiseSummary::calibCountHF_, HcalNoiseSummary::calibCountTS45_, calibdigiHBHEthreshold_, calibdigiHBHEtimeslices_, calibdigiHFthreshold_, calibdigiHFtimeslices_, HcalCalibDetId::calibFlavor(), HcalCalibDetId::cboxChannel(), RecoTauCleanerPlugins::charge, submitPVValidationJobs::conditions, counter, isotrackApplyRegressor::df, digiCollName_, HcalSeverityLevelComputer::dropChannel(), edm::RefVector< C, T, F >::end(), Exception, fillLaserMonitor_, reco::HcalNoiseRBXArray::findHPD(), reco::HcalNoiseRBXArray::findRBX(), edm::Event::getByToken(), edm::EventSetup::getHandle(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), patZpeak::handle, hbhedigi_token_, HcalBarrel, hcalcalibdigi_token_, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary, HcalEndcap, HcalForward, HcalCalibDetId::HOCrosstalk, mps_fire::i, HBHEDataFrame::id(), HcalChannelStatus::isBitSet(), isum, edm::HandleBase::isValid(), laserMonCBoxList_, lasermondigi_token_, laserMonIEtaList_, laserMonIPhiList_, laserMonitorSamples_, laserMonitorTSEnd_, laserMonitorTSStart_, HcalNoiseSummary::lasmonCharge_, reco::HcalNoiseHPD::maxZeros_, convertSQLiteXML::ok, HcalCalibrations::pedestal(), edm::ESHandle< class >::product(), edm::errors::ProductNotFound, quality_token_, DetId::rawId(), HI_PhotonSkim_cff::rechits, reco::HcalNoiseHPD::rechits_, service_token_, severitycomputer_token_, CaloSamples::size(), totalCalibCharge, totalLasmonCharge, reco::HcalNoiseHPD::totalZeros_, and HBHEDataFrame::zsMarkAndPass().

Referenced by produce().

863  {
864  // Some initialization
865  totalCalibCharge = 0;
866  totalLasmonCharge = 0;
867 
868  // Starting with this version (updated by Jeff Temple on Dec. 6, 2012), the "TS45" names in the variables are mis-nomers. The actual time slices used are determined from the digiTimeSlices_ variable, which may not be limited to only time slices 4 and 5. For now, "TS45" name kept, because that is what is used in HcalNoiseSummary object (in GetCalibCountTS45, etc.). Likewise, the charge value in 'gt15' is now configurable, though the name remains the same. For HBHE, we track both the number of calibration channels (NcalibTS45) and the number of calibration channels above threshold (NcalibTS45gt15). For HF, we track only the number of channels above the given threshold in the given time window (NcalibHFgtX). Default for HF in 2012 is to use the full time sample with effectively no threshold (threshold=-999)
869  int NcalibTS45 = 0;
870  int NcalibTS45gt15 = 0;
871  int NcalibHFgtX = 0;
872 
873  double chargecalibTS45 = 0;
874  double chargecalibgt15TS45 = 0;
875 
876  // get the conditions and channel quality
877  edm::ESHandle<HcalDbService> conditions = iSetup.getHandle(service_token_);
878  edm::ESHandle<HcalChannelQuality> qualhandle = iSetup.getHandle(quality_token_);
879  const HcalChannelQuality* myqual = qualhandle.product();
880  edm::ESHandle<HcalSeverityLevelComputer> mycomputer = iSetup.getHandle(severitycomputer_token_);
881  const HcalSeverityLevelComputer* mySeverity = mycomputer.product();
882 
883  // get the digis
884  edm::Handle<HBHEDigiCollection> handle;
885  // iEvent.getByLabel(digiCollName_, handle);
886  iEvent.getByToken(hbhedigi_token_, handle);
887 
888  if (!handle.isValid()) {
889  throw edm::Exception(edm::errors::ProductNotFound)
890  << " could not find HBHEDigiCollection named " << digiCollName_ << "\n.";
891  return;
892  }
893 
894  // loop over all of the digi information
895  for (HBHEDigiCollection::const_iterator it = handle->begin(); it != handle->end(); ++it) {
896  const HBHEDataFrame& digi = (*it);
897  HcalDetId cell = digi.id();
898  DetId detcell = (DetId)cell;
899 
900  // check on cells to be ignored and dropped
901  const HcalChannelStatus* mydigistatus = myqual->getValues(detcell.rawId());
902  if (mySeverity->dropChannel(mydigistatus->getValue()))
903  continue;
904  if (digi.zsMarkAndPass())
905  continue;
906  // Drop if exclude bit set
907  if (mydigistatus->isBitSet(HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary))
908  continue;
909 
910  // get the calibrations and coder
911  const HcalCalibrations& calibrations = conditions->getHcalCalibrations(cell);
912  const HcalQIECoder* channelCoder = conditions->getHcalCoder(cell);
913  const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
914  HcalCoderDb coder(*channelCoder, *shape);
915 
916  // match the digi to an rbx and hpd
917  HcalNoiseRBX& rbx = (*array.findRBX(digi));
918  HcalNoiseHPD& hpd = (*array.findHPD(digi));
919 
920  // determine if the digi is one the highest energy hits in the HPD
921  // this works because the rechits are sorted by energy (see fillrechits() below)
922  bool isBig = false, isBig5 = false, isRBX = false;
923  int counter = 0;
924  edm::RefVector<HBHERecHitCollection>& rechits = hpd.rechits_;
925  for (edm::RefVector<HBHERecHitCollection>::const_iterator rit = rechits.begin(); rit != rechits.end();
926  ++rit, ++counter) {
927  const HcalDetId& detid = (*rit)->idFront();
928  if (DetId(detid) == digi.id()) {
929  if (counter == 0)
930  isBig = isBig5 = true; // digi is also the highest energy rechit
931  if (counter < 5)
932  isBig5 = true; // digi is one of 5 highest energy rechits
933  isRBX = true;
934  }
935  }
936 
937  // loop over each of the digi's time slices
938  int totalzeros = 0;
939  CaloSamples tool;
940  coder.adc2fC(digi, tool);
941  for (int ts = 0; ts < tool.size(); ++ts) {
942  // count zero's
943  if (digi[ts].adc() == 0) {
944  ++hpd.totalZeros_;
945  ++totalzeros;
946  }
947 
948  // get the fC's
949  double corrfc = tool[ts] - calibrations.pedestal(digi[ts].capid());
950 
951  // fill the relevant digi arrays
952  if (isBig)
953  hpd.bigCharge_[ts] += corrfc;
954  if (isBig5)
955  hpd.big5Charge_[ts] += corrfc;
956  if (isRBX)
957  rbx.allCharge_[ts] += corrfc;
958  }
959 
960  // record the maximum number of zero's found
961  if (totalzeros > hpd.maxZeros_)
962  hpd.maxZeros_ = totalzeros;
963  }
964 
965  // get the calibration digis
966  edm::Handle<HcalCalibDigiCollection> hCalib;
967  // iEvent.getByLabel("hcalDigis", hCalib);
968  iEvent.getByToken(hcalcalibdigi_token_, hCalib);
969 
970  // get the lasermon digis
971  edm::Handle<QIE10DigiCollection> hLasermon;
972  iEvent.getByToken(lasermondigi_token_, hLasermon);
973 
974  // get total charge in calibration channels
975  if (hCalib.isValid() == true) {
976  for (HcalCalibDigiCollection::const_iterator digi = hCalib->begin(); digi != hCalib->end(); digi++) {
977  if (digi->id().hcalSubdet() == 0)
978  continue;
979 
980  for (unsigned i = 0; i < (unsigned)digi->size(); i++)
981  totalCalibCharge = totalCalibCharge + adc2fC[digi->sample(i).adc() & 0xff];
982 
983  HcalCalibDetId myid = (HcalCalibDetId)digi->id();
984  if (myid.calibFlavor() == HcalCalibDetId::HOCrosstalk)
985  continue; // ignore HOCrosstalk channels
986  if (digi->zsMarkAndPass())
987  continue; // skip "mark-and-pass" channels when computing charge in calib channels
988 
989  if (digi->id().hcalSubdet() == HcalForward) // check HF
990  {
991  double sumChargeHF = 0;
992  for (unsigned int i = 0; i < calibdigiHFtimeslices_.size(); ++i) {
993  // skip unphysical time slices
994  if (calibdigiHFtimeslices_[i] < 0 || calibdigiHFtimeslices_[i] > digi->size())
995  continue;
996  sumChargeHF += adc2fC[digi->sample(calibdigiHFtimeslices_[i]).adc() & 0xff];
997  }
998  if (sumChargeHF > calibdigiHFthreshold_)
999  ++NcalibHFgtX;
1000  } // end of HF check
1001  else if (digi->id().hcalSubdet() == HcalBarrel || digi->id().hcalSubdet() == HcalEndcap) // now check HBHE
1002  {
1003  double sumChargeHBHE = 0;
1004  for (unsigned int i = 0; i < calibdigiHBHEtimeslices_.size(); ++i) {
1005  // skip unphysical time slices
1006  if (calibdigiHBHEtimeslices_[i] < 0 || calibdigiHBHEtimeslices_[i] > digi->size())
1007  continue;
1008  sumChargeHBHE += adc2fC[digi->sample(calibdigiHBHEtimeslices_[i]).adc() & 0xff];
1009  }
1010  ++NcalibTS45;
1011  chargecalibTS45 += sumChargeHBHE;
1012  if (sumChargeHBHE > calibdigiHBHEthreshold_) {
1013  ++NcalibTS45gt15;
1014  chargecalibgt15TS45 += sumChargeHBHE;
1015  }
1016  } // end of HBHE check
1017  } // loop on HcalCalibDigiCollection
1018 
1019  } // if (hCalib.isValid()==true)
1020  if (fillLaserMonitor_ && (hLasermon.isValid() == true)) {
1021  int icombts = -1;
1022  float max_charge = 0;
1023  int max_ts = -1;
1024  std::vector<float> comb_charge;
1025 
1026  unsigned nch = laserMonCBoxList_.size();
1027  // loop over channels in the order provided
1028  for (unsigned ich = 0; ich < nch; ++ich) {
1029  int cboxch = laserMonCBoxList_[ich];
1030  int iphi = laserMonIPhiList_[ich];
1031  int ieta = laserMonIEtaList_[ich];
1032 
1033  // loop over digis, find the digi that matches this channel
1034  for (const QIE10DataFrame df : (*hLasermon)) {
1035  HcalCalibDetId calibId(df.id());
1036 
1037  int ch_cboxch = calibId.cboxChannel();
1038  int ch_iphi = calibId.iphi();
1039  int ch_ieta = calibId.ieta();
1040 
1041  if (cboxch == ch_cboxch && iphi == ch_iphi && ieta == ch_ieta) {
1042  unsigned ts_size = df.samples();
1043 
1044  // loop over time slices
1045  for (unsigned its = 0; its < ts_size; ++its) {
1046  // if we are on the last channel, use all the data
1047  // otherwise only take the unique samples
1048  if (((ich + 1) < nch) && its >= laserMonitorSamples_)
1049  continue;
1050 
1051  bool ok = df[its].ok();
1052  int adc = df[its].adc();
1053 
1054  icombts++;
1055  // apply integration limits
1056  if (icombts < laserMonitorTSStart_)
1057  continue;
1058  if (laserMonitorTSEnd_ > 0 && icombts > laserMonitorTSEnd_)
1059  continue;
1060 
1061  if (ok && adc >= 0) { // protection against QIE reset or bad ADC values
1062 
1063  float charge = adc2fCHF[adc];
1064  if (charge > max_charge) {
1065  max_charge = charge;
1066  max_ts = icombts;
1067  }
1068 
1069  comb_charge.push_back(charge);
1070  } // if( ok && adc >= 0 )
1071  } // loop over time slices
1072  } // if( cboxch == ch_cboxch && iphi == ch_iphi && ieta == ch_ieta )
1073  } // loop over digi collection
1074  } // loop over channel list
1075 
1076  // do not continue with the calculation
1077  // if the vector was not filled
1078  if (comb_charge.empty()) {
1079  totalLasmonCharge = -1;
1080  } else {
1081  // integrate from +- 3 TS around the time sample
1082  // having the maximum charge
1083  int start_ts = max_ts - 3;
1084  int end_ts = max_ts + 3;
1085 
1086  // Change the integration limits
1087  // if outside of the range
1088  if (start_ts < 0)
1089  start_ts = 0;
1090  if (end_ts >= int(comb_charge.size()))
1091  end_ts = comb_charge.size() - 1;
1092 
1093  for (int isum = start_ts; isum <= end_ts; ++isum) {
1094  totalLasmonCharge += comb_charge[isum];
1095  }
1096  }
1097  } // if( fillLaserMonitor_ && (hLasermon.isValid() == true) )
1098 
1099  summary.calibCharge_ = totalCalibCharge;
1100  summary.lasmonCharge_ = totalLasmonCharge;
1101  summary.calibCountTS45_ = NcalibTS45;
1102  summary.calibCountgt15TS45_ = NcalibTS45gt15;
1103  summary.calibChargeTS45_ = chargecalibTS45;
1104  summary.calibChargegt15TS45_ = chargecalibgt15TS45;
1105  summary.calibCountHF_ = NcalibHFgtX;
1106 
1107  return;
1108 }
reco::HcalNoiseInfoProducer::lasermondigi_token_
edm::EDGetTokenT< QIE10DigiCollection > lasermondigi_token_
Definition: HcalNoiseInfoProducer.h:116
reco::HcalNoiseInfoProducer::calibdigiHFtimeslices_
std::vector< int > calibdigiHFtimeslices_
Definition: HcalNoiseInfoProducer.h:140
mps_fire.i
i
Definition: mps_fire.py:428
HcalCalibDetId::calibFlavor
CalibDetType calibFlavor() const
get the flavor of this calibration detid
Definition: HcalCalibDetId.h:77
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
reco::HcalNoiseInfoProducer::severitycomputer_token_
edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd > severitycomputer_token_
Definition: HcalNoiseInfoProducer.h:124
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
reco::HcalNoiseHPD::bigCharge_
std::vector< float > bigCharge_
Definition: HcalNoiseHPD.h:137
reco::HcalNoiseHPD::big5Charge_
std::vector< float > big5Charge_
Definition: HcalNoiseHPD.h:138
HBHEDataFrame::id
constexpr const HcalDetId & id() const
Definition: HBHEDataFrame.h:23
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
isum
double isum
Definition: MuScleFitUtils.cc:77
edm::RefVector::end
const_iterator end() const
Termination of iteration.
Definition: RefVector.h:228
edm::RefVector::begin
const_iterator begin() const
Initialize an iterator over the RefVector.
Definition: RefVector.h:223
reco::HcalNoiseHPD::rechits_
edm::RefVector< HBHERecHitCollection > rechits_
Definition: HcalNoiseHPD.h:141
reco::HcalNoiseRBXArray::findRBX
HcalNoiseRBXArray::iterator findRBX(int rbxindex)
Definition: HcalNoiseRBXArray.cc:87
reco::HcalNoiseInfoProducer::quality_token_
edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd > quality_token_
Definition: HcalNoiseInfoProducer.h:123
HBHEDataFrame::zsMarkAndPass
constexpr bool zsMarkAndPass() const
was ZS MarkAndPass?
Definition: HBHEDataFrame.h:31
reco::HcalNoiseInfoProducer::hcalcalibdigi_token_
edm::EDGetTokenT< HcalCalibDigiCollection > hcalcalibdigi_token_
Definition: HcalNoiseInfoProducer.h:115
patZpeak.handle
tuple handle
Definition: patZpeak.py:25
reco::HcalNoiseRBXArray::findHPD
std::vector< HcalNoiseHPD >::iterator findHPD(int hpdindex)
Definition: HcalNoiseRBXArray.cc:46
HcalSeverityLevelComputer::dropChannel
bool dropChannel(const uint32_t &mystatus) const
Definition: HcalSeverityLevelComputer.cc:395
HBHEDataFrame::sample
constexpr HcalQIESample const & sample(int i) const
access a sample
Definition: HBHEDataFrame.h:40
reco::HcalNoiseInfoProducer::hbhedigi_token_
edm::EDGetTokenT< HBHEDigiCollection > hbhedigi_token_
Definition: HcalNoiseInfoProducer.h:114
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HcalQIESample::adc
constexpr int adc() const
get the ADC sample
Definition: HcalQIESample.h:43
HBHEDataFrame::size
constexpr int size() const
total number of samples in the digi
Definition: HBHEDataFrame.h:27
reco::HcalNoiseRBX::allCharge_
std::vector< float > allCharge_
Definition: HcalNoiseRBX.h:108
DetId
Definition: DetId.h:17
HcalCalibrations::pedestal
constexpr double pedestal(int fCapId) const
get pedestal for capid=0..3
Definition: HcalCalibrations.h:46
CaloSamples::size
int size() const
get the size
Definition: CaloSamples.h:24
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
reco::HcalNoiseInfoProducer::service_token_
edm::ESGetToken< HcalDbService, HcalDbRecord > service_token_
Definition: HcalNoiseInfoProducer.h:122
HcalChannelStatus::isBitSet
bool isBitSet(unsigned int bitnumber) const
Definition: HcalChannelStatus.h:53
HcalCalibDetId::cboxChannel
int cboxChannel() const
get the calibration box channel (if relevant)
Definition: HcalCalibDetId.cc:86
HI_PhotonSkim_cff.rechits
tuple rechits
Definition: HI_PhotonSkim_cff.py:76
counter
static std::atomic< unsigned int > counter
Definition: SharedResourceNames.cc:18
reco::HcalNoiseInfoProducer::calibdigiHBHEtimeslices_
std::vector< int > calibdigiHBHEtimeslices_
Definition: HcalNoiseInfoProducer.h:137
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:151
HcalChannelStatus::getValue
uint32_t getValue() const
Definition: HcalChannelStatus.h:60
gpuClustering::adc
uint16_t *__restrict__ uint16_t const *__restrict__ adc
Definition: gpuClusterChargeCut.h:21
void HcalNoiseInfoProducer::filljetinfo ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
HcalNoiseSummary &  summary 
) const
private

Definition at line 1327 of file HcalNoiseInfoProducer.cc.

References edm::Event::getByToken(), HcalNoiseSummary::goodJetFoundInLowBVRegion_, edm::HandleBase::isValid(), metsig::jet, jet_token_, jetCollName_, maxjetindex_, maxNHF_, and reco::btau::neutralHadronEnergyFraction.

Referenced by produce().

1329  {
1330  bool goodJetFoundInLowBVRegion = false; // checks whether a jet is in
1331  // a low BV region, where false
1332  // noise flagging rate is higher.
1333  if (!jetCollName_.empty()) {
1334  edm::Handle<reco::PFJetCollection> pfjet_h;
1335  iEvent.getByToken(jet_token_, pfjet_h);
1336 
1337  if (pfjet_h.isValid()) {
1338  int jetindex = 0;
1339  for (reco::PFJetCollection::const_iterator jet = pfjet_h->begin(); jet != pfjet_h->end(); ++jet) {
1340  if (jetindex > maxjetindex_)
1341  break; // only look at jets with
1342  // indices up to maxjetindex_
1343 
1344  // Check whether jet is in low-BV region (0<eta<1.4, -1.8<phi<-1.4)
1345  if (jet->eta() > 0.0 && jet->eta() < 1.4 && jet->phi() > -1.8 && jet->phi() < -1.4) {
1346  // Look for a good jet in low BV region;
1347  // if found, we will keep event
1348  if (maxNHF_ < 0.0 || jet->neutralHadronEnergyFraction() < maxNHF_) {
1349  goodJetFoundInLowBVRegion = true;
1350  break;
1351  }
1352  }
1353  ++jetindex;
1354  }
1355  }
1356  }
1357 
1358  summary.goodJetFoundInLowBVRegion_ = goodJetFoundInLowBVRegion;
1359 }
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
reco::HcalNoiseInfoProducer::jet_token_
edm::EDGetTokenT< reco::PFJetCollection > jet_token_
Definition: HcalNoiseInfoProducer.h:120
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
void HcalNoiseInfoProducer::fillOtherSummaryVariables ( HcalNoiseSummary &  summary,
const CommonHcalNoiseRBXData &  data 
) const
private

Definition at line 756 of file HcalNoiseInfoProducer.cc.

References algo_, HcalNoiseSummary::cnthit10_, HcalNoiseSummary::cnthit25_, CommonHcalNoiseRBXData::e10ts(), CommonHcalNoiseRBXData::e2ts(), HcalNoiseSummary::filterstatus_, HcalNoiseSummary::hasBadRBXRechitR45Loose_, HcalNoiseSummary::hasBadRBXRechitR45Tight_, HcalNoiseSummary::hasBadRBXTS4TS5_, CommonHcalNoiseRBXData::highEHitTimeSqrd(), HcalNoiseSummary::hlnoisetwrs_, CommonHcalNoiseRBXData::HPDEMF(), if(), join(), HcalNoiseSummary::loosenoisetwrs_, CommonHcalNoiseRBXData::lowEHitTimeSqrd(), HcalNoiseSummary::max10_, HcalNoiseSummary::max10GeVHitTime(), HcalNoiseSummary::max25_, HcalNoiseSummary::max25GeVHitTime(), HcalNoiseSummary::maxe10ts_, HcalNoiseSummary::maxE2Over10TS(), HcalNoiseSummary::maxe2ts_, CommonHcalNoiseRBXData::maxHighEHitTime(), HcalNoiseSummary::maxHPDHits(), HcalNoiseSummary::maxhpdhits_, HcalNoiseSummary::maxhpdhitsnoother_, HcalNoiseSummary::maxHPDNoOtherHits(), CommonHcalNoiseRBXData::maxLowEHitTime(), HcalNoiseSummary::maxRBXHits(), HcalNoiseSummary::maxrbxhits_, HcalNoiseSummary::maxZeros(), HcalNoiseSummary::maxzeros_, HcalNoiseSummary::min10_, HcalNoiseSummary::min10GeVHitTime(), HcalNoiseSummary::min25_, HcalNoiseSummary::min25GeVHitTime(), HcalNoiseSummary::mine10ts_, HcalNoiseSummary::minE2Over10TS(), HcalNoiseSummary::mine2ts_, CommonHcalNoiseRBXData::minHighEHitTime(), HcalNoiseSummary::minHPDEMF(), HcalNoiseSummary::minhpdemf_, CommonHcalNoiseRBXData::minLowEHitTime(), HcalNoiseSummary::minRBXEMF(), HcalNoiseSummary::minrbxemf_, CommonHcalNoiseRBXData::numHighEHits(), CommonHcalNoiseRBXData::numHPDHits(), CommonHcalNoiseRBXData::numHPDNoOtherHits(), CommonHcalNoiseRBXData::numLowEHits(), CommonHcalNoiseRBXData::numRBXHits(), CommonHcalNoiseRBXData::numZeros(), HcalNoiseAlgo::passEMFThreshold(), HcalNoiseAlgo::passHighLevelNoiseFilter(), HcalNoiseAlgo::passLooseHits(), HcalNoiseAlgo::passLooseNoiseFilter(), HcalNoiseAlgo::passLooseRatio(), HcalNoiseAlgo::passLooseRBXRechitR45(), HcalNoiseAlgo::passLooseTiming(), HcalNoiseAlgo::passLooseZeros(), HcalNoiseAlgo::passRatioThreshold(), HcalNoiseAlgo::passTightHits(), HcalNoiseAlgo::passTightNoiseFilter(), HcalNoiseAlgo::passTightRatio(), HcalNoiseAlgo::passTightRBXRechitR45(), HcalNoiseAlgo::passTightTiming(), HcalNoiseAlgo::passTightZeros(), CommonHcalNoiseRBXData::PassTS4TS5(), HcalNoiseAlgo::passZerosThreshold(), CommonHcalNoiseRBXData::ratio(), CommonHcalNoiseRBXData::RBXEMF(), CommonHcalNoiseRBXData::rbxTowers(), HcalNoiseSummary::rms10_, HcalNoiseSummary::rms25_, HcalNoiseSummary::tightnoisetwrs_, and CommonHcalNoiseRBXData::validRatio().

Referenced by produce().

757  {
758  // charge ratio
759  if (algo_.passRatioThreshold(data) && data.validRatio()) {
760  if (data.ratio() < summary.minE2Over10TS()) {
761  summary.mine2ts_ = data.e2ts();
762  summary.mine10ts_ = data.e10ts();
763  }
764  if (data.ratio() > summary.maxE2Over10TS()) {
765  summary.maxe2ts_ = data.e2ts();
766  summary.maxe10ts_ = data.e10ts();
767  }
768  }
769 
770  // ADC zeros
771  if (algo_.passZerosThreshold(data)) {
772  if (data.numZeros() > summary.maxZeros()) {
773  summary.maxzeros_ = data.numZeros();
774  }
775  }
776 
777  // hits count
778  if (data.numHPDHits() > summary.maxHPDHits()) {
779  summary.maxhpdhits_ = data.numHPDHits();
780  }
781  if (data.numRBXHits() > summary.maxRBXHits()) {
782  summary.maxrbxhits_ = data.numRBXHits();
783  }
784  if (data.numHPDNoOtherHits() > summary.maxHPDNoOtherHits()) {
785  summary.maxhpdhitsnoother_ = data.numHPDNoOtherHits();
786  }
787 
788  // TS4TS5
789  if (data.PassTS4TS5() == false)
790  summary.hasBadRBXTS4TS5_ = true;
791 
792  if (algo_.passLooseRBXRechitR45(data) == false)
793  summary.hasBadRBXRechitR45Loose_ = true;
794  if (algo_.passTightRBXRechitR45(data) == false)
795  summary.hasBadRBXRechitR45Tight_ = true;
796 
797  // hit timing
798  if (data.minLowEHitTime() < summary.min10GeVHitTime()) {
799  summary.min10_ = data.minLowEHitTime();
800  }
801  if (data.maxLowEHitTime() > summary.max10GeVHitTime()) {
802  summary.max10_ = data.maxLowEHitTime();
803  }
804  summary.rms10_ += data.lowEHitTimeSqrd();
805  summary.cnthit10_ += data.numLowEHits();
806  if (data.minHighEHitTime() < summary.min25GeVHitTime()) {
807  summary.min25_ = data.minHighEHitTime();
808  }
809  if (data.maxHighEHitTime() > summary.max25GeVHitTime()) {
810  summary.max25_ = data.maxHighEHitTime();
811  }
812  summary.rms25_ += data.highEHitTimeSqrd();
813  summary.cnthit25_ += data.numHighEHits();
814 
815  // EMF
816  if (algo_.passEMFThreshold(data)) {
817  if (summary.minHPDEMF() > data.HPDEMF()) {
818  summary.minhpdemf_ = data.HPDEMF();
819  }
820  if (summary.minRBXEMF() > data.RBXEMF()) {
821  summary.minrbxemf_ = data.RBXEMF();
822  }
823  }
824 
825  // summary flag
826  if (!algo_.passLooseRatio(data))
827  summary.filterstatus_ |= 0x1;
828  if (!algo_.passLooseHits(data))
829  summary.filterstatus_ |= 0x2;
830  if (!algo_.passLooseZeros(data))
831  summary.filterstatus_ |= 0x4;
832  if (!algo_.passLooseTiming(data))
833  summary.filterstatus_ |= 0x8;
834 
835  if (!algo_.passTightRatio(data))
836  summary.filterstatus_ |= 0x100;
837  if (!algo_.passTightHits(data))
838  summary.filterstatus_ |= 0x200;
839  if (!algo_.passTightZeros(data))
840  summary.filterstatus_ |= 0x400;
841  if (!algo_.passTightTiming(data))
842  summary.filterstatus_ |= 0x800;
843 
844  if (!algo_.passHighLevelNoiseFilter(data))
845  summary.filterstatus_ |= 0x10000;
846 
847  // summary refvectors
848  JoinCaloTowerRefVectorsWithoutDuplicates join;
849  if (!algo_.passLooseNoiseFilter(data))
850  join(summary.loosenoisetwrs_, data.rbxTowers());
851  if (!algo_.passTightNoiseFilter(data))
852  join(summary.tightnoisetwrs_, data.rbxTowers());
853  if (!algo_.passHighLevelNoiseFilter(data))
854  join(summary.hlnoisetwrs_, data.rbxTowers());
855 
856  return;
857 }
CommonHcalNoiseRBXData::e10ts
double e10ts(void) const
Definition: HcalNoiseAlgo.h:27
CommonHcalNoiseRBXData::HPDEMF
double HPDEMF(void) const
Definition: HcalNoiseAlgo.h:42
CommonHcalNoiseRBXData::e2ts
double e2ts(void) const
Definition: HcalNoiseAlgo.h:26
HcalNoiseSummary::hlnoisetwrs_
edm::RefVector< CaloTowerCollection > hlnoisetwrs_
Definition: HcalNoiseSummary.h:220
CommonHcalNoiseRBXData::lowEHitTimeSqrd
double lowEHitTimeSqrd(void) const
Definition: HcalNoiseAlgo.h:35
HcalNoiseAlgo::passTightRatio
bool passTightRatio(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:299
CommonHcalNoiseRBXData::minLowEHitTime
double minLowEHitTime(void) const
Definition: HcalNoiseAlgo.h:33
HcalNoiseAlgo::passLooseRatio
bool passLooseRatio(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:244
HcalNoiseAlgo::passTightRBXRechitR45
bool passTightRBXRechitR45(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:339
CommonHcalNoiseRBXData::maxHighEHitTime
double maxHighEHitTime(void) const
Definition: HcalNoiseAlgo.h:38
HcalNoiseAlgo::passTightTiming
bool passTightTiming(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:327
CommonHcalNoiseRBXData::maxLowEHitTime
double maxLowEHitTime(void) const
Definition: HcalNoiseAlgo.h:34
CommonHcalNoiseRBXData::numLowEHits
int numLowEHits(void) const
Definition: HcalNoiseAlgo.h:36
CommonHcalNoiseRBXData::ratio
double ratio(void) const
Definition: HcalNoiseAlgo.h:25
HcalNoiseSummary::maxE2Over10TS
float maxE2Over10TS(void) const
Definition: HcalNoiseSummary.cc:142
HcalNoiseSummary::loosenoisetwrs_
edm::RefVector< CaloTowerCollection > loosenoisetwrs_
Definition: HcalNoiseSummary.h:218
HcalNoiseAlgo::passLooseRBXRechitR45
bool passLooseRBXRechitR45(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:284
HcalNoiseSummary::min25GeVHitTime
float min25GeVHitTime(void) const
Definition: HcalNoiseSummary.cc:122
HcalNoiseAlgo::passHighLevelNoiseFilter
bool passHighLevelNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:234
CommonHcalNoiseRBXData::RBXEMF
double RBXEMF(void) const
Definition: HcalNoiseAlgo.h:41
CommonHcalNoiseRBXData::numRBXHits
int numRBXHits(void) const
Definition: HcalNoiseAlgo.h:30
HcalNoiseSummary::min10GeVHitTime
float min10GeVHitTime(void) const
Definition: HcalNoiseSummary.cc:116
CommonHcalNoiseRBXData::highEHitTimeSqrd
double highEHitTimeSqrd(void) const
Definition: HcalNoiseAlgo.h:39
HcalNoiseSummary::minE2Over10TS
float minE2Over10TS(void) const
Definition: HcalNoiseSummary.cc:136
HcalNoiseAlgo::passLooseZeros
bool passLooseZeros(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:264
HcalNoiseSummary::maxRBXHits
int maxRBXHits(void) const
Definition: HcalNoiseSummary.cc:150
HcalNoiseAlgo::passTightZeros
bool passTightZeros(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:319
if
if(conf_.getParameter< bool >("UseStripCablingDB"))
Definition: SiStripLAProfileBooker.cc:72
CommonHcalNoiseRBXData::numHighEHits
int numHighEHits(void) const
Definition: HcalNoiseAlgo.h:40
HcalNoiseAlgo::passRatioThreshold
bool passRatioThreshold(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:354
CommonHcalNoiseRBXData::numHPDNoOtherHits
int numHPDNoOtherHits(void) const
Definition: HcalNoiseAlgo.h:31
CommonHcalNoiseRBXData::numZeros
int numZeros(void) const
Definition: HcalNoiseAlgo.h:32
HcalNoiseAlgo::passZerosThreshold
bool passZerosThreshold(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:358
HcalNoiseSummary::minHPDEMF
float minHPDEMF(void) const
Definition: HcalNoiseSummary.cc:152
HcalNoiseAlgo::passLooseNoiseFilter
bool passLooseNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:224
CommonHcalNoiseRBXData::rbxTowers
edm::RefVector< CaloTowerCollection > rbxTowers(void) const
Definition: HcalNoiseAlgo.h:44
HcalNoiseAlgo::passLooseHits
bool passLooseHits(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:254
join
static std::string join(char **cmd)
Definition: RemoteFile.cc:19
HcalNoiseSummary::maxZeros
int maxZeros(void) const
Definition: HcalNoiseSummary.cc:144
HcalNoiseSummary::tightnoisetwrs_
edm::RefVector< CaloTowerCollection > tightnoisetwrs_
Definition: HcalNoiseSummary.h:219
HcalNoiseAlgo::passLooseTiming
bool passLooseTiming(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:272
HcalNoiseAlgo::passTightNoiseFilter
bool passTightNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:229
HcalNoiseSummary::maxHPDNoOtherHits
int maxHPDNoOtherHits(void) const
Definition: HcalNoiseSummary.cc:148
CommonHcalNoiseRBXData::minHighEHitTime
double minHighEHitTime(void) const
Definition: HcalNoiseAlgo.h:37
HcalNoiseSummary::minRBXEMF
float minRBXEMF(void) const
Definition: HcalNoiseSummary.cc:154
CommonHcalNoiseRBXData::PassTS4TS5
bool PassTS4TS5(void) const
Definition: HcalNoiseAlgo.h:43
HcalNoiseAlgo::passEMFThreshold
bool passEMFThreshold(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:362
HcalNoiseSummary::maxHPDHits
int maxHPDHits(void) const
Definition: HcalNoiseSummary.cc:146
CommonHcalNoiseRBXData::numHPDHits
int numHPDHits(void) const
Definition: HcalNoiseAlgo.h:29
HcalNoiseSummary::max25GeVHitTime
float max25GeVHitTime(void) const
Definition: HcalNoiseSummary.cc:124
HcalNoiseSummary::max10GeVHitTime
float max10GeVHitTime(void) const
Definition: HcalNoiseSummary.cc:118
CommonHcalNoiseRBXData::validRatio
bool validRatio(void) const
Definition: HcalNoiseAlgo.h:28
HcalNoiseAlgo::passTightHits
bool passTightHits(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:309
void HcalNoiseInfoProducer::fillrechits ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
HcalNoiseRBXArray &  array,
HcalNoiseSummary &  summary 
) const
private

Definition at line 1111 of file HcalNoiseInfoProducer.cc.

References HBHERecHit::auxPhase1(), calogeometry_token_, HcalNoiseSummary::energyInLaserRegion_, HcalNoiseSummary::energyInNonLaserRegion_, HBHERecHit::eraw(), Exception, reco::HcalNoiseRBXArray::findHPD(), CaloRecHit::flags(), HcalNoiseSummary::flatnoisee_, HcalNoiseSummary::flatnoiseet_, CaloRecHitAuxSetter::getBit(), edm::Event::getByToken(), edm::EventSetup::getHandle(), CaloGeometry::getPosition(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), runTauDisplay::gp, patZpeak::handle, HcalCaloFlagLabels::HBHEFlatNoise, HcalCaloFlagLabels::HBHEIsolatedNoise, HcalCaloFlagLabels::HBHENegativeNoise, hbherechit_token_, HcalCaloFlagLabels::HBHESpikeNoise, HcalCaloFlagLabels::HBHETriangleNoise, HcalCaloFlagLabels::HBHETS4TS5Noise, HcalAcceptSeverityLevel_, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummaryR45, HcalRecHitFlagsToBeExcluded_, HcalNoiseSummary::hitsInLaserRegion_, HcalNoiseSummary::hitsInNonLaserRegion_, mps_fire::i, HBHERecHit::id(), HBHERecHit::idFront(), HcalChannelStatus::isBitSet(), HcalNoiseSummary::isolnoisee_, HcalNoiseSummary::isolnoiseet_, edm::HandleBase::isValid(), HcalNoiseSummary::negativenoisee_, HcalNoiseSummary::negativenoiseet_, HcalNoiseSummary::nflatnoise_, HcalNoiseSummary::nisolnoise_, HcalNoiseSummary::nnegativenoise_, HcalNoiseSummary::nspikenoise_, HcalNoiseSummary::ntrianglenoise_, HcalNoiseSummary::nts4ts5noise_, HBHERecHitAuxSetter::OFF_COMBINED, HBHERecHitAuxSetter::OFF_TDC_TIME, edm::ESHandle< class >::product(), edm::errors::ProductNotFound, quality_token_, recHitCollName_, HcalNoiseSummary::rechitCount15_, HcalNoiseSummary::rechitCount_, HcalNoiseSummary::rechitEnergy15_, HcalNoiseSummary::rechitEnergy_, HcalSeverityLevelComputer::recoveredRecHit(), reco::HcalNoiseHPD::refrechitset_, severitycomputer_token_, HcalNoiseSummary::spikenoisee_, HcalNoiseSummary::spikenoiseet_, HcalNoiseSummary::trianglenoisee_, HcalNoiseSummary::trianglenoiseet_, HcalNoiseSummary::ts4ts5noisee_, and HcalNoiseSummary::ts4ts5noiseet_.

Referenced by produce().

1114  {
1115  // get the HCAL channel status map
1116  edm::ESHandle<HcalChannelQuality> hcalChStatus = iSetup.getHandle(quality_token_);
1117  const HcalChannelQuality* dbHcalChStatus = hcalChStatus.product();
1118 
1119  // get the severity level computer
1120  edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl = iSetup.getHandle(severitycomputer_token_);
1121  const HcalSeverityLevelComputer* hcalSevLvlComputer = hcalSevLvlComputerHndl.product();
1122 
1123  // get the calo geometry
1124  edm::ESHandle<CaloGeometry> pG = iSetup.getHandle(calogeometry_token_);
1125  const CaloGeometry* geo = pG.product();
1126 
1127  // get the rechits
1128  edm::Handle<HBHERecHitCollection> handle;
1129  // iEvent.getByLabel(recHitCollName_, handle);
1130  iEvent.getByToken(hbherechit_token_, handle);
1131 
1132  if (!handle.isValid()) {
1133  throw edm::Exception(edm::errors::ProductNotFound)
1134  << " could not find HBHERecHitCollection named " << recHitCollName_ << "\n.";
1135  return;
1136  }
1137 
1138  summary.rechitCount_ = 0;
1139  summary.rechitCount15_ = 0;
1140  summary.rechitEnergy_ = 0;
1141  summary.rechitEnergy15_ = 0;
1142 
1143  summary.hitsInLaserRegion_ = 0;
1144  summary.hitsInNonLaserRegion_ = 0;
1145  summary.energyInLaserRegion_ = 0;
1146  summary.energyInNonLaserRegion_ = 0;
1147 
1148  // loop over all of the rechit information
1149  for (HBHERecHitCollection::const_iterator it = handle->begin(); it != handle->end(); ++it) {
1150  const HBHERecHit& rechit = (*it);
1151 
1152  // skip bad rechits (other than those flagged by the isolated noise, triangle, flat, and spike algorithms)
1153  const DetId id = rechit.idFront();
1154 
1155  uint32_t recHitFlag = rechit.flags();
1156  uint32_t isolbitset = (1 << HcalCaloFlagLabels::HBHEIsolatedNoise);
1157  uint32_t flatbitset = (1 << HcalCaloFlagLabels::HBHEFlatNoise);
1158  uint32_t spikebitset = (1 << HcalCaloFlagLabels::HBHESpikeNoise);
1159  uint32_t trianglebitset = (1 << HcalCaloFlagLabels::HBHETriangleNoise);
1160  uint32_t ts4ts5bitset = (1 << HcalCaloFlagLabels::HBHETS4TS5Noise);
1161  uint32_t negativebitset = (1 << HcalCaloFlagLabels::HBHENegativeNoise);
1162  for (unsigned int i = 0; i < HcalRecHitFlagsToBeExcluded_.size(); i++) {
1163  uint32_t bitset = (1 << HcalRecHitFlagsToBeExcluded_[i]);
1164  recHitFlag = (recHitFlag & bitset) ? recHitFlag - bitset : recHitFlag;
1165  }
1166  const HcalChannelStatus* dbStatus = dbHcalChStatus->getValues(id);
1167 
1168  // Ignore rechit if exclude bit set, regardless of severity of other bits
1169  if (dbStatus->isBitSet(HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary))
1170  continue;
1171 
1172  int severityLevel = hcalSevLvlComputer->getSeverityLevel(id, recHitFlag, dbStatus->getValue());
1173  bool isRecovered = hcalSevLvlComputer->recoveredRecHit(id, recHitFlag);
1174  if (severityLevel != 0 && !isRecovered && severityLevel > static_cast<int>(HcalAcceptSeverityLevel_))
1175  continue;
1176 
1177  // do some rechit counting and energies
1178  summary.rechitCount_ = summary.rechitCount_ + 1;
1179  summary.rechitEnergy_ = summary.rechitEnergy_ + rechit.eraw();
1180  if (dbStatus->isBitSet(
1181  HcalChannelStatus::
1182  HcalBadLaserSignal)) // hit comes from a region where no laser calibration pulse is normally seen
1183  {
1184  ++summary.hitsInNonLaserRegion_;
1185  summary.energyInNonLaserRegion_ += rechit.eraw();
1186  } else // hit comes from region where laser calibration pulse is seen
1187  {
1188  ++summary.hitsInLaserRegion_;
1189  summary.energyInLaserRegion_ += rechit.eraw();
1190  }
1191 
1192  if (rechit.eraw() > 1.5) {
1193  summary.rechitCount15_ = summary.rechitCount15_ + 1;
1194  summary.rechitEnergy15_ = summary.rechitEnergy15_ + rechit.eraw();
1195  }
1196 
1197  // Exclude uncollapsed QIE11 channels
1198  if (CaloRecHitAuxSetter::getBit(rechit.auxPhase1(), HBHERecHitAuxSetter::OFF_TDC_TIME) &&
1199  !CaloRecHitAuxSetter::getBit(rechit.auxPhase1(), HBHERecHitAuxSetter::OFF_COMBINED))
1200  continue;
1201 
1202  // if it was ID'd as isolated noise, update the summary object
1203  if (rechit.flags() & isolbitset) {
1204  summary.nisolnoise_++;
1205  summary.isolnoisee_ += rechit.eraw();
1206  GlobalPoint gp = geo->getPosition(rechit.id());
1207  double et = rechit.eraw() * gp.perp() / gp.mag();
1208  summary.isolnoiseet_ += et;
1209  }
1210 
1211  if (rechit.flags() & flatbitset) {
1212  summary.nflatnoise_++;
1213  summary.flatnoisee_ += rechit.eraw();
1214  GlobalPoint gp = geo->getPosition(rechit.id());
1215  double et = rechit.eraw() * gp.perp() / gp.mag();
1216  summary.flatnoiseet_ += et;
1217  }
1218 
1219  if (rechit.flags() & spikebitset) {
1220  summary.nspikenoise_++;
1221  summary.spikenoisee_ += rechit.eraw();
1222  GlobalPoint gp = geo->getPosition(rechit.id());
1223  double et = rechit.eraw() * gp.perp() / gp.mag();
1224  summary.spikenoiseet_ += et;
1225  }
1226 
1227  if (rechit.flags() & trianglebitset) {
1228  summary.ntrianglenoise_++;
1229  summary.trianglenoisee_ += rechit.eraw();
1230  GlobalPoint gp = geo->getPosition(rechit.id());
1231  double et = rechit.eraw() * gp.perp() / gp.mag();
1232  summary.trianglenoiseet_ += et;
1233  }
1234 
1235  if (rechit.flags() & ts4ts5bitset) {
1236  // only add to TS4TS5 if the bit is not marked as "HcalCellExcludeFromHBHENoiseSummaryR45"
1237  if (not dbStatus->isBitSet(HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummaryR45)) {
1238  summary.nts4ts5noise_++;
1239  summary.ts4ts5noisee_ += rechit.eraw();
1240  GlobalPoint gp = geo->getPosition(rechit.id());
1241  double et = rechit.eraw() * gp.perp() / gp.mag();
1242  summary.ts4ts5noiseet_ += et;
1243  }
1244  }
1245 
1246  if (rechit.flags() & negativebitset) {
1247  summary.nnegativenoise_++;
1248  summary.negativenoisee_ += rechit.eraw();
1249  GlobalPoint gp = geo->getPosition(rechit.id());
1250  double et = rechit.eraw() * gp.perp() / gp.mag();
1251  summary.negativenoiseet_ += et;
1252  }
1253 
1254  // find the hpd that the rechit is in
1255  HcalNoiseHPD& hpd = (*array.findHPD(rechit));
1256 
1257  // create a persistent reference to the rechit
1258  edm::Ref<HBHERecHitCollection> myRef(handle, it - handle->begin());
1259 
1260  // store it in a place so that it remains sorted by energy
1261  hpd.refrechitset_.insert(myRef);
1262 
1263  } // end loop over rechits
1264 
1265  // loop over all HPDs and transfer the information from refrechitset_ to rechits_;
1266  for (HcalNoiseRBXArray::iterator rbxit = array.begin(); rbxit != array.end(); ++rbxit) {
1267  for (std::vector<HcalNoiseHPD>::iterator hpdit = rbxit->hpds_.begin(); hpdit != rbxit->hpds_.end(); ++hpdit) {
1268  // loop over all of the entries in the set and add them to rechits_
1269  for (std::set<edm::Ref<HBHERecHitCollection>, RefHBHERecHitEnergyComparison>::const_iterator it =
1270  hpdit->refrechitset_.begin();
1271  it != hpdit->refrechitset_.end();
1272  ++it) {
1273  hpdit->rechits_.push_back(*it);
1274  }
1275  }
1276  }
1277  // now the rechits in all the HPDs are sorted by energy!
1278 
1279  return;
1280 }
reco::HcalNoiseInfoProducer::HcalRecHitFlagsToBeExcluded_
std::vector< int > HcalRecHitFlagsToBeExcluded_
Definition: HcalNoiseInfoProducer.h:147
mps_fire.i
i
Definition: mps_fire.py:428
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
HBHERecHit::auxPhase1
constexpr uint32_t auxPhase1() const
Definition: HBHERecHit.h:56
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
reco::HcalNoiseInfoProducer::severitycomputer_token_
edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd > severitycomputer_token_
Definition: HcalNoiseInfoProducer.h:124
reco::HcalNoiseHPD::refrechitset_
std::set< edm::Ref< HBHERecHitCollection >, RefHBHERecHitEnergyComparison > refrechitset_
Definition: HcalNoiseHPD.h:145
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
HBHERecHit::idFront
HcalDetId idFront() const
Definition: HBHERecHit.cc:5
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
HBHERecHitAuxSetter::OFF_TDC_TIME
static const unsigned OFF_TDC_TIME
Definition: HBHERecHitAuxSetter.h:35
reco::HcalNoiseInfoProducer::quality_token_
edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd > quality_token_
Definition: HcalNoiseInfoProducer.h:123
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:64
HBHERecHit::id
constexpr HcalDetId id() const
get the id
Definition: HBHERecHit.h:41
HcalSeverityLevelComputer::recoveredRecHit
bool recoveredRecHit(const DetId &myid, const uint32_t &myflag) const
Definition: HcalSeverityLevelComputer.cc:363
CaloGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
patZpeak.handle
tuple handle
Definition: patZpeak.py:25
reco::HcalNoiseRBXArray::findHPD
std::vector< HcalNoiseHPD >::iterator findHPD(int hpdindex)
Definition: HcalNoiseRBXArray.cc:46
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
reco::HcalNoiseInfoProducer::calogeometry_token_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > calogeometry_token_
Definition: HcalNoiseInfoProducer.h:125
CaloRecHitAuxSetter::getBit
constexpr bool getBit(const uint32_t u, const unsigned bitnum)
Definition: CaloRecHitAuxSetter.h:37
DetId
Definition: DetId.h:17
reco::HcalNoiseInfoProducer::hbherechit_token_
edm::EDGetTokenT< HBHERecHitCollection > hbherechit_token_
Definition: HcalNoiseInfoProducer.h:117
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
HcalSeverityLevelComputer::getSeverityLevel
int getSeverityLevel(const DetId &myid, const uint32_t &myflag, const uint32_t &mystatus) const
Definition: HcalSeverityLevelComputer.cc:304
HcalChannelStatus::isBitSet
bool isBitSet(unsigned int bitnumber) const
Definition: HcalChannelStatus.h:53
HBHERecHit::eraw
constexpr float eraw() const
Definition: HBHERecHit.h:47
HcalNoiseSummary::energyInNonLaserRegion_
double energyInNonLaserRegion_
Definition: HcalNoiseSummary.h:214
HBHERecHitAuxSetter::OFF_COMBINED
static const unsigned OFF_COMBINED
Definition: HBHERecHitAuxSetter.h:41
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:151
HcalChannelStatus::getValue
uint32_t getValue() const
Definition: HcalChannelStatus.h:60
CaloRecHit::flags
constexpr uint32_t flags() const
Definition: CaloRecHit.h:34
void HcalNoiseInfoProducer::filltracks ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup,
HcalNoiseSummary &  summary 
) const
private

Definition at line 1362 of file HcalNoiseInfoProducer.cc.

References reco::TrackBase::eta(), edm::Event::getByToken(), patZpeak::handle, edm::HandleBase::isValid(), maxTrackEta_, minTrackPt_, reco::TrackBase::p(), reco::TrackBase::pt(), track_token_, and HcalNoiseSummary::trackenergy_.

Referenced by produce().

1364  {
1365  edm::Handle<reco::TrackCollection> handle;
1366  // iEvent.getByLabel(trackCollName_, handle);
1367  iEvent.getByToken(track_token_, handle);
1368 
1369  // don't throw exception, just return quietly
1370  if (!handle.isValid()) {
1371  // throw edm::Exception(edm::errors::ProductNotFound)
1372  // << " could not find trackCollection named " << trackCollName_ << "\n.";
1373  return;
1374  }
1375 
1376  summary.trackenergy_ = 0.0;
1377  for (reco::TrackCollection::const_iterator iTrack = handle->begin(); iTrack != handle->end(); ++iTrack) {
1378  reco::Track trk = *iTrack;
1379  if (trk.pt() < minTrackPt_ || fabs(trk.eta()) > maxTrackEta_)
1380  continue;
1381 
1382  summary.trackenergy_ += trk.p();
1383  }
1384 
1385  return;
1386 }
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:631
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
patZpeak.handle
tuple handle
Definition: patZpeak.py:25
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
reco::HcalNoiseInfoProducer::track_token_
edm::EDGetTokenT< reco::TrackCollection > track_token_
Definition: HcalNoiseInfoProducer.h:119
void HcalNoiseInfoProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivate

Definition at line 674 of file HcalNoiseInfoProducer.cc.

References algo_, data, CommonHcalNoiseRBXData::energy(), fillCaloTowers_, fillcalotwrs(), filldigis(), fillDigis_, filljetinfo(), fillOtherSummaryVariables(), fillrechits(), fillRecHits_, filltracks(), fillTracks_, edm::EventSetup::getHandle(), hcaltopo_token_, HcalNoiseAlgo::isProblematic(), maxProblemRBXs_, minHighHitE_, minLowHitE_, minR45HitE_, minRecHitE_, eostools::move(), HcalNoiseSummary::nproblemRBXs_, HcalNoiseAlgo::passHighLevelNoiseFilter(), HcalNoiseAlgo::passLooseNoiseFilter(), HcalNoiseAlgo::passTightNoiseFilter(), edm::ESHandle< class >::product(), edm::Event::put(), edmLumisInFiles::summary, theHcalTopology_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

674  {
675  // this is what we're going to actually write to the EDM
676  auto result1 = std::make_unique<HcalNoiseRBXCollection>();
677  auto result2 = std::make_unique<HcalNoiseSummary>();
678 
679  // define an empty HcalNoiseRBXArray that we're going to fill
680  HcalNoiseRBXArray rbxarray;
681  HcalNoiseSummary& summary = *result2;
682 
683  // Get topology class to use later
684  edm::ESHandle<HcalTopology> topo = iSetup.getHandle(hcaltopo_token_);
685  theHcalTopology_ = topo.product();
686 
687  // fill them with the various components
688  // digi assumes that rechit information is available
689  if (fillRecHits_)
690  fillrechits(iEvent, iSetup, rbxarray, summary);
691  if (fillDigis_)
692  filldigis(iEvent, iSetup, rbxarray, summary);
693  if (fillCaloTowers_)
694  fillcalotwrs(iEvent, iSetup, rbxarray, summary);
695  if (fillTracks_)
696  filltracks(iEvent, iSetup, summary);
697 
698  filljetinfo(iEvent, iSetup, summary);
699 
700  // select those RBXs which are interesting
701  // also look for the highest energy RBX
702  HcalNoiseRBXArray::iterator maxit = rbxarray.begin();
703  double maxenergy = -999;
704  bool maxwritten = false;
705  for (HcalNoiseRBXArray::iterator rit = rbxarray.begin(); rit != rbxarray.end(); ++rit) {
706  HcalNoiseRBX& rbx = (*rit);
707  CommonHcalNoiseRBXData data(rbx,
708  minRecHitE_,
709  minLowHitE_,
710  minHighHitE_,
711  TS4TS5EnergyThreshold_,
712  TS4TS5UpperCut_,
713  TS4TS5LowerCut_,
714  minR45HitE_);
715 
716  // find the highest energy rbx
717  if (data.energy() > maxenergy) {
718  maxenergy = data.energy();
719  maxit = rit;
720  maxwritten = false;
721  }
722 
723  // find out if the rbx is problematic/noisy/etc.
724  bool writerbx = algo_.isProblematic(data) || !algo_.passLooseNoiseFilter(data) ||
725  !algo_.passTightNoiseFilter(data) || !algo_.passHighLevelNoiseFilter(data);
726 
727  // fill variables in the summary object not filled elsewhere
728  fillOtherSummaryVariables(summary, data);
729 
730  if (writerbx) {
731  summary.nproblemRBXs_++;
732  if (summary.nproblemRBXs_ <= maxProblemRBXs_) {
733  result1->push_back(rbx);
734  if (maxit == rit)
735  maxwritten = true;
736  }
737  }
738  } // end loop over rbxs
739 
740  // if we still haven't written the maximum energy rbx, write it now
741  if (!maxwritten) {
742  HcalNoiseRBX& rbx = (*maxit);
743 
744  // add the RBX to the event
745  result1->push_back(rbx);
746  }
747 
748  // put the rbxcollection and summary into the EDM
749  iEvent.put(std::move(result1));
750  iEvent.put(std::move(result2));
751 
752  return;
753 }
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
HcalNoiseAlgo::isProblematic
bool isProblematic(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:187
reco::HcalNoiseInfoProducer::fillrechits
void fillrechits(edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const
Definition: HcalNoiseInfoProducer.cc:1111
reco::HcalNoiseInfoProducer::TS4TS5LowerCut_
std::vector< std::pair< double, double > > TS4TS5LowerCut_
Definition: HcalNoiseInfoProducer.h:144
reco::HcalNoiseInfoProducer::fillcalotwrs
void fillcalotwrs(edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const
Definition: HcalNoiseInfoProducer.cc:1283
reco::HcalNoiseInfoProducer::hcaltopo_token_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > hcaltopo_token_
Definition: HcalNoiseInfoProducer.h:121
HcalNoiseAlgo::passHighLevelNoiseFilter
bool passHighLevelNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:234
reco::HcalNoiseInfoProducer::filljetinfo
void filljetinfo(edm::Event &, const edm::EventSetup &, HcalNoiseSummary &) const
Definition: HcalNoiseInfoProducer.cc:1327
eostools.move
def move
Definition: eostools.py:511
reco::HcalNoiseInfoProducer::TS4TS5UpperCut_
std::vector< std::pair< double, double > > TS4TS5UpperCut_
Definition: HcalNoiseInfoProducer.h:143
HcalNoiseAlgo::passLooseNoiseFilter
bool passLooseNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:224
reco::HcalNoiseInfoProducer::filldigis
void filldigis(edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &)
Definition: HcalNoiseInfoProducer.cc:860
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
HcalNoiseAlgo::passTightNoiseFilter
bool passTightNoiseFilter(const CommonHcalNoiseRBXData &) const
Definition: HcalNoiseAlgo.cc:229
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:79
reco::HcalNoiseInfoProducer::theHcalTopology_
const HcalTopology * theHcalTopology_
Definition: HcalNoiseInfoProducer.h:106
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:151
reco::HcalNoiseInfoProducer::filltracks
void filltracks(edm::Event &, const edm::EventSetup &, HcalNoiseSummary &) const
Definition: HcalNoiseInfoProducer.cc:1362
reco::HcalNoiseInfoProducer::fillOtherSummaryVariables
void fillOtherSummaryVariables(HcalNoiseSummary &summary, const CommonHcalNoiseRBXData &data) const
Definition: HcalNoiseInfoProducer.cc:756

Member Data Documentation

std::vector<float> reco::HcalNoiseInfoProducer::adc2fC
private

Definition at line 157 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<float> reco::HcalNoiseInfoProducer::adc2fCHF
private

Definition at line 158 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

HcalNoiseAlgo reco::HcalNoiseInfoProducer::algo_
private

Definition at line 131 of file HcalNoiseInfoProducer.h.

Referenced by fillOtherSummaryVariables(), and produce().

double reco::HcalNoiseInfoProducer::calibdigiHBHEthreshold_
private

Definition at line 136 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::calibdigiHBHEtimeslices_
private

Definition at line 137 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

double reco::HcalNoiseInfoProducer::calibdigiHFthreshold_
private

Definition at line 139 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::calibdigiHFtimeslices_
private

Definition at line 140 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> reco::HcalNoiseInfoProducer::calogeometry_token_
private

Definition at line 125 of file HcalNoiseInfoProducer.h.

Referenced by fillrechits(), and HcalNoiseInfoProducer().

edm::EDGetTokenT<CaloTowerCollection> reco::HcalNoiseInfoProducer::calotower_token_
private

Definition at line 118 of file HcalNoiseInfoProducer.h.

Referenced by fillcalotwrs(), and HcalNoiseInfoProducer().

std::string reco::HcalNoiseInfoProducer::caloTowerCollName_
private

Definition at line 110 of file HcalNoiseInfoProducer.h.

Referenced by fillcalotwrs(), and HcalNoiseInfoProducer().

std::string reco::HcalNoiseInfoProducer::digiCollName_
private

Definition at line 108 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

bool reco::HcalNoiseInfoProducer::fillCaloTowers_
private

Definition at line 92 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

bool reco::HcalNoiseInfoProducer::fillDigis_
private

Definition at line 90 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

bool reco::HcalNoiseInfoProducer::fillLaserMonitor_
private

Definition at line 94 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

bool reco::HcalNoiseInfoProducer::fillRecHits_
private

Definition at line 91 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

bool reco::HcalNoiseInfoProducer::fillTracks_
private

Definition at line 93 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

edm::EDGetTokenT<HBHEDigiCollection> reco::HcalNoiseInfoProducer::hbhedigi_token_
private

Definition at line 114 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

edm::EDGetTokenT<HBHERecHitCollection> reco::HcalNoiseInfoProducer::hbherechit_token_
private

Definition at line 117 of file HcalNoiseInfoProducer.h.

Referenced by fillrechits(), and HcalNoiseInfoProducer().

uint32_t reco::HcalNoiseInfoProducer::HcalAcceptSeverityLevel_
private

Definition at line 146 of file HcalNoiseInfoProducer.h.

Referenced by fillrechits(), and HcalNoiseInfoProducer().

edm::EDGetTokenT<HcalCalibDigiCollection> reco::HcalNoiseInfoProducer::hcalcalibdigi_token_
private

Definition at line 115 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::HcalRecHitFlagsToBeExcluded_
private

Definition at line 147 of file HcalNoiseInfoProducer.h.

Referenced by fillrechits(), and HcalNoiseInfoProducer().

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> reco::HcalNoiseInfoProducer::hcaltopo_token_
private

Definition at line 121 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

edm::EDGetTokenT<reco::PFJetCollection> reco::HcalNoiseInfoProducer::jet_token_
private

Definition at line 120 of file HcalNoiseInfoProducer.h.

Referenced by filljetinfo(), and HcalNoiseInfoProducer().

std::string reco::HcalNoiseInfoProducer::jetCollName_
private

Definition at line 112 of file HcalNoiseInfoProducer.h.

Referenced by filljetinfo(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::laserMonCBoxList_
private

Definition at line 149 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

edm::EDGetTokenT<QIE10DigiCollection> reco::HcalNoiseInfoProducer::lasermondigi_token_
private

Definition at line 116 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::laserMonIEtaList_
private

Definition at line 151 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

std::vector<int> reco::HcalNoiseInfoProducer::laserMonIPhiList_
private

Definition at line 150 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

unsigned reco::HcalNoiseInfoProducer::laserMonitorSamples_
private

Definition at line 155 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

int reco::HcalNoiseInfoProducer::laserMonitorTSEnd_
private

Definition at line 154 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

int reco::HcalNoiseInfoProducer::laserMonitorTSStart_
private

Definition at line 153 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

int reco::HcalNoiseInfoProducer::maxCaloTowerIEta_
private

Definition at line 100 of file HcalNoiseInfoProducer.h.

Referenced by fillcalotwrs(), and HcalNoiseInfoProducer().

int reco::HcalNoiseInfoProducer::maxjetindex_
private

Definition at line 104 of file HcalNoiseInfoProducer.h.

Referenced by filljetinfo(), and HcalNoiseInfoProducer().

double reco::HcalNoiseInfoProducer::maxNHF_
private

Definition at line 103 of file HcalNoiseInfoProducer.h.

Referenced by filljetinfo(), and HcalNoiseInfoProducer().

int reco::HcalNoiseInfoProducer::maxProblemRBXs_
private

Definition at line 97 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

double reco::HcalNoiseInfoProducer::maxTrackEta_
private

Definition at line 101 of file HcalNoiseInfoProducer.h.

Referenced by filltracks(), and HcalNoiseInfoProducer().

double reco::HcalNoiseInfoProducer::minHighHitE_
private

Definition at line 130 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

double reco::HcalNoiseInfoProducer::minLowHitE_
private

Definition at line 130 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

double reco::HcalNoiseInfoProducer::minR45HitE_
private

Definition at line 130 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

double reco::HcalNoiseInfoProducer::minRecHitE_
private

Definition at line 130 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

double reco::HcalNoiseInfoProducer::minTrackPt_
private

Definition at line 102 of file HcalNoiseInfoProducer.h.

Referenced by filltracks(), and HcalNoiseInfoProducer().

edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> reco::HcalNoiseInfoProducer::quality_token_
private

Definition at line 123 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), fillrechits(), and HcalNoiseInfoProducer().

std::string reco::HcalNoiseInfoProducer::recHitCollName_
private

Definition at line 109 of file HcalNoiseInfoProducer.h.

Referenced by fillrechits(), and HcalNoiseInfoProducer().

edm::ESGetToken<HcalDbService, HcalDbRecord> reco::HcalNoiseInfoProducer::service_token_
private

Definition at line 122 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), and HcalNoiseInfoProducer().

edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> reco::HcalNoiseInfoProducer::severitycomputer_token_
private

Definition at line 124 of file HcalNoiseInfoProducer.h.

Referenced by filldigis(), fillrechits(), and HcalNoiseInfoProducer().

const HcalTopology* reco::HcalNoiseInfoProducer::theHcalTopology_
private

Definition at line 106 of file HcalNoiseInfoProducer.h.

Referenced by produce().

double reco::HcalNoiseInfoProducer::totalCalibCharge
private

Definition at line 127 of file HcalNoiseInfoProducer.h.

Referenced by filldigis().

double reco::HcalNoiseInfoProducer::totalLasmonCharge
private

Definition at line 128 of file HcalNoiseInfoProducer.h.

Referenced by filldigis().

edm::EDGetTokenT<reco::TrackCollection> reco::HcalNoiseInfoProducer::track_token_
private

Definition at line 119 of file HcalNoiseInfoProducer.h.

Referenced by filltracks(), and HcalNoiseInfoProducer().

std::string reco::HcalNoiseInfoProducer::trackCollName_
private

Definition at line 111 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer().

double reco::HcalNoiseInfoProducer::TS4TS5EnergyThreshold_
private

Definition at line 142 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

std::vector<std::pair<double, double> > reco::HcalNoiseInfoProducer::TS4TS5LowerCut_
private

Definition at line 144 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

std::vector<std::pair<double, double> > reco::HcalNoiseInfoProducer::TS4TS5UpperCut_
private

Definition at line 143 of file HcalNoiseInfoProducer.h.

Referenced by HcalNoiseInfoProducer(), and produce().

bool reco::HcalNoiseInfoProducer::useCalibDigi_
private

Definition at line 133 of file HcalNoiseInfoProducer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5