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

#include <HcalHotCellMonitor.h>

Inheritance diagram for HcalHotCellMonitor:
HcalBaseDQMonitor DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (edm::Event const &e, edm::EventSetup const &s)
 
void beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
 
void bookHistograms (DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c)
 
void done ()
 
void endJob ()
 
void endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
 
void endRun (const edm::Run &run, const edm::EventSetup &c)
 
 HcalHotCellMonitor (const edm::ParameterSet &ps)
 
void periodicReset ()
 
void processEvent (const HBHERecHitCollection &hbHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, const HcalTopology &topology)
 
void processEvent_rechitenergy (const HBHERecHitCollection &hbheHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, const HcalTopology &topology)
 
template<class R , class C >
void processHit_rechitNeighbors (R &rechititer, C &collection, hotNeighborParams &params, const HcalTopology &topology)
 
template<class RECHIT , class RECHITCOLLECTION >
void processHit_rechitNeighbors (RECHIT &rechit, RECHITCOLLECTION &coll, hotNeighborParams &params, const HcalTopology &topology)
 
void reset ()
 
void setup (DQMStore::IBooker &)
 
 ~HcalHotCellMonitor ()
 
- Public Member Functions inherited from HcalBaseDQMonitor
 HcalBaseDQMonitor (const edm::ParameterSet &ps)
 
 HcalBaseDQMonitor ()
 
virtual ~HcalBaseDQMonitor ()
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

void fillNevents_energy (const HcalTopology &)
 
void fillNevents_neighbor (const HcalTopology &)
 
void fillNevents_persistentenergy (const HcalTopology &)
 
void fillNevents_problemCells (const HcalTopology &)
 
void zeroCounters ()
 

Private Attributes

int aboveenergy [85][72][4]
 
EtaPhiHists AboveEnergyThresholdCellsByDepth
 
int aboveet [85][72][4]
 
EtaPhiHists AboveETThresholdCellsByDepth
 
int aboveneighbors [85][72][4]
 
EtaPhiHists AboveNeighborsHotCellsByDepth
 
int abovepersistent [85][72][4]
 
int abovepersistentET [85][72][4]
 
EtaPhiHists AbovePersistentETThresholdCellsByDepth
 
EtaPhiHists AbovePersistentThresholdCellsByDepth
 
MonitorElement * d_HBenergyVsNeighbor
 
MonitorElement * d_HEenergyVsNeighbor
 
MonitorElement * d_HFenergyVsNeighbor
 
MonitorElement * d_HOenergyVsNeighbor
 
double energyThreshold_
 
double ETThreshold_
 
bool excludeHORing2_
 
double HBenergyThreshold_
 
double HBETThreshold_
 
hotNeighborParams HBHENeighborParams_
 
edm::InputTag hbheRechitLabel_
 
double HBnsigma_
 
double HBpersistentETThreshold_
 
double HBpersistentThreshold_
 
int hbVsNeighbor [500]
 
double HEenergyThreshold_
 
double HEETThreshold_
 
double HEnsigma_
 
double HEpersistentETThreshold_
 
double HEpersistentThreshold_
 
int heVsNeighbor [500]
 
double HFenergyThreshold_
 
double HFETThreshold_
 
double HFfarfwdScale_
 
hotNeighborParams HFNeighborParams_
 
double HFnsigma_
 
double HFpersistentETThreshold_
 
double HFpersistentThreshold_
 
edm::InputTag hfRechitLabel_
 
int hfVsNeighbor [500]
 
double HOenergyThreshold_
 
double HOETThreshold_
 
hotNeighborParams HONeighborParams_
 
double HOnsigma_
 
double HOpersistentETThreshold_
 
double HOpersistentThreshold_
 
edm::InputTag hoRechitLabel_
 
int hoVsNeighbor [500]
 
double minErrorFlag_
 
int minEvents_
 
double nsigma_
 
double persistentETThreshold_
 
double persistentThreshold_
 
int rechit_occupancy_sum [85][72][4]
 
bool setupDone_
 
double SiPMscale_
 
bool test_energy_
 
bool test_et_
 
bool test_neighbor_
 
bool test_persistent_
 
edm::EDGetTokenT
< HBHERecHitCollection > 		tok_hbhe_
 
edm::EDGetTokenT
< HFRecHitCollection > 		tok_hf_
 
edm::EDGetTokenT
< HORecHitCollection > 		tok_ho_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...> CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T...> HasAbility
 
typedef
CacheTypes::LuminosityBlockCache 		LuminosityBlockCache
 
typedef
LuminosityBlockContextT
< LuminosityBlockCache,
RunCache, GlobalCache > 		LuminosityBlockContext
 
typedef
CacheTypes::LuminosityBlockSummaryCache 		LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache,
GlobalCache > 		RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr
< dqmDetails::NoCache > 		globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from HcalBaseDQMonitor
virtual void cleanup (void)
 
virtual void dqmBeginRun (const edm::Run &run, const edm::EventSetup &c)
 
void getLogicalMap (const edm::EventSetup &c)
 
bool IsAllowedCalibType ()
 
bool LumiInOrder (int lumisec)
 
void SetupEtaPhiHists (DQMStore::IBooker &ib, EtaPhiHists &hh, std::string Name, std::string Units)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
- Protected Attributes inherited from HcalBaseDQMonitor
std::vector< int > AllowedCalibTypes_
 
int badChannelStatusMask_
 
int currentLS
 
int currenttype_
 
int debug_
 
bool enableCleanup_
 
bool eventAllowed_
 
bool HBpresent_
 
bool HEpresent_
 
bool HFpresent_
 
bool HOpresent_
 
int ievt_
 
std::map< unsigned int, int > KnownBadCells_
 
int levt_
 
HcalLogicalMap * logicalMap_
 
bool makeDiagnostics_
 
MonitorElement * meIevt_
 
MonitorElement * meLevt_
 
bool mergeRuns_
 
MonitorElement * meTevt_
 
MonitorElement * meTevtHist_
 
bool needLogicalMap_
 
int NLumiBlocks_
 
bool Online_
 
std::string prefixME_
 
MonitorElement * ProblemsCurrentLB
 
MonitorElement * ProblemsVsLB
 
MonitorElement * ProblemsVsLB_HB
 
MonitorElement * ProblemsVsLB_HBHEHF
 
MonitorElement * ProblemsVsLB_HE
 
MonitorElement * ProblemsVsLB_HF
 
MonitorElement * ProblemsVsLB_HO
 
bool skipOutOfOrderLS_
 
std::string subdir_
 
int tevt_
 

Detailed Description

Author
J. Temple - Univ. of Maryland

Definition at line 29 of file HcalHotCellMonitor.h.

Constructor & Destructor Documentation

HcalHotCellMonitor::HcalHotCellMonitor ( const edm::ParameterSet &  ps)

Definition at line 5 of file HcalHotCellMonitor.cc.

References HcalBaseDQMonitor::AllowedCalibTypes_, HcalBaseDQMonitor::debug_, hotNeighborParams::DeltaDepth, hotNeighborParams::DeltaIeta, hotNeighborParams::DeltaIphi, HcalBaseDQMonitor::enableCleanup_, energyThreshold_, ETThreshold_, excludeHORing2_, edm::ParameterSet::getUntrackedParameter(), HBenergyThreshold_, HBETThreshold_, HBHENeighborParams_, hbheRechitLabel_, HBpersistentETThreshold_, HBpersistentThreshold_, HEenergyThreshold_, HEETThreshold_, HEpersistentETThreshold_, HEpersistentThreshold_, HFenergyThreshold_, HFETThreshold_, HFfarfwdScale_, HFNeighborParams_, HFpersistentETThreshold_, HFpersistentThreshold_, hfRechitLabel_, HOenergyThreshold_, HOETThreshold_, HONeighborParams_, HOpersistentETThreshold_, HOpersistentThreshold_, hoRechitLabel_, hotNeighborParams::HotEnergyFrac, HcalBaseDQMonitor::makeDiagnostics_, hotNeighborParams::maxEnergy, HcalBaseDQMonitor::mergeRuns_, hotNeighborParams::minCellEnergy, minErrorFlag_, minEvents_, hotNeighborParams::minNeighborEnergy, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::Online_, persistentETThreshold_, persistentThreshold_, HcalBaseDQMonitor::prefixME_, setupDone_, SiPMscale_, HcalBaseDQMonitor::skipOutOfOrderLS_, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, test_energy_, test_et_, test_neighbor_, test_persistent_, tok_hbhe_, tok_hf_, and tok_ho_.

5  :HcalBaseDQMonitor(ps) {
6  // Standard information, inherited from base class
7  Online_ = ps.getUntrackedParameter<bool>("online",false);
8  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
9  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
10  debug_ = ps.getUntrackedParameter<int>("debug",0);
11  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
12  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
13  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
14  prefixME_.append("/");
15  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","HotCellMonitor_Hcal/"); // HotCellMonitor_Hcal
16  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
17  subdir_.append("/");
18  subdir_=prefixME_+subdir_;
19  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
20  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",true);
21  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
22 
23  // Collection type info
24  hbheRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hbheRechitLabel");
25  hoRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hoRechitLabel");
26  hfRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hfRechitLabel");
27 
28  // register for data access
29  tok_hbhe_ = consumes<HBHERecHitCollection>(hbheRechitLabel_);
30  tok_ho_ = consumes<HORecHitCollection>(hoRechitLabel_);
31  tok_hf_ = consumes<HFRecHitCollection>(hfRechitLabel_);
32 
33  // Hot Cell-specific tests
34  minEvents_ = ps.getUntrackedParameter<int>("minEvents");
35  minErrorFlag_ = ps.getUntrackedParameter<double>("minErrorFlag",1);
36 
37  // exclude HO ring 2
38  excludeHORing2_ = ps.getUntrackedParameter<bool>("excludeHORing2",false);
39 
40 
41  // Set which hot cell checks will be performed
42  test_persistent_ = ps.getUntrackedParameter<bool>("test_persistent"); // true by default
43  test_neighbor_ = ps.getUntrackedParameter<bool>("test_neighbor"); // false by default; test disabled
44  test_energy_ = ps.getUntrackedParameter<bool>("test_energy"); // true by default
45  test_et_ = ps.getUntrackedParameter<bool>("test_et"); // true by default
46 
47 
48  // rechit energy test -- cell must be above threshold value for a number of consecutive events to be considered hot
49  energyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold");
50  ETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold");
51 
52  HBenergyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold_HB",energyThreshold_);
53  HEenergyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold_HE",energyThreshold_);
54  HOenergyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold_HO",energyThreshold_);
55  HFenergyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold_HF",energyThreshold_);
56 
57  HBETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold_HB",ETThreshold_);
58  HEETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold_HE",ETThreshold_);
59  HOETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold_HO",ETThreshold_);
60  HFETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold_HF",ETThreshold_);
61 
62  // rechit event-by-event energy test -- cell must be above threshold to be considered hot
63  persistentThreshold_ = ps.getUntrackedParameter<double>("persistentThreshold");
64 
65  HBpersistentThreshold_ = ps.getUntrackedParameter<double>("persistentThreshold_HB",persistentThreshold_);
66  HEpersistentThreshold_ = ps.getUntrackedParameter<double>("persistentThreshold_HE",persistentThreshold_);
67  HOpersistentThreshold_ = ps.getUntrackedParameter<double>("persistentThreshold_HO",persistentThreshold_);
68  HFpersistentThreshold_ = ps.getUntrackedParameter<double>("persistentThreshold_HF",persistentThreshold_);
69 
70  persistentETThreshold_ = ps.getUntrackedParameter<double>("persistentETThreshold");
71 
72  HBpersistentETThreshold_ = ps.getUntrackedParameter<double>("persistentETThreshold_HB",persistentETThreshold_);
73  HEpersistentETThreshold_ = ps.getUntrackedParameter<double>("persistentETThreshold_HE",persistentETThreshold_);
74  HOpersistentETThreshold_ = ps.getUntrackedParameter<double>("persistentETThreshold_HO",persistentETThreshold_);
75  HFpersistentETThreshold_ = ps.getUntrackedParameter<double>("persistentETThreshold_HF",persistentETThreshold_);
76 
77  HFfarfwdScale_ = ps.getUntrackedParameter<double>("HFfwdScale",2.);
78  SiPMscale_ = ps.getUntrackedParameter<double>("HO_SiPMscalefactor",1.); // default scale factor of 4?
79 
80  // neighboring-cell tests
81  HBHENeighborParams_.DeltaIphi = ps.getUntrackedParameter<int>("HBHE_neighbor_deltaIphi", 1);
82  HBHENeighborParams_.DeltaIeta = ps.getUntrackedParameter<int>("HBHE_neighbor_deltaIeta", 1);
83  HBHENeighborParams_.DeltaDepth = ps.getUntrackedParameter<int>("HBHE_neighbor_deltaDepth", 2);
84  HBHENeighborParams_.minCellEnergy = ps.getUntrackedParameter<double>("HBHE_neighbor_minCellEnergy",3.);
85  HBHENeighborParams_.minNeighborEnergy = ps.getUntrackedParameter<double>("HBHE_neighbor_minNeighborEnergy",0.);
86  HBHENeighborParams_.maxEnergy = ps.getUntrackedParameter<double>("HBHE_neighbor_maxEnergy",100);
87  HBHENeighborParams_.HotEnergyFrac = ps.getUntrackedParameter<double>("HBHE_neighbor_HotEnergyFrac",0.05);
88 
89  HONeighborParams_.DeltaIphi = ps.getUntrackedParameter<int>("HO_neighbor_deltaIphi", 1);
90  HONeighborParams_.DeltaIeta = ps.getUntrackedParameter<int>("HO_neighbor_deltaIeta", 1);
91  HONeighborParams_.DeltaDepth = ps.getUntrackedParameter<int>("HO_neighbor_deltaDepth", 0);
92  HONeighborParams_.minCellEnergy = ps.getUntrackedParameter<double>("HO_neighbor_minCellEnergy",10.);
93  HONeighborParams_.minNeighborEnergy = ps.getUntrackedParameter<double>("HO_neighbor_minNeighborEnergy",0.);
94  HONeighborParams_.maxEnergy = ps.getUntrackedParameter<double>("HO_neighbor_maxEnergy",100);
95  HONeighborParams_.HotEnergyFrac = ps.getUntrackedParameter<double>("HO_neighbor_HotEnergyFrac",0.01);
96 
97  HFNeighborParams_.DeltaIphi = ps.getUntrackedParameter<int>("HF_neighbor_deltaIphi", 1);
98  HFNeighborParams_.DeltaIeta = ps.getUntrackedParameter<int>("HF_neighbor_deltaIeta", 1);
99  HFNeighborParams_.DeltaDepth = ps.getUntrackedParameter<int>("HF_neighbor_deltaDepth", 1);
100  HFNeighborParams_.minCellEnergy = ps.getUntrackedParameter<double>("HF_neighbor_minCellEnergy",10.);
101  HFNeighborParams_.minNeighborEnergy = ps.getUntrackedParameter<double>("HF_neighbor_minNeighborEnergy",0.);
102  HFNeighborParams_.maxEnergy = ps.getUntrackedParameter<double>("HF_neighbor_maxEnergy",100);
103  HFNeighborParams_.HotEnergyFrac = ps.getUntrackedParameter<double>("HF_neighbor_HotEnergyFrac",0.01);
104  setupDone_=false;
105 } //constructor
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalHotCellMonitor::hoRechitLabel_
edm::InputTag hoRechitLabel_
Definition: HcalHotCellMonitor.h:130
HcalHotCellMonitor::HONeighborParams_
hotNeighborParams HONeighborParams_
Definition: HcalHotCellMonitor.h:128
HcalBaseDQMonitor::AllowedCalibTypes_
std::vector< int > AllowedCalibTypes_
Definition: HcalBaseDQMonitor.h:84
HcalHotCellMonitor::tok_ho_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
Definition: HcalHotCellMonitor.h:132
HcalHotCellMonitor::tok_hf_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
Definition: HcalHotCellMonitor.h:133
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HcalHotCellMonitor::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalHotCellMonitor.h:131
HcalHotCellMonitor::HBHENeighborParams_
hotNeighborParams HBHENeighborParams_
Definition: HcalHotCellMonitor.h:128
HcalHotCellMonitor::hbheRechitLabel_
edm::InputTag hbheRechitLabel_
Definition: HcalHotCellMonitor.h:130
HcalHotCellMonitor::hfRechitLabel_
edm::InputTag hfRechitLabel_
Definition: HcalHotCellMonitor.h:130
HcalHotCellMonitor::HFNeighborParams_
hotNeighborParams HFNeighborParams_
Definition: HcalHotCellMonitor.h:128
HcalHotCellMonitor::~HcalHotCellMonitor ( )

Definition at line 107 of file HcalHotCellMonitor.cc.

107 { } //destructor

Member Function Documentation

void HcalHotCellMonitor::analyze ( edm::Event const &  e,
edm::EventSetup const &  s 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 374 of file HcalHotCellMonitor.cc.

References HcalBaseDQMonitor::analyze(), gather_cfg::cout, HcalBaseDQMonitor::debug_, edm::EventSetup::get(), edm::Event::getByToken(), hbheRechitLabel_, hfRechitLabel_, hoRechitLabel_, HcalBaseDQMonitor::ievt_, HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), processEvent(), tok_hbhe_, tok_hf_, and tok_ho_.

375 {
376  HcalBaseDQMonitor::analyze(e,s);
377  if (!IsAllowedCalibType()) return;
378  if (LumiInOrder(e.luminosityBlock())==false) return;
379 
380  // try to get rechits
381  edm::Handle<HBHERecHitCollection> hbhe_rechit;
382  edm::Handle<HORecHitCollection> ho_rechit;
383  edm::Handle<HFRecHitCollection> hf_rechit;
384 
385  if (!(e.getByToken(tok_hbhe_,hbhe_rechit)))
386  {
387  edm::LogWarning("HcalHotCellMonitor")<< hbheRechitLabel_<<" hbhe_rechit not available";
388  return;
389  }
390 
391  if (!(e.getByToken(tok_hf_,hf_rechit)))
392  {
393  edm::LogWarning("HcalHotCellMonitor")<< hfRechitLabel_<<" hf_rechit not available";
394  return;
395  }
396  if (!(e.getByToken(tok_ho_,ho_rechit)))
397  {
398  edm::LogWarning("HcalHotCellMonitor")<< hoRechitLabel_<<" ho_rechit not available";
399  return;
400  }
401 
402  // Good event found; increment counter (via base class analyze method)
403  edm::ESHandle<HcalTopology> topo;
404  s.get<HcalRecNumberingRecord>().get(topo);
405 
406  // HcalBaseDQMonitor::analyze(e,s);
407  if (debug_>1) std::cout <<"\t<HcalHotCellMonitor::analyze> Processing good event! event # = "<<ievt_<<std::endl;
408 
409  processEvent(*hbhe_rechit, *ho_rechit, *hf_rechit, *topo);
410 
411 } // void HcalHotCellMonitor::analyze(...)
HcalHotCellMonitor::processEvent
void processEvent(const HBHERecHitCollection &hbHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, const HcalTopology &topology)
Definition: HcalHotCellMonitor.cc:419
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:203
HcalHotCellMonitor::hoRechitLabel_
edm::InputTag hoRechitLabel_
Definition: HcalHotCellMonitor.h:130
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:244
HcalHotCellMonitor::tok_ho_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
Definition: HcalHotCellMonitor.h:132
HcalHotCellMonitor::tok_hf_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
Definition: HcalHotCellMonitor.h:133
HcalHotCellMonitor::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalHotCellMonitor.h:131
HcalHotCellMonitor::hbheRechitLabel_
edm::InputTag hbheRechitLabel_
Definition: HcalHotCellMonitor.h:130
HcalHotCellMonitor::hfRechitLabel_
edm::InputTag hfRechitLabel_
Definition: HcalHotCellMonitor.h:130
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::beginLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 328 of file HcalHotCellMonitor.cc.

References HcalBaseDQMonitor::beginLuminosityBlock(), HcalBaseDQMonitor::LumiInOrder(), edm::LuminosityBlockBase::luminosityBlock(), HcalBaseDQMonitor::ProblemsCurrentLB, MonitorElement::Reset(), and zeroCounters().

330 {
331  if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
332  HcalBaseDQMonitor::beginLuminosityBlock(lumiSeg,c);
333  zeroCounters(); // zero hot cell counters at the start of each luminosity block
334  ProblemsCurrentLB->Reset();
335  return;
336 } // beginLuminosityBlock(...)
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:203
HcalBaseDQMonitor::ProblemsCurrentLB
MonitorElement * ProblemsCurrentLB
Definition: HcalBaseDQMonitor.h:112
edm::LuminosityBlockBase::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockBase.h:40
HcalBaseDQMonitor::beginLuminosityBlock
virtual void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:184
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:453
void HcalHotCellMonitor::bookHistograms ( DQMStore::IBooker &  ib,
const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 282 of file HcalHotCellMonitor.cc.

References HcalBaseDQMonitor::bookHistograms(), gather_cfg::cout, HcalBaseDQMonitor::debug_, HcalBaseDQMonitor::mergeRuns_, reset(), setup(), and HcalBaseDQMonitor::tevt_.

283 {
284  if (debug_>1) std::cout <<"HcalHotCellMonitor::bookHistograms"<<std::endl;
285  HcalBaseDQMonitor::bookHistograms(ib,run,c);
286 
287  if (tevt_==0) this->setup(ib); // set up histograms if they have not been created before
288  if (mergeRuns_==false)
289  this->reset();
290 
291  return;
292 } //void HcalHotCellMonitor::bookHistograms(...)
HcalBaseDQMonitor::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: HcalBaseDQMonitor.cc:107
HcalHotCellMonitor::setup
void setup(DQMStore::IBooker &)
Definition: HcalHotCellMonitor.cc:112
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::done ( )

Definition at line 366 of file HcalHotCellMonitor.cc.

367 {
368  // moved database dumps to client; we want to be able to sum over results in offline
369  return;
370 
371 } // void HcalHotCellMonitor::done()
void HcalHotCellMonitor::endJob ( void  )

Definition at line 1088 of file HcalHotCellMonitor.cc.

References HcalBaseDQMonitor::cleanup(), gather_cfg::cout, HcalBaseDQMonitor::debug_, and HcalBaseDQMonitor::enableCleanup_.

1089 {
1090  if (debug_>0) std::cout <<"HcalHotCellMonitor::endJob()"<<std::endl;
1091  if (enableCleanup_) cleanup(); // when do we force cleanup?
1092 }
HcalBaseDQMonitor::cleanup
virtual void cleanup(void)
Definition: HcalBaseDQMonitor.cc:146
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::endLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 339 of file HcalHotCellMonitor.cc.

References fillNevents_energy(), fillNevents_neighbor(), fillNevents_persistentenergy(), fillNevents_problemCells(), edm::EventSetup::get(), HcalBaseDQMonitor::LumiInOrder(), edm::LuminosityBlockBase::luminosityBlock(), HcalBaseDQMonitor::makeDiagnostics_, test_energy_, test_et_, test_neighbor_, and test_persistent_.

341 {
342  //FIX with check on whether LB already processed
343 
344  if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
345 
346  edm::ESHandle<HcalTopology> topo;
347  c.get<HcalRecNumberingRecord>().get(topo);
348 
349  if (test_neighbor_ || makeDiagnostics_)
350  fillNevents_neighbor(*topo);
351 
352  if (test_energy_ || test_et_)
353  fillNevents_energy(*topo);
354 
355  if (test_persistent_)
356  fillNevents_persistentenergy(*topo);
357 
358  fillNevents_problemCells(*topo);
359  return;
360 } //endLuminosityBlock(...)
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:203
HcalHotCellMonitor::fillNevents_problemCells
void fillNevents_problemCells(const HcalTopology &)
Definition: HcalHotCellMonitor.cc:922
HcalHotCellMonitor::fillNevents_energy
void fillNevents_energy(const HcalTopology &)
Definition: HcalHotCellMonitor.cc:778
edm::LuminosityBlockBase::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockBase.h:40
HcalHotCellMonitor::fillNevents_neighbor
void fillNevents_neighbor(const HcalTopology &)
Definition: HcalHotCellMonitor.cc:860
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
HcalHotCellMonitor::fillNevents_persistentenergy
void fillNevents_persistentenergy(const HcalTopology &)
Definition: HcalHotCellMonitor.cc:682
void HcalHotCellMonitor::endRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 1083 of file HcalHotCellMonitor.cc.

1084 {
1085  // Anything to do here?
1086 }
void HcalHotCellMonitor::fillNevents_energy ( const HcalTopology &  topology)
private

Definition at line 778 of file HcalHotCellMonitor.cc.

References aboveenergy, AboveEnergyThresholdCellsByDepth, aboveet, AboveETThresholdCellsByDepth, HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, HLT_FULL_cff::depth, eta, HcalObjRepresent::FillUnphysicalHEHFBins(), h, HcalForward, HcalBaseDQMonitor::ievt_, phi, test_energy_, test_et_, and HcalTopology::validDetId().

Referenced by endLuminosityBlock().

778  {
779  // Fill Histograms showing rec hits that are above some energy value
780  // (Fill for each instance when cell is above energy; don't require it to be hot for a number of consecutive events)
781 
782  if (debug_>0)
783  std::cout <<"<HcalHotCellMonitor::fillNevents_energy> ABOVE-ENERGY-THRESHOLD PLOTS"<<std::endl;
784 
785  if (test_energy_) {
786  for (unsigned int h=0;h<AboveEnergyThresholdCellsByDepth.depth.size();++h)
787  AboveEnergyThresholdCellsByDepth.depth[h]->setBinContent(0,0,ievt_);
788  }
789  if (test_et_) {
790  for (unsigned int h=0;h<AboveETThresholdCellsByDepth.depth.size();++h)
791  AboveETThresholdCellsByDepth.depth[h]->setBinContent(0,0,ievt_);
792  }
793 
794  int ieta=0;
795  int iphi=0;
796  int etabins=0;
797  int phibins=0;
798  unsigned int maxdepth=0;
799 
800  if (test_energy_)
801  maxdepth = AboveEnergyThresholdCellsByDepth.depth.size();
802  if (maxdepth==0 && test_et_)
803  maxdepth = AboveETThresholdCellsByDepth.depth.size();
804  for (unsigned int depth=0;depth<maxdepth;++depth) {
805  if (test_energy_) {
806  etabins=AboveEnergyThresholdCellsByDepth.depth[depth]->getNbinsX();
807  phibins=AboveEnergyThresholdCellsByDepth.depth[depth]->getNbinsY();
808  }
809  if (test_et_) {
810  etabins=AboveETThresholdCellsByDepth.depth[depth]->getNbinsX();
811  phibins=AboveETThresholdCellsByDepth.depth[depth]->getNbinsY();
812  }
813  for (int eta=0;eta<etabins;++eta) {
814  for (int phi=0;phi<phibins;++phi) {
815  iphi=phi+1;
816  for (int subdet=1;subdet<=4;++subdet) {
817  ieta=CalcIeta((HcalSubdetector)subdet,eta,depth+1); //converts bin to ieta
818  if (ieta==-9999) continue;
819  if (!(topology.validDetId((HcalSubdetector)subdet, ieta, iphi, depth+1)))
820  continue;
821  if (subdet==HcalForward) // shift HcalForward ieta by 1 for filling purposes
822  ieta<0 ? ieta-- : ieta++;
823 
824  if (test_energy_) {
825  if (aboveenergy[eta][phi][depth]>0) {
826  if (debug_>2)
827  std::cout <<"HOT CELL; ABOVE ENERGY THRESHOLD at subdet = "<<subdet<<", eta = "<<ieta<<", phi = "<<iphi<<" depth = "<<depth+1<<" ABOVE THRESHOLD IN "<<aboveenergy[eta][phi][depth]<<" EVENTS"<<std::endl;
828  AboveEnergyThresholdCellsByDepth.depth[depth]->Fill(ieta,iphi, aboveenergy[eta][phi][depth]);
829  aboveenergy[eta][phi][depth]=0;
830  } // if (aboveenergy[eta][phi][depth])
831  } // if (test_energy_)
832  if (test_et_) {
833  if (aboveet[eta][phi][depth]>0) {
834  if (debug_>2)
835  std::cout <<"HOT CELL; ABOVE ET THRESHOLD at subdet = "<<subdet<<", eta = "<<ieta<<", phi = "<<iphi<<" depth = "<<depth+1<<" ABOVE THRESHOLD IN "<<aboveet[eta][phi][depth]<<" EVENTS"<<std::endl;
836  AboveETThresholdCellsByDepth.depth[depth]->Fill(ieta,iphi, aboveet[eta][phi][depth]);
837  aboveet[eta][phi][depth]=0;
838  } // if (aboveet[eta][phi][depth])
839  } // if (test_et_)
840  } // for (int subdet=0)
841  } // for (int phi=0;...)
842  } // for (int eta=0;...)
843  } // for (int depth=0;...)
844 
845  if (test_energy_)
846  FillUnphysicalHEHFBins(AboveEnergyThresholdCellsByDepth);
847 
848  if (test_et_)
849  FillUnphysicalHEHFBins(AboveETThresholdCellsByDepth);
850 
851  return;
852 
853 
854 } // void HcalHotCellMonitor::fillNevents_energy(const HcalTopology&)
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
HcalTopology::validDetId
bool validDetId(HcalSubdetector subdet, int ieta, int iphi, int depth) const
Definition: HcalTopology.cc:195
HcalHotCellMonitor::aboveenergy
int aboveenergy[85][72][4]
Definition: HcalHotCellMonitor.h:108
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
HcalHotCellMonitor::aboveet
int aboveet[85][72][4]
Definition: HcalHotCellMonitor.h:109
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::fillNevents_neighbor ( const HcalTopology &  topology)
private

Definition at line 860 of file HcalHotCellMonitor.cc.

References aboveneighbors, AboveNeighborsHotCellsByDepth, HcalObjRepresent::CalcIeta(), gather_cfg::cout, d_HBenergyVsNeighbor, d_HEenergyVsNeighbor, d_HFenergyVsNeighbor, d_HOenergyVsNeighbor, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, HLT_FULL_cff::depth, eta, MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), h, hbVsNeighbor, HcalForward, heVsNeighbor, hfVsNeighbor, hoVsNeighbor, i, HcalBaseDQMonitor::ievt_, HcalBaseDQMonitor::makeDiagnostics_, phi, and HcalTopology::validDetId().

Referenced by endLuminosityBlock().

860  {
861  // Fill Histograms showing rec hits with energy much less than neighbors' average
862 
863  if (debug_>0)
864  std::cout <<"<HcalHotCellMonitor::fillNevents_neighbor> FILLING ABOVE-NEIGHBOR-ENERGY PLOTS"<<std::endl;
865 
866  for (unsigned int h=0;h<AboveNeighborsHotCellsByDepth.depth.size();++h)
867  AboveNeighborsHotCellsByDepth.depth[h]->setBinContent(0,0,ievt_);
868 
869  int ieta=0;
870  int iphi=0;
871  int etabins=0;
872  int phibins=0;
873 
874  for (unsigned int depth=0;depth<AboveNeighborsHotCellsByDepth.depth.size();++depth) {
875  etabins=AboveNeighborsHotCellsByDepth.depth[depth]->getNbinsX();
876  phibins=AboveNeighborsHotCellsByDepth.depth[depth]->getNbinsY();
877 
878  for (int eta=0;eta<etabins;++eta) {
879  for (int phi=0;phi<phibins;++phi) {
880  iphi=phi+1;
881  for (int subdet=1;subdet<=4;++subdet) {
882  ieta=CalcIeta((HcalSubdetector)subdet,eta,depth+1); //converts bin to ieta
883  if (ieta==-9999) continue;
884  if (!(topology.validDetId((HcalSubdetector)subdet, ieta, iphi, depth+1)))
885  continue;
886  if (subdet==HcalForward) // shift HcalForward ieta by 1 for filling purposes
887  ieta<0 ? ieta-- : ieta++;
888 
889  if (aboveneighbors[eta][phi][depth]>0) {
890  if (debug_>2) std::cout <<"HOT CELL; ABOVE NEIGHBORS at eta = "<<ieta<<", phi = "<<iphi<<" depth = "<<(depth>4 ? depth+1 : depth-3)<<std::endl;
891  AboveNeighborsHotCellsByDepth.depth[depth]->Fill(ieta,iphi,aboveneighbors[eta][phi][depth]);
892  //reset counter
893  aboveneighbors[eta][phi][depth]=0;
894  } // if (aboveneighbors[eta][phi][mydepth]>0)
895  } // for (int subdet=1;...)
896  } // for (int phi=0;...)
897  } // for (int eta=0;...)
898  } // for (unsigned int depth=0;...)
899  FillUnphysicalHEHFBins(AboveNeighborsHotCellsByDepth);
900 
901  if (!makeDiagnostics_) return;
902  for (int i=0;i<500;++i) {
903  d_HBenergyVsNeighbor->Fill(i/50.,hbVsNeighbor[i]);
904  hbVsNeighbor[i]=0;
905  d_HEenergyVsNeighbor->Fill(i/50.,heVsNeighbor[i]);
906  heVsNeighbor[i]=0;
907  d_HOenergyVsNeighbor->Fill(i/50.,hoVsNeighbor[i]);
908  hoVsNeighbor[i]=0;
909  d_HFenergyVsNeighbor->Fill(i/50.,hfVsNeighbor[i]);
910  hfVsNeighbor[i]=0;
911  }
912 
913  return;
914 
915 } // void HcalHotCellMonitor::fillNevents_neighbor(const HcalTopology&)
i
int i
Definition: DBlmapReader.cc:9
HcalHotCellMonitor::aboveneighbors
int aboveneighbors[85][72][4]
Definition: HcalHotCellMonitor.h:107
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
HcalTopology::validDetId
bool validDetId(HcalSubdetector subdet, int ieta, int iphi, int depth) const
Definition: HcalTopology.cc:195
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalHotCellMonitor::d_HFenergyVsNeighbor
MonitorElement * d_HFenergyVsNeighbor
Definition: HcalHotCellMonitor.h:121
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
HcalHotCellMonitor::d_HBenergyVsNeighbor
MonitorElement * d_HBenergyVsNeighbor
Definition: HcalHotCellMonitor.h:118
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
HcalHotCellMonitor::d_HOenergyVsNeighbor
MonitorElement * d_HOenergyVsNeighbor
Definition: HcalHotCellMonitor.h:120
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
HcalHotCellMonitor::d_HEenergyVsNeighbor
MonitorElement * d_HEenergyVsNeighbor
Definition: HcalHotCellMonitor.h:119
void HcalHotCellMonitor::fillNevents_persistentenergy ( const HcalTopology &  topology)
private

Definition at line 682 of file HcalHotCellMonitor.cc.

References abovepersistent, abovepersistentET, AbovePersistentETThresholdCellsByDepth, AbovePersistentThresholdCellsByDepth, HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, HLT_FULL_cff::depth, eta, HcalObjRepresent::FillUnphysicalHEHFBins(), h, HcalForward, HcalBaseDQMonitor::ievt_, HcalBaseDQMonitor::levt_, minEvents_, phi, test_energy_, test_et_, and HcalTopology::validDetId().

Referenced by endLuminosityBlock().

682  {
683  // Fill Histograms showing rechits with energies > some threshold for N consecutive events
684 
685  if (levt_<minEvents_) return;
686 
687  if (debug_>0)
688  std::cout <<"<HcalHotCellMonitor::fillNevents_persistentenergy> FILLING PERSISTENT ENERGY PLOTS"<<std::endl;
689 
690  if (test_energy_) {
691  for (unsigned int h=0;h<AbovePersistentThresholdCellsByDepth.depth.size();++h)
692  AbovePersistentThresholdCellsByDepth.depth[h]->setBinContent(0,0,ievt_);
693 
694  int ieta=0;
695  int iphi=0;
696  int etabins=0;
697  int phibins=0;
698 
699  for (unsigned int depth=0;depth<AbovePersistentThresholdCellsByDepth.depth.size();++depth) {
700  etabins=AbovePersistentThresholdCellsByDepth.depth[depth]->getNbinsX();
701  phibins=AbovePersistentThresholdCellsByDepth.depth[depth]->getNbinsY();
702 
703  for (int eta=0;eta<etabins;++eta) {
704  for (int phi=0;phi<phibins;++phi) {
705  iphi=phi+1;
706  for (int subdet=1;subdet<=4;++subdet) {
707  ieta=CalcIeta((HcalSubdetector)subdet,eta,depth+1); //converts bin to ieta
708  if (ieta==-9999) continue;
709  if (!(topology.validDetId((HcalSubdetector)subdet, ieta, iphi, depth+1)))
710  continue;
711  if (subdet==HcalForward) // shift HcalForward ieta by 1 for filling purposes
712  ieta<0 ? ieta-- : ieta++;
713 
714  // MUST BE ABOVE ENERGY THRESHOLD FOR ALL N EVENTS in a luminosity block
715  if (abovepersistent[eta][phi][depth]<levt_) {
716  abovepersistent[eta][phi][depth]=0;
717  continue;
718  }
719  if (debug_>0) std::cout <<"HOT CELL; PERSISTENT ENERGY at subdet = "<<subdet<<", eta = "<<ieta<<", phi = "<<iphi<<" depth = "<<depth<<std::endl;
720  AbovePersistentThresholdCellsByDepth.depth[depth]->Fill(ieta,iphi,abovepersistent[eta][phi][depth]);
721  AbovePersistentThresholdCellsByDepth.depth[depth]->setBinContent(0,0,ievt_);
722  abovepersistent[eta][phi][depth]=0; // reset counter
723  } // for (int subdet=1; subdet<=4;++subdet)
724  } // for (int phi=0;...)
725  } // for (int eta=0;...)
726  } // for (unsigned int depth=0;...)
727  FillUnphysicalHEHFBins(AbovePersistentThresholdCellsByDepth);
728  } // if (test_energy_)
729 
730  if (test_et_) {
731  for (unsigned int h=0;h<AbovePersistentETThresholdCellsByDepth.depth.size();++h)
732  AbovePersistentETThresholdCellsByDepth.depth[h]->setBinContent(0,0,ievt_);
733 
734  int ieta=0;
735  int iphi=0;
736  int etabins=0;
737  int phibins=0;
738 
739  for (unsigned int depth=0;depth<AbovePersistentETThresholdCellsByDepth.depth.size();++depth) {
740  etabins=AbovePersistentETThresholdCellsByDepth.depth[depth]->getNbinsX();
741  phibins=AbovePersistentETThresholdCellsByDepth.depth[depth]->getNbinsY();
742 
743  for (int eta=0;eta<etabins;++eta) {
744  for (int phi=0;phi<phibins;++phi) {
745  iphi=phi+1;
746  for (int subdet=1;subdet<=4;++subdet) {
747  ieta=CalcIeta((HcalSubdetector)subdet,eta,depth+1); //converts bin to ieta
748  if (ieta==-9999) continue;
749  if (!(topology.validDetId((HcalSubdetector)subdet, ieta, iphi, depth+1)))
750  continue;
751  if (subdet==HcalForward) // shift HcalForward ieta by 1 for filling purposes
752  ieta<0 ? ieta-- : ieta++;
753 
754  // MUST BE ABOVE ET THRESHOLD FOR ALL N EVENTS in a luminosity block
755  if (abovepersistentET[eta][phi][depth]<levt_) {
756  abovepersistentET[eta][phi][depth]=0;
757  continue;
758  }
759  if (debug_>0) std::cout <<"HOT CELL; PERSISTENT ENERGY at subdet = "<<subdet<<", eta = "<<ieta<<", phi = "<<iphi<<" depth = "<<depth<<std::endl;
760  AbovePersistentETThresholdCellsByDepth.depth[depth]->Fill(ieta,iphi,abovepersistentET[eta][phi][depth]);
761  AbovePersistentETThresholdCellsByDepth.depth[depth]->setBinContent(0,0,ievt_);
762  abovepersistentET[eta][phi][depth]=0; // reset counter
763  } // for (int subdet=1; subdet<=4;++subdet)
764  } // for (int phi=0;...)
765  } // for (int eta=0;...)
766  } // for (unsigned int depth=0;...)
767  FillUnphysicalHEHFBins(AbovePersistentETThresholdCellsByDepth);
768 
769  } // if (test_et_)
770  // Add test_ET
771  return;
772 } // void HcalHotCellMonitor::fillNevents_persistentenergy(const HcalTopology&)
HcalHotCellMonitor::abovepersistent
int abovepersistent[85][72][4]
Definition: HcalHotCellMonitor.h:110
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
HcalTopology::validDetId
bool validDetId(HcalSubdetector subdet, int ieta, int iphi, int depth) const
Definition: HcalTopology.cc:195
HcalHotCellMonitor::abovepersistentET
int abovepersistentET[85][72][4]
Definition: HcalHotCellMonitor.h:111
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
EtaPhiHists AbovePersistentETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:104
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::fillNevents_problemCells ( const HcalTopology &  topology)
private

Definition at line 922 of file HcalHotCellMonitor.cc.

References AboveEnergyThresholdCellsByDepth, AboveETThresholdCellsByDepth, AboveNeighborsHotCellsByDepth, AbovePersistentETThresholdCellsByDepth, AbovePersistentThresholdCellsByDepth, funct::abs(), HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQMonitor::currentLS, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, HLT_FULL_cff::depth, eta, MonitorElement::Fill(), HcalBaseDQMonitor::ievt_, isHB(), isHE(), isHF(), isHO(), HcalBaseDQMonitor::levt_, minErrorFlag_, phi, HcalBaseDQMonitor::ProblemsCurrentLB, HcalBaseDQMonitor::ProblemsVsLB, HcalBaseDQMonitor::ProblemsVsLB_HB, HcalBaseDQMonitor::ProblemsVsLB_HBHEHF, HcalBaseDQMonitor::ProblemsVsLB_HE, HcalBaseDQMonitor::ProblemsVsLB_HF, HcalBaseDQMonitor::ProblemsVsLB_HO, test_energy_, test_et_, test_neighbor_, and test_persistent_.

Referenced by endLuminosityBlock().

922  {
923  if (debug_>0)
924  std::cout <<"<HcalHotCellMonitor::fillNevents_problemCells> FILLING PROBLEM CELL PLOTS"<<std::endl;
925 
926  if (ievt_==0) return; // no events; no need to bother with this
927 
928  int ieta=0;
929  int etabins=0;
930  int phibins=0;
931  bool problemvalue=false;
932 
933  // Count problem cells in each subdetector
934  int NumBadHB=0;
935  int NumBadHE=0;
936  int NumBadHO=0;
937  int NumBadHF=0;
938  int NumBadHO0=0;
939  int NumBadHO12=0;
940  int NumBadHFLUMI=0;
941 
942  unsigned int DEPTH = 0;
943 
944  if (test_persistent_)
945  {
946  if (test_energy_)
947  DEPTH = AbovePersistentThresholdCellsByDepth.depth.size();
948  else if (test_et_)
949  DEPTH = AbovePersistentETThresholdCellsByDepth.depth.size();
950  }
951  else if (test_energy_ && DEPTH==0) DEPTH = AboveEnergyThresholdCellsByDepth.depth.size();
952  else if (test_et_ && DEPTH==0) DEPTH = AboveETThresholdCellsByDepth.depth.size();
953  else if (test_neighbor_ && DEPTH==0) DEPTH = AboveNeighborsHotCellsByDepth.depth.size();
954 
955  if (DEPTH==0) return;
956 
957  for (unsigned int depth=0;depth<DEPTH;++depth)
958  {
959  if (test_persistent_)
960  {
961  if (test_energy_)
962  {
963  etabins=AbovePersistentThresholdCellsByDepth.depth[depth]->getNbinsX();
964  phibins=AbovePersistentThresholdCellsByDepth.depth[depth]->getNbinsY();
965  }
966  else if (test_et_)
967  {
968  etabins=AbovePersistentETThresholdCellsByDepth.depth[depth]->getNbinsX();
969  phibins=AbovePersistentETThresholdCellsByDepth.depth[depth]->getNbinsY();
970  }
971  }
972 
973  if (test_neighbor_ && (etabins==0 || phibins==0))
974  {
975  etabins=AboveNeighborsHotCellsByDepth.depth[depth]->getNbinsX();
976  phibins=AboveNeighborsHotCellsByDepth.depth[depth]->getNbinsY();
977  }
978 
979  if (test_energy_ && (etabins==0 || phibins==0))
980  {
981  etabins=AboveEnergyThresholdCellsByDepth.depth[depth]->getNbinsX();
982  phibins=AboveEnergyThresholdCellsByDepth.depth[depth]->getNbinsY();
983  }
984 
985  if (test_et_ && (etabins==0 || phibins==0))
986  {
987  etabins=AboveETThresholdCellsByDepth.depth[depth]->getNbinsX();
988  phibins=AboveETThresholdCellsByDepth.depth[depth]->getNbinsY();
989  }
990 
991  for (int eta=0;eta<etabins;++eta)
992  {
993  ieta=CalcIeta(eta,depth+1);
994  if (ieta==-9999) continue;
995  for (int phi=0;phi<phibins;++phi)
996  {
997  if (abs(ieta)>20 && phi%2==1) continue; //skip non-physical cells
998  else if (abs(ieta)>39 && (phi+1)%4!=3) continue;
999  // find problem rate for particular cell
1000  problemvalue=false;
1001  if (test_energy_ && test_persistent_ && AbovePersistentThresholdCellsByDepth.depth[depth]->getBinContent(eta+1,phi+1)>minErrorFlag_*ievt_)
1002  problemvalue=true;
1003  if (test_neighbor_ && AboveNeighborsHotCellsByDepth.depth[depth]->getBinContent(eta+1,phi+1)>minErrorFlag_*ievt_)
1004  problemvalue=true;
1005  if (test_energy_ && AboveEnergyThresholdCellsByDepth.depth[depth]->getBinContent(eta+1,phi+1)>minErrorFlag_*ievt_)
1006  problemvalue=true;
1007  if (test_et_ && AboveETThresholdCellsByDepth.depth[depth]->getBinContent(eta+1,phi+1)>minErrorFlag_*ievt_)
1008  problemvalue=true;
1009  if (test_et_ && test_persistent_ && AbovePersistentETThresholdCellsByDepth.depth[depth]->getBinContent(eta+1,phi+1)>minErrorFlag_*ievt_)
1010  problemvalue=true;
1011  if (problemvalue==false) continue;
1012  if (isHB(eta,depth+1)) ++NumBadHB;
1013  else if (isHE(eta,depth+1))
1014  ++NumBadHE;
1015  else if (isHO(eta,depth+1))
1016  {
1017  ++NumBadHO;
1018  if (abs(ieta)<5) ++NumBadHO0;
1019  else ++NumBadHO12;
1020  }
1021  else if (isHF(eta,depth+1))
1022  {
1023  ++NumBadHF;
1024  if (depth+1==1 && (abs(ieta)==33 || abs(ieta)==34)) ++NumBadHFLUMI;
1025  else if (depth+1==2 && (abs(ieta)==35 || abs(ieta)==36)) ++NumBadHFLUMI;
1026  }
1027  } // for (int phi=0;...)
1028  } //for (int eta=0;...)
1029  } // for (int depth=0;...)
1030 
1031  if (debug_>2) std::cout <<"<HcalHotCellMonitor::fillNevents_problemCells> Num Bad HB = "<<NumBadHB<<" Num Bad HE = "<<NumBadHE<<" Num Bad HO = "<<NumBadHO<<" Num Bad HF = "<<NumBadHF<<" CURRENT LS = "<<currentLS<<std::endl;
1032  // Fill number of problem cells
1033  ProblemsVsLB_HB->Fill(currentLS,NumBadHB);
1034  ProblemsVsLB_HE->Fill(currentLS,NumBadHE);
1035  ProblemsVsLB_HO->Fill(currentLS,NumBadHO);
1036  ProblemsVsLB_HF->Fill(currentLS,NumBadHF);
1037  ProblemsVsLB_HBHEHF->Fill(currentLS,NumBadHB+NumBadHE+NumBadHF);
1038  ProblemsVsLB->Fill(currentLS,NumBadHB+NumBadHE+NumBadHO+NumBadHF);
1039 
1040  ProblemsCurrentLB->Fill(-1,-1,levt_);
1041  ProblemsCurrentLB->Fill(0,0,NumBadHB);
1042  ProblemsCurrentLB->Fill(1,0,NumBadHE);
1043  ProblemsCurrentLB->Fill(2,0,NumBadHO);
1044  ProblemsCurrentLB->Fill(3,0,NumBadHF);
1045  ProblemsCurrentLB->Fill(4,0,NumBadHO0);
1046  ProblemsCurrentLB->Fill(5,0,NumBadHO12);
1047  ProblemsCurrentLB->Fill(6,0,NumBadHFLUMI);
1048 
1049 } // void HcalHotCellMonitor::fillNevents_problemCells(const HcalTopology&)
isHO
bool isHO(int etabin, int depth)
Definition: HcalEtaPhiHists.h:479
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
HcalBaseDQMonitor::ProblemsCurrentLB
MonitorElement * ProblemsCurrentLB
Definition: HcalBaseDQMonitor.h:112
HcalBaseDQMonitor::ProblemsVsLB_HF
MonitorElement * ProblemsVsLB_HF
Definition: HcalBaseDQMonitor.h:111
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
EtaPhiHists AbovePersistentETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:104
isHE
bool isHE(int etabin, int depth)
Definition: HcalEtaPhiHists.h:429
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalBaseDQMonitor::ProblemsVsLB_HBHEHF
MonitorElement * ProblemsVsLB_HBHEHF
Definition: HcalBaseDQMonitor.h:111
HcalBaseDQMonitor::ProblemsVsLB_HB
MonitorElement * ProblemsVsLB_HB
Definition: HcalBaseDQMonitor.h:110
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
isHB
bool isHB(int etabin, int depth)
Definition: HcalEtaPhiHists.h:407
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
isHF
bool isHF(int etabin, int depth)
Definition: HcalEtaPhiHists.h:455
HcalBaseDQMonitor::ProblemsVsLB_HE
MonitorElement * ProblemsVsLB_HE
Definition: HcalBaseDQMonitor.h:110
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
HcalBaseDQMonitor::ProblemsVsLB_HO
MonitorElement * ProblemsVsLB_HO
Definition: HcalBaseDQMonitor.h:111
HcalBaseDQMonitor::ProblemsVsLB
MonitorElement * ProblemsVsLB
Definition: HcalBaseDQMonitor.h:109
void HcalHotCellMonitor::periodicReset ( )

Definition at line 1114 of file HcalHotCellMonitor.cc.

References AboveEnergyThresholdCellsByDepth, AboveETThresholdCellsByDepth, AboveNeighborsHotCellsByDepth, AbovePersistentETThresholdCellsByDepth, AbovePersistentThresholdCellsByDepth, HcalBaseDQMonitor::makeDiagnostics_, EtaPhiHists::Reset(), test_energy_, test_et_, test_neighbor_, test_persistent_, and zeroCounters().

1115 {
1116 
1117  // first reset base class objects
1118  //FIX HcalBaseMonitor::periodicReset();
1119 
1120  // then reset the temporary histograms
1121  zeroCounters();
1122 
1123  // now reset all the MonitorElements
1124 
1125  // resetting eta-phi histograms
1126  if (test_neighbor_)
1127  AboveNeighborsHotCellsByDepth.Reset();
1128  if (test_energy_ || makeDiagnostics_)
1129  AboveEnergyThresholdCellsByDepth.Reset();
1130  if (test_et_ || makeDiagnostics_)
1131  AboveETThresholdCellsByDepth.Reset();
1132  if (test_persistent_)
1133  {
1134  AbovePersistentThresholdCellsByDepth.Reset();
1135  if (test_et_)
1136  AbovePersistentETThresholdCellsByDepth.Reset();
1137  }
1138  return;
1139 }
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
EtaPhiHists::Reset
void Reset(void)
Definition: HcalEtaPhiHists.h:179
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
EtaPhiHists AbovePersistentETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:104
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
void HcalHotCellMonitor::processEvent ( const HBHERecHitCollection &  hbHits,
const HORecHitCollection &  hoHits,
const HFRecHitCollection &  hfHits,
const HcalTopology &  topology 
)

Definition at line 419 of file HcalHotCellMonitor.cc.

References gather_cfg::cout, HcalBaseDQMonitor::debug_, processEvent_rechitenergy(), test_energy_, test_et_, and test_persistent_.

Referenced by analyze().

422  {
423 
424  if (debug_>1) std::cout <<"<HcalHotCellMonitor::processEvent> Processing event..."<<std::endl;
425 
426  // Search for hot cells above a certain energy
427  if (test_energy_ || test_et_ || test_persistent_) {
428  processEvent_rechitenergy(hbHits, hoHits,hfHits, topology);
429  }
430 
431  return;
432 } // void HcalHotCellMonitor::processEvent(...)
HcalHotCellMonitor::processEvent_rechitenergy
void processEvent_rechitenergy(const HBHERecHitCollection &hbheHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, const HcalTopology &topology)
Definition: HcalHotCellMonitor.cc:438
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
void HcalHotCellMonitor::processEvent_rechitenergy ( const HBHERecHitCollection &  hbheHits,
const HORecHitCollection &  hoHits,
const HFRecHitCollection &  hfHits,
const HcalTopology &  topology 
)

Definition at line 438 of file HcalHotCellMonitor.cc.

References aboveenergy, AboveEnergyThresholdCellsByDepth, aboveet, AboveETThresholdCellsByDepth, AboveNeighborsHotCellsByDepth, abovepersistent, abovepersistentET, AbovePersistentThresholdCellsByDepth, funct::abs(), edm::SortedCollection< T, SORT >::begin(), CalcEtaBin(), gather_cfg::cout, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, HLT_FULL_cff::depth, edm::SortedCollection< T, SORT >::end(), HcalTopology::etaRange(), EnergyCorrector::etas, excludeHORing2_, HBenergyThreshold_, HBETThreshold_, HBHENeighborParams_, HBpersistentETThreshold_, HBpersistentThreshold_, HcalBarrel, HcalEndcap, HEenergyThreshold_, HEETThreshold_, HEpersistentETThreshold_, HEpersistentThreshold_, HFenergyThreshold_, HFETThreshold_, HFfarfwdScale_, HFNeighborParams_, HFpersistentETThreshold_, HFpersistentThreshold_, HOenergyThreshold_, HOETThreshold_, HONeighborParams_, HOpersistentETThreshold_, HOpersistentThreshold_, i, HcalDetId::ieta(), isSiPM(), HcalBaseDQMonitor::makeDiagnostics_, processHit_rechitNeighbors(), SiPMscale_, test_energy_, test_et_, test_neighbor_, and dtDQMClient_cfg::threshold.

Referenced by processEvent().

441  {
442 
443  // Looks at rechits of cells and compares to threshold energies.
444  // Cells above thresholds get marked as hot candidates
445 
446  if (debug_>1) std::cout <<"<HcalHotCellMonitor::processEvent_rechitenergy> Processing rechits..."<<std::endl;
447 
448  // loop over HBHE
449  for (HBHERecHitCollection::const_iterator HBHEiter=hbheHits.begin(); HBHEiter!=hbheHits.end(); ++HBHEiter) { // loop over all hits
450  float en = HBHEiter->energy();
451  //float ti = HBHEiter->time();
452 
453  HcalDetId id(HBHEiter->detid().rawId());
454  int ieta = id.ieta();
455  int iphi = id.iphi();
456  int depth = id.depth();
457  std::pair<double,double> etas = topology.etaRange(id.subdet(),abs(ieta));
458  double fEta=fabs(0.5*(etas.first+etas.second));
459  float et = en/cosh(fEta);
460 
461  if (test_neighbor_ || makeDiagnostics_) {
462  processHit_rechitNeighbors(HBHEiter, hbheHits, HBHENeighborParams_, topology);
463  }
464  if (id.subdet()==HcalBarrel) {
465  if (en>=HBenergyThreshold_)
466  ++aboveenergy[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
467  if (et>=HBETThreshold_)
468  ++aboveet[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
469  if (test_energy_ && en>=HBpersistentThreshold_)
470  ++abovepersistent[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
471  if (test_et_ && et>=HBpersistentETThreshold_)
472  ++abovepersistentET[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
473  } else if (id.subdet()==HcalEndcap) {
474  if (en>=HEenergyThreshold_)
475  ++aboveenergy[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
476  if (et>=HEETThreshold_)
477  ++aboveet[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
478  if (test_energy_)
479  if (en>=HEpersistentThreshold_)
480  ++abovepersistent[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
481  if (test_et_)
482  if (et>=HEpersistentETThreshold_)
483  ++abovepersistentET[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
484  }
485  } //for (HBHERecHitCollection::const_iterator HBHEiter=...)
486 
487  // loop over HO
488  for (HORecHitCollection::const_iterator HOiter=hoHits.begin(); HOiter!=hoHits.end(); ++HOiter) { // loop over all hits
489  float en = HOiter->energy();
490 
491  HcalDetId id(HOiter->detid().rawId());
492  int ieta = id.ieta();
493  int iphi = id.iphi();
494  int depth = id.depth();
495  std::pair<double,double> etas = topology.etaRange(id.subdet(),abs(ieta));
496  double fEta=fabs(0.5*(etas.first+etas.second));
497  float et = en/cosh(fEta);
498 
499  if (test_neighbor_ || makeDiagnostics_)
500  processHit_rechitNeighbors(HOiter, hoHits, HONeighborParams_, topology);
501 
502  if (isSiPM(ieta,iphi,depth)) {
503  if (en>=HOenergyThreshold_*SiPMscale_)
504  ++aboveenergy[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
505  if (et>=HOETThreshold_*SiPMscale_)
506  ++aboveet[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
507  if (test_energy_)
508  if (en>=HOpersistentThreshold_*SiPMscale_)
509  ++abovepersistent[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
510  if (test_et_)
511  if (et>=HOpersistentETThreshold_*SiPMscale_)
512  ++abovepersistentET[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
513  } else {
514  // Skip HO ring 2 when required
515  if (abs(ieta)>10 && excludeHORing2_==true)
516  continue;
517 
518  if (en>=HOenergyThreshold_)
519  ++aboveenergy[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
520  if (et>=HOETThreshold_)
521  ++aboveet[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
522  if (test_energy_)
523  if (en>=HOpersistentThreshold_)
524  ++abovepersistent[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
525  if (test_et_)
526  if (en>=HOpersistentETThreshold_)
527  ++abovepersistentET[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
528  }
529  }
530 
531  // loop over HF
532  for (HFRecHitCollection::const_iterator HFiter=hfHits.begin(); HFiter!=hfHits.end(); ++HFiter) { // loop over all hits
533  float en = HFiter->energy();
534  float threshold=HFenergyThreshold_;
535  float threshold_pers = HFpersistentThreshold_;
536  float etthreshold=HFETThreshold_;
537  HcalDetId id(HFiter->detid().rawId());
538  int ieta = id.ieta();
539  int iphi = id.iphi();
540  int depth = id.depth();
541  std::pair<double,double> etas = topology.etaRange(id.subdet(),abs(ieta));
542  double fEta=fabs(0.5*(etas.first+etas.second));
543  float et = en/cosh(fEta);
544 
545  if (test_neighbor_ || makeDiagnostics_)
546  processHit_rechitNeighbors(HFiter, hfHits, HFNeighborParams_, topology);
547 
548  if (abs(ieta)>39) { // increase the thresholds in far-forward part of HF
549  threshold*=HFfarfwdScale_;
550  threshold_pers*=HFfarfwdScale_;
551  }
552 
553  if (en>=threshold)
554  ++aboveenergy[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
555  if (et>=etthreshold)
556  ++aboveet[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
557  if (test_energy_) {
558  if (en>=threshold_pers)
559  ++abovepersistent[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
560  }
561  if (test_et_) {
562  if (et>=HFpersistentETThreshold_)
563  ++abovepersistentET[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1];
564  }
565  }
566 
567  // call update every event -- still necessary?
568 
569  for (unsigned int i=0;i<AbovePersistentThresholdCellsByDepth.depth.size();++i)
570  AbovePersistentThresholdCellsByDepth.depth[i]->update();
571  for (unsigned int i=0;i<AboveEnergyThresholdCellsByDepth.depth.size();++i)
572  AboveEnergyThresholdCellsByDepth.depth[i]->update();
573  for (unsigned int i=0;i<AboveETThresholdCellsByDepth.depth.size();++i)
574  AboveETThresholdCellsByDepth.depth[i]->update();
575  for (unsigned int i=0;i<AboveNeighborsHotCellsByDepth.depth.size();++i)
576  AboveNeighborsHotCellsByDepth.depth[i]->update();
577 
578  return;
579 } // void HcalHotCellMonitor::processEvent_rechitenergy
HcalHotCellMonitor::abovepersistent
int abovepersistent[85][72][4]
Definition: HcalHotCellMonitor.h:110
i
int i
Definition: DBlmapReader.cc:9
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
HcalHotCellMonitor::HONeighborParams_
hotNeighborParams HONeighborParams_
Definition: HcalHotCellMonitor.h:128
HcalHotCellMonitor::abovepersistentET
int abovepersistentET[85][72][4]
Definition: HcalHotCellMonitor.h:111
edm::SortedCollection< HBHERecHit >::const_iterator
std::vector< HBHERecHit >::const_iterator const_iterator
Definition: SortedCollection.h:82
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::aboveenergy
int aboveenergy[85][72][4]
Definition: HcalHotCellMonitor.h:108
isSiPM
bool isSiPM(int ieta, int iphi, int depth)
Definition: HcalEtaPhiHists.h:489
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HcalHotCellMonitor::HBHENeighborParams_
hotNeighborParams HBHENeighborParams_
Definition: HcalHotCellMonitor.h:128
dtDQMClient_cfg.threshold
tuple threshold
Definition: dtDQMClient_cfg.py:16
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
HcalHotCellMonitor::processHit_rechitNeighbors
void processHit_rechitNeighbors(R &rechititer, C &collection, hotNeighborParams &params, const HcalTopology &topology)
HcalHotCellMonitor::aboveet
int aboveet[85][72][4]
Definition: HcalHotCellMonitor.h:109
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
HcalTopology::etaRange
std::pair< double, double > etaRange(HcalSubdetector subdet, int ieta) const
Definition: HcalTopology.cc:892
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
HcalHotCellMonitor::HFNeighborParams_
hotNeighborParams HFNeighborParams_
Definition: HcalHotCellMonitor.h:128
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
template<class R , class C >
void HcalHotCellMonitor::processHit_rechitNeighbors ( R &  rechititer,
C &  collection,
hotNeighborParams &  params,
const HcalTopology &  topology 
)

Referenced by processEvent_rechitenergy().

template<class RECHIT , class RECHITCOLLECTION >
void HcalHotCellMonitor::processHit_rechitNeighbors ( RECHIT &  rechit,
RECHITCOLLECTION &  coll,
hotNeighborParams &  params,
const HcalTopology &  topology 
)

Definition at line 585 of file HcalHotCellMonitor.cc.

References aboveneighbors, funct::abs(), CalcEtaBin(), d_HBenergyVsNeighbor, d_HEenergyVsNeighbor, d_HFenergyVsNeighbor, d_HOenergyVsNeighbor, hotNeighborParams::DeltaDepth, hotNeighborParams::DeltaIeta, hotNeighborParams::DeltaIphi, HLT_FULL_cff::depth, HcalTopology::etaRange(), EnergyCorrector::etas, MonitorElement::Fill(), hbVsNeighbor, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, heVsNeighbor, hfVsNeighbor, hotNeighborParams::HotEnergyFrac, hoVsNeighbor, HcalDetId::ieta(), HcalBaseDQMonitor::makeDiagnostics_, hotNeighborParams::maxEnergy, hotNeighborParams::minCellEnergy, hotNeighborParams::minNeighborEnergy, and pftools::RECHIT.

588  {
589  // Compares energy to energy of neighboring cells.
590  // This is a slightly simplified version of D0's NADA algorithm
591  // 17 June 2009 -- this needs major work. I'm not sure I have the [eta][phi][depth] array mapping correct everywhere.
592  // Maybe even tear it apart and start again?
593 
594  int ieta, iphi, depth;
595  float en;
596 
597  int neighborsfound=0;
598  float enNeighbor=0;
599 
600  en = rechit->energy();
601  HcalDetId id(rechit->detid().rawId());
602  ieta = id.ieta();
603  iphi = id.iphi();
604  depth = id.depth();
605 
606  std::pair<double,double> etas = topology.etaRange(id.subdet(),abs(ieta));
607  double fEta=fabs(0.5*(etas.first+etas.second));
608 
609  float et = en/cosh(fEta);
610 
611  // Case 0: ET too low to trigger hot cell check
612  if (et<=params.minCellEnergy) return;
613 
614  // Case 1: above threshold energy; always count as hot
615  if (et>params.maxEnergy) {
616  if (makeDiagnostics_) {
617  // fill overflow bin when energy > max threshold
618  if (id.subdet()==HcalBarrel) d_HBenergyVsNeighbor->Fill(1000);
619  else if (id.subdet()==HcalEndcap) d_HEenergyVsNeighbor->Fill(1000);
620  else if (id.subdet()==HcalOuter) d_HOenergyVsNeighbor->Fill(1000);
621  else if (id.subdet()==HcalForward) d_HFenergyVsNeighbor->Fill(1000);
622  }
623  return;
624  }
625 
626  // Case 2: Search keys for neighboring cells
627 
628  neighborsfound=0;
629  enNeighbor=0;
630 
631  int mydeltaphi=params.DeltaIphi;
632  // scale appropriately for larger cells at higher eta values
633  if (abs(ieta)>39) mydeltaphi*=4;
634  else if (abs(ieta)>20) mydeltaphi*=2;
635 
636  for (int nD=-1*params.DeltaDepth;nD<=params.DeltaDepth;++nD) {
637  for (int nP =-1*mydeltaphi;nP<=mydeltaphi;++nP) {
638  for (int nE =-1*params.DeltaIeta;nE<=params.DeltaIeta;++nE) {
639  if (nD==0 && nE==0 && nP==0)
640  continue; // don't count the cell itself
641  int myphi=(nP+iphi)%72;
642  HcalDetId myid((HcalSubdetector)(1), nE+ieta, myphi, nD+depth); // HB
643  RECHIT part=coll.find(myid);
644  if (part==coll.end())
645  continue;
646  if (part->energy()<params.minNeighborEnergy)
647  continue;
648  ++neighborsfound;
649  enNeighbor+=part->energy();
650  } // loop over nE (neighbor eta)
651  } // loop over nP (neighbor phi)
652  } // loop over nD depths
653 
654  // Case 2a: Not enough good neighbors found -- do we want to implement this?
655  //if (neighborsfound==0)
656  // return;
657 
658  // Case 2b: (avg. neighbor energy)/energy too large for cell to be considered hot
659  if (makeDiagnostics_) {
660  int myval=(int)(enNeighbor/en*50);
661  if (myval<0) myval=0;
662  if (myval>499) myval=499;
663  if (enNeighbor/en<0 || enNeighbor/en>=10) return;
664  if (id.subdet()==HcalBarrel) ++hbVsNeighbor[myval];
665  else if (id.subdet()==HcalEndcap) ++heVsNeighbor[myval];
666  else if (id.subdet()==HcalOuter) ++hoVsNeighbor[myval];
667  else if (id.subdet()==HcalForward) ++hfVsNeighbor[myval];
668  }
669  if ((1.*enNeighbor/en)>params.HotEnergyFrac && en>0 && enNeighbor>0)
670  return;
671 
672  // Case 2c: Tests passed; cell marked as hot
673  aboveneighbors[CalcEtaBin(id.subdet(),ieta,depth)][iphi-1][depth-1]++;
674 
675  return;
676 } // void HcalHotCellMonitor::processEvent_rechitneighbor
HcalHotCellMonitor::aboveneighbors
int aboveneighbors[85][72][4]
Definition: HcalHotCellMonitor.h:107
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
HcalHotCellMonitor::d_HFenergyVsNeighbor
MonitorElement * d_HFenergyVsNeighbor
Definition: HcalHotCellMonitor.h:121
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HLT_FULL_cff.depth
tuple depth
Definition: HLT_FULL_cff.py:11144
HcalHotCellMonitor::d_HBenergyVsNeighbor
MonitorElement * d_HBenergyVsNeighbor
Definition: HcalHotCellMonitor.h:118
coll
JetCorrectorParametersCollection coll
Definition: classes.h:10
part
part
Definition: HCALResponse.h:20
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
HcalHotCellMonitor::d_HOenergyVsNeighbor
MonitorElement * d_HOenergyVsNeighbor
Definition: HcalHotCellMonitor.h:120
HcalTopology::etaRange
std::pair< double, double > etaRange(HcalSubdetector subdet, int ieta) const
Definition: HcalTopology.cc:892
HcalHotCellMonitor::d_HEenergyVsNeighbor
MonitorElement * d_HEenergyVsNeighbor
Definition: HcalHotCellMonitor.h:119
void HcalHotCellMonitor::reset ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 297 of file HcalHotCellMonitor.cc.

References AboveEnergyThresholdCellsByDepth, AboveETThresholdCellsByDepth, AboveNeighborsHotCellsByDepth, AbovePersistentETThresholdCellsByDepth, AbovePersistentThresholdCellsByDepth, d_HBenergyVsNeighbor, d_HEenergyVsNeighbor, d_HFenergyVsNeighbor, d_HOenergyVsNeighbor, HcalBaseDQMonitor::makeDiagnostics_, HcalBaseDQMonitor::reset(), EtaPhiHists::Reset(), MonitorElement::Reset(), test_energy_, test_et_, test_neighbor_, test_persistent_, and zeroCounters().

Referenced by bookHistograms(), and setup().

298 {
299  HcalBaseDQMonitor::reset();
300  zeroCounters();
301 
302  // now reset all the MonitorElements
303 
304  // resetting eta-phi histograms
305  if (test_neighbor_ || makeDiagnostics_)
306  AboveNeighborsHotCellsByDepth.Reset();
307 
308  if (test_energy_ )
309  AboveEnergyThresholdCellsByDepth.Reset();
310 
311  if ( test_et_ )
312  AboveETThresholdCellsByDepth.Reset();
313 
314  if (test_persistent_)
315  {
316  if (test_energy_) AbovePersistentThresholdCellsByDepth.Reset();
317  if (test_et_) AbovePersistentETThresholdCellsByDepth.Reset();
318  }
319  if (makeDiagnostics_)
320  {
321  d_HBenergyVsNeighbor->Reset();
322  d_HEenergyVsNeighbor->Reset();
323  d_HOenergyVsNeighbor->Reset();
324  d_HFenergyVsNeighbor->Reset();
325  }
326 }
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
EtaPhiHists::Reset
void Reset(void)
Definition: HcalEtaPhiHists.h:179
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
EtaPhiHists AbovePersistentETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:104
HcalHotCellMonitor::d_HFenergyVsNeighbor
MonitorElement * d_HFenergyVsNeighbor
Definition: HcalHotCellMonitor.h:121
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
HcalHotCellMonitor::d_HBenergyVsNeighbor
MonitorElement * d_HBenergyVsNeighbor
Definition: HcalHotCellMonitor.h:118
HcalBaseDQMonitor::reset
virtual void reset(void)
Definition: HcalBaseDQMonitor.cc:127
HcalHotCellMonitor::d_HOenergyVsNeighbor
MonitorElement * d_HOenergyVsNeighbor
Definition: HcalHotCellMonitor.h:120
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:453
HcalHotCellMonitor::d_HEenergyVsNeighbor
MonitorElement * d_HEenergyVsNeighbor
Definition: HcalHotCellMonitor.h:119
void HcalHotCellMonitor::setup ( DQMStore::IBooker &  ib)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 112 of file HcalHotCellMonitor.cc.

References AboveEnergyThresholdCellsByDepth, AboveETThresholdCellsByDepth, AboveNeighborsHotCellsByDepth, AbovePersistentETThresholdCellsByDepth, AbovePersistentThresholdCellsByDepth, DQMStore::IBooker::book1D(), DQMStore::IBooker::bookFloat(), DQMStore::IBooker::bookInt(), DQMStore::IBooker::bookProfile(), gather_cfg::cout, d_HBenergyVsNeighbor, d_HEenergyVsNeighbor, d_HFenergyVsNeighbor, d_HOenergyVsNeighbor, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, MonitorElement::Fill(), MonitorElement::getTProfile(), HBenergyThreshold_, HBETThreshold_, HBpersistentETThreshold_, HBpersistentThreshold_, HEenergyThreshold_, HEETThreshold_, HEpersistentETThreshold_, HEpersistentThreshold_, HFenergyThreshold_, HFETThreshold_, HFpersistentETThreshold_, HFpersistentThreshold_, HOenergyThreshold_, HOETThreshold_, HOpersistentETThreshold_, HOpersistentThreshold_, HcalBaseDQMonitor::makeDiagnostics_, minErrorFlag_, minEvents_, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::ProblemsVsLB, HcalBaseDQMonitor::ProblemsVsLB_HB, HcalBaseDQMonitor::ProblemsVsLB_HBHEHF, HcalBaseDQMonitor::ProblemsVsLB_HE, HcalBaseDQMonitor::ProblemsVsLB_HF, HcalBaseDQMonitor::ProblemsVsLB_HO, reset(), DQMStore::IBooker::setCurrentFolder(), HcalBaseDQMonitor::setup(), setupDone_, HcalBaseDQMonitor::SetupEtaPhiHists(), HcalBaseDQMonitor::subdir_, test_energy_, test_et_, test_neighbor_, test_persistent_, and units().

Referenced by bookHistograms().

113 {
114  if (setupDone_)
115  return;
116  setupDone_ = true;
117  // Call base class setup
118  HcalBaseDQMonitor::setup(ib);
119 
120  if (debug_>1)
121  std::cout <<"<HcalHotCellMonitor::setup> Setting up histograms"<<std::endl;
122 
123  ib.setCurrentFolder(subdir_);
124 
125  MonitorElement* me;
126  me=ib.bookFloat("minErrorFractionPerLumiSection");
127  me->Fill(minErrorFlag_);
128  // Create plots of problems vs LB
129 
130  // 1D plots count number of bad cells vs. luminosity block
131  ProblemsVsLB=ib.bookProfile("TotalHotCells_HCAL_vs_LS",
132  "Total Number of Hot Hcal Cells vs lumi section",
133  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,10000);
134 
135  ProblemsVsLB_HB=ib.bookProfile("TotalHotCells_HB_vs_LS",
136  "Total Number of Hot HB Cells vs lumi section",
137  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
138  ProblemsVsLB_HE=ib.bookProfile("TotalHotCells_HE_vs_LS",
139  "Total Number of Hot HE Cells vs lumi section",
140  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
141  ProblemsVsLB_HO=ib.bookProfile("TotalHotCells_HO_vs_LS",
142  "Total Number of Hot HO Cells vs lumi section",
143  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,3000);
144  ProblemsVsLB_HF=ib.bookProfile("TotalHotCells_HF_vs_LS",
145  "Total Number of Hot HF Cells vs lumi section",
146  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,2000);
147  ProblemsVsLB_HBHEHF=ib.bookProfile("TotalHotCells_HBHEHF_vs_LS",
148  "Total Number of Hot HBHEHF Cells vs lumi section",
149  NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,2000);
150 
151  ProblemsVsLB->getTProfile()->SetMarkerStyle(20);
152  ProblemsVsLB_HB->getTProfile()->SetMarkerStyle(20);
153  ProblemsVsLB_HE->getTProfile()->SetMarkerStyle(20);
154  ProblemsVsLB_HO->getTProfile()->SetMarkerStyle(20);
155  ProblemsVsLB_HF->getTProfile()->SetMarkerStyle(20);
156  ProblemsVsLB_HBHEHF->getTProfile()->SetMarkerStyle(20);
157 
158  // Set up plots for each failure mode of hot cells
159  std::stringstream units; // We'll need to set the titles individually, rather than passing units to SetupEtaPhiHists (since this also would affect the name of the histograms)
160 
161  ib.setCurrentFolder(subdir_+"hot_rechit_above_threshold");
162  me=ib.bookInt("HotCellAboveThresholdTestEnabled");
163  me->Fill(0);
164 
165  if (test_energy_)
166  {
167  me->Fill(1);
168  SetupEtaPhiHists(ib,AboveEnergyThresholdCellsByDepth,
169  "Hot Cells Above Energy Threshold","");
170  //setMinMaxHists2D(AboveEnergyThresholdCellsByDepth,0.,1.);
171 
172  // set more descriptive titles for plots
173  units.str("");
174  units<<"Hot Cells: Depth 1 -- HB > "<<HBenergyThreshold_<<" GeV, HE > "<<HEenergyThreshold_<<" GeV, HF > "<<HFenergyThreshold_<<" GeV";
175  AboveEnergyThresholdCellsByDepth.depth[0]->setTitle(units.str().c_str());
176  units.str("");
177  units<<"Hot Cells: Depth 2 -- HB > "<<HBenergyThreshold_<<" GeV, HE > "<<HEenergyThreshold_<<" GeV, HF > "<<HFenergyThreshold_<<" GeV";
178  AboveEnergyThresholdCellsByDepth.depth[1]->setTitle(units.str().c_str());
179  units.str("");
180  units<<"Hot Cells: Depth 3 -- HE > "<<HEenergyThreshold_<<" GeV";
181  AboveEnergyThresholdCellsByDepth.depth[2]->setTitle(units.str().c_str());
182  units.str("");
183  units<<"Hot Cells: HO > "<<HOenergyThreshold_<<" GeV";
184  AboveEnergyThresholdCellsByDepth.depth[3]->setTitle(units.str().c_str());
185  units.str("");
186  }
187  if (test_et_)
188  {
189  me->Fill(1);
190  SetupEtaPhiHists(ib,AboveETThresholdCellsByDepth,
191  "Hot Cells Above ET Threshold","");
192  //setMinMaxHists2D(AboveETThresholdCellsByDepth,0.,1.);
193 
194  // set more descriptive titles for plots
195  units.str("");
196  units<<"Hot Cells: Depth 1 -- HB > "<<HBETThreshold_<<" GeV (ET), HE > "<<HEETThreshold_<<" GeV (ET), HF > "<<HFETThreshold_<<" GeV (ET)";
197  AboveETThresholdCellsByDepth.depth[0]->setTitle(units.str().c_str());
198  units.str("");
199  units<<"Hot Cells: Depth 2 -- HB > "<<HBETThreshold_<<" GeV (ET), HE > "<<HEETThreshold_<<" GeV (ET), HF > "<<HFETThreshold_<<" GeV (ET)";
200  AboveETThresholdCellsByDepth.depth[1]->setTitle(units.str().c_str());
201  units.str("");
202  units<<"Hot Cells: Depth 3 -- HE > "<<HEETThreshold_<<" GeV (ET)";
203  AboveETThresholdCellsByDepth.depth[2]->setTitle(units.str().c_str());
204  units.str("");
205  units<<"Hot Cells: HO > "<<HOETThreshold_<<" GeV (ET)";
206  AboveETThresholdCellsByDepth.depth[3]->setTitle(units.str().c_str());
207  units.str("");
208  }
209 
210  ib.setCurrentFolder(subdir_+"hot_rechit_always_above_threshold");
211  me=ib.bookInt("PersistentHotCellTestEnabled");
212  me->Fill(0);
213  if (test_persistent_)
214  {
215  me->Fill(1);
216  me=ib.bookInt("minEventsPerLS");
217  me->Fill(minEvents_);
218 
219  if (test_energy_) {
220  SetupEtaPhiHists(ib,AbovePersistentThresholdCellsByDepth,
221  "Hot Cells Persistently Above Energy Threshold","");
222  //setMinMaxHists2D(AbovePersistentThresholdCellsByDepth,0.,1.);
223 
224  // set more descriptive titles for plots
225  units.str("");
226  units<<"Hot Cells: Depth 1 -- HB > "<<HBpersistentThreshold_<<" GeV, HE > "<<HEpersistentThreshold_<<", HF > "<<HFpersistentThreshold_<<" GeV for 1 full Lumi Block";
227  AbovePersistentThresholdCellsByDepth.depth[0]->setTitle(units.str().c_str());
228  units.str("");
229  units<<"Hot Cells: Depth 2 -- HB > "<<HBpersistentThreshold_<<" GeV, HE > "<<HEpersistentThreshold_<<", HF > "<<HFpersistentThreshold_<<" GeV for 1 full Lumi Block";
230  AbovePersistentThresholdCellsByDepth.depth[1]->setTitle(units.str().c_str());
231  units.str("");
232  units<<"Hot Cells: Depth 3 -- HE > "<<HEpersistentThreshold_<<" GeV for 1 full Lumi Block";
233  AbovePersistentThresholdCellsByDepth.depth[2]->setTitle(units.str().c_str());
234  units.str("");
235  units<<"Hot Cells: HO > "<<HOpersistentThreshold_<<" GeV for 1 full Lumi Block";
236  AbovePersistentThresholdCellsByDepth.depth[3]->setTitle(units.str().c_str());
237  units.str("");
238  }
239 
240  if (test_et_) {
241  SetupEtaPhiHists(ib,AbovePersistentETThresholdCellsByDepth,
242  "Hot Cells Persistently Above ET Threshold","");
243  //setMinMaxHists2D(AbovePersistentThresholdCellsByDepth,0.,1.);
244 
245  // set more descriptive titles for plots
246  units.str("");
247  units<<"Hot Cells: Depth 1 -- HB > "<<HBpersistentETThreshold_<<" GeV (ET), HE > "<<HEpersistentETThreshold_<<" GeV (ET), HF > "<<HFpersistentETThreshold_<<" GeV (ET) for 1 full Lumi Block";
248  AbovePersistentETThresholdCellsByDepth.depth[0]->setTitle(units.str().c_str());
249  units.str("");
250  units<<"Hot Cells: Depth 2 -- HB > "<<HBpersistentETThreshold_<<" GeV (ET), HE > "<<HEpersistentETThreshold_<<" GeV (ET), HF > "<<HFpersistentETThreshold_<<" GeV (ET) for 1 full Lumi Block";
251  AbovePersistentETThresholdCellsByDepth.depth[1]->setTitle(units.str().c_str());
252  units.str("");
253  units<<"Hot Cells: Depth 3 -- HE > "<<HEpersistentETThreshold_<<" GeV (ET) for 1 full Lumi Block";
254  AbovePersistentETThresholdCellsByDepth.depth[2]->setTitle(units.str().c_str());
255  units.str("");
256  units<<"Hot Cells: HO > "<<HOpersistentETThreshold_<<" GeV (ET) for 1 full Lumi Block";
257  AbovePersistentETThresholdCellsByDepth.depth[3]->setTitle(units.str().c_str());
258  units.str("");
259  }
260  }
261 
262  ib.setCurrentFolder(subdir_+"hot_neighbortest");
263  me=ib.bookInt("NeighborTestEnabled");
264  me->Fill(0);
265  if (test_neighbor_)
266  me->Fill(1);
267  if (test_neighbor_ || makeDiagnostics_)
268  {
269  SetupEtaPhiHists(ib,AboveNeighborsHotCellsByDepth,"Hot Cells Failing Neighbor Test","");
270  if (makeDiagnostics_)
271  {
272  d_HBenergyVsNeighbor=ib.book1D("NeighborSumOverEnergyHB","HB Neighbor Sum Energy/Cell Energy;sum(neighbors)/E_cell",500,0,10);
273  d_HEenergyVsNeighbor=ib.book1D("NeighborSumOverEnergyHE","HE Neighbor Sum Energy/Cell Energy;sum(neighbors)/E_cell",500,0,10);
274  d_HOenergyVsNeighbor=ib.book1D("NeighborSumOverEnergyHO","HO Neighbor Sum Energy/Cell Energy;sum(neighbors)/E_cell",500,0,10);
275  d_HFenergyVsNeighbor=ib.book1D("NeighborSumOverEnergyHF","HF Neighbor Sum Energy/Cell Energy;sum(neighbors)/E_cell",500,0,10);
276  }
277  } // if (test_neighbor_ || makeDiagnostics_)
278 
279  this->reset();
280 } // void HcalHotCellMonitor::setup(...)
HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
EtaPhiHists AboveEnergyThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:101
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
HcalBaseDQMonitor::ProblemsVsLB_HF
MonitorElement * ProblemsVsLB_HF
Definition: HcalBaseDQMonitor.h:111
DQMStore::IBooker::bookInt
MonitorElement * bookInt(Args &&...args)
Definition: DQMStore.h:103
HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
EtaPhiHists AbovePersistentThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:103
HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
EtaPhiHists AbovePersistentETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:104
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
HcalBaseDQMonitor::SetupEtaPhiHists
void SetupEtaPhiHists(DQMStore::IBooker &ib, EtaPhiHists &hh, std::string Name, std::string Units)
Definition: HcalBaseDQMonitor.h:74
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalHotCellMonitor::d_HFenergyVsNeighbor
MonitorElement * d_HFenergyVsNeighbor
Definition: HcalHotCellMonitor.h:121
HcalBaseDQMonitor::ProblemsVsLB_HBHEHF
MonitorElement * ProblemsVsLB_HBHEHF
Definition: HcalBaseDQMonitor.h:111
HcalBaseDQMonitor::ProblemsVsLB_HB
MonitorElement * ProblemsVsLB_HB
Definition: HcalBaseDQMonitor.h:110
HcalHotCellMonitor::AboveETThresholdCellsByDepth
EtaPhiHists AboveETThresholdCellsByDepth
Definition: HcalHotCellMonitor.h:102
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
HcalHotCellMonitor::d_HBenergyVsNeighbor
MonitorElement * d_HBenergyVsNeighbor
Definition: HcalHotCellMonitor.h:118
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
HcalBaseDQMonitor::ProblemsVsLB_HE
MonitorElement * ProblemsVsLB_HE
Definition: HcalBaseDQMonitor.h:110
HcalHotCellMonitor::d_HOenergyVsNeighbor
MonitorElement * d_HOenergyVsNeighbor
Definition: HcalHotCellMonitor.h:120
HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
EtaPhiHists AboveNeighborsHotCellsByDepth
Definition: HcalHotCellMonitor.h:100
units
TString units(TString variable, Char_t axis)
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1513
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
HcalBaseDQMonitor::ProblemsVsLB_HO
MonitorElement * ProblemsVsLB_HO
Definition: HcalBaseDQMonitor.h:111
DQMStore::IBooker::bookFloat
MonitorElement * bookFloat(Args &&...args)
Definition: DQMStore.h:109
HcalBaseDQMonitor::ProblemsVsLB
MonitorElement * ProblemsVsLB
Definition: HcalBaseDQMonitor.h:109
HcalBaseDQMonitor::setup
virtual void setup(DQMStore::IBooker &)
Definition: HcalBaseDQMonitor.cc:151
HcalHotCellMonitor::d_HEenergyVsNeighbor
MonitorElement * d_HEenergyVsNeighbor
Definition: HcalHotCellMonitor.h:119
void HcalHotCellMonitor::zeroCounters ( void  )
private

Definition at line 1052 of file HcalHotCellMonitor.cc.

References aboveenergy, aboveet, aboveneighbors, abovepersistent, abovepersistentET, hbVsNeighbor, heVsNeighbor, hfVsNeighbor, hoVsNeighbor, i, j, relval_2017::k, and rechit_occupancy_sum.

Referenced by beginLuminosityBlock(), periodicReset(), and reset().

1053 {
1054 
1055  // zero all counters
1056  for (int i=0;i<85;++i)
1057  {
1058  for (int j=0;j<72;++j)
1059  {
1060  for (int k=0;k<4;++k)
1061  {
1062  abovepersistent[i][j][k]=0;
1063  abovepersistentET[i][j][k]=0;
1064  aboveneighbors[i][j][k]=0;
1065  aboveenergy[i][j][k]=0;
1066  aboveet[i][j][k]=0;
1067  rechit_occupancy_sum[i][j][k]=0;
1068  }
1069  }
1070  }
1071 
1072  for (int i=0;i<500;++i)
1073  {
1074  hbVsNeighbor[i]=0;
1075  heVsNeighbor[i]=0;
1076  hoVsNeighbor[i]=0;
1077  hfVsNeighbor[i]=0;
1078  }
1079  return;
1080 
1081 } // void HcalHotCellMonitor::zeroCounters()
HcalHotCellMonitor::abovepersistent
int abovepersistent[85][72][4]
Definition: HcalHotCellMonitor.h:110
i
int i
Definition: DBlmapReader.cc:9
HcalHotCellMonitor::aboveneighbors
int aboveneighbors[85][72][4]
Definition: HcalHotCellMonitor.h:107
relval_2017.k
list k
Definition: relval_2017.py:23
HcalHotCellMonitor::abovepersistentET
int abovepersistentET[85][72][4]
Definition: HcalHotCellMonitor.h:111
HcalHotCellMonitor::rechit_occupancy_sum
int rechit_occupancy_sum[85][72][4]
Definition: HcalHotCellMonitor.h:112
HcalHotCellMonitor::aboveenergy
int aboveenergy[85][72][4]
Definition: HcalHotCellMonitor.h:108
j
int j
Definition: DBlmapReader.cc:9
HcalHotCellMonitor::aboveet
int aboveet[85][72][4]
Definition: HcalHotCellMonitor.h:109

Member Data Documentation

int HcalHotCellMonitor::aboveenergy[85][72][4]
private

Definition at line 108 of file HcalHotCellMonitor.h.

Referenced by fillNevents_energy(), processEvent_rechitenergy(), and zeroCounters().

EtaPhiHists HcalHotCellMonitor::AboveEnergyThresholdCellsByDepth
private

Definition at line 101 of file HcalHotCellMonitor.h.

Referenced by fillNevents_energy(), fillNevents_problemCells(), periodicReset(), processEvent_rechitenergy(), reset(), and setup().

int HcalHotCellMonitor::aboveet[85][72][4]
private

Definition at line 109 of file HcalHotCellMonitor.h.

Referenced by fillNevents_energy(), processEvent_rechitenergy(), and zeroCounters().

EtaPhiHists HcalHotCellMonitor::AboveETThresholdCellsByDepth
private

Definition at line 102 of file HcalHotCellMonitor.h.

Referenced by fillNevents_energy(), fillNevents_problemCells(), periodicReset(), processEvent_rechitenergy(), reset(), and setup().

int HcalHotCellMonitor::aboveneighbors[85][72][4]
private

Definition at line 107 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), and zeroCounters().

EtaPhiHists HcalHotCellMonitor::AboveNeighborsHotCellsByDepth
private

Definition at line 100 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), fillNevents_problemCells(), periodicReset(), processEvent_rechitenergy(), reset(), and setup().

int HcalHotCellMonitor::abovepersistent[85][72][4]
private

Definition at line 110 of file HcalHotCellMonitor.h.

Referenced by fillNevents_persistentenergy(), processEvent_rechitenergy(), and zeroCounters().

int HcalHotCellMonitor::abovepersistentET[85][72][4]
private

Definition at line 111 of file HcalHotCellMonitor.h.

Referenced by fillNevents_persistentenergy(), processEvent_rechitenergy(), and zeroCounters().

EtaPhiHists HcalHotCellMonitor::AbovePersistentETThresholdCellsByDepth
private

Definition at line 104 of file HcalHotCellMonitor.h.

Referenced by fillNevents_persistentenergy(), fillNevents_problemCells(), periodicReset(), reset(), and setup().

EtaPhiHists HcalHotCellMonitor::AbovePersistentThresholdCellsByDepth
private

Definition at line 103 of file HcalHotCellMonitor.h.

Referenced by fillNevents_persistentenergy(), fillNevents_problemCells(), periodicReset(), processEvent_rechitenergy(), reset(), and setup().

MonitorElement* HcalHotCellMonitor::d_HBenergyVsNeighbor
private

Definition at line 118 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), reset(), and setup().

MonitorElement* HcalHotCellMonitor::d_HEenergyVsNeighbor
private

Definition at line 119 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), reset(), and setup().

MonitorElement* HcalHotCellMonitor::d_HFenergyVsNeighbor
private

Definition at line 121 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), reset(), and setup().

MonitorElement* HcalHotCellMonitor::d_HOenergyVsNeighbor
private

Definition at line 120 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), reset(), and setup().

double HcalHotCellMonitor::energyThreshold_
private

Definition at line 89 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor().

double HcalHotCellMonitor::ETThreshold_
private

Definition at line 90 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor().

bool HcalHotCellMonitor::excludeHORing2_
private

Definition at line 114 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

double HcalHotCellMonitor::HBenergyThreshold_
private

Definition at line 89 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HBETThreshold_
private

Definition at line 90 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

hotNeighborParams HcalHotCellMonitor::HBHENeighborParams_
private

Definition at line 128 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

edm::InputTag HcalHotCellMonitor::hbheRechitLabel_
private

Definition at line 130 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

double HcalHotCellMonitor::HBnsigma_
private

Definition at line 99 of file HcalHotCellMonitor.h.

double HcalHotCellMonitor::HBpersistentETThreshold_
private

Definition at line 92 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HBpersistentThreshold_
private

Definition at line 91 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

int HcalHotCellMonitor::hbVsNeighbor[500]
private

Definition at line 123 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), and zeroCounters().

double HcalHotCellMonitor::HEenergyThreshold_
private

Definition at line 89 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HEETThreshold_
private

Definition at line 90 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HEnsigma_
private

Definition at line 99 of file HcalHotCellMonitor.h.

double HcalHotCellMonitor::HEpersistentETThreshold_
private

Definition at line 92 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HEpersistentThreshold_
private

Definition at line 91 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

int HcalHotCellMonitor::heVsNeighbor[500]
private

Definition at line 124 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), and zeroCounters().

double HcalHotCellMonitor::HFenergyThreshold_
private

Definition at line 89 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HFETThreshold_
private

Definition at line 90 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HFfarfwdScale_
private

Definition at line 94 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

hotNeighborParams HcalHotCellMonitor::HFNeighborParams_
private

Definition at line 128 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

double HcalHotCellMonitor::HFnsigma_
private

Definition at line 99 of file HcalHotCellMonitor.h.

double HcalHotCellMonitor::HFpersistentETThreshold_
private

Definition at line 92 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HFpersistentThreshold_
private

Definition at line 91 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

edm::InputTag HcalHotCellMonitor::hfRechitLabel_
private

Definition at line 130 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

int HcalHotCellMonitor::hfVsNeighbor[500]
private

Definition at line 126 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), and zeroCounters().

double HcalHotCellMonitor::HOenergyThreshold_
private

Definition at line 89 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HOETThreshold_
private

Definition at line 90 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

hotNeighborParams HcalHotCellMonitor::HONeighborParams_
private

Definition at line 128 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

double HcalHotCellMonitor::HOnsigma_
private

Definition at line 99 of file HcalHotCellMonitor.h.

double HcalHotCellMonitor::HOpersistentETThreshold_
private

Definition at line 92 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

double HcalHotCellMonitor::HOpersistentThreshold_
private

Definition at line 91 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), processEvent_rechitenergy(), and setup().

edm::InputTag HcalHotCellMonitor::hoRechitLabel_
private

Definition at line 130 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

int HcalHotCellMonitor::hoVsNeighbor[500]
private

Definition at line 125 of file HcalHotCellMonitor.h.

Referenced by fillNevents_neighbor(), processHit_rechitNeighbors(), and zeroCounters().

double HcalHotCellMonitor::minErrorFlag_
private

Definition at line 96 of file HcalHotCellMonitor.h.

Referenced by fillNevents_problemCells(), HcalHotCellMonitor(), and setup().

int HcalHotCellMonitor::minEvents_
private

Definition at line 81 of file HcalHotCellMonitor.h.

Referenced by fillNevents_persistentenergy(), HcalHotCellMonitor(), and setup().

double HcalHotCellMonitor::nsigma_
private

Definition at line 98 of file HcalHotCellMonitor.h.

double HcalHotCellMonitor::persistentETThreshold_
private

Definition at line 92 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor().

double HcalHotCellMonitor::persistentThreshold_
private

Definition at line 91 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor().

int HcalHotCellMonitor::rechit_occupancy_sum[85][72][4]
private

Definition at line 112 of file HcalHotCellMonitor.h.

Referenced by zeroCounters().

bool HcalHotCellMonitor::setupDone_
private

Definition at line 87 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and setup().

double HcalHotCellMonitor::SiPMscale_
private

Definition at line 106 of file HcalHotCellMonitor.h.

Referenced by HcalHotCellMonitor(), and processEvent_rechitenergy().

bool HcalHotCellMonitor::test_energy_
private

Definition at line 84 of file HcalHotCellMonitor.h.

Referenced by endLuminosityBlock(), fillNevents_energy(), fillNevents_persistentenergy(), fillNevents_problemCells(), HcalHotCellMonitor(), periodicReset(), processEvent(), processEvent_rechitenergy(), reset(), and setup().

bool HcalHotCellMonitor::test_et_
private

Definition at line 85 of file HcalHotCellMonitor.h.

Referenced by endLuminosityBlock(), fillNevents_energy(), fillNevents_persistentenergy(), fillNevents_problemCells(), HcalHotCellMonitor(), periodicReset(), processEvent(), processEvent_rechitenergy(), reset(), and setup().

bool HcalHotCellMonitor::test_neighbor_
private

Definition at line 83 of file HcalHotCellMonitor.h.

Referenced by endLuminosityBlock(), fillNevents_problemCells(), HcalHotCellMonitor(), periodicReset(), processEvent_rechitenergy(), reset(), and setup().

bool HcalHotCellMonitor::test_persistent_
private

Definition at line 86 of file HcalHotCellMonitor.h.

Referenced by endLuminosityBlock(), fillNevents_problemCells(), HcalHotCellMonitor(), periodicReset(), processEvent(), reset(), and setup().

edm::EDGetTokenT<HBHERecHitCollection> HcalHotCellMonitor::tok_hbhe_
private

Definition at line 131 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

edm::EDGetTokenT<HFRecHitCollection> HcalHotCellMonitor::tok_hf_
private

Definition at line 133 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

edm::EDGetTokenT<HORecHitCollection> HcalHotCellMonitor::tok_ho_
private

Definition at line 132 of file HcalHotCellMonitor.h.

Referenced by analyze(), and HcalHotCellMonitor().

Generated for CMSSW Reference Manual by doxygen 1.8.5