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 Attributes
HcalRecHitClient Class Reference

#include <HcalRecHitClient.h>

Inheritance diagram for HcalRecHitClient:
HcalBaseDQClient

Public Member Functions

void analyze (void)
 
void beginJob (void)
 
void beginRun (void)
 
void calculateProblems (void)
 
void cleanup (void)
 
void endJob (void)
 
void endRun (void)
 
bool hasErrors_Temp (void)
 
bool hasOther_Temp (void)
 
bool hasWarnings_Temp (void)
 
 HcalRecHitClient ()
 Constructors. More...
 
 HcalRecHitClient (std::string myname)
 
 HcalRecHitClient (std::string myname, const edm::ParameterSet &ps)
 
void setup (void)
 
bool test_enabled (void)
 
void updateChannelStatus (std::map< HcalDetId, unsigned int > &myqual)
 
 ~HcalRecHitClient ()
 Destructor. More...
 
- Public Member Functions inherited from HcalBaseDQClient
void getLogicalMap (const edm::EventSetup &es)
 
 HcalBaseDQClient ()
 
 HcalBaseDQClient (std::string s, const edm::ParameterSet &ps)
 
virtual void htmlOutput (std::string htmlDir)
 
std::string name ()
 
void setEventSetup (const edm::EventSetup &es)
 
virtual void setStatusMap (std::map< HcalDetId, unsigned int > &map)
 
virtual bool validHtmlOutput ()
 
virtual ~HcalBaseDQClient (void)
 

Private Attributes

EtaPhiHists * meEnergyByDepth
 
EtaPhiHists * meEnergyThreshByDepth
 
MonitorElement * meHBEnergy_1D
 
MonitorElement * meHBEnergyRMS_1D
 
MonitorElement * meHBEnergyRMSThresh_1D
 
MonitorElement * meHBEnergyThresh_1D
 
MonitorElement * meHEEnergy_1D
 
MonitorElement * meHEEnergyRMS_1D
 
MonitorElement * meHEEnergyRMSThresh_1D
 
MonitorElement * meHEEnergyThresh_1D
 
MonitorElement * meHFEnergy_1D
 
MonitorElement * meHFEnergyRMS_1D
 
MonitorElement * meHFEnergyRMSThresh_1D
 
MonitorElement * meHFEnergyThresh_1D
 
MonitorElement * meHOEnergy_1D
 
MonitorElement * meHOEnergyRMS_1D
 
MonitorElement * meHOEnergyRMSThresh_1D
 
MonitorElement * meHOEnergyThresh_1D
 
EtaPhiHists * meSqrtSumEnergy2ByDepth
 
EtaPhiHists * meSqrtSumEnergy2ThreshByDepth
 
EtaPhiHists * meTimeByDepth
 
EtaPhiHists * meTimeThreshByDepth
 
int nevts_
 

Additional Inherited Members

- Public Attributes inherited from HcalBaseDQClient
int badChannelStatusMask_
 
std::map< HcalDetId, unsigned int > badstatusmap
 
const edm::EventSetup * c
 
bool cloneME_
 
int debug_
 
DQMStore * dqmStore_
 
bool enableCleanup_
 
bool enoughevents_
 
HcalLogicalMap * logicalMap_
 
double minerrorrate_
 
int minevents_
 
std::string name_
 
bool needLogicalMap_
 
bool Online_
 
std::string prefixME_
 
MonitorElement * ProblemCells
 
EtaPhiHists * ProblemCellsByDepth
 
std::vector< std::string > problemnames_
 
std::string subdir_
 
bool testenabled_
 
bool validHtmlOutput_
 

Detailed Description

Definition at line 10 of file HcalRecHitClient.h.

Constructor & Destructor Documentation

HcalRecHitClient::HcalRecHitClient ( )
inline

Constructors.

Definition at line 15 of file HcalRecHitClient.h.

References HcalBaseDQClient::name_.

15 {name_="";};
HcalRecHitClient::HcalRecHitClient ( std::string  myname)

Definition at line 18 of file HcalRecHitClient.cc.

References HcalBaseDQClient::name_.

19 {
20  name_=myname;
21 }
HcalRecHitClient::HcalRecHitClient ( std::string  myname,
const edm::ParameterSet &  ps 
)

Definition at line 23 of file HcalRecHitClient.cc.

References HcalBaseDQClient::badChannelStatusMask_, HcalBaseDQClient::cloneME_, HcalBaseDQClient::debug_, HcalBaseDQClient::enableCleanup_, HcalBaseDQClient::enoughevents_, edm::ParameterSet::getUntrackedParameter(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name_, HcalBaseDQClient::Online_, HcalBaseDQClient::prefixME_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and HcalBaseDQClient::validHtmlOutput_.

24 {
25  name_=myname;
26  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
27  debug_ = ps.getUntrackedParameter<int>("debug",0);
28  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
29  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
30  prefixME_.append("/");
31  subdir_ = ps.getUntrackedParameter<std::string>("RecHitFolder","RecHitMonitor_Hcal/"); // RecHitMonitor_Hcal
32  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
33  subdir_.append("/");
34  subdir_=prefixME_+subdir_;
35 
36  validHtmlOutput_ = ps.getUntrackedParameter<bool>("RecHit_validHtmlOutput",true);
37  cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
38  badChannelStatusMask_ = ps.getUntrackedParameter<int>("RecHit_BadChannelStatusMask",
39  ps.getUntrackedParameter<int>("BadChannelStatusMask",
40  0)); // identify channel status values to mask
41 
42  minerrorrate_ = ps.getUntrackedParameter<double>("RecHit_minerrorrate",
43  ps.getUntrackedParameter<double>("minerrorrate",0.01));
44  minevents_ = ps.getUntrackedParameter<int>("RecHit_minevents",
45  ps.getUntrackedParameter<int>("minevents",1));
46  enoughevents_=false;
47  Online_ = ps.getUntrackedParameter<bool>("online",false);
48 
49  ProblemCells=0;
50  ProblemCellsByDepth=0;
51 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalBaseDQClient::ProblemCells
MonitorElement * ProblemCells
Definition: HcalBaseDQClient.h:75
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HcalBaseDQClient::ProblemCellsByDepth
EtaPhiHists * ProblemCellsByDepth
Definition: HcalBaseDQClient.h:76
HcalBaseDQClient::subdir_
std::string subdir_
Definition: HcalBaseDQClient.h:62
HcalBaseDQClient::prefixME_
std::string prefixME_
Definition: HcalBaseDQClient.h:61
HcalRecHitClient::~HcalRecHitClient ( )

Destructor.

Definition at line 502 of file HcalRecHitClient.cc.

503 {}

Member Function Documentation

void HcalRecHitClient::analyze ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 53 of file HcalRecHitClient.cc.

References HcalObjRepresent::CalcIeta(), calculateProblems(), HcalBaseDQClient::cloneME_, gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), DQMStore::get(), HcalBarrel, HcalEndcap, HcalEtaPhiHistNames(), HcalForward, HcalOuter, i, isHB(), isHE(), isHF(), isHO(), meEnergyByDepth, meEnergyThreshByDepth, meHBEnergy_1D, meHBEnergyRMS_1D, meHBEnergyRMSThresh_1D, meHBEnergyThresh_1D, meHEEnergy_1D, meHEEnergyRMS_1D, meHEEnergyRMSThresh_1D, meHEEnergyThresh_1D, meHFEnergy_1D, meHFEnergyRMS_1D, meHFEnergyRMSThresh_1D, meHFEnergyThresh_1D, meHOEnergy_1D, meHOEnergyRMS_1D, meHOEnergyRMSThresh_1D, meHOEnergyThresh_1D, meTimeByDepth, meTimeThreshByDepth, HcalBaseDQClient::name(), phi, funct::pow(), MonitorElement::Reset(), RMS, alignCSCRings::s, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQClient::subdir_, and validDetId().

Referenced by endRun().

54 {
55  if (debug_>2) std::cout <<"\tHcalRecHitClient::analyze()"<<std::endl;
56 
57  TH2F* OccupancyByDepth[4];
58  TH2F* SumEnergyByDepth[4];
59  TH2F* SumTimeByDepth[4];
60  TH2F* SqrtSumEnergy2ByDepth[4];
61 
62  TH2F* OccupancyThreshByDepth[4];
63  TH2F* SumEnergyThreshByDepth[4];
64  TH2F* SumTimeThreshByDepth[4];
65  TH2F* SqrtSumEnergy2ThreshByDepth[4];
66 
67  std::vector<std::string> name = HcalEtaPhiHistNames();
68 
69  MonitorElement* me;
70  bool gotHistos=true;
71 
72  for (int i=0;i<4;++i)
73  {
74  std::string s=subdir_+"Distributions_AllRecHits/"+name[i]+"RecHit Occupancy";
75  me=dqmStore_->get(s.c_str());
76  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
77  OccupancyByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, OccupancyByDepth[i], debug_);
78  s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Summed Energy GeV";
79  me=dqmStore_->get(s.c_str());
80  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
81  SumEnergyByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumEnergyByDepth[i], debug_);
82  s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Summed Time nS";
83  me=dqmStore_->get(s.c_str());
84  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
85  SumTimeByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumTimeByDepth[i], debug_);
86  s=subdir_+"Distributions_AllRecHits/sumplots/"+name[i]+"RecHit Sqrt Summed Energy2 GeV";
87  me=dqmStore_->get(s.c_str());
88  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
89  SqrtSumEnergy2ByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SqrtSumEnergy2ByDepth[i], debug_);
90 
91  // Threshold histograms
92  s=subdir_+"Distributions_PassedMinBias/"+name[i]+"Above Threshold RecHit Occupancy";
93  me=dqmStore_->get(s.c_str());
94  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
95  OccupancyThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, OccupancyThreshByDepth[i], debug_);
96  s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Summed Energy GeV";
97  me=dqmStore_->get(s.c_str());
98  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
99  SumEnergyThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumEnergyThreshByDepth[i], debug_);
100  s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Summed Time nS";
101  me=dqmStore_->get(s.c_str());
102  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
103  SumTimeThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SumTimeThreshByDepth[i], debug_);
104  s=subdir_+"Distributions_PassedMinBias/sumplots/"+name[i]+"Above Threshold RecHit Sqrt Summed Energy2 GeV";
105  me=dqmStore_->get(s.c_str());
106  if (me==0) {if (debug_>0) std::cout <<"Could not get histogram "<<s<<std::endl; gotHistos=false; break;}
107  SqrtSumEnergy2ThreshByDepth[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, SqrtSumEnergy2ThreshByDepth[i], debug_);
108  }
109  if (gotHistos==false)
110  {
111  if (debug_>0) std::cout <<"<HcalRecHitClient::calculateProblems()> Not all histograms could be found; skipping normalization"<<std::endl;
112  return;
113  }
114 
115  // Clear histograms before re-filling
116  meHBEnergy_1D->Reset();
117  meHBEnergyRMS_1D->Reset();
118  meHEEnergy_1D->Reset();
119  meHEEnergyRMS_1D->Reset();
120  meHOEnergy_1D->Reset();
121  meHOEnergyRMS_1D->Reset();
122  meHFEnergy_1D->Reset();
123  meHFEnergyRMS_1D->Reset();
124  meHBEnergyThresh_1D->Reset();
125  meHBEnergyRMSThresh_1D->Reset();
126  meHEEnergyThresh_1D->Reset();
127  meHEEnergyRMSThresh_1D->Reset();
128  meHOEnergyThresh_1D->Reset();
129  meHOEnergyRMSThresh_1D->Reset();
130  meHFEnergyThresh_1D->Reset();
131  meHFEnergyRMSThresh_1D->Reset();
132 
133  for (int mydepth=0;mydepth<4;++mydepth)
134  {
135  for (int eta=0;eta<OccupancyByDepth[mydepth]->GetNbinsX();++eta)
136  {
137  // eta+1=1: ieta = -42
138  // eta+1=13: ieta = -29
139 
140  for (int phi=0;phi<72;++phi)
141  {
142  if (OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)>0)
143  {
144  // fill 1D plots
145  if (isHB(eta,mydepth+1))
146  {
147  if (validDetId(HcalBarrel, CalcIeta(HcalBarrel, eta, mydepth+1), phi+1, mydepth+1))
148  {
149  meHBEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
150  meHBEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
151  }
152  }
153  else if (isHE(eta,mydepth+1))
154  {
155  if (validDetId(HcalEndcap, CalcIeta(HcalEndcap, eta, mydepth+1), phi+1, mydepth+1))
156  {
157 
158  meHEEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
159  meHEEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
160  }
161  }
162  else if (isHO(eta,mydepth+1))
163  {
164  if (validDetId(HcalOuter, CalcIeta(HcalOuter, eta, mydepth+1), phi+1, mydepth+1))
165  {
166  meHOEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
167  meHOEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
168  }
169  }
170  else if (isHF(eta,mydepth+1))
171  {
172  if (validDetId(HcalForward, CalcIeta(HcalForward, eta, mydepth+1), phi+1, mydepth+1))
173  {
174  meHFEnergy_1D->Fill(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
175  meHFEnergyRMS_1D->Fill(sqrt(pow(SqrtSumEnergy2ByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1)-pow(SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)));
176  }
177  }
178  // normalize 2D plots by number of events
179 
180  meEnergyByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumEnergyByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
181  meTimeByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumTimeByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyByDepth[mydepth]->GetBinContent(eta+1, phi+1));
182  } // (OccupancyByDepth[mydepth]->GetBinContent(eta+1,phi+1)>0)
183 
184  if (OccupancyThreshByDepth[mydepth]==0) continue;
185  if (OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)>0)
186  {
187  // fill 1D plots
188  if (isHB(eta,mydepth+1))
189  {
190  if (validDetId(HcalBarrel, CalcIeta(HcalBarrel, eta, mydepth+1), phi+1, mydepth+1))
191  {
192  meHBEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
193  double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
194  RMS=pow(fabs(RMS),0.5);
195  meHBEnergyRMSThresh_1D->Fill(RMS);
196  }
197  }
198  else if (isHE(eta,mydepth+1))
199  {
200  if (validDetId(HcalEndcap, CalcIeta(HcalEndcap, eta, mydepth+1), phi+1, mydepth+1))
201  {
202 
203  meHEEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
204  double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
205  RMS=pow(fabs(RMS),0.5);
206  meHEEnergyRMSThresh_1D->Fill(RMS);
207  }
208  }
209  else if (isHO(eta,mydepth+1))
210  {
211  if (validDetId(HcalOuter, CalcIeta(HcalOuter, eta, mydepth+1), phi+1, mydepth+1))
212  {
213  meHOEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
214  double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
215  RMS=pow(fabs(RMS),0.5);
216  meHOEnergyRMSThresh_1D->Fill(RMS);
217  }
218  }
219  else if (isHF(eta,mydepth+1))
220  {
221  if (validDetId(HcalForward, CalcIeta(HcalForward, eta, mydepth+1), phi+1, mydepth+1))
222  {
223  meHFEnergyThresh_1D->Fill(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
224  double RMS=pow(SqrtSumEnergy2ThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1,phi+1)-pow(SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1),2);
225  RMS=pow(fabs(RMS),0.5);
226  meHFEnergyRMSThresh_1D->Fill(RMS);
227  }
228  }
229  // fill 2D plots
230  meEnergyThreshByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumEnergyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
231  meTimeThreshByDepth->depth[mydepth]->setBinContent(eta+1, phi+1, SumTimeThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1)/OccupancyThreshByDepth[mydepth]->GetBinContent(eta+1, phi+1));
232  }
233  } // for (int phi=0;phi<72;++phi)
234  } // for (int eta=0;eta<OccupancyByDepth->..;++eta)
235  } // for (int mydepth=0;...)
236 
237  FillUnphysicalHEHFBins(*meEnergyByDepth);
238  FillUnphysicalHEHFBins(*meTimeByDepth);
239  FillUnphysicalHEHFBins(*meEnergyThreshByDepth);
240  FillUnphysicalHEHFBins(*meTimeThreshByDepth);
241 
242  calculateProblems();
243 }
isHO
bool isHO(int etabin, int depth)
Definition: HcalEtaPhiHists.h:479
HcalRecHitClient::meHFEnergyThresh_1D
MonitorElement * meHFEnergyThresh_1D
Definition: HcalRecHitClient.h:60
i
int i
Definition: DBlmapReader.cc:9
HcalBaseDQClient::name
std::string name()
Definition: HcalBaseDQClient.h:58
HcalRecHitClient::meHOEnergyRMS_1D
MonitorElement * meHOEnergyRMS_1D
Definition: HcalRecHitClient.h:54
HcalRecHitClient::meHEEnergyRMS_1D
MonitorElement * meHEEnergyRMS_1D
Definition: HcalRecHitClient.h:53
HcalRecHitClient::meHFEnergyRMS_1D
MonitorElement * meHFEnergyRMS_1D
Definition: HcalRecHitClient.h:55
HcalRecHitClient::meHEEnergyRMSThresh_1D
MonitorElement * meHEEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:63
HcalEtaPhiHistNames
std::vector< std::string > HcalEtaPhiHistNames()
Definition: HcalEtaPhiHists.h:396
HcalRecHitClient::meHBEnergyRMS_1D
MonitorElement * meHBEnergyRMS_1D
Definition: HcalRecHitClient.h:52
isHE
bool isHE(int etabin, int depth)
Definition: HcalEtaPhiHists.h:429
eta
T eta() const
Definition: Basic3DVectorLD.h:177
HcalRecHitClient::meTimeThreshByDepth
EtaPhiHists * meTimeThreshByDepth
Definition: HcalRecHitClient.h:43
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HcalRecHitClient::meHFEnergyRMSThresh_1D
MonitorElement * meHFEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:65
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
HcalRecHitClient::meHEEnergy_1D
MonitorElement * meHEEnergy_1D
Definition: HcalRecHitClient.h:48
HcalRecHitClient::meHBEnergy_1D
MonitorElement * meHBEnergy_1D
Definition: HcalRecHitClient.h:47
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalRecHitClient::meHOEnergyRMSThresh_1D
MonitorElement * meHOEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:64
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
HcalRecHitClient::meEnergyThreshByDepth
EtaPhiHists * meEnergyThreshByDepth
Definition: HcalRecHitClient.h:41
HcalRecHitClient::meHOEnergyThresh_1D
MonitorElement * meHOEnergyThresh_1D
Definition: HcalRecHitClient.h:59
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1708
isHB
bool isHB(int etabin, int depth)
Definition: HcalEtaPhiHists.h:407
HcalRecHitClient::meHOEnergy_1D
MonitorElement * meHOEnergy_1D
Definition: HcalRecHitClient.h:49
HcalRecHitClient::meHEEnergyThresh_1D
MonitorElement * meHEEnergyThresh_1D
Definition: HcalRecHitClient.h:58
isHF
bool isHF(int etabin, int depth)
Definition: HcalEtaPhiHists.h:455
HcalRecHitClient::meHFEnergy_1D
MonitorElement * meHFEnergy_1D
Definition: HcalRecHitClient.h:50
HcalRecHitClient::meEnergyByDepth
EtaPhiHists * meEnergyByDepth
Definition: HcalRecHitClient.h:40
HcalBaseDQClient::subdir_
std::string subdir_
Definition: HcalBaseDQClient.h:62
HcalRecHitClient::meHBEnergyThresh_1D
MonitorElement * meHBEnergyThresh_1D
Definition: HcalRecHitClient.h:57
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
HcalRecHitClient::calculateProblems
void calculateProblems(void)
Definition: HcalRecHitClient.cc:245
HcalRecHitClient::meTimeByDepth
EtaPhiHists * meTimeByDepth
Definition: HcalRecHitClient.h:42
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
HcalRecHitClient::meHBEnergyRMSThresh_1D
MonitorElement * meHBEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:62
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:475
HcalBaseDQClient::dqmStore_
DQMStore * dqmStore_
Definition: HcalBaseDQClient.h:81
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
validDetId
bool validDetId(HcalSubdetector sd, int ies, int ip, int dp)
Definition: HcalEtaPhiHists.h:502
phi
Definition: DDAxes.h:10
void HcalRecHitClient::beginJob ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 353 of file HcalRecHitClient.cc.

References gather_cfg::cout, HcalBaseDQClient::debug_, HcalBaseDQClient::dqmStore_, cppFunctionSkipper::operator, and DQMStore::showDirStructure().

354 {
355  dqmStore_ = edm::Service<DQMStore>().operator->();
356  if (debug_>0)
357  {
358  std::cout <<"<HcalRecHitClient::beginJob()> Displaying dqmStore directory structure:"<<std::endl;
359  dqmStore_->showDirStructure();
360  }
361 }
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:3332
HcalBaseDQClient::dqmStore_
DQMStore * dqmStore_
Definition: HcalBaseDQClient.h:81
void HcalRecHitClient::beginRun ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 364 of file HcalRecHitClient.cc.

References DQMStore::book1D(), DQMStore::book2D(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, MonitorElement::getName(), i, meEnergyByDepth, meEnergyThreshByDepth, meHBEnergy_1D, meHBEnergyRMS_1D, meHBEnergyRMSThresh_1D, meHBEnergyThresh_1D, meHEEnergy_1D, meHEEnergyRMS_1D, meHEEnergyRMSThresh_1D, meHEEnergyThresh_1D, meHFEnergy_1D, meHFEnergyRMS_1D, meHFEnergyRMSThresh_1D, meHFEnergyThresh_1D, meHOEnergy_1D, meHOEnergyRMS_1D, meHOEnergyRMSThresh_1D, meHOEnergyThresh_1D, meTimeByDepth, meTimeThreshByDepth, nevts_, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, HcalBaseDQClient::problemnames_, DQMStore::setCurrentFolder(), EtaPhiHists::setup(), HcalBaseDQClient::subdir_, x, and detailsBasic3DVector::y.

365 {
366  if (debug_>1) std::cout <<"<HcalRecHitClient::endRun>"<<std::endl;
367 
368  if (!dqmStore_)
369  {
370  if (debug_>0) std::cout <<"<HcalRecHitClient::beginRun> dqmStore does not exist!"<<std::endl;
371  return;
372  }
373  dqmStore_->setCurrentFolder(subdir_);
374  problemnames_.clear();
375  ProblemCells=dqmStore_->book2D(" ProblemRecHits",
376  "Problem RecHit Rate for all HCAL;ieta;iphi",
377  85,-42.5,42.5,
378  72,0.5,72.5);
379  problemnames_.push_back(ProblemCells->getName());
380  if (debug_>1)
381  std::cout << "Tried to create ProblemCells Monitor Element in directory "<<subdir_<<" \t Failed? "<<(ProblemCells==0)<<std::endl;
382  dqmStore_->setCurrentFolder(subdir_+"problem_rechits");
383  ProblemCellsByDepth=new EtaPhiHists();
384  ProblemCellsByDepth->setup(dqmStore_," Problem RecHit Rate");
385  for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
386  problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());
387 
388  nevts_=0;
389 
390  dqmStore_->setCurrentFolder(subdir_+"Distributions_AllRecHits");
391  meEnergyByDepth = new EtaPhiHists();
392  meEnergyByDepth->setup(dqmStore_,"RecHit Average Energy","GeV");
393  meTimeByDepth = new EtaPhiHists();
394  meTimeByDepth->setup(dqmStore_,"RecHit Average Time","nS");
395  // set all average times to -1000 by default (so that they don't show up on plots
396  for (unsigned int i=0;i<meTimeByDepth->depth.size();++i)
397  {
398  (meTimeByDepth->depth[i]->getTH2F())->SetMinimum(-150);
399  (meTimeByDepth->depth[i]->getTH2F())->SetMaximum(150);
400  int etabins=(meTimeByDepth->depth[i]->getTH2F())->GetNbinsX();
401  int phibins=(meTimeByDepth->depth[i]->getTH2F())->GetNbinsY();
402  for (int x=1;x<=etabins;++x)
403  for (int y=1;y<=phibins;++y)
404  meTimeByDepth->depth[i]->setBinContent(x,y,-1000);
405  }
406 
407  dqmStore_->setCurrentFolder(subdir_+"Distributions_PassedMinBias");
408  meEnergyThreshByDepth = new EtaPhiHists();
409  meEnergyThreshByDepth->setup(dqmStore_,"Above Threshold RecHit Average Energy","GeV");
410  meTimeThreshByDepth = new EtaPhiHists();
411  meTimeThreshByDepth->setup(dqmStore_,"Above Threshold RecHit Average Time","nS");
412  // set all average times to -1000 by default (so that they don't show up on plots
413  for (unsigned int i=0;i<meTimeThreshByDepth->depth.size();++i)
414  {
415  (meTimeThreshByDepth->depth[i]->getTH2F())->SetMinimum(-150);
416  (meTimeThreshByDepth->depth[i]->getTH2F())->SetMaximum(150);
417  int etabins=(meTimeThreshByDepth->depth[i]->getTH2F())->GetNbinsX();
418  int phibins=(meTimeThreshByDepth->depth[i]->getTH2F())->GetNbinsY();
419  for (int x=1;x<=etabins;++x)
420  for (int y=1;y<=phibins;++y)
421  meTimeThreshByDepth->depth[i]->setBinContent(x,y,-1000);
422  }
423 
424  dqmStore_->setCurrentFolder(subdir_+"Distributions_AllRecHits/rechit_1D_plots/");
425  meHBEnergy_1D=dqmStore_->book1D("HB_energy_1D","HB Average Energy Per RecHit;Energy (GeV)",200,-5,15);
426  meHEEnergy_1D=dqmStore_->book1D("HE_energy_1D","HE Average Energy Per RecHit;Energy (GeV)",200,-5,15);
427  meHOEnergy_1D=dqmStore_->book1D("HO_energy_1D","HO Average Energy Per RecHit;Energy (GeV)",200,-10,20);
428  meHFEnergy_1D=dqmStore_->book1D("HF_energy_1D","HF Average Energy Per RecHit;Energy (GeV)",200,-5,15);
429 
430  meHBEnergyRMS_1D=dqmStore_->book1D("HB_energy_RMS_1D","HB Energy RMS Per RecHit;Energy (GeV)",250,0,5);
431  meHEEnergyRMS_1D=dqmStore_->book1D("HE_energy_RMS_1D","HE Energy RMS Per RecHit;Energy (GeV)",250,0,5);
432  meHOEnergyRMS_1D=dqmStore_->book1D("HO_energy_RMS_1D","HO Energy RMS Per RecHit;Energy (GeV)",250,0,5);
433  meHFEnergyRMS_1D=dqmStore_->book1D("HF_energy_RMS_1D","HF Energy RMS Per RecHit;Energy (GeV)",250,0,5);
434 
435  dqmStore_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/rechit_1D_plots/");
436  meHBEnergyThresh_1D=dqmStore_->book1D("HB_energyThresh_1D","HB Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,35);
437  meHEEnergyThresh_1D=dqmStore_->book1D("HE_energyThresh_1D","HE Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,35);
438  meHOEnergyThresh_1D=dqmStore_->book1D("HO_energyThresh_1D","HO Average Energy Per RecHit Above Threshold;Energy (GeV)",300,-10,50);
439  meHFEnergyThresh_1D=dqmStore_->book1D("HF_energyThresh_1D","HF Average Energy Per RecHit Above Threshold;Energy (GeV)",200,-5,95);
440 
441  meHBEnergyRMSThresh_1D=dqmStore_->book1D("HB_energy_RMSThresh_1D","HB Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
442  meHEEnergyRMSThresh_1D=dqmStore_->book1D("HE_energy_RMSThresh_1D","HE Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
443  meHOEnergyRMSThresh_1D=dqmStore_->book1D("HO_energy_RMSThresh_1D","HO Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,10);
444  meHFEnergyRMSThresh_1D=dqmStore_->book1D("HF_energy_RMSThresh_1D","HF Energy RMS Per RecHit Above Threshold;Energy (GeV)",200,0,20);
445 }
HcalRecHitClient::meHFEnergyThresh_1D
MonitorElement * meHFEnergyThresh_1D
Definition: HcalRecHitClient.h:60
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:113
i
int i
Definition: DBlmapReader.cc:9
HcalBaseDQClient::ProblemCells
MonitorElement * ProblemCells
Definition: HcalBaseDQClient.h:75
HcalRecHitClient::meHOEnergyRMS_1D
MonitorElement * meHOEnergyRMS_1D
Definition: HcalRecHitClient.h:54
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:954
HcalRecHitClient::meHEEnergyRMS_1D
MonitorElement * meHEEnergyRMS_1D
Definition: HcalRecHitClient.h:53
HcalRecHitClient::meHFEnergyRMS_1D
MonitorElement * meHFEnergyRMS_1D
Definition: HcalRecHitClient.h:55
HcalRecHitClient::meHEEnergyRMSThresh_1D
MonitorElement * meHEEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:63
EtaPhiHists::setup
void setup(DQMStore *&m_dbe, std::string Name, std::string Units="")
Definition: HcalEtaPhiHists.h:24
HcalRecHitClient::meHBEnergyRMS_1D
MonitorElement * meHBEnergyRMS_1D
Definition: HcalRecHitClient.h:52
HcalRecHitClient::meTimeThreshByDepth
EtaPhiHists * meTimeThreshByDepth
Definition: HcalRecHitClient.h:43
HcalRecHitClient::meHFEnergyRMSThresh_1D
MonitorElement * meHFEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:65
HcalRecHitClient::meHEEnergy_1D
MonitorElement * meHEEnergy_1D
Definition: HcalRecHitClient.h:48
HcalRecHitClient::meHBEnergy_1D
MonitorElement * meHBEnergy_1D
Definition: HcalRecHitClient.h:47
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalRecHitClient::meHOEnergyRMSThresh_1D
MonitorElement * meHOEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:64
HcalBaseDQClient::problemnames_
std::vector< std::string > problemnames_
Definition: HcalBaseDQClient.h:78
HcalBaseDQClient::ProblemCellsByDepth
EtaPhiHists * ProblemCellsByDepth
Definition: HcalBaseDQClient.h:76
HcalRecHitClient::meEnergyThreshByDepth
EtaPhiHists * meEnergyThreshByDepth
Definition: HcalRecHitClient.h:41
HcalRecHitClient::meHOEnergyThresh_1D
MonitorElement * meHOEnergyThresh_1D
Definition: HcalRecHitClient.h:59
HcalRecHitClient::meHOEnergy_1D
MonitorElement * meHOEnergy_1D
Definition: HcalRecHitClient.h:49
HcalRecHitClient::meHEEnergyThresh_1D
MonitorElement * meHEEnergyThresh_1D
Definition: HcalRecHitClient.h:58
HcalRecHitClient::meHFEnergy_1D
MonitorElement * meHFEnergy_1D
Definition: HcalRecHitClient.h:50
HcalRecHitClient::meEnergyByDepth
EtaPhiHists * meEnergyByDepth
Definition: HcalRecHitClient.h:40
HcalBaseDQClient::subdir_
std::string subdir_
Definition: HcalBaseDQClient.h:62
HcalRecHitClient::meHBEnergyThresh_1D
MonitorElement * meHBEnergyThresh_1D
Definition: HcalRecHitClient.h:57
HcalRecHitClient::meTimeByDepth
EtaPhiHists * meTimeByDepth
Definition: HcalRecHitClient.h:42
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
HcalRecHitClient::meHBEnergyRMSThresh_1D
MonitorElement * meHBEnergyRMSThresh_1D
Definition: HcalRecHitClient.h:62
x
Definition: DDAxes.h:10
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:1082
HcalBaseDQClient::dqmStore_
DQMStore * dqmStore_
Definition: HcalBaseDQClient.h:81
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:667
void HcalRecHitClient::calculateProblems ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 245 of file HcalRecHitClient.cc.

References HcalBaseDQClient::badstatusmap, HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, HcalBaseDQClient::dqmStore_, HcalBaseDQClient::enoughevents_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), MonitorElement::getBinContent(), MonitorElement::getTH2F(), HcalBarrel, HcalEmpty, HcalEndcap, HcalEtaPhiHistNames(), HcalForward, HcalOuter, isHB(), isHE(), isHF(), isHO(), min(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::minevents_, HcalBaseDQClient::name(), phi, HcalBaseDQClient::ProblemCells, HcalBaseDQClient::ProblemCellsByDepth, MonitorElement::Reset(), MonitorElement::setBinContent(), and ecaldqm::zside().

Referenced by analyze().

246 {
247  if (debug_>2) std::cout <<"\t\tHcalRecHitClient::calculateProblems()"<<std::endl;
248  if(!dqmStore_) return;
249  double totalevents=0;
250  int etabins=0, phibins=0, zside=0;
251  double problemvalue=0;
252  std::vector<std::string> name = HcalEtaPhiHistNames(); // use this to get EtaPhiHistograms that feed problem calculation, once they exist (see analyze function for example usage of HcalEtaPhiHistNames())
253 
254  // Get histograms used in determining rechit problem rate,
255  // and get totalevents from their underflow bins.
256  // No such problem histograms are defined so far.
257 
258  enoughevents_=true;
259 
260  if (totalevents==0)
261  return;
262 
263  // Clear away old problems
264  if (ProblemCells!=0)
265  {
266  ProblemCells->Reset();
267  (ProblemCells->getTH2F())->SetMaximum(1.05);
268  (ProblemCells->getTH2F())->SetMinimum(0.);
269  }
270  for (unsigned int d=0;d<ProblemCellsByDepth->depth.size();++d)
271  {
272  if (ProblemCellsByDepth->depth[d]!=0)
273  {
274  ProblemCellsByDepth->depth[d]->Reset();
275  (ProblemCellsByDepth->depth[d]->getTH2F())->SetMaximum(1.05);
276  (ProblemCellsByDepth->depth[d]->getTH2F())->SetMinimum(0.);
277  }
278  }
279 
280  for (unsigned int d=0;d<ProblemCellsByDepth->depth.size();++d)
281  {
282  if (ProblemCellsByDepth->depth[d]==0) continue;
283 
284  // Get total events from some future histogram
285  //totalevents=DigiPresentByDepth[d]->GetBinContent(0); // get totalevents from each depth, in case they differ
286  if (totalevents==0 || totalevents<minevents_) continue;
287  etabins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsX();
288  phibins=(ProblemCellsByDepth->depth[d]->getTH2F())->GetNbinsY();
289  problemvalue=0;
290  for (int eta=0;eta<etabins;++eta)
291  {
292  int ieta=CalcIeta(eta,d+1);
293  if (ieta==-9999) continue;
294  for (int phi=0;phi<phibins;++phi)
295  {
296  problemvalue=0;
297  // Add histograms here when rechit testing decided upon
298  /*if (DigiPresentByDepth[d]!=0 && DigiPresentByDepth[d]->GetBinContent(eta+1,phi+1)>0)
299  problemvalue=totalevents; */
300 
301  if (problemvalue==0) continue;
302  problemvalue/=totalevents; // problem value is a rate; should be between 0 and 1
303  problemvalue = std::min(1.,problemvalue);
304 
305  zside=0;
306  if (isHF(eta,d+1)) // shift ieta by 1 for HF
307  ieta<0 ? zside = -1 : zside = 1;
308 
309  // For problem cells that exceed our allowed rate,
310  // set the values to 999 if the cells are already marked in the status database
311  if (problemvalue>minerrorrate_)
312  {
313  HcalSubdetector subdet=HcalEmpty;
314  if (isHB(eta,d+1))subdet=HcalBarrel;
315  else if (isHE(eta,d+1)) subdet=HcalEndcap;
316  else if (isHF(eta,d+1)) subdet=HcalForward;
317  else if (isHO(eta,d+1)) subdet=HcalOuter;
318  HcalDetId hcalid(subdet, ieta, phi+1, (int)(d+1));
319  if (badstatusmap.find(hcalid)!=badstatusmap.end())
320  problemvalue=999;
321  }
322 
323 
324  ProblemCellsByDepth->depth[d]->setBinContent(eta+1,phi+1,problemvalue);
325  if (ProblemCells!=0) ProblemCells->Fill(ieta+zside,phi+1,problemvalue);
326  } // loop on phi
327  } // loop on eta
328  } // loop on depth
329 
330  if (ProblemCells==0)
331  {
332  if (debug_>0) std::cout <<"<HcalRecHitClient::analyze> ProblemCells histogram does not exist!"<<std::endl;
333  return;
334  }
335 
336  // Normalization of ProblemCell plot, in the case where there are errors in multiple depths
337  etabins=(ProblemCells->getTH2F())->GetNbinsX();
338  phibins=(ProblemCells->getTH2F())->GetNbinsY();
339  for (int eta=0;eta<etabins;++eta)
340  {
341  for (int phi=0;phi<phibins;++phi)
342  {
343  if (ProblemCells->getBinContent(eta+1,phi+1)>1. && ProblemCells->getBinContent(eta+1,phi+1)<999)
344  ProblemCells->setBinContent(eta+1,phi+1,1.);
345  }
346  }
347  FillUnphysicalHEHFBins(*ProblemCellsByDepth);
348  FillUnphysicalHEHFBins(ProblemCells);
349  return;
350 
351 }
isHO
bool isHO(int etabin, int depth)
Definition: HcalEtaPhiHists.h:479
HcalBaseDQClient::ProblemCells
MonitorElement * ProblemCells
Definition: HcalBaseDQClient.h:75
HcalBaseDQClient::name
std::string name()
Definition: HcalBaseDQClient.h:58
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:866
HcalEtaPhiHistNames
std::vector< std::string > HcalEtaPhiHistNames()
Definition: HcalEtaPhiHists.h:396
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:210
isHE
bool isHE(int etabin, int depth)
Definition: HcalEtaPhiHists.h:429
eta
T eta() const
Definition: Basic3DVectorLD.h:177
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalBaseDQClient::ProblemCellsByDepth
EtaPhiHists * ProblemCellsByDepth
Definition: HcalBaseDQClient.h:76
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
min
T min(T a, T b)
Definition: MathUtil.h:58
isHB
bool isHB(int etabin, int depth)
Definition: HcalEtaPhiHists.h:407
isHF
bool isHF(int etabin, int depth)
Definition: HcalEtaPhiHists.h:455
HcalBaseDQClient::badstatusmap
std::map< HcalDetId, unsigned int > badstatusmap
Definition: HcalBaseDQClient.h:80
HcalObjRepresent::FillUnphysicalHEHFBins
void FillUnphysicalHEHFBins(std::vector< TH2F > &hh)
Definition: HcalObjRepresent.h:528
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:765
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1563
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:475
HcalBaseDQClient::dqmStore_
DQMStore * dqmStore_
Definition: HcalBaseDQClient.h:81
phi
Definition: DDAxes.h:10
void HcalRecHitClient::cleanup ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 454 of file HcalRecHitClient.cc.

454 {}
void HcalRecHitClient::endJob ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 362 of file HcalRecHitClient.cc.

362 {}
void HcalRecHitClient::endRun ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 447 of file HcalRecHitClient.cc.

References analyze(), gather_cfg::cout, and HcalBaseDQClient::debug_.

448 {
449  if (debug_>1) std::cout <<"<HcalRecHitClient::endRun>"<<std::endl;
450  analyze();
451 }
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
bool HcalRecHitClient::hasErrors_Temp ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 456 of file HcalRecHitClient.cc.

References HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQClient::debug_, EtaPhiHists::depth, MonitorElement::getTH2F(), HcalBaseDQClient::minerrorrate_, HcalBaseDQClient::ProblemCells, and HcalBaseDQClient::ProblemCellsByDepth.

457 {
458  if (!ProblemCells)
459  {
460  if (debug_>1) std::cout <<"<HcalRecHitClient::hasErrors_Temp> ProblemCells histogram does not exist!"<<std::endl;
461  return false;
462  }
463  int problemcount=0;
464  int ieta=-9999;
465 
466  for (int depth=0;depth<4; ++depth)
467  {
468  int etabins = (ProblemCells->getTH2F())->GetNbinsX();
469  int phibins = (ProblemCells->getTH2F())->GetNbinsY();
470  for (int hist_eta=0;hist_eta<etabins;++hist_eta)
471  {
472  for (int hist_phi=0; hist_phi<phibins;++hist_phi)
473  {
474  ieta=CalcIeta(hist_eta,depth+1);
475  if (ieta==-9999) continue;
476  if (ProblemCellsByDepth->depth[depth]==0)
477  continue;
478  if (ProblemCellsByDepth->depth[depth]->getBinContent(hist_eta,hist_phi)>minerrorrate_)
479  ++problemcount;
480  } // for (int hist_phi=1;...)
481  } // for (int hist_eta=1;...)
482  } // for (int depth=0;...)
483 
484  if (problemcount>0) return true;
485  return false;
486 }
HcalBaseDQClient::ProblemCells
MonitorElement * ProblemCells
Definition: HcalBaseDQClient.h:75
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::CalcIeta
int CalcIeta(int subdet, int eta, int depth)
Definition: HcalObjRepresent.h:200
HcalBaseDQClient::ProblemCellsByDepth
EtaPhiHists * ProblemCellsByDepth
Definition: HcalBaseDQClient.h:76
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1563
bool HcalRecHitClient::hasOther_Temp ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 489 of file HcalRecHitClient.cc.

489 {return false;}
bool HcalRecHitClient::hasWarnings_Temp ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 488 of file HcalRecHitClient.cc.

488 {return false;}
void HcalRecHitClient::setup ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 453 of file HcalRecHitClient.cc.

453 {}
bool HcalRecHitClient::test_enabled ( void  )
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 490 of file HcalRecHitClient.cc.

490 {return true;}
void HcalRecHitClient::updateChannelStatus ( std::map< HcalDetId, unsigned int > &  myqual)
virtual

Reimplemented from HcalBaseDQClient.

Definition at line 494 of file HcalRecHitClient.cc.

495 {
496  // rechit quality not used to update channel status yet; see dead cell client for example
497 
498 
499 
500 } //void HcalRecHitClient::updateChannelStatus

Member Data Documentation

EtaPhiHists* HcalRecHitClient::meEnergyByDepth
private

Definition at line 40 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

EtaPhiHists* HcalRecHitClient::meEnergyThreshByDepth
private

Definition at line 41 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHBEnergy_1D
private

Definition at line 47 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHBEnergyRMS_1D
private

Definition at line 52 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHBEnergyRMSThresh_1D
private

Definition at line 62 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHBEnergyThresh_1D
private

Definition at line 57 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHEEnergy_1D
private

Definition at line 48 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHEEnergyRMS_1D
private

Definition at line 53 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHEEnergyRMSThresh_1D
private

Definition at line 63 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHEEnergyThresh_1D
private

Definition at line 58 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHFEnergy_1D
private

Definition at line 50 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHFEnergyRMS_1D
private

Definition at line 55 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHFEnergyRMSThresh_1D
private

Definition at line 65 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHFEnergyThresh_1D
private

Definition at line 60 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHOEnergy_1D
private

Definition at line 49 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHOEnergyRMS_1D
private

Definition at line 54 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHOEnergyRMSThresh_1D
private

Definition at line 64 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

MonitorElement* HcalRecHitClient::meHOEnergyThresh_1D
private

Definition at line 59 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

EtaPhiHists* HcalRecHitClient::meSqrtSumEnergy2ByDepth
private

Definition at line 44 of file HcalRecHitClient.h.

EtaPhiHists* HcalRecHitClient::meSqrtSumEnergy2ThreshByDepth
private

Definition at line 45 of file HcalRecHitClient.h.

EtaPhiHists* HcalRecHitClient::meTimeByDepth
private

Definition at line 42 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

EtaPhiHists* HcalRecHitClient::meTimeThreshByDepth
private

Definition at line 43 of file HcalRecHitClient.h.

Referenced by analyze(), and beginRun().

int HcalRecHitClient::nevts_
private

Definition at line 38 of file HcalRecHitClient.h.

Referenced by beginRun().

Generated for CMSSW Reference Manual by doxygen 1.8.5