12 if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!=
"/")
13 prefixME_.append(
"/");
73 std::cout <<
"<HcalDeadCellMonitor::setup> Setting up histograms"<<std::endl;
82 Nevents =
dbe_->
book1D(
"NumberOfDeadCellEvents",
"Number of Events Seen by DeadCellMonitor",2,0,2);
87 "Total Number of Dead Hcal Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
91 "Total Number of Dead HB Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
95 "Total Number of Dead HE Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
98 "Total Number of Dead HO Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
101 "Total Number of Dead HF Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
104 "Total Number of Dead HBHEHF Cells (excluding known problems) vs LS;Lumi Section;Dead Cells",
117 me=
dbe_->
bookInt(
"Test_DigiMissing_Periodic_Lumi_Check");
122 me=
dbe_->
bookInt(
"Min_Events_Required_Periodic_Lumi_Check");
134 me=
dbe_->
bookInt(
"Test_RecHitsMissing_Periodic_Lumi_Check");
140 std::stringstream units;
141 std::stringstream
name;
147 "Digi Present At Least Once",
"");
150 "Total Number of Never-Present Hcal Cells vs LS;Lumi Section;Dead Cells",
154 "Total Number of Never-Present HB Cells vs LS;Lumi Section;Dead Cells",
158 "Total Number of Never-Present HE Cells vs LS;Lumi Section;Dead Cells",
162 "Total Number of Never-Present HO Cells vs LS;Lumi Section;Dead Cells",
166 "Total Number of Never-Present HF Cells vs LS;Lumi Section;Dead Cells",
180 name<<
"Dead Cells with No Digis";
185 name<<
"HB HE HF Depth 1 Dead Cells with No Digis for at least 1 Full Luminosity Block";
189 name<<
"HB HE HF Depth 2 Dead Cells with No Digis for at least 1 Full Luminosity Block";
193 name<<
"HE Depth 3 Dead Cells with No Digis for at least 1 Full Luminosity Block";
197 name<<
"HO Depth 4 Dead Cells with No Digis for at least 1 Full Luminosity Block";
202 name<<
"Total Number of Hcal Digis Unoccupied for at least 1 Full Luminosity Block";
207 name<<
"Total Number of HB Digis Unoccupied for at least 1 Full LS vs LS;Lumi Section; Dead Cells";
212 name<<
"Total Number of HE Digis Unoccupied for at least 1 Full LS vs LS;Lumi Section; Dead Cells";
217 name<<
"Total Number of HO Digis Unoccupied for at least 1 Full LS vs LS;Lumi Section; Dead Cells";
222 name<<
"Total Number of HF Digis Unoccupied for at least 1 Full LS vs LS;Lumi Section; Dead Cells";
247 units<<
"Cells Above Energy Threshold At Least Once: Depth 3 -- HE >="<<
HEenergyThreshold_<<
" GeV";
250 units<<
"Cells Above Energy Threshold At Least Once: Depth 4 -- HO >="<<
HOenergyThreshold_<<
" GeV";
256 "Total Number of Hcal Rechits with Low Energy;Lumi Section;Dead Cells",
259 name<<
"Total Number of HB RecHits with Energy Never >= "<<
HBenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
264 name<<
"Total Number of HE RecHits with Energy Never >= "<<HEenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
269 name<<
"Total Number of HO RecHits with Energy Never >= "<<HOenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
274 name<<
"Total Number of HF RecHits with Energy Never >= "<<HFenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
288 units<<
"RecHits with Consistent Low Energy Depth 1 -- HB <"<<HBenergyThreshold_<<
" GeV, HE < "<<HEenergyThreshold_<<
", HF <"<<HFenergyThreshold_<<
" GeV";
291 units<<
"RecHits with Consistent Low Energy Depth 2 -- HB <"<<HBenergyThreshold_<<
" GeV, HE < "<<HEenergyThreshold_<<
", HF <"<<HFenergyThreshold_<<
" GeV";
294 units<<
"RecHits with Consistent Low Energy Depth 3 -- HE <"<<HEenergyThreshold_<<
" GeV";
297 units<<
"RecHits with Consistent Low Energy Depth 4 -- HO <"<<HOenergyThreshold_<<
" GeV";
303 name<<
"Total Number of Hcal RecHits with Consistent Low Energy;Lumi Section;Dead Cells";
308 name<<
"Total Number of HB RecHits with Consistent Low Energy < "<<HBenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
313 name<<
"Total Number of HE RecHits with Consistent Low Energy < "<<HEenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
318 name<<
"Total Number of HO RecHits with Consistent Low Energy < "<<HOenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
323 name<<
"Total Number of HF RecHits with Consistent Low Energy < "<<HFenergyThreshold_<<
" GeV;Lumi Section;Dead Cells";
367 for (std::vector<DetId>::const_iterator
i = mydetids.begin();
398 for (
int ieta=11;ieta<=15;++ieta)
399 for (
int iphi=1;iphi<=72;++iphi)
402 if (
isSiPM(ieta,iphi,4)==
false)
434 for (
int ieta=11;ieta<=15;++ieta)
435 for (
int iphi=1;iphi<=72;++iphi)
573 if (
debug_>1)
std::cout <<
"\t<HcalDeadCellMonitor::analyze> Processing good event! event # = "<<
ievt_<<std::endl;
580 processEvent(*hbhe_rechit, *ho_rechit, *hf_rechit, *hbhe_digi, *ho_digi, *hf_digi);
594 if (
debug_>1)
std::cout <<
"<HcalDeadCellMonitor::processEvent> Processing event..."<<std::endl;
622 j!=hbhedigi.
end(); ++
j)
629 j!=hodigi.
end(); ++
j)
635 j!=hfdigi.
end(); ++
j)
656 j!=hbHits.
end(); ++
j)
660 k!=hoHits.
end(); ++
k)
664 j!=hfHits.
end(); ++
j)
677 for (
int d=0;d<4;++d)
684 if (d==3) ++hopresent;
685 else if (d==2) ++hepresent;
690 if (
abs(ieta)>29) ++hfpresent;
691 else if (
abs(ieta)<17) ++hbpresent;
697 if (
abs(ieta)>29) ++hfpresent;
698 else if (
abs(ieta)<17) ++hbpresent;
728 int ieta=digi.id().ieta();
729 int iphi=digi.id().iphi();
730 int depth=digi.id().depth();
749 template<
class RECHIT>
752 float en = rechit->energy();
754 int ieta =
id.ieta();
755 int iphi =
id.iphi();
756 int depth =
id.depth();
809 std::cout <<
"<HcalDeadCellMonitor::fillNevents_recentdigis> CHECKING FOR RECENT MISSING DIGIS evtcount = "<<
deadevt_<<std::endl;
824 for (
int subdet=1;subdet<=4;++subdet)
827 if (ieta==-9999)
continue;
844 ieta<0 ? zside=-1 : zside=+1;
850 std::cout <<
"DEAD CELL; NO RECENT OCCUPANCY: subdet = "<<subdet<<
", ieta = "<<ieta<<
", iphi = "<<iphi<<
" depth = "<<depth+1<<std::endl;
851 std::cout <<
"\t RAW COORDINATES: eta = "<<
eta<<
" phi = "<<
phi<<
" depth = "<<depth<<std::endl;
892 std::cout <<
"<HcalDeadCellMonitor::fillNevents_energy> BELOW-ENERGY-THRESHOLD PLOTS"<<std::endl;
906 for (
int subdet=1;subdet<=4;++subdet)
909 if (ieta==-9999)
continue;
928 ieta<0 ? zside=-1 : zside=+1;
932 std::cout <<
"DEAD CELL; BELOW ENERGY THRESHOLD; subdet = "<<subdet<<
" ieta = "<<ieta<<
", phi = "<<iphi<<
" depth = "<<depth+1<<std::endl;
954 std::cout <<
"<HcalDeadCellMonitor::fillNevents_problemCells> FILLING PROBLEM CELL PLOTS"<<std::endl;
973 int knownBadHFLUMI=0;
978 unsigned int neverpresentHB=0;
979 unsigned int neverpresentHE=0;
980 unsigned int neverpresentHO=0;
981 unsigned int neverpresentHF=0;
983 unsigned int unoccupiedHB=0;
984 unsigned int unoccupiedHE=0;
985 unsigned int unoccupiedHO=0;
986 unsigned int unoccupiedHF=0;
988 unsigned int belowenergyHB=0;
989 unsigned int belowenergyHE=0;
990 unsigned int belowenergyHO=0;
991 unsigned int belowenergyHF=0;
993 unsigned int energyneverpresentHB=0;
994 unsigned int energyneverpresentHE=0;
995 unsigned int energyneverpresentHO=0;
996 unsigned int energyneverpresentHF=0;
1018 for (
int subdet=1;subdet<=4;++subdet)
1021 if (ieta==-9999)
continue;
1073 if (
abs(ieta)<5) ++knownBadHO0;
1074 else ++knownBadHO12;
1086 if (depth==1 && (
abs(ieta)==33 ||
abs(ieta)==34))
1088 else if (depth==2 && (
abs(ieta)==35 ||
abs(ieta)==36))
1093 if (depth==1 && (
abs(ieta)==33 ||
abs(ieta)==34))
1095 else if (depth==2 && (
abs(ieta)==35 ||
abs(ieta)==36))
1103 else if (subdet==
HcalEndcap) ++neverpresentHE;
1130 if (subdet==
HcalBarrel) ++energyneverpresentHB;
1131 else if (subdet==
HcalEndcap) ++energyneverpresentHE;
1134 ++energyneverpresentHO;
1136 --energyneverpresentHO;
1138 else if (subdet==
HcalForward) ++energyneverpresentHF;
1143 else if (subdet==
HcalEndcap) ++belowenergyHE;
1211 for (
unsigned int i=0;
i<85;++
i)
1213 for (
unsigned int j=0;
j<72;++
j)
1215 for (
unsigned int k=0;
k<4;++
k)
edm::InputTag hfRechitLabel_
void beginRun(const edm::Run &run, const edm::EventSetup &c)
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement * NumberOfNeverPresentRecHitsHO
bool LumiInOrder(int lumisec)
void fillNevents_problemCells()
MonitorElement * NumberOfRecentMissingRecHitsHF
MonitorElement * NumberOfNeverPresentRecHits
void setBinContent(int binx, double content)
set content of bin (1-D)
MonitorElement * ProblemsCurrentLB
MonitorElement * NumberOfNeverPresentDigisHB
void analyze(edm::Event const &e, edm::EventSetup const &s)
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
MonitorElement * HEDeadVsEvent
MonitorElement * NumberOfRecentMissingRecHits
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
void zeroCounters(bool resetpresent=false)
bool present_digi[85][72][4]
edm::InputTag hbheRechitLabel_
#define DEFINE_FWK_MODULE(type)
std::map< unsigned int, int > KnownBadCells_
MonitorElement * NumberOfNeverPresentDigisHE
std::vector< int > AllowedCalibTypes_
void endRun(const edm::Run &run, const edm::EventSetup &c)
std::vector< T >::const_iterator const_iterator
HcalDeadCellMonitor(const edm::ParameterSet &ps)
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)
edm::LuminosityBlockNumber_t luminosityBlock() const
MonitorElement * NumberOfRecentMissingDigisHB
void update(void)
Mark the object updated.
void process_RecHit(T &rechit)
void fillNevents_recentrechits()
MonitorElement * NumberOfRecentMissingDigisHE
bool present_rechit[85][72][4]
MonitorElement * bookFloat(const char *name)
Book float.
LuminosityBlockNumber_t luminosityBlock() const
uint32_t rawId() const
get the raw id
double HFenergyThreshold_
bool isSiPM(int ieta, int iphi, int depth)
EtaPhiHists RecHitPresentByDepth
std::vector< MonitorElement * > depth
virtual void beginRun(const edm::Run &run, const edm::EventSetup &c)
MonitorElement * ProblemsVsLB
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
bool IsAllowedCalibType()
MonitorElement * NumberOfNeverPresentDigisHF
MonitorElement * ProblemsVsLB_HBHEHF
edm::InputTag hoRechitLabel_
double HEenergyThreshold_
std::vector< DetId > getAllChannels() const
MonitorElement * NumberOfNeverPresentDigisHO
void processEvent_HBHEdigi(const HBHEDataFrame digi)
MonitorElement * NumberOfNeverPresentRecHitsHB
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
unsigned int recentoccupancy_rechit[85][72][4]
MonitorElement * NumberOfRecentMissingRecHitsHE
MonitorElement * ProblemsVsLB_HF
MonitorElement * NumberOfRecentMissingDigisHO
MonitorElement * NumberOfNeverPresentRecHitsHF
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
int CalcIeta(int subdet, int eta, int depth)
MonitorElement * ProblemsVsLB_HB
const_iterator end() const
unsigned int recentoccupancy_digi[85][72][4]
void processEvent(const HBHERecHitCollection &hbHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, const HBHEDigiCollection &hbhedigi, const HODigiCollection &hodigi, const HFDigiCollection &hfdigi)
int badChannelStatusMask_
EtaPhiHists RecentMissingDigisByDepth
MonitorElement * NumberOfRecentMissingRecHitsHO
T const * product() const
EtaPhiHists DigiPresentByDepth
void FillUnphysicalHEHFBins(EtaPhiHists &hh)
EtaPhiHists RecentMissingRecHitsByDepth
int CalcEtaBin(int subdet, int ieta, int depth)
void SetupEtaPhiHists(EtaPhiHists &hh, std::string Name, std::string Units)
double HOenergyThreshold_
MonitorElement * HBDeadVsEvent
MonitorElement * NumberOfNeverPresentRecHitsHE
MonitorElement * ProblemsVsLB_HE
TProfile * getTProfile(void) const
MonitorElement * NumberOfNeverPresentDigis
double HBenergyThreshold_
MonitorElement * NumberOfRecentMissingDigis
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
MonitorElement * NumberOfRecentMissingDigisHF
void process_Digi(T &digi)
void fillNevents_recentdigis()
bool deadmon_test_rechits_
MonitorElement * bookInt(const char *name)
Book int.
const Item * getValues(DetId fId) const
MonitorElement * HFDeadVsEvent
void Reset(void)
reset ME (ie. contents, errors, etc)
void setCurrentFolder(const std::string &fullpath)
MonitorElement * HODeadVsEvent
const_iterator begin() const
bool validDetId(HcalSubdetector sd, int ies, int ip, int dp)
MonitorElement * NumberOfRecentMissingRecHitsHB
MonitorElement * ProblemsVsLB_HO