CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalRecHitMonitor.cc
Go to the documentation of this file.
1 #include "DQM/HcalMonitorTasks/interface/HcalRecHitMonitor.h"
2 #include "FWCore/Framework/interface/LuminosityBlock.h"
3 #include "FWCore/Common/interface/TriggerNames.h"
4 #include "CondFormats/HcalObjects/interface/HcalChannelStatus.h"
5 #include "CondFormats/HcalObjects/interface/HcalChannelQuality.h"
6 
7 #include "CalibCalorimetry/HcalAlgos/interface/HcalLogicalMapGenerator.h"
8 #include "CondFormats/HcalObjects/interface/HcalLogicalMap.h"
9 #include "RecoLocalCalo/HcalRecAlgos/interface/HcalCaloFlagLabels.h"
10 #include "DataFormats/Common/interface/TriggerResults.h"
11 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
12 #include "Geometry/HcalTowerAlgo/src/HcalHardcodeGeometryData.h" // for eta bounds
13 #include "FWCore/MessageLogger/interface/MessageLogger.h"
14 #include <cmath>
15 
16 HcalRecHitMonitor::HcalRecHitMonitor(const edm::ParameterSet& ps)
17 {
18  // Common Base Class parameters
19  Online_ = ps.getUntrackedParameter<bool>("online",false);
20  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
21  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
22  debug_ = ps.getUntrackedParameter<int>("debug",0);
23  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
24  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
25  prefixME_.append("/");
26  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","RecHitMonitor_Hcal/");
27  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
28  subdir_.append("/");
29  subdir_=prefixME_+subdir_;
30  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
31  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",false);
32  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
33  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
34 
35  hbheRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hbheRechitLabel");
36  hoRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hoRechitLabel");
37  hfRechitLabel_ = ps.getUntrackedParameter<edm::InputTag>("hfRechitLabel");
38  l1gtLabel_ = ps.getUntrackedParameter<edm::InputTag>("L1GTLabel"); // should be l1GtUnpack
39 
40  hltresultsLabel_ = ps.getUntrackedParameter<edm::InputTag>("HLTResultsLabel");
41  HcalHLTBits_ = ps.getUntrackedParameter<std::vector<std::string> >("HcalHLTBits");
42  MinBiasHLTBits_ = ps.getUntrackedParameter<std::vector<std::string> >("MinBiasHLTBits");
43 
44  // energy/ET threshold plots also require that at least one MinBias trigger bit fires
45  energyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold",2);
46  HBenergyThreshold_ = ps.getUntrackedParameter<double>("HB_energyThreshold",energyThreshold_);
47  HEenergyThreshold_ = ps.getUntrackedParameter<double>("HE_energyThreshold",energyThreshold_);
48  HOenergyThreshold_ = ps.getUntrackedParameter<double>("HO_energyThreshold",energyThreshold_);
49  HFenergyThreshold_ = ps.getUntrackedParameter<double>("HF_energyThreshold",energyThreshold_);
50 
51  ETThreshold_ = ps.getUntrackedParameter<double>("ETThreshold",0);
52  HBETThreshold_ = ps.getUntrackedParameter<double>("HB_ETThreshold",ETThreshold_);
53  HEETThreshold_ = ps.getUntrackedParameter<double>("HE_ETThreshold",ETThreshold_);
54  HOETThreshold_ = ps.getUntrackedParameter<double>("HO_ETThreshold",ETThreshold_);
55  HFETThreshold_ = ps.getUntrackedParameter<double>("HF_ETThreshold",ETThreshold_);
56 
57  timediffThresh_ = ps.getUntrackedParameter<double>("collisiontimediffThresh",10.);
58 
59 } //constructor
60 
61 HcalRecHitMonitor::~HcalRecHitMonitor()
62 {
63 } //destructor
64 
65 
66 /* ------------------------------------ */
67 
68 
69 void HcalRecHitMonitor::setup()
70 {
71  HcalBaseDQMonitor::setup();
72 
73 
74  if (debug_>0)
75  std::cout <<"<HcalRecHitMonitor::setup> Setting up histograms"<<std::endl;
76 
77  // Can we include this just in the setup, or do we need to get a new logical map with every run?
78  HcalLogicalMapGenerator gen;
79  logicalMap=new HcalLogicalMap(gen.createMap());
80 
81  // RecHit Monitor - specific cfg variables
82 
83  if (debug_>1)
84  std::cout <<"<HcalRecHitMonitor::setup> Creating Histograms"<<std::endl;
85 
86  dbe_->setCurrentFolder(subdir_);
87  h_TriggeredEvents=dbe_->book1D("EventTriggers","EventTriggers",3,-0.5,2.5);
88  h_TriggeredEvents->setBinLabel(1,"AllEvents");
89  h_TriggeredEvents->setBinLabel(2,"HLT_Minbias");
90  h_TriggeredEvents->setBinLabel(3,"HLT_Hcal");
91 
92  dbe_->setCurrentFolder(subdir_+"rechit_parameters");
93  MonitorElement* THR;
94  dbe_->setCurrentFolder(subdir_+"rechit_parameters/thresholds");
95  THR=dbe_->bookFloat("HB_Rechit_Energy_Threshold");
96  THR->Fill(HBenergyThreshold_);
97  THR=dbe_->bookFloat("HE_Rechit_Energy_Threshold");
98  THR->Fill(HEenergyThreshold_);
99  THR=dbe_->bookFloat("HO_Rechit_Energy_Threshold");
100  THR->Fill(HOenergyThreshold_);
101  THR=dbe_->bookFloat("HF_Rechit_Energy_Threshold");
102  THR->Fill(HFenergyThreshold_);
103  THR=dbe_->bookFloat("HB_Rechit_ET_Threshold");
104  THR->Fill(HBETThreshold_);
105  THR=dbe_->bookFloat("HE_Rechit_ET_Threshold");
106  THR->Fill(HEETThreshold_);
107  THR=dbe_->bookFloat("HO_Rechit_ET_Threshold");
108  THR->Fill(HOETThreshold_);
109  THR=dbe_->bookFloat("HF_Rechit_ET_Threshold");
110  THR->Fill(HFETThreshold_);
111  THR=dbe_->bookFloat("Maximum_HFM_HFP_time_difference_for_luminosityplots");
112  THR->Fill(timediffThresh_);
113 
114 
115  // Set up histograms that are filled by all rechits
116  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits");
117  SetupEtaPhiHists(OccupancyByDepth,"RecHit Occupancy","");
118  h_rechitieta = dbe_->book1D("HcalRecHitIeta",
119  "Hcal RecHit ieta",
120  83,-41.5,41.5);
121  h_rechitiphi = dbe_->book1D("HcalRecHitIphi",
122  "Hcal RecHit iphi",
123  72,0.5,72.5);
124 
125  h_rechitieta_05 = dbe_->book1D("HcalRecHitIeta05",
126  "Hcal RecHit ieta E>0.5 GeV",
127  83,-41.5,41.5);
128  h_rechitiphi_05 = dbe_->book1D("HcalRecHitIphi05",
129  "Hcal RecHit iphi E>0.5 GeV",
130  72,0.5,72.5);
131  h_rechitieta_10 = dbe_->book1D("HcalRecHitIeta10",
132  "Hcal RecHit ieta E>1.0 GeV",
133  83,-41.5,41.5);
134  h_rechitiphi_10 = dbe_->book1D("HcalRecHitIphi10",
135  "Hcal RecHit iphi E>1.0 GeV",
136  72,0.5,72.5);
137  h_rechitieta_25 = dbe_->book1D("HcalRecHitIeta25",
138  "Hcal RecHit ieta E>2.5 GeV",
139  83,-41.5,41.5);
140  h_rechitiphi_25 = dbe_->book1D("HcalRecHitIphi25",
141  "Hcal RecHit iphi E>2.5 GeV",
142  72,0.5,72.5);
143  h_rechitieta_100 = dbe_->book1D("HcalRecHitIeta100",
144  "Hcal RecHit ieta E>10.0 GeV",
145  83,-41.5,41.5);
146  h_rechitiphi_100 = dbe_->book1D("HcalRecHitIphi100",
147  "Hcal RecHit iphi E>10.0 GeV",
148  72,0.5,72.5);
149 
150 
151 
152  h_LumiPlot_LS_allevents = dbe_->book1D("AllEventsPerLS",
153  "LS # of all events",
154  NLumiBlocks_,0.5,NLumiBlocks_+0.5);
155  h_LumiPlot_BX_allevents = dbe_->book1D("BX_allevents",
156  "BX # of all events",
157  3600,0,3600);
158  h_LumiPlot_MinTime_vs_MinHT = dbe_->book2D("MinTime_vs_MinSumET",
159  "Energy-Weighted Time vs Min (HF+,HF-) Scalar Sum ET;min Sum ET(GeV);time(ns)",
160  100,0,10,80,-40,40);
161 
162  h_LumiPlot_timeHT_HFM = dbe_->book2D("HFM_Time_vs_SumET",
163  "Energy-Weighted Time vs HFMinus Scalar Sum ET;Sum ET(GeV);time(ns)",
164  100,0,10,80,-40,40);
165 
166  h_LumiPlot_timeHT_HFP = dbe_->book2D("HFP_Time_vs_SumET",
167  "Energy-Weighted Time vs HFPlus Scalar Sum ET;Sum ET(GeV);time(ns)",
168  100,0,10,80,-40,40);
169 
170 
171  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits/sumplots");
172  SetupEtaPhiHists(SumEnergyByDepth,"RecHit Summed Energy","GeV");
173  SetupEtaPhiHists(SqrtSumEnergy2ByDepth,"RecHit Sqrt Summed Energy2","GeV");
174  SetupEtaPhiHists(SumTimeByDepth,"RecHit Summed Time","nS");
175 
176  // Histograms for events that passed MinBias triggers
177  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias");
178 
179  h_HBP_weightedTime = dbe_->book1D("WeightedTime_HBP","Weighted Time for HBP",
180  300,-150,150);
181  h_HBM_weightedTime = dbe_->book1D("WeightedTime_HBM","Weighted Time for HBM",
182  300,-150,150);
183  h_HEP_weightedTime = dbe_->book1D("WeightedTime_HEP","Weighted Time for HEP",
184  300,-150,150);
185  h_HEM_weightedTime = dbe_->book1D("WeightedTime_HEM","Weighted Time for HEM",
186  300,-150,150);
187  h_HFP_weightedTime = dbe_->book1D("WeightedTime_HFP","Weighted Time for HFP",
188  300,-150,150);
189  h_HFM_weightedTime = dbe_->book1D("WeightedTime_HFM","Weighted Time for HFM",
190  300,-150,150);
191 
192  h_HFtimedifference = dbe_->book1D("HFweightedtimeDifference",
193  "Energy-Weighted time difference between HF+ and HF- passing MinBias (no HT cut)",
194  251,-250.5,250.5);
195  h_HEtimedifference = dbe_->book1D("HEweightedtimeDifference",
196  "Energy-Weighted time difference between HE+ and HE- passing MinBias (no HT cut)",
197  251,-250.5,250.5);
198 
199  HFP_HFM_Energy = dbe_->book2D("HFP_HFM_Energy",
200  "HFP VS HFM Energy; Total Energy in HFMinus (TeV); Total Energy in HFPlus (TeV)",
201  100,0,100, 100,0,100);
202 
203  // Would these work better as 2D plots?
204  h_HFenergydifference = dbe_->book1D("HFenergyDifference",
205  "Sum(E_HFPlus - E_HFMinus)/Sum(E_HFPlus + E_HFMinus)",
206  200,-1,1);
207  h_HEenergydifference = dbe_->book1D("HEenergyDifference",
208  "Sum(E_HEPlus - E_HEMinus)/Sum(E_HEPlus + E_HEMinus)",
209  200,-1,1);
210 
211  h_LumiPlot_LS_MinBiasEvents=dbe_->book1D("MinBiasEventsPerLS",
212  "Number of MinBias Events vs LS (HT cut and HFM-HFP time cut)",
213  NLumiBlocks_/10,0.5,NLumiBlocks_+0.5);
214  h_LumiPlot_LS_MinBiasEvents_notimecut=dbe_->book1D("MinBiasEventsPerLS_notimecut",
215  "Number of Events with MinBias vs LS (HFM,HFP HT>1,no time cut)",
216  NLumiBlocks_/10,0.5,NLumiBlocks_+0.5);
217 
218  h_LumiPlot_SumHT_HFPlus_vs_HFMinus = dbe_->book2D("SumHT_plus_minus",
219  "HF+ Sum HT vs HF- Sum HT",60,0,30,60,0,30);
220  h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus = dbe_->book2D("SumEnergy_plus_minus",
221  "HF+ Sum Energy vs HF- Sum Energy",
222  60,0,150,60,0,150);
223  h_LumiPlot_timeHFPlus_vs_timeHFMinus = dbe_->book2D("timeHFplus_vs_timeHFminus",
224  "Energy-weighted time average of HF+ vs HF-",
225  60,-60,60,60,-60,60);
226  h_LumiPlot_BX_MinBiasEvents = dbe_->book1D("BX_MinBias_Events_TimeCut",
227  "BX # of MinBias events (HFM & HFP HT>1 & HFM-HFP time cut)",
228  3600,0,3600);
229  h_LumiPlot_BX_MinBiasEvents_notimecut = dbe_->book1D("BX_MinBias_Events_notimecut",
230  "BX # of MinBias events (HFM,HFP HT>1, no time cut)",
231  3600,0,3600);
232  // threshold plots must pass MinBias Trigger
233  SetupEtaPhiHists(OccupancyThreshByDepth,"Above Threshold RecHit Occupancy","");
234  h_rechitieta_thresh = dbe_->book1D("HcalRecHitIeta_thresh",
235  "Hcal RecHit ieta above energy and ET threshold",
236  83,-41.5,41.5);
237  h_rechitiphi_thresh = dbe_->book1D("HcalRecHitIphi_thresh",
238  "Hcal RecHit iphi above energy and ET threshold",
239  72,0.5,72.5);
240 
241  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/sumplots");
242  SetupEtaPhiHists(SumEnergyThreshByDepth,"Above Threshold RecHit Summed Energy","GeV");
243  SetupEtaPhiHists(SumTimeThreshByDepth,"Above Threshold RecHit Summed Time","nS");
244  SetupEtaPhiHists(SqrtSumEnergy2ThreshByDepth,"Above Threshold RecHit Sqrt Summed Energy2","GeV");
245 
246  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/rechit_1D_plots");
247  h_HBThreshTime=dbe_->book1D("HB_time_thresh",
248  "HB RecHit Time Above Threshold",
249  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
250  h_HBThreshOccupancy=dbe_->book1D("HB_occupancy_thresh",
251  "HB RecHit Occupancy Above Threshold",260,-0.5,2599.5);
252  h_HEThreshTime=dbe_->book1D("HE_time_thresh",
253  "HE RecHit Time Above Threshold",
254  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
255  h_HEThreshOccupancy=dbe_->book1D("HE_occupancy_thresh",
256  "HE RecHit Occupancy Above Threshold",260,-0.5,2599.5);
257  h_HOThreshTime=dbe_->book1D("HO_time_thresh",
258  "HO RecHit Time Above Threshold",
259  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
260  h_HOThreshOccupancy=dbe_->book1D("HO_occupancy_thresh",
261  "HO RecHit Occupancy Above Threshold",217,-0.5,2169.5);
262  h_HFThreshTime=dbe_->book1D("HF_time_thresh",
263  "HF RecHit Time Above Threshold",
264  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
265  h_HFThreshOccupancy=dbe_->book1D("HF_occupancy_thresh",
266  "HF RecHit Occupancy Above Threshold",
267  173,-0.5,1729.5);
268 
269  // Histograms for events that did passed Hcal-specfied HLT triggers
270  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers");
271 
272  h_LumiPlot_BX_HcalHLTEvents = dbe_->book1D("BX_HcalHLT_Events_TimeCut",
273  "BX # of HcalHLT events (HFM & HFP HT>1 & HFM-HFP time cut)",
274  3600,0,3600);
275  h_LumiPlot_BX_HcalHLTEvents_notimecut = dbe_->book1D("BX_HcalHLT_Events_notimecut",
276  "BX # of HcalHLT events (HFM,HFP HT>1, no time cut)",
277  3600,0,3600);
278  h_LumiPlot_LS_HcalHLTEvents=dbe_->book1D("HcalHLTEventsPerLS",
279  "Number of HcalHLT Events vs LS (HT cut and HFM-HFP time cut)",
280  NLumiBlocks_/10,0.5,NLumiBlocks_+0.5);
281  h_LumiPlot_LS_HcalHLTEvents_notimecut=dbe_->book1D("HcalHLTEventsPerLS_notimecut",
282  "Number of Events with HcalHLT vs LS (HFM,HFP HT>1,no time cut)",
283  NLumiBlocks_/10,0.5,NLumiBlocks_+0.5);
284 
285 
286  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers/");
287  h_HF_HcalHLT_weightedtimedifference = dbe_->book1D("HF_HcalHLT_weightedtimeDifference",
288  "Energy-Weighted time difference between HF+ and HF- Hcal HLT",
289  251,-250.5,250.5);
290  h_HE_HcalHLT_weightedtimedifference = dbe_->book1D("HE_HcalHLT_weightedtimeDifference",
291  "Energy-Weighted time difference between HE+ and HE- Hcal HLT",
292  251,-250.5,250.5);
293  h_HF_HcalHLT_energydifference = dbe_->book1D("HF_HcalHLT_energyDifference",
294  "Sum(E_HFPlus - E_HFMinus)/Sum(E_HFPlus + E_HFMinus)",
295  200,-1,1);
296  h_HE_HcalHLT_energydifference = dbe_->book1D("HE_HcalHLT_energyDifference",
297  "Sum(E_HEPlus - E_HEMinus)/Sum(E_HEPlus + E_HEMinus)",
298  200,-1,1);
299 
300  // Do we want separate directories for Minbias, other flags at some point?
301  dbe_->setCurrentFolder(subdir_+"AnomalousCellFlags");// HB Flag Histograms
302 
303 
304  h_HFLongShort_vs_LS=dbe_->book1D("HFLongShort_vs_LS",
305  "HFLongShort Flags vs Lumi Section",
306  NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
307  h_HFDigiTime_vs_LS=dbe_->book1D("HFDigiTime_vs_LS",
308  "HFDigiTime Flags vs Lumi Section",
309  NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
310  h_HBHEHPDMult_vs_LS=dbe_->book1D("HBHEHPDMult_vs_LS",
311  "HBHEHPDMult Flags vs Lumi Section",
312  NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
313  h_HBHEPulseShape_vs_LS=dbe_->book1D("HBHEPulseShape_vs_LS",
314  "HBHEPulseShape Flags vs Lumi Section",
315  NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
316 
317  h_HF_FlagCorr=dbe_->book2D("HF_FlagCorrelation",
318  "HF LongShort vs. DigiTime flags; DigiTime; LongShort",
319  2,-0.5,1.5,2,-0.5,1.5);
320  h_HF_FlagCorr->setBinLabel(1,"OFF",1);
321  h_HF_FlagCorr->setBinLabel(2,"ON",1);
322  h_HF_FlagCorr->setBinLabel(1,"OFF",2);
323  h_HF_FlagCorr->setBinLabel(2,"ON",2);
324 
325  h_HBHE_FlagCorr=dbe_->book2D("HBHE_FlagCorrelation",
326  "HBHE HpdHitMultiplicity vs. PulseShape flags; PulseShape; HpdHitMultiplicity",
327  2,-0.5,1.5,2,-0.5,1.5);
328  h_HBHE_FlagCorr->setBinLabel(1,"OFF",1);
329  h_HBHE_FlagCorr->setBinLabel(2,"ON",1);
330  h_HBHE_FlagCorr->setBinLabel(1,"OFF",2);
331  h_HBHE_FlagCorr->setBinLabel(2,"ON",2);
332 
333  h_FlagMap_HPDMULT=dbe_->book2D("FlagMap_HPDMULT",
334  "RBX Map of HBHEHpdHitMultiplicity Flags;RBX;RM",
335  72,-0.5,71.5,4,0.5,4.5);
336  h_FlagMap_PULSESHAPE=dbe_->book2D("FlagMap_PULSESHAPE",
337  "RBX Map of HBHEPulseShape Flags;RBX;RM",
338  72,-0.5,71.5,4,0.5,4.5);
339  h_FlagMap_DIGITIME=dbe_->book2D("FlagMap_DIGITIME",
340  "RBX Map of HFDigiTime Flags;RBX;RM",
341  24,131.5,155.5,4,0.5,4.5);
342  h_FlagMap_LONGSHORT=dbe_->book2D("FlagMap_LONGSHORT",
343  "RBX Map of HFLongShort Flags;RBX;RM",
344  24,131.5,155.5,4,0.5,4.5);
345 
346  h_FlagMap_TIMEADD=dbe_->book2D("FlagMap_TIMEADD",
347  "RBX Map of Timing Added Flags;RBX;RM",
348  156,-0.5,155.5,4,0.5,4.5);
349  h_FlagMap_TIMESUBTRACT=dbe_->book2D("FlagMap_TIMESUBTRACT",
350  "RBX Map of Timing Subtracted Flags;RBX;RM",
351  156,-0.5,155.5,4,0.5,4.5);
352  h_FlagMap_TIMEERROR=dbe_->book2D("FlagMap_TIMEERROR",
353  "RBX Map of Timing Error Flags;RBX;RM",
354  156,-0.5,155.5,4,0.5,4.5);
355 
356  h_HBflagcounter=dbe_->book1D("HBflags","HB flags",32,-0.5,31.5);
357  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEHpdHitMultiplicity, "HpdHitMult",1);
358  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEPulseShape, "PulseShape",1);
359  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_R1R2, "HSCP R1R2",1);
360  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_FracLeader, "HSCP FracLeader",1);
361  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_OuterEnergy, "HSCP OuterEnergy",1);
362  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_ExpFit, "HSCP ExpFit",1);
363  // 2-bit timing counter
364  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust1",1);
365  h_HBflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust2",1);
366  //3-bit timing shape cut
367  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape1",1);
368  h_HBflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape2",1);
369  h_HBflagcounter->setBinLabel(3+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape3",1);
370 
371  // common flags
372  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
373  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
374  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
375  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
376  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);
377 
378  // HE Flag Histograms
379  h_HEflagcounter=dbe_->book1D("HEflags","HE flags",32,-0.5,31.5);
380  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEHpdHitMultiplicity, "HpdHitMult",1);
381  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEPulseShape, "PulseShape",1);
382  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_R1R2, "HSCP R1R2",1);
383  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_FracLeader, "HSCP FracLeader",1);
384  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_OuterEnergy, "HSCP OuterEnergy",1);
385  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_ExpFit, "HSCP ExpFit",1);
386  // 2-bit timing counter
387  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust1",1);
388  h_HEflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust2",1);
389  //3-bit timing shape cut
390  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape1",1);
391  h_HEflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape2",1);
392  h_HEflagcounter->setBinLabel(3+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape3",1);
393  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
394  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
395  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
396  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
397  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);
398 
399  // HO Flag Histograms
400  h_HOflagcounter=dbe_->book1D("HOflags","HO flags",32,-0.5,31.5);
401  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
402  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
403  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
404  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
405  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);
406 
407  // HF Flag Histograms
408  h_HFflagcounter=dbe_->book1D("HFflags","HF flags",32,-0.5,31.5);
409  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFLongShort, "LongShort",1);
410  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFDigiTime, "DigiTime",1);
411  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFTimingTrustBits,"TimingTrust1",1);
412  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
413  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
414  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
415  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
416  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);
417 
418  h_HBflagcounter->getTH1F()->LabelsOption("v");
419  h_HEflagcounter->getTH1F()->LabelsOption("v");
420  h_HOflagcounter->getTH1F()->LabelsOption("v");
421  h_HFflagcounter->getTH1F()->LabelsOption("v");
422 
423 
424  // Diagnostic plots are currently filled for all rechits (no trigger/threshold requirement)
425  // hb
426  dbe_->setCurrentFolder(subdir_+"diagnostics/hb");
427 
428  h_HBTimeVsEnergy=dbe_->book2D("HBTimeVsEnergy","HB Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);
429 
430  h_HBsizeVsLS=dbe_->bookProfile("HBRecHitsVsLB","HB RecHits vs Luminosity Block",
431  NLumiBlocks_,0.5,NLumiBlocks_+0.5,
432  100,0,10000);
433 
434  h_HBTime=dbe_->book1D("HB_time","HB RecHit Time",
435  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
436  h_HBOccupancy=dbe_->book1D("HB_occupancy",
437  "HB RecHit Occupancy",260,-0.5,2599.5);
438 
439  //he
440  dbe_->setCurrentFolder(subdir_+"diagnostics/he");
441 
442  h_HETimeVsEnergy=dbe_->book2D("HETimeVsEnergy","HE Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);
443 
444  h_HEsizeVsLS=dbe_->bookProfile("HERecHitsVsLB","HE RecHits vs Luminosity Block",
445  NLumiBlocks_,0.5,NLumiBlocks_+0.5,
446  100,0,10000);
447 
448  h_HETime=dbe_->book1D("HE_time","HE RecHit Time",
449  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
450  h_HEOccupancy=dbe_->book1D("HE_occupancy","HE RecHit Occupancy",260,-0.5,2599.5);
451 
452  // ho
453  dbe_->setCurrentFolder(subdir_+"diagnostics/ho");
454 
455  h_HOTimeVsEnergy=dbe_->book2D("HOTimeVsEnergy","HO Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);
456 
457  h_HOsizeVsLS=dbe_->bookProfile("HORecHitsVsLB","HO RecHits vs Luminosity Block",
458  NLumiBlocks_,0.5,NLumiBlocks_+0.5,
459  100,0,10000);
460  h_HOTime=dbe_->book1D("HO_time",
461  "HO RecHit Time",
462  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
463  h_HOOccupancy=dbe_->book1D("HO_occupancy",
464  "HO RecHit Occupancy",217,-0.5,2169.5);
465 
466  // hf
467  dbe_->setCurrentFolder(subdir_+"diagnostics/hf");
468 
469  h_HFTimeVsEnergy=dbe_->book2D("HFTimeVsEnergy","HF Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);
470 
471  h_HFsizeVsLS=dbe_->bookProfile("HFRecHitsVsLB",
472  "HF RecHits vs Luminosity Block",
473  NLumiBlocks_,0.5,NLumiBlocks_+0.5,
474  100, 0,10000);
475  h_HFTime=dbe_->book1D("HF_time","HF RecHit Time",
476  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
477  h_HFOccupancy=dbe_->book1D("HF_occupancy","HF RecHit Occupancy",173,-0.5,1729.5);
478 
479  // clear all counters, reset histograms
480  this->zeroCounters();
481  this->reset();
482 
483  return;
484 } //void HcalRecHitMonitor::setup(...)
485 
486 void HcalRecHitMonitor::beginRun(const edm::Run& run, const edm::EventSetup& c)
487 {
488 
489  if (debug_>0) std::cout <<"HcalRecHitMonitor::beginRun(): task = '"<<subdir_<<"'"<<std::endl;
490  HcalBaseDQMonitor::beginRun(run, c);
491  if (tevt_==0) // create histograms, if they haven't been created already
492  this->setup();
493  // Clear histograms at the start of each run if not merging runs
494  if (mergeRuns_==false)
495  this->reset();
496 
497  if (tevt_!=0) return;
498  // create histograms displaying trigger parameters? Specify names?
499  dbe_->setCurrentFolder(subdir_+"rechit_parameters");
500  MonitorElement* trig;
501  std::string tnames="";
502  if (HcalHLTBits_.size()>0)
503  tnames=HcalHLTBits_[0];
504  for (unsigned int i=1;i<HcalHLTBits_.size();++i)
505  tnames=tnames + " OR " + HcalHLTBits_[i];
506  trig=dbe_->bookString("HcalHLTriggerRequirements",tnames);
507  tnames="";
508  if (MinBiasHLTBits_.size()>0)
509  tnames=MinBiasHLTBits_[0];
510  for (unsigned int i=1;i<MinBiasHLTBits_.size();++i)
511  tnames=tnames + " OR " + MinBiasHLTBits_[i];
512  trig=dbe_->bookString("MinBiasHLTriggerRequirements",tnames);
513  return;
514 
515 } //void HcalRecHitMonitor::beginRun(...)
516 
517 
518 void HcalRecHitMonitor::endRun(const edm::Run& run, const edm::EventSetup& c)
519 {
520  if (debug_>0) std::cout <<"HcalRecHitMonitor::endRun(): task = '"<<subdir_<<"'"<<std::endl;
521 
522  //Any special fill calls needed? Shouldn't be necessary; last endLuminosityBlock should do necessary fills
523 } // void HcalRecHitMonitor::endrun(...)
524 
525 
526 /* --------------------------- */
527 
528 void HcalRecHitMonitor::reset()
529 {
530  std::vector<MonitorElement*> hists = dbe_->getAllContents(subdir_);
531  for (unsigned int i=0;i<hists.size();++i)
532  {
533  if (hists[i]->kind()==MonitorElement::DQM_KIND_TH1F ||
534  hists[i]->kind()==MonitorElement::DQM_KIND_TH2F ||
535  hists[i]->kind()==MonitorElement::DQM_KIND_TPROFILE)
536  hists[i]->Reset();
537  }
538 }
539 
540 void HcalRecHitMonitor::endJob()
541 {
542  if (!enableCleanup_) return;
543  HcalBaseDQMonitor::cleanup();
544  this->cleanup();
545 }
546 
547 
548 /* --------------------------------- */
549 
550 void HcalRecHitMonitor::cleanup()
551 {
552  //Add code to clean out subdirectories
553  if (!enableCleanup_) return;
554  if (dbe_)
555  {
556  dbe_->setCurrentFolder(subdir_); dbe_->removeContents();
557  dbe_->setCurrentFolder(subdir_+"rechit_parameters"); dbe_->removeContents();
558  dbe_->setCurrentFolder(subdir_+"rechit_parameters/thresholds"); dbe_->removeContents();
559  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits"); dbe_->removeContents();
560  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits/sumplots"); dbe_->removeContents();
561  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias"); dbe_->removeContents();
562  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/sumplots"); dbe_->removeContents();
563  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers"); dbe_->removeContents();
564  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers/passedTechTriggers/"); dbe_->removeContents();
565 
566  dbe_->setCurrentFolder(subdir_+"AnomalousCellFlags"); dbe_->removeContents();
567  dbe_->setCurrentFolder(subdir_+"diagnostics/hb"); dbe_->removeContents();
568  dbe_->setCurrentFolder(subdir_+"diagnostics/he"); dbe_->removeContents();
569  dbe_->setCurrentFolder(subdir_+"diagnostics/ho"); dbe_->removeContents();
570  dbe_->setCurrentFolder(subdir_+"diagnostics/hf"); dbe_->removeContents();
571  }
572  return;
573 } // void HcalRecHitMonitor::cleanup()
574 
575 /* -------------------------------- */
576 
577 void HcalRecHitMonitor::analyze(edm::Event const&e, edm::EventSetup const&s)
578 {
579  if (debug_>0) std::cout <<"HcalRecHitMonitor::analyze; debug = "<<debug_<<std::endl;
580 
581  if (!IsAllowedCalibType()) return;
582  if (LumiInOrder(e.luminosityBlock())==false) return;
583 
584  // Get objects
585  edm::Handle<HBHERecHitCollection> hbhe_rechit;
586  edm::Handle<HORecHitCollection> ho_rechit;
587  edm::Handle<HFRecHitCollection> hf_rechit;
588 
589  if (!(e.getByLabel(hbheRechitLabel_,hbhe_rechit)))
590  {
591  edm::LogWarning("HcalHotCellMonitor")<< hbheRechitLabel_<<" hbhe_rechit not available";
592  return;
593  }
594 
595  if (!(e.getByLabel(hfRechitLabel_,hf_rechit)))
596  {
597  edm::LogWarning("HcalHotCellMonitor")<< hfRechitLabel_<<" hf_rechit not available";
598  return;
599  }
600 
601  if (!(e.getByLabel(hoRechitLabel_,ho_rechit)))
602  {
603  edm::LogWarning("HcalHotCellMonitor")<< hoRechitLabel_<<" ho_rechit not available";
604  return;
605  }
606 
607 
608  h_LumiPlot_LS_allevents->Fill(currentLS);
609  h_LumiPlot_BX_allevents->Fill(e.bunchCrossing());
610  processEvent(*hbhe_rechit, *ho_rechit, *hf_rechit, e.bunchCrossing(), e);
611 
612  HcalBaseDQMonitor::analyze(e,s);
613 } // void HcalRecHitMonitor::analyze()
614 
615 
616 void HcalRecHitMonitor::processEvent(const HBHERecHitCollection& hbHits,
617  const HORecHitCollection& hoHits,
618  const HFRecHitCollection& hfHits,
619  int BCN,
620  const edm::Event & iEvent
621  )
622 {
623 
624 
625  if (debug_>1) std::cout <<"<HcalRecHitMonitor::processEvent> Processing event..."<<std::endl;
626 
627 
628  bool passedHcalHLT=false;
629  bool passedMinBiasHLT=false;
630 
631  edm::Handle<edm::TriggerResults> hltRes;
632  if (!(iEvent.getByLabel(hltresultsLabel_,hltRes)))
633  {
634  if (debug_>0) edm::LogWarning("HcalRecHitMonitor")<<" Could not get HLT results with tag "<<hltresultsLabel_<<std::endl;
635  }
636  else
637  {
638  const edm::TriggerNames & triggerNames = iEvent.triggerNames(*hltRes);
639  const unsigned int nTrig(triggerNames.size());
640  for (unsigned int i=0;i<nTrig;++i)
641  {
642  // trigger decision is based on 'OR' of any specified trigger names
643  for (unsigned int k=0;k<HcalHLTBits_.size();++k)
644  {
645  // if (triggerNames.triggerName(i)==HcalHLTBits_[k] && hltRes->accept(i))
646  if (triggerNames.triggerName(i).find(HcalHLTBits_[k])!=std::string::npos && hltRes->accept(i))
647  {
648  passedHcalHLT=true;
649  break;
650  }
651  }
652  // repeat for minbias triggers
653  for (unsigned int k=0;k<MinBiasHLTBits_.size();++k)
654  {
655  // if (triggerNames.triggerName(i)==MinBiasHLTBits_[k] && hltRes->accept(i))
656  if (triggerNames.triggerName(i).find(MinBiasHLTBits_[k])!=std::string::npos && hltRes->accept(i))
657  {
658  passedMinBiasHLT=true;
659  break;
660  }
661  }
662  }
663  } //else
664 
665  if (debug_>2 && passedHcalHLT) std::cout <<"\t<HcalRecHitMonitor::processEvent> Passed Hcal HLT trigger "<<std::endl;
666  if (debug_>2 && passedMinBiasHLT) std::cout <<"\t<HcalRecHitMonitor::processEvent> Passed MinBias HLT trigger "<<std::endl;
667 
668  h_TriggeredEvents->Fill(0); // all events
669  if (passedMinBiasHLT) h_TriggeredEvents->Fill(1); // Minbias;
670  if (passedHcalHLT) h_TriggeredEvents->Fill(2); // hcal HLT
671  processEvent_rechit(hbHits, hoHits, hfHits,passedHcalHLT,passedMinBiasHLT,BCN);
672 
673  return;
674 } // void HcalRecHitMonitor::processEvent(...)
675 
676 
677 /* --------------------------------------- */
678 
679 
680 void HcalRecHitMonitor::processEvent_rechit( const HBHERecHitCollection& hbheHits,
681  const HORecHitCollection& hoHits,
682  const HFRecHitCollection& hfHits,
683  bool passedHcalHLT,
684  bool passedMinBiasHLT,
685  int BCN)
686 {
687  // Gather rechit info
688 
689  //const float area[]={0.111,0.175,0.175,0.175,0.175,0.175,0.174,0.178,0.172,0.175,0.178,0.346,0.604};
690 
691  if (debug_>1) std::cout <<"<HcalRecHitMonitor::processEvent_rechitenergy> Processing rechits..."<<std::endl;
692 
693  // loop over HBHE
694 
695  int hbocc=0;
696  int heocc=0;
697  int hboccthresh=0;
698  int heoccthresh=0;
699 
700  double HtPlus =0, HtMinus=0;
701  double HFePlus=0, HFeMinus=0;
702  double HBePlus=0, HBeMinus=0;
703  double HEePlus=0, HEeMinus=0;
704  double HFtPlus=0, HFtMinus=0;
705  double HBtPlus=0, HBtMinus=0;
706  double HEtPlus=0, HEtMinus=0;
707 
708  int hbpocc=0, hbmocc=0, hepocc=0, hemocc=0, hfpocc=0, hfmocc=0;
709 
710  for (unsigned int i=0;i<4;++i)
711  {
712  OccupancyByDepth.depth[i]->update();
713  OccupancyThreshByDepth.depth[i]->update();
714  SumEnergyByDepth.depth[i]->update();
715  SqrtSumEnergy2ByDepth.depth[i]->update();
716  SumTimeByDepth.depth[i]->update();
717  }
718 
719  h_HBflagcounter->update();
720  h_HEflagcounter->update();
721  h_HFflagcounter->update();
722  h_HOflagcounter->update();
723 
724 
725  for (HBHERecHitCollection::const_iterator HBHEiter=hbheHits.begin(); HBHEiter!=hbheHits.end(); ++HBHEiter)
726  { // loop over all hits
727  float en = HBHEiter->energy();
728  float ti = HBHEiter->time();
729  HcalDetId id(HBHEiter->detid().rawId());
730  int ieta = id.ieta();
731  int iphi = id.iphi();
732  int depth = id.depth();
733 
734  if (en>0.5)
735  {
736  h_rechitieta_05->Fill(ieta);
737  h_rechitiphi_05->Fill(iphi);
738  if (en>1.)
739  {
740  h_rechitieta_10->Fill(ieta);
741  h_rechitiphi_10->Fill(iphi);
742  if (en>2.5)
743  {
744  h_rechitieta_25->Fill(ieta);
745  h_rechitiphi_25->Fill(iphi);
746  if (en>10.)
747  {
748  h_rechitieta_100->Fill(ieta);
749  h_rechitiphi_100->Fill(iphi);
750  }
751  }
752  }
753  }
754 
755 
756 
757  HcalSubdetector subdet = id.subdet();
758  double fEta=fabs(0.5*(theHBHEEtaBounds[abs(ieta)-1]+theHBHEEtaBounds[abs(ieta)]));
759 
760  int calcEta = CalcEtaBin(subdet,ieta,depth);
761  int rbxindex=logicalMap->getHcalFrontEndId(HBHEiter->detid()).rbxIndex();
762  int rm= logicalMap->getHcalFrontEndId(HBHEiter->detid()).rm();
763 
764  // Fill HBHE flag plots
765  h_HBHE_FlagCorr->Fill(HBHEiter->flagField(HcalCaloFlagLabels::HBHEPulseShape),
766  HBHEiter->flagField(HcalCaloFlagLabels::HBHEHpdHitMultiplicity));
767 
768  if (HBHEiter->flagField(HcalCaloFlagLabels::HBHEHpdHitMultiplicity))
769  {
770  h_FlagMap_HPDMULT->Fill(rbxindex,rm);
771  h_HBHEHPDMult_vs_LS->Fill(currentLS,1);
772  }
773  if (HBHEiter->flagField(HcalCaloFlagLabels::HBHEPulseShape))
774  {
775  h_FlagMap_PULSESHAPE->Fill(rbxindex,rm);
776  h_HBHEPulseShape_vs_LS->Fill(currentLS,1);
777  }
778  if (HBHEiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
779  h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
780  else if (HBHEiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
781  h_FlagMap_TIMEADD->Fill(rbxindex,rm);
782  else if (HBHEiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
783  h_FlagMap_TIMEERROR->Fill(rbxindex,rm);
784 
785  if (subdet==HcalBarrel)
786  {
787  if (en>HBenergyThreshold_)
788  h_HBTimeVsEnergy->Fill(en,ti);
789  //Looping over HB searching for flags --- cris
790  for (int f=0;f<32;f++)
791  {
792  // Let's display HSCP just to see if these bits are set
793  /*
794  if (f == HcalCaloFlagLabels::HSCP_R1R2) continue;
795  if (f == HcalCaloFlagLabels::HSCP_FracLeader) continue;
796  if (f == HcalCaloFlagLabels::HSCP_OuterEnergy) continue;
797  if (f == HcalCaloFlagLabels::HSCP_ExpFit) continue;
798  */
799  if (HBHEiter->flagField(f))
800  ++HBflagcounter_[f];
801  }
802  ++occupancy_[calcEta][iphi-1][depth-1];
803  energy_[calcEta][iphi-1][depth-1]+=en;
804  energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
805  time_[calcEta][iphi-1][depth-1]+=ti;
806  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
807  h_HBTime->Fill(ti);
808  else
809  ++HBtime_[int(ti-RECHITMON_TIME_MIN)];
810  ++hbocc;
811 
812  // Threshold plots; require E> threshold and minbias trigger
813  if (
814  en>=HBenergyThreshold_ &&
815  en/cosh(fEta)>=HBETThreshold_
816  )
817  {
818  if (passedMinBiasHLT==true)
819  {
820  ++occupancy_thresh_[calcEta][iphi-1][depth-1];
821  energy_thresh_[calcEta][iphi-1][depth-1]+=en;
822  energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
823  time_thresh_[calcEta][iphi-1][depth-1]+=ti;
824 
825  ++hboccthresh;
826  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
827  h_HBThreshTime->Fill(ti);
828  else
829  ++HBtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
830  }
831 
832  if (ieta>0)
833  {
834  HBePlus+=en;
835  HBtPlus+=ti*en;
836  hbpocc++;
837  }
838  else
839  {
840  HBeMinus+=en;
841  HBtMinus+=ti*en;
842  hbmocc++;
843  }
844  } // if (HB en>thresh, ET>thresh)
845  } // if (id.subdet()==HcalBarrel)
846 
847  else if (subdet==HcalEndcap)
848  {
849  if (en>HEenergyThreshold_)
850  h_HETimeVsEnergy->Fill(en,ti);
851  //Looping over HE searching for flags --- cris
852  for (int f=0;f<32;f++)
853  {
854  if (HBHEiter->flagField(f))
855  ++HEflagcounter_[f];
856  }
857 
858  ++occupancy_[calcEta][iphi-1][depth-1];
859  energy_[calcEta][iphi-1][depth-1]+=en;
860  energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
861  time_[calcEta][iphi-1][depth-1]+=ti;
862 
863  ++heocc;
864  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
865  h_HETime->Fill(ti);
866  else
867  ++HEtime_[int(ti-RECHITMON_TIME_MIN)];
868 
869  // Threshold plots require e>E_thresh, ET>ET_thresh
870  if (en>=HEenergyThreshold_
871  && en/cosh(fEta)>=HEETThreshold_
872  )
873  {
874  // occupancy plots also require passedMinBiasHLT
875  if (passedMinBiasHLT==true)
876  {
877  ++occupancy_thresh_[calcEta][iphi-1][depth-1];
878  energy_thresh_[calcEta][iphi-1][depth-1]+=en;
879  energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
880  time_thresh_[calcEta][iphi-1][depth-1]+=ti;
881  ++heoccthresh;
882  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
883  h_HEThreshTime->Fill(ti);
884  else
885  ++HEtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
886  }
887  // ePlus, tPlus calculated regardless of trigger
888  if (ieta>0)
889  {
890  HEePlus+=en;
891  HEtPlus+=ti*en;
892  hepocc++;
893  }
894  else
895  {
896  HEeMinus+=en;
897  HEtMinus+=ti*en;
898  hemocc++;
899  }
900  } // if (en>=HEenergyThreshold_ && ET>threshold)
901 
902  } // else if (id.subdet()==HcalEndcap)
903 
904  } //for (HBHERecHitCollection::const_iterator HBHEiter=...)
905 
906  // Calculate normalized time
907  HEePlus>0 ? HEtPlus/=HEePlus : HEtPlus=10000;
908  HEeMinus>0 ? HEtMinus/=HEeMinus : HEtMinus=-10000;
909  HBePlus>0 ? HBtPlus/=HBePlus : HBtPlus=10000;
910  HBeMinus>0 ? HBtMinus/=HBeMinus : HBtMinus=-10000;
911 
912  ++HB_occupancy_[hbocc/10];
913  ++HE_occupancy_[heocc/10];
914  ++HB_occupancy_thresh_[hboccthresh/10];
915  ++HE_occupancy_thresh_[heoccthresh/10];
916  h_HBsizeVsLS->Fill(currentLS,hbocc);
917  h_HEsizeVsLS->Fill(currentLS,heocc);
918 
919  // loop over HO
920 
921  h_HOsizeVsLS->Fill(currentLS,hoHits.size());
922  int hoocc=0;
923  int hooccthresh=0;
924  for (HORecHitCollection::const_iterator HOiter=hoHits.begin(); HOiter!=hoHits.end(); ++HOiter)
925  { // loop over all hits
926  float en = HOiter->energy();
927  float ti = HOiter->time();
928  if (en>HOenergyThreshold_)
929  h_HOTimeVsEnergy->Fill(en,ti);
930 
931  HcalDetId id(HOiter->detid().rawId());
932  int ieta = id.ieta();
933  int iphi = id.iphi();
934  int depth = id.depth();
935 
936  if (en>0.5)
937  {
938  h_rechitieta_05->Fill(ieta);
939  h_rechitiphi_05->Fill(iphi);
940  if (en>1.)
941  {
942  h_rechitieta_10->Fill(ieta);
943  h_rechitiphi_10->Fill(iphi);
944  if (en>2.5)
945  {
946  h_rechitieta_25->Fill(ieta);
947  h_rechitiphi_25->Fill(iphi);
948  if (en>10.)
949  {
950  h_rechitieta_100->Fill(ieta);
951  h_rechitiphi_100->Fill(iphi);
952  }
953  }
954  }
955  }
956 
957 
958 
959  int calcEta = CalcEtaBin(HcalOuter,ieta,depth);
960  double fEta=fabs(0.5*(theHBHEEtaBounds[abs(ieta)-1]+theHBHEEtaBounds[abs(ieta)]));
961 
962  int rbxindex=logicalMap->getHcalFrontEndId(HOiter->detid()).rbxIndex();
963  int rm= logicalMap->getHcalFrontEndId(HOiter->detid()).rm();
964 
965  if (HOiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
966  h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
967  else if (HOiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
968  h_FlagMap_TIMEADD->Fill(rbxindex,rm);
969  else if (HOiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
970  h_FlagMap_TIMEERROR->Fill(rbxindex,rm);
971 
972 
973  //Looping over HO searching for flags --- cris
974  for (int f=0;f<32;f++)
975  {
976  if (HOiter->flagField(f))
977  HOflagcounter_[f]++;
978  }
979 
980  ++occupancy_[calcEta][iphi-1][depth-1];
981  energy_[calcEta][iphi-1][depth-1]+=en;
982  energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
983  time_[calcEta][iphi-1][depth-1]+=ti;
984  ++hoocc;
985  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
986  h_HOTime->Fill(ti);
987  else
988  ++HOtime_[int(ti-RECHITMON_TIME_MIN)];
989 
990  // We don't calculate HOplus/HOminus values (independent of trigger), so require min bias trigger
991  // along with E, ET thresholds directly in this HO loop:
992 
993  if (en>=HOenergyThreshold_
994  && en/cosh(fEta)>=HOETThreshold_
995  && passedMinBiasHLT==true
996  )
997  {
998  ++occupancy_thresh_[calcEta][iphi-1][depth-1];
999  energy_thresh_[calcEta][iphi-1][depth-1]+=en;
1000  energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
1001  time_thresh_[calcEta][iphi-1][depth-1]+=ti;
1002 
1003  ++hooccthresh;
1004  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
1005  h_HOThreshTime->Fill(ti);
1006  else
1007  ++HOtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
1008  }
1009  } // loop over all HO hits
1010 
1011  ++HO_occupancy_[hoocc/10];
1012  ++HO_occupancy_thresh_[hooccthresh/10];
1013 
1014  // loop over HF
1015  h_HFsizeVsLS->Fill(currentLS,hfHits.size());
1016 
1017  HtPlus=0; HtMinus=0;
1018 
1019  int hfocc=0;
1020  int hfoccthresh=0;
1021  for (HFRecHitCollection::const_iterator HFiter=hfHits.begin(); HFiter!=hfHits.end(); ++HFiter)
1022  { // loop over all hits
1023  float en = HFiter->energy();
1024  float ti = HFiter->time();
1025  if (en> HFenergyThreshold_)
1026  h_HFTimeVsEnergy->Fill(en,ti);
1027 
1028  HcalDetId id(HFiter->detid().rawId());
1029  int ieta = id.ieta();
1030  int iphi = id.iphi();
1031  int depth = id.depth();
1032 
1033  if (en>0.5)
1034  {
1035  h_rechitieta_05->Fill(ieta);
1036  h_rechitiphi_05->Fill(iphi);
1037  if (en>1.)
1038  {
1039  h_rechitieta_10->Fill(ieta);
1040  h_rechitiphi_10->Fill(iphi);
1041  if (en>2.5)
1042  {
1043  h_rechitieta_25->Fill(ieta);
1044  h_rechitiphi_25->Fill(iphi);
1045  if (en>10.)
1046  {
1047  h_rechitieta_100->Fill(ieta);
1048  h_rechitiphi_100->Fill(iphi);
1049  }
1050  }
1051  }
1052  }
1053 
1054  double fEta=fabs(0.5*(theHFEtaBounds[abs(ieta)-29]+theHFEtaBounds[abs(ieta)-28]));
1055  int calcEta = CalcEtaBin(HcalForward,ieta,depth);
1056 
1057  int rbxindex=logicalMap->getHcalFrontEndId(HFiter->detid()).rbxIndex();
1058  int rm= logicalMap->getHcalFrontEndId(HFiter->detid()).rm();
1059 
1060  h_HF_FlagCorr->Fill(HFiter->flagField(HcalCaloFlagLabels::HFDigiTime),HFiter->flagField(HcalCaloFlagLabels::HFLongShort));
1061  if (HFiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
1062  h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
1063  else if (HFiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
1064  h_FlagMap_TIMEADD->Fill(rbxindex,rm);
1065  else if (HFiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
1066  h_FlagMap_TIMEERROR->Fill(rbxindex,rm);
1067 
1068  if (HFiter->flagField(HcalCaloFlagLabels::HFDigiTime))
1069  {
1070  h_FlagMap_DIGITIME->Fill(rbxindex,rm);
1071  h_HFDigiTime_vs_LS->Fill(currentLS,1);
1072  }
1073  if (HFiter->flagField(HcalCaloFlagLabels::HFLongShort))
1074  {
1075  h_FlagMap_LONGSHORT->Fill(rbxindex,rm);
1076  h_HFLongShort_vs_LS->Fill(currentLS,1);
1077  }
1078  //Looping over HF searching for flags --- cris
1079  for (int f=0;f<32;f++)
1080  {
1081  if (HFiter->flagField(f))
1082  HFflagcounter_[f]++;
1083  }
1084 
1085  // Occupancy plots, without threshold
1086  ++occupancy_[calcEta][iphi-1][depth-1];
1087  energy_[calcEta][iphi-1][depth-1]+=en;
1088  energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
1089  time_[calcEta][iphi-1][depth-1]+=ti;
1090  ++hfocc;
1091  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
1092  h_HFTime->Fill(ti);
1093  else
1094  ++HFtime_[int(ti-RECHITMON_TIME_MIN)];
1095 
1096  ieta>0 ? HtPlus+=en/cosh(fEta) : HtMinus+=en/cosh(fEta); // add energy from all cells, or only those > threshold?
1097 
1098  if (en>=HFenergyThreshold_ &&
1099  en/cosh(fEta)>=HFETThreshold_
1100  )
1101  {
1102  // Occupancy plots require min bias trigger, along with thresholds exceeded
1103  if (passedMinBiasHLT)
1104  {
1105  ++occupancy_thresh_[calcEta][iphi-1][depth-1];
1106  energy_thresh_[calcEta][iphi-1][depth-1]+=en;
1107  energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
1108  time_thresh_[calcEta][iphi-1][depth-1]+=ti;
1109 
1110  ++hfoccthresh;
1111  if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
1112  h_HFThreshTime->Fill(ti);
1113  else
1114  ++HFtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
1115  }
1116 
1117  if (ieta>0)
1118  {
1119  HFtPlus+=en*ti;
1120  HFePlus+=en;
1121  hfpocc++;
1122  }
1123  else if (ieta<0)
1124  {
1125  HFtMinus+=en*ti;
1126  HFeMinus+=en;
1127  hfmocc++;
1128  }
1129  } // if (en>thresh, ET>thresh)
1130  } // loop over all HF hits
1131 
1132  ++HF_occupancy_[hfocc/10];
1133  ++HF_occupancy_thresh_[hfoccthresh/10];
1134 
1135 
1136  // Form event-wide variables (time averages, etc.), and plot them
1137 
1138  // Calculate weighted times. (Set tPlus, tMinus to overflow in case where total energy < 0)
1139  HFePlus>0 ? HFtPlus/=HFePlus : HFtPlus = 10000;
1140  HFeMinus>0 ? HFtMinus/=HFeMinus : HFtMinus = -10000;
1141 
1142  double mintime=99; // used to be min(tPlus,tMinus);
1143  double minHT=std::min(HtMinus,HtPlus);
1144  minHT==HtMinus ? mintime=HFtMinus : mintime = HFtPlus;
1145  //mintime = min(HFtPlus,HFtMinus); // I think we might want to use this value for mintime?
1146 
1147 
1148  h_LumiPlot_MinTime_vs_MinHT->Fill(minHT, mintime);
1149  h_LumiPlot_timeHT_HFM->Fill(HtMinus,HFtMinus);
1150  h_LumiPlot_timeHT_HFP->Fill(HtPlus,HFtPlus);
1151 
1152  if (passedMinBiasHLT==true)
1153  {
1154  h_LumiPlot_SumHT_HFPlus_vs_HFMinus->Fill(HtMinus,HtPlus);
1155  // HtMinus, HtPlus require no energy cuts for their contributing cells
1156  // HFeMinus, HFePlus require that cells be > threshold cut
1157 
1158  if (HtMinus>1 && HtPlus > 1) // is this the condition we want, or do we want hfmocc>0 && hfpocc >0?
1159  {
1160  h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus->Fill(HFeMinus,HFePlus);
1161  h_LumiPlot_timeHFPlus_vs_timeHFMinus->Fill(HFtMinus,HFtPlus);
1162 
1163  h_HFP_weightedTime->Fill(HFtPlus);
1164  h_HFM_weightedTime->Fill(HFtMinus);
1165  h_HBP_weightedTime->Fill(HBtPlus);
1166  h_HBM_weightedTime->Fill(HBtMinus);
1167 
1168  h_HEP_weightedTime->Fill(HEtPlus);
1169  h_HEM_weightedTime->Fill(HEtMinus);
1170 
1171  if (hepocc>0 && hemocc>0)
1172  {
1173  h_HEtimedifference->Fill(HEtPlus-HEtMinus);
1174  if (HEePlus-HEeMinus!=0) h_HEenergydifference->Fill((HEePlus-HEeMinus)/(HEePlus+HEeMinus));
1175  }
1176  if (hfpocc>0 && hfmocc>0) // which condition do we want?
1177  {
1178  h_HFtimedifference->Fill((HFtPlus)-(HFtMinus));
1179  if (HFePlus+HFeMinus!=0) h_HFenergydifference->Fill((HFePlus-HFeMinus)/(HFePlus+HFeMinus));
1180  }
1181 
1182  h_LumiPlot_LS_MinBiasEvents_notimecut->Fill(currentLS);
1183  h_LumiPlot_BX_MinBiasEvents_notimecut->Fill(BCN);
1184  if (fabs(HFtPlus-HFtMinus)<timediffThresh_)
1185  {
1186  h_LumiPlot_LS_MinBiasEvents->Fill(currentLS);
1187  h_LumiPlot_BX_MinBiasEvents->Fill(BCN);
1188  }
1189 
1190  HFP_HFM_Energy->Fill(HFeMinus/1000., HFePlus/1000.);
1191  }
1192 
1193  if (debug_>1) std::cout <<"\t<HcalRecHitMonitor:: HF averages> TPLUS = "<<HFtPlus<<" EPLUS = "<<HFePlus<<" TMINUS = "<<HFtMinus<<" EMINUS = "<<HFeMinus<<" Weighted Time Diff = "<<((HFtPlus)-(HFtMinus))<<std::endl;
1194 
1195 
1196  } // if (passedMinBiasHLT)
1197 
1198  if (passedHcalHLT && HtMinus>1 && HtPlus> 1 )
1199  {
1200  if (hfpocc>0 && hfmocc>0)
1201  {
1202  h_HF_HcalHLT_weightedtimedifference->Fill(HFtPlus-HFtMinus);
1203  if (HFePlus+HFeMinus!=0) h_HF_HcalHLT_energydifference->Fill((HFePlus-HFeMinus)/(HFePlus+HFeMinus));
1204  }
1205  if (hepocc>0 && hemocc>0)
1206  {
1207  h_HE_HcalHLT_weightedtimedifference->Fill(HEtPlus-HEtMinus);
1208  if (HEePlus-HEeMinus!=0) h_HE_HcalHLT_energydifference->Fill((HEePlus-HEeMinus)/(HEePlus+HEeMinus));
1209  }
1210 
1211  h_LumiPlot_LS_HcalHLTEvents_notimecut->Fill(currentLS);
1212  h_LumiPlot_BX_HcalHLTEvents_notimecut->Fill(BCN);
1213  if (fabs(HFtPlus-HFtMinus)<timediffThresh_)
1214  {
1215  h_LumiPlot_LS_HcalHLTEvents->Fill(currentLS);
1216  h_LumiPlot_BX_HcalHLTEvents->Fill(BCN);
1217  }
1218  } // passsed Hcal HLT
1219 
1220  return;
1221 } // void HcalRecHitMonitor::processEvent_rechitenergy
1222 
1223 /* --------------------------------------- */
1224 
1225 
1226 void HcalRecHitMonitor::endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,
1227  const edm::EventSetup& c)
1228 
1229 {
1230  // don't fill lumi block information if it's already been filled
1231  if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
1232  fill_Nevents();
1233  return;
1234 } //endLuminosityBlock
1235 
1236 
1237 void HcalRecHitMonitor::fill_Nevents(void)
1238 {
1239  // looking at the contents of HbFlagcounters
1240  if (debug_>0)
1241  {
1242  for (int k = 0; k < 32; k++){
1243  std::cout << "<HcalRecHitMonitor::fill_Nevents> HF Flag counter: Bin #" << k+1 << " = "<< HFflagcounter_[k] << std::endl;
1244  }
1245  }
1246 
1247  for (int i=0;i<32;i++)
1248  {
1249  h_HBflagcounter->Fill(i,HBflagcounter_[i]);
1250  h_HEflagcounter->Fill(i,HEflagcounter_[i]);
1251  h_HOflagcounter->Fill(i,HOflagcounter_[i]);
1252  h_HFflagcounter->Fill(i,HFflagcounter_[i]);
1253  HBflagcounter_[i]=0;
1254  HEflagcounter_[i]=0;
1255  HOflagcounter_[i]=0;
1256  HFflagcounter_[i]=0;
1257  }
1258 
1259  // Fill Occupancy & Sum Energy, Time plots
1260  int myieta=-1;
1261  if (ievt_>0)
1262  {
1263  for (int mydepth=0;mydepth<4;++mydepth)
1264  {
1265  for (int eta=0;eta<OccupancyByDepth.depth[mydepth]->getNbinsX();++eta)
1266  {
1267  myieta=CalcIeta(eta,mydepth+1);
1268 
1269  for (int phi=0;phi<72;++phi)
1270  {
1271  if (occupancy_[eta][phi][mydepth]>0)
1272  {
1273  h_rechitieta->Fill(myieta,occupancy_[eta][phi][mydepth]);
1274  h_rechitiphi->Fill(phi+1,occupancy_[eta][phi][mydepth]);
1275  }
1276  if (occupancy_thresh_[eta][phi][mydepth]>0)
1277  {
1278  h_rechitieta_thresh->Fill(myieta,occupancy_thresh_[eta][phi][mydepth]);
1279  h_rechitiphi_thresh->Fill(phi+1,occupancy_thresh_[eta][phi][mydepth]);
1280  }
1281  OccupancyByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,occupancy_[eta][phi][mydepth]);
1282  SumEnergyByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,energy_[eta][phi][mydepth]);
1283  SqrtSumEnergy2ByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,sqrt(energy2_[eta][phi][mydepth]));
1284  SumTimeByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,time_[eta][phi][mydepth]);
1285 
1286  OccupancyThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,occupancy_thresh_[eta][phi][mydepth]);
1287  SumEnergyThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,energy_thresh_[eta][phi][mydepth]);
1288  SqrtSumEnergy2ThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,sqrt(energy2_thresh_[eta][phi][mydepth]));
1289  SumTimeThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,time_thresh_[eta][phi][mydepth]);
1290  } // for (int phi=0;phi<72;++phi)
1291  } // for (int eta=0;eta<OccupancyByDepth...;++eta)
1292  } // for (int mydepth=0;...)
1293 
1294  FillUnphysicalHEHFBins(OccupancyByDepth);
1295  FillUnphysicalHEHFBins(OccupancyThreshByDepth);
1296  FillUnphysicalHEHFBins(SumEnergyByDepth);
1297  FillUnphysicalHEHFBins(SqrtSumEnergy2ByDepth);
1298  FillUnphysicalHEHFBins(SumEnergyThreshByDepth);
1299  FillUnphysicalHEHFBins(SqrtSumEnergy2ThreshByDepth);
1300  FillUnphysicalHEHFBins(SumTimeByDepth);
1301  FillUnphysicalHEHFBins(SumTimeThreshByDepth);
1302 
1303  } // if (ievt_>0)
1304 
1305  // Fill subdet plots
1306 
1307  for (int i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)
1308  {
1309  if (HBtime_[i]!=0)
1310  {
1311  h_HBTime->setBinContent(i+1,HBtime_[i]);
1312  }
1313  if (HBtime_thresh_[i]!=0)
1314  {
1315  h_HBThreshTime->setBinContent(i+1,HBtime_thresh_[i]);
1316  }
1317  if (HEtime_[i]!=0)
1318  {
1319 
1320  h_HETime->setBinContent(i+1,HEtime_[i]);
1321  }
1322  if (HEtime_thresh_[i]!=0)
1323  {
1324  h_HEThreshTime->setBinContent(i+1,HEtime_thresh_[i]);
1325  }
1326  if (HOtime_[i]!=0)
1327  {
1328  h_HOTime->setBinContent(i+1,HOtime_[i]);
1329  }
1330  if (HOtime_thresh_[i]!=0)
1331  {
1332  h_HOThreshTime->setBinContent(i+1,HOtime_thresh_[i]);
1333  }
1334  if (HFtime_[i]!=0)
1335  {
1336  h_HFTime->setBinContent(i+1,HFtime_[i]);
1337  }
1338  if (HFtime_thresh_[i]!=0)
1339  {
1340  h_HFThreshTime->setBinContent(i+1,HFtime_thresh_[i]);
1341  }
1342  } // for (int i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)
1343 
1344  for (int i=0;i<260;++i)
1345  {
1346  if (HB_occupancy_[i]>0)
1347  {
1348  h_HBOccupancy->setBinContent(i+1,HB_occupancy_[i]);
1349  }
1350  if (HB_occupancy_thresh_[i]>0)
1351  {
1352  h_HBThreshOccupancy->setBinContent(i+1,HB_occupancy_thresh_[i]);
1353  }
1354  if (HE_occupancy_[i]>0)
1355  {
1356  h_HEOccupancy->setBinContent(i+1,HE_occupancy_[i]);
1357  }
1358  if (HE_occupancy_thresh_[i]>0)
1359  {
1360  h_HEThreshOccupancy->setBinContent(i+1,HE_occupancy_thresh_[i]);
1361  }
1362  }//for (int i=0;i<260;++i)
1363 
1364  for (int i=0;i<217;++i)
1365  {
1366  if (HO_occupancy_[i]>0)
1367  {
1368  h_HOOccupancy->setBinContent(i+1,HO_occupancy_[i]);
1369  }
1370  if (HO_occupancy_thresh_[i]>0)
1371  {
1372  h_HOThreshOccupancy->setBinContent(i+1,HO_occupancy_thresh_[i]);
1373  }
1374  }// for (int i=0;i<217;++i)
1375 
1376  for (int i=0;i<173;++i)
1377  {
1378  if (HF_occupancy_[i]>0)
1379  {
1380  h_HFOccupancy->setBinContent(i+1,HF_occupancy_[i]);
1381  }
1382  if (HF_occupancy_thresh_[i]>0)
1383  {
1384  h_HFThreshOccupancy->setBinContent(i+1,HF_occupancy_thresh_[i]);
1385  }
1386  }// for (int i=0;i<173;++i)
1387 
1388  //zeroCounters();
1389 
1390  if (debug_>0)
1391  std::cout <<"<HcalRecHitMonitor::fill_Nevents> FILLED REC HIT CELL PLOTS"<<std::endl;
1392 
1393 } // void HcalRecHitMonitor::fill_Nevents(void)
1394 
1395 
1396 void HcalRecHitMonitor::zeroCounters(void)
1397 {
1398  // Set all histogram counters back to zero
1399 
1400  for (int i=0;i<32;++i)
1401  {
1402  HBflagcounter_[i]=0;
1403  HEflagcounter_[i]=0;
1404  HOflagcounter_[i]=0;
1405  HFflagcounter_[i]=0;
1406 
1407  }
1408  // TH2F counters
1409  for (int i=0;i<85;++i)
1410  {
1411  for (int j=0;j<72;++j)
1412  {
1413  for (int k=0;k<4;++k)
1414  {
1415  occupancy_[i][j][k]=0;
1416  occupancy_thresh_[i][j][k]=0;
1417  energy_[i][j][k]=0;
1418  energy2_[i][j][k]=0;
1419  energy_thresh_[i][j][k]=0;
1420  energy2_thresh_[i][j][k]=0;
1421  time_[i][j][k]=0;
1422  time_thresh_[i][j][k]=0;
1423  }
1424  } // for (int j=0;j<PHIBINS;++j)
1425  } // for (int i=0;i<87;++i)
1426 
1427  // TH1F counters
1428 
1429  for (int i=0;i<200;++i)
1430  {
1431  HFenergyLong_[i]=0;
1432  HFenergyLong_thresh_[i]=0;
1433  HFenergyShort_[i]=0;
1434  HFenergyShort_thresh_[i]=0;
1435  }
1436 
1437  // time
1438  for (int i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)
1439  {
1440  HBtime_[i]=0;
1441  HBtime_thresh_[i]=0;
1442  HEtime_[i]=0;
1443  HEtime_thresh_[i]=0;
1444  HOtime_[i]=0;
1445  HOtime_thresh_[i]=0;
1446  HFtime_[i]=0;
1447  HFtime_thresh_[i]=0;
1448  HFtimeLong_[i]=0;
1449  HFtimeLong_thresh_[i]=0;
1450  HFtimeShort_[i]=0;
1451  HFtimeShort_thresh_[i]=0;
1452  }
1453 
1454  // occupancy
1455  for (int i=0;i<865;++i)
1456  {
1457  if (i<260)
1458  {
1459  HB_occupancy_[i]=0;
1460  HE_occupancy_[i]=0;
1461  HB_occupancy_thresh_[i]=0;
1462  HE_occupancy_thresh_[i]=0;
1463  }
1464  if (i<218)
1465  {
1466  HO_occupancy_[i]=0;
1467  HO_occupancy_thresh_[i]=0;
1468  }
1469  if (i<174)
1470  {
1471  HF_occupancy_[i]=0;
1472  HF_occupancy_thresh_[i]=0;
1473  }
1474 
1475  HFlong_occupancy_[i] =0;
1476  HFshort_occupancy_[i]=0;
1477  HFlong_occupancy_thresh_[i] =0;
1478  HFshort_occupancy_thresh_[i]=0;
1479  } // for (int i=0;i<865;++i)
1480 
1481  return;
1482 } //void HcalRecHitMonitor::zeroCounters(void)
1483 
1484 DEFINE_FWK_MODULE(HcalRecHitMonitor);
HcalRecHitMonitor::HFenergyShort_
double HFenergyShort_[200]
Definition: HcalRecHitMonitor.h:106
HcalRecHitMonitor::h_rechitieta_100
MonitorElement * h_rechitieta_100
Definition: HcalRecHitMonitor.h:130
HcalRecHitMonitor::time_
double time_[85][72][4]
Definition: HcalRecHitMonitor.h:85
HcalRecHitMonitor::HB_occupancy_
double HB_occupancy_[260]
Definition: HcalRecHitMonitor.h:90
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
HcalRecHitMonitor::h_HOThreshTime
MonitorElement * h_HOThreshTime
Definition: HcalRecHitMonitor.h:155
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:169
HcalRecHitMonitor::h_HFP_weightedTime
MonitorElement * h_HFP_weightedTime
Definition: HcalRecHitMonitor.h:219
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:208
HcalRecHitMonitor::h_HFLongShort_vs_LS
MonitorElement * h_HFLongShort_vs_LS
Definition: HcalRecHitMonitor.h:177
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:800
HcalRecHitMonitor::h_HFtimedifference
MonitorElement * h_HFtimedifference
Definition: HcalRecHitMonitor.h:187
HcalRecHitMonitor::h_LumiPlot_timeHT_HFP
MonitorElement * h_LumiPlot_timeHT_HFP
Definition: HcalRecHitMonitor.h:216
HcalRecHitMonitor::h_HBHE_FlagCorr
MonitorElement * h_HBHE_FlagCorr
Definition: HcalRecHitMonitor.h:183
HcalRecHitMonitor::h_HBsizeVsLS
MonitorElement * h_HBsizeVsLS
Definition: HcalRecHitMonitor.h:139
HcalRecHitMonitor::h_LumiPlot_MinTime_vs_MinHT
MonitorElement * h_LumiPlot_MinTime_vs_MinHT
Definition: HcalRecHitMonitor.h:214
HcalRecHitMonitor::HFshort_occupancy_
double HFshort_occupancy_[865]
Definition: HcalRecHitMonitor.h:114
HcalRecHitMonitor::OccupancyThreshByDepth
EtaPhiHists OccupancyThreshByDepth
Definition: HcalRecHitMonitor.h:70
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:206
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:514
HcalRecHitMonitor::h_HF_HcalHLT_weightedtimedifference
MonitorElement * h_HF_HcalHLT_weightedtimedifference
Definition: HcalRecHitMonitor.h:192
HcalRecHitMonitor::HF_occupancy_thresh_
double HF_occupancy_thresh_[174]
Definition: HcalRecHitMonitor.h:111
HcalRecHitMonitor::h_HBThreshTime
MonitorElement * h_HBThreshTime
Definition: HcalRecHitMonitor.h:145
HcalRecHitMonitor::processEvent
void processEvent(const HBHERecHitCollection &hbHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, int BCN, const edm::Event &iEvent)
Definition: HcalRecHitMonitor.cc:616
HcalRecHitMonitor::HFenergyLong_
double HFenergyLong_[200]
Definition: HcalRecHitMonitor.h:102
HcalRecHitMonitor::h_FlagMap_PULSESHAPE
MonitorElement * h_FlagMap_PULSESHAPE
Definition: HcalRecHitMonitor.h:170
HcalRecHitMonitor::HFshort_occupancy_thresh_
double HFshort_occupancy_thresh_[865]
Definition: HcalRecHitMonitor.h:115
HcalRecHitMonitor::h_HBHEPulseShape_vs_LS
MonitorElement * h_HBHEPulseShape_vs_LS
Definition: HcalRecHitMonitor.h:180
HcalRecHitMonitor::h_HF_FlagCorr
MonitorElement * h_HF_FlagCorr
Definition: HcalRecHitMonitor.h:182
HcalRecHitMonitor::h_HOflagcounter
MonitorElement * h_HOflagcounter
Definition: HcalRecHitMonitor.h:166
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
HcalRecHitMonitor::h_rechitiphi_05
MonitorElement * h_rechitiphi_05
Definition: HcalRecHitMonitor.h:131
HcalBaseDQMonitor::AllowedCalibTypes_
std::vector< int > AllowedCalibTypes_
Definition: HcalBaseDQMonitor.h:82
HcalRecHitMonitor::HEtime_
double HEtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:92
HcalRecHitMonitor::h_rechitieta_10
MonitorElement * h_rechitieta_10
Definition: HcalRecHitMonitor.h:128
HcalRecHitMonitor::HO_occupancy_thresh_
double HO_occupancy_thresh_[218]
Definition: HcalRecHitMonitor.h:99
HcalRecHitMonitor::beginRun
void beginRun(const edm::Run &run, const edm::EventSetup &c)
Definition: HcalRecHitMonitor.cc:486
HcalRecHitMonitor::h_LumiPlot_timeHFPlus_vs_timeHFMinus
MonitorElement * h_LumiPlot_timeHFPlus_vs_timeHFMinus
Definition: HcalRecHitMonitor.h:204
HcalRecHitMonitor::HBtime_thresh_
double HBtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:89
HcalRecHitMonitor::h_HFThreshTime
MonitorElement * h_HFThreshTime
Definition: HcalRecHitMonitor.h:160
HcalRecHitMonitor::h_rechitieta_25
MonitorElement * h_rechitieta_25
Definition: HcalRecHitMonitor.h:129
HcalBaseDQMonitor::cleanup
virtual void cleanup(void)
Definition: HcalBaseDQMonitor.cc:115
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:81
HcalRecHitMonitor::h_rechitiphi_10
MonitorElement * h_rechitiphi_10
Definition: HcalRecHitMonitor.h:132
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:62
HcalRecHitMonitor::h_HFTime
MonitorElement * h_HFTime
Definition: HcalRecHitMonitor.h:159
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:876
HcalRecHitMonitor::h_HFThreshOccupancy
MonitorElement * h_HFThreshOccupancy
Definition: HcalRecHitMonitor.h:162
HcalRecHitMonitor::h_FlagMap_TIMEADD
MonitorElement * h_FlagMap_TIMEADD
Definition: HcalRecHitMonitor.h:173
abs
#define abs(x)
Definition: mlp_lapack.h:159
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:59
HcalRecHitMonitor::h_LumiPlot_LS_MinBiasEvents
MonitorElement * h_LumiPlot_LS_MinBiasEvents
Definition: HcalRecHitMonitor.h:198
HcalRecHitMonitor::h_rechitieta
MonitorElement * h_rechitieta
Definition: HcalRecHitMonitor.h:124
HcalRecHitMonitor::h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus
MonitorElement * h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus
Definition: HcalRecHitMonitor.h:206
HcalRecHitMonitor::HF_occupancy_
double HF_occupancy_[174]
Definition: HcalRecHitMonitor.h:110
min
#define min(a, b)
Definition: mlp_lapack.h:161
HcalRecHitMonitor::l1gtLabel_
edm::InputTag l1gtLabel_
Definition: HcalRecHitMonitor.h:235
HcalRecHitMonitor::h_HFflagcounter
MonitorElement * h_HFflagcounter
Definition: HcalRecHitMonitor.h:167
edm::TriggerNames::size
Strings::size_type size() const
Definition: TriggerNames.cc:39
HcalRecHitMonitor::h_HFM_weightedTime
MonitorElement * h_HFM_weightedTime
Definition: HcalRecHitMonitor.h:220
MonitorElement::update
void update(void)
Mark the object updated.
Definition: MonitorElement.h:123
HcalRecHitMonitor::h_HBTimeVsEnergy
MonitorElement * h_HBTimeVsEnergy
Definition: HcalRecHitMonitor.h:226
eta
T eta() const
Definition: Basic3DVectorLD.h:168
HcalRecHitMonitor::hfRechitLabel_
edm::InputTag hfRechitLabel_
Definition: HcalRecHitMonitor.h:234
HcalRecHitMonitor::h_HEOccupancy
MonitorElement * h_HEOccupancy
Definition: HcalRecHitMonitor.h:151
HcalRecHitMonitor::h_LumiPlot_LS_allevents
MonitorElement * h_LumiPlot_LS_allevents
Definition: HcalRecHitMonitor.h:197
HcalRecHitMonitor::h_HEtimedifference
MonitorElement * h_HEtimedifference
Definition: HcalRecHitMonitor.h:189
HcalRecHitMonitor::h_HEM_weightedTime
MonitorElement * h_HEM_weightedTime
Definition: HcalRecHitMonitor.h:222
HcalRecHitMonitor::h_HE_HcalHLT_weightedtimedifference
MonitorElement * h_HE_HcalHLT_weightedtimedifference
Definition: HcalRecHitMonitor.h:194
HcalRecHitMonitor::h_FlagMap_TIMEERROR
MonitorElement * h_FlagMap_TIMEERROR
Definition: HcalRecHitMonitor.h:175
HcalRecHitMonitor::h_HF_HcalHLT_energydifference
MonitorElement * h_HF_HcalHLT_energydifference
Definition: HcalRecHitMonitor.h:193
HcalRecHitMonitor::SumEnergyByDepth
EtaPhiHists SumEnergyByDepth
Definition: HcalRecHitMonitor.h:72
HcalRecHitMonitor::h_HOThreshOccupancy
MonitorElement * h_HOThreshOccupancy
Definition: HcalRecHitMonitor.h:157
HcalRecHitMonitor::HE_occupancy_thresh_
double HE_occupancy_thresh_[260]
Definition: HcalRecHitMonitor.h:95
DQMStore::bookFloat
MonitorElement * bookFloat(const char *name)
Book float.
Definition: DQMStore.cc:451
HcalRecHitMonitor::occupancy_
unsigned int occupancy_[85][72][4]
Definition: HcalRecHitMonitor.h:79
DQMStore::getAllContents
std::vector< MonitorElement * > getAllContents(const std::string &path) const
Definition: DQMStore.cc:1469
HcalRecHitMonitor::h_rechitiphi_100
MonitorElement * h_rechitiphi_100
Definition: HcalRecHitMonitor.h:134
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
HcalRecHitMonitor::h_HETime
MonitorElement * h_HETime
Definition: HcalRecHitMonitor.h:149
HcalRecHitMonitor::h_rechitieta_thresh
MonitorElement * h_rechitieta_thresh
Definition: HcalRecHitMonitor.h:136
edm::LuminosityBlockBase::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockBase.h:44
HcalRecHitMonitor::OccupancyByDepth
EtaPhiHists OccupancyByDepth
Definition: HcalRecHitMonitor.h:69
HcalRecHitMonitor::h_HOOccupancy
MonitorElement * h_HOOccupancy
Definition: HcalRecHitMonitor.h:156
HcalRecHitMonitor::h_HFenergydifference
MonitorElement * h_HFenergydifference
Definition: HcalRecHitMonitor.h:188
HcalRecHitMonitor::time_thresh_
double time_thresh_[85][72][4]
Definition: HcalRecHitMonitor.h:86
HcalRecHitMonitor::h_HEsizeVsLS
MonitorElement * h_HEsizeVsLS
Definition: HcalRecHitMonitor.h:140
HcalRecHitMonitor::h_LumiPlot_BX_MinBiasEvents_notimecut
MonitorElement * h_LumiPlot_BX_MinBiasEvents_notimecut
Definition: HcalRecHitMonitor.h:210
HcalRecHitMonitor::SumTimeThreshByDepth
EtaPhiHists SumTimeThreshByDepth
Definition: HcalRecHitMonitor.h:77
HcalRecHitMonitor::h_HOsizeVsLS
MonitorElement * h_HOsizeVsLS
Definition: HcalRecHitMonitor.h:141
HcalLogicalMapGenerator::createMap
HcalLogicalMap createMap(unsigned int mapIOV=4)
Definition: HcalLogicalMapGenerator.cc:38
iEvent
int iEvent
Definition: GenABIO.cc:243
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalBaseDQMonitor::beginRun
virtual void beginRun(const edm::Run &run, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:83
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
DQMStore::removeContents
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
Definition: DQMStore.cc:2345
HcalRecHitMonitor::h_LumiPlot_BX_HcalHLTEvents_notimecut
MonitorElement * h_LumiPlot_BX_HcalHLTEvents_notimecut
Definition: HcalRecHitMonitor.h:212
HcalRecHitMonitor::SqrtSumEnergy2ThreshByDepth
EtaPhiHists SqrtSumEnergy2ThreshByDepth
Definition: HcalRecHitMonitor.h:75
HcalRecHitMonitor::SumTimeByDepth
EtaPhiHists SumTimeByDepth
Definition: HcalRecHitMonitor.h:76
HcalRecHitMonitor::HFtimeLong_
double HFtimeLong_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:104
HcalRecHitMonitor::processEvent_rechit
void processEvent_rechit(const HBHERecHitCollection &hbheHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, bool passedHcalHLT, bool passedMinBiasHLT, int BCN)
Definition: HcalRecHitMonitor.cc:680
DQMStore::bookString
MonitorElement * bookString(const char *name, const char *value)
Book string.
Definition: DQMStore.cc:480
HcalRecHitMonitor::occupancy_thresh_
unsigned int occupancy_thresh_[85][72][4]
Definition: HcalRecHitMonitor.h:80
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
HcalRecHitMonitor::h_HETimeVsEnergy
MonitorElement * h_HETimeVsEnergy
Definition: HcalRecHitMonitor.h:227
HcalRecHitMonitor::HEtime_thresh_
double HEtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:93
HcalRecHitMonitor::h_HEenergydifference
MonitorElement * h_HEenergydifference
Definition: HcalRecHitMonitor.h:190
submit.rm
string rm
Definition: submit.py:76
HcalRecHitMonitor::h_HFTimeVsEnergy
MonitorElement * h_HFTimeVsEnergy
Definition: HcalRecHitMonitor.h:229
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:38
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:32
HcalRecHitMonitor::h_LumiPlot_BX_HcalHLTEvents
MonitorElement * h_LumiPlot_BX_HcalHLTEvents
Definition: HcalRecHitMonitor.h:211
HcalRecHitMonitor::HFtimeShort_
double HFtimeShort_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:108
j
int j
Definition: DBlmapReader.cc:9
HcalRecHitMonitor::HFlong_occupancy_thresh_
double HFlong_occupancy_thresh_[865]
Definition: HcalRecHitMonitor.h:113
HcalRecHitMonitor::h_HBM_weightedTime
MonitorElement * h_HBM_weightedTime
Definition: HcalRecHitMonitor.h:224
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
HcalRecHitMonitor::HcalHLTBits_
std::vector< std::string > HcalHLTBits_
Definition: HcalRecHitMonitor.h:237
HcalRecHitMonitor::h_HFsizeVsLS
MonitorElement * h_HFsizeVsLS
Definition: HcalRecHitMonitor.h:142
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:828
HcalRecHitMonitor::SqrtSumEnergy2ByDepth
EtaPhiHists SqrtSumEnergy2ByDepth
Definition: HcalRecHitMonitor.h:73
HcalRecHitMonitor::h_HBOccupancy
MonitorElement * h_HBOccupancy
Definition: HcalRecHitMonitor.h:146
HcalRecHitMonitor::h_HFDigiTime_vs_LS
MonitorElement * h_HFDigiTime_vs_LS
Definition: HcalRecHitMonitor.h:178
HcalRecHitMonitor::analyze
void analyze(const edm::Event &, const edm::EventSetup &)
Definition: HcalRecHitMonitor.cc:577
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:355
HcalRecHitMonitor::h_HBHEHPDMult_vs_LS
MonitorElement * h_HBHEHPDMult_vs_LS
Definition: HcalRecHitMonitor.h:179
HcalRecHitMonitor::energy2_
double energy2_[85][72][4]
Definition: HcalRecHitMonitor.h:82
HcalRecHitMonitor::h_HOTime
MonitorElement * h_HOTime
Definition: HcalRecHitMonitor.h:154
gen::k
int k[5][pyjets_maxn]
Definition: Pythia6Hadronizer.cc:86
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:302
HcalRecHitMonitor::h_rechitiphi_25
MonitorElement * h_rechitiphi_25
Definition: HcalRecHitMonitor.h:133
HcalRecHitMonitor::h_FlagMap_TIMESUBTRACT
MonitorElement * h_FlagMap_TIMESUBTRACT
Definition: HcalRecHitMonitor.h:174
HcalLogicalMap::getHcalFrontEndId
const HcalFrontEndId getHcalFrontEndId(const DetId &)
Definition: HcalLogicalMap.cc:262
HcalRecHitMonitor::HOtime_
double HOtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:96
HcalRecHitMonitor::energy2_thresh_
double energy2_thresh_[85][72][4]
Definition: HcalRecHitMonitor.h:84
HcalRecHitMonitor::h_HEP_weightedTime
MonitorElement * h_HEP_weightedTime
Definition: HcalRecHitMonitor.h:221
HcalRecHitMonitor::h_FlagMap_LONGSHORT
MonitorElement * h_FlagMap_LONGSHORT
Definition: HcalRecHitMonitor.h:172
HcalRecHitMonitor::h_LumiPlot_BX_allevents
MonitorElement * h_LumiPlot_BX_allevents
Definition: HcalRecHitMonitor.h:208
HcalRecHitMonitor::hltresultsLabel_
edm::InputTag hltresultsLabel_
Definition: HcalRecHitMonitor.h:236
HcalRecHitMonitor::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
Definition: HcalRecHitMonitor.cc:1226
HcalRecHitMonitor::energy_thresh_
double energy_thresh_[85][72][4]
Definition: HcalRecHitMonitor.h:83
theHFEtaBounds
static const double theHFEtaBounds[]
Definition: HcalHardcodeGeometryData.h:12
HcalRecHitMonitor::MinBiasHLTBits_
std::vector< std::string > MinBiasHLTBits_
Definition: HcalRecHitMonitor.h:238
HcalRecHitMonitor::h_HBThreshOccupancy
MonitorElement * h_HBThreshOccupancy
Definition: HcalRecHitMonitor.h:147
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
HcalRecHitMonitor::endRun
void endRun(const edm::Run &run, const edm::EventSetup &c)
Definition: HcalRecHitMonitor.cc:518
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1473
HcalRecHitMonitor::h_HE_HcalHLT_energydifference
MonitorElement * h_HE_HcalHLT_energydifference
Definition: HcalRecHitMonitor.h:195
HcalRecHitMonitor::h_rechitiphi
MonitorElement * h_rechitiphi
Definition: HcalRecHitMonitor.h:125
HcalRecHitMonitor::h_LumiPlot_LS_MinBiasEvents_notimecut
MonitorElement * h_LumiPlot_LS_MinBiasEvents_notimecut
Definition: HcalRecHitMonitor.h:199
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
HcalRecHitMonitor::h_HFOccupancy
MonitorElement * h_HFOccupancy
Definition: HcalRecHitMonitor.h:161
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
HcalRecHitMonitor::HFtime_
double HFtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:100
HcalRecHitMonitor::h_rechitieta_05
MonitorElement * h_rechitieta_05
Definition: HcalRecHitMonitor.h:127
HcalRecHitMonitor::HFtime_thresh_
double HFtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:101
HcalBaseDQMonitor::SetupEtaPhiHists
void SetupEtaPhiHists(EtaPhiHists &hh, std::string Name, std::string Units)
Definition: HcalBaseDQMonitor.h:72
RECHITMON_TIME_MIN
#define RECHITMON_TIME_MIN
Definition: HcalRecHitMonitor.h:8
HcalRecHitMonitor::h_HEflagcounter
MonitorElement * h_HEflagcounter
Definition: HcalRecHitMonitor.h:165
HcalRecHitMonitor::hbheRechitLabel_
edm::InputTag hbheRechitLabel_
Definition: HcalRecHitMonitor.h:234
HcalRecHitMonitor::h_HBTime
MonitorElement * h_HBTime
Definition: HcalRecHitMonitor.h:144
HcalRecHitMonitor::h_rechitiphi_thresh
MonitorElement * h_rechitiphi_thresh
Definition: HcalRecHitMonitor.h:137
HcalRecHitMonitor::hoRechitLabel_
edm::InputTag hoRechitLabel_
Definition: HcalRecHitMonitor.h:234
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:228
HcalRecHitMonitor::h_HBP_weightedTime
MonitorElement * h_HBP_weightedTime
Definition: HcalRecHitMonitor.h:223
HcalRecHitMonitor::HB_occupancy_thresh_
double HB_occupancy_thresh_[260]
Definition: HcalRecHitMonitor.h:91
HcalRecHitMonitor::h_LumiPlot_BX_MinBiasEvents
MonitorElement * h_LumiPlot_BX_MinBiasEvents
Definition: HcalRecHitMonitor.h:209
HcalRecHitMonitor::h_FlagMap_HPDMULT
MonitorElement * h_FlagMap_HPDMULT
Definition: HcalRecHitMonitor.h:169
RECHITMON_TIME_MAX
#define RECHITMON_TIME_MAX
Definition: HcalRecHitMonitor.h:9
HcalRecHitMonitor::h_HOTimeVsEnergy
MonitorElement * h_HOTimeVsEnergy
Definition: HcalRecHitMonitor.h:228
HcalRecHitMonitor::h_HEThreshTime
MonitorElement * h_HEThreshTime
Definition: HcalRecHitMonitor.h:150
HcalRecHitMonitor::HcalRecHitMonitor
HcalRecHitMonitor(const edm::ParameterSet &ps)
Definition: HcalRecHitMonitor.cc:16
HcalRecHitMonitor::HFtimeShort_thresh_
double HFtimeShort_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:109
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:41
HcalRecHitMonitor::SumEnergyThreshByDepth
EtaPhiHists SumEnergyThreshByDepth
Definition: HcalRecHitMonitor.h:74
HcalObjRepresent::Reset
void Reset(std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:620
asciidump.s
string s
Definition: asciidump.py:422
HcalRecHitMonitor::h_HEThreshOccupancy
MonitorElement * h_HEThreshOccupancy
Definition: HcalRecHitMonitor.h:152
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:642
HcalBaseDQMonitor::setup
virtual void setup(void)
Definition: HcalBaseDQMonitor.cc:120
HcalRecHitMonitor::HE_occupancy_
double HE_occupancy_[260]
Definition: HcalRecHitMonitor.h:94
HcalRecHitMonitor::h_LumiPlot_LS_HcalHLTEvents_notimecut
MonitorElement * h_LumiPlot_LS_HcalHLTEvents_notimecut
Definition: HcalRecHitMonitor.h:201
HcalRecHitMonitor::HFlong_occupancy_
double HFlong_occupancy_[865]
Definition: HcalRecHitMonitor.h:112
HcalRecHitMonitor::h_LumiPlot_LS_HcalHLTEvents
MonitorElement * h_LumiPlot_LS_HcalHLTEvents
Definition: HcalRecHitMonitor.h:200
theHBHEEtaBounds
static const double theHBHEEtaBounds[]
Definition: HcalHardcodeGeometryData.h:5
HcalRecHitMonitor::HBtime_
double HBtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:88
HcalRecHitMonitor::HFP_HFM_Energy
MonitorElement * HFP_HFM_Energy
Definition: HcalRecHitMonitor.h:230
HcalRecHitMonitor::HO_occupancy_
double HO_occupancy_[218]
Definition: HcalRecHitMonitor.h:98
HcalRecHitMonitor::HFtimeLong_thresh_
double HFtimeLong_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:105
HcalRecHitMonitor::h_FlagMap_DIGITIME
MonitorElement * h_FlagMap_DIGITIME
Definition: HcalRecHitMonitor.h:171
HcalRecHitMonitor::h_LumiPlot_timeHT_HFM
MonitorElement * h_LumiPlot_timeHT_HFM
Definition: HcalRecHitMonitor.h:215
HcalRecHitMonitor::HFenergyShort_thresh_
double HFenergyShort_thresh_[200]
Definition: HcalRecHitMonitor.h:107
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
HcalRecHitMonitor::HFenergyLong_thresh_
double HFenergyLong_thresh_[200]
Definition: HcalRecHitMonitor.h:103
HcalRecHitMonitor::h_HBflagcounter
MonitorElement * h_HBflagcounter
Definition: HcalRecHitMonitor.h:164
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:232
HcalRecHitMonitor::HOtime_thresh_
double HOtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
Definition: HcalRecHitMonitor.h:97
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:295
edm::Run
Definition: Run.h:32
HcalRecHitMonitor::h_TriggeredEvents
MonitorElement * h_TriggeredEvents
Definition: HcalRecHitMonitor.h:218
HcalRecHitMonitor::h_LumiPlot_SumHT_HFPlus_vs_HFMinus
MonitorElement * h_LumiPlot_SumHT_HFPlus_vs_HFMinus
Definition: HcalRecHitMonitor.h:203
HcalRecHitMonitor::logicalMap
HcalLogicalMap * logicalMap
Definition: HcalRecHitMonitor.h:66
HcalRecHitMonitor::energy_
double energy_[85][72][4]
Definition: HcalRecHitMonitor.h:81
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5