1 #include "../interface/SummaryClient.h"
13 towerBadFraction_(0.),
28 std::vector<std::string> sourceList(_params.
getUntrackedParameter<std::vector<std::string> >(
"activeSources"));
31 if(
std::find(sourceList.begin(), sourceList.end(),
"TriggerPrimitives") == sourceList.end())
sources_.erase(
std::string(
"TriggerPrimitives"));
40 MESet& meReportSummaryContents(
MEs_.at(
"ReportSummaryContents"));
41 MESet& meReportSummary(
MEs_.at(
"ReportSummary"));
42 MESet& meReportSummaryMap(
MEs_.at(
"ReportSummaryMap"));
44 for(
unsigned iDCC(0); iDCC <
nDCC; ++iDCC){
46 meReportSummaryContents.fill(dccid, -1.);
49 meReportSummary.fill(-1.);
51 meReportSummaryMap.reset(-1.);
58 MESet& meReportSummaryContents(
MEs_.at(
"ReportSummaryContents"));
59 MESet& meReportSummary(
MEs_.at(
"ReportSummary"));
61 for(
unsigned iDCC(0); iDCC <
nDCC; ++iDCC){
63 meReportSummaryContents.fill(dccid, -1.);
65 meReportSummary.fill(-1.);
67 MESet const& sIntegrityByLumi(
sources_.at(
"IntegrityByLumi"));
71 double integrityByLumi[
nDCC];
72 double rawDataByLumi[
nDCC];
73 for(
unsigned iDCC(0); iDCC <
nDCC; ++iDCC){
74 integrityByLumi[iDCC] = sIntegrityByLumi.
getBinContent(iDCC + 1);
75 rawDataByLumi[iDCC] = sDesyncByLumi.getBinContent(iDCC + 1) + sFEByLumi.getBinContent(iDCC + 1);
78 MESet& meQualitySummary(
MEs_.at(
"QualitySummary"));
79 MESet& meReportSummaryMap(
MEs_.at(
"ReportSummaryMap"));
85 MESet const* sTriggerPrimitives(
using_(
"TriggerPrimitives") ? &
sources_.at(
"TriggerPrimitives") : 0);
88 float totalChannels(0.);
89 float totalGood(0.), totalGoodRaw(0);
91 double dccChannels[
nDCC];
92 std::fill_n(dccChannels, nDCC, 0.);
93 double dccGood[
nDCC], dccGoodRaw[
nDCC];
94 std::fill_n(dccGood, nDCC, 0.);
95 std::fill_n(dccGoodRaw, nDCC, 0.);
97 std::map<uint32_t, int> badChannelsCount;
103 std::vector<bool> hasMismatchDCC(nDCC,
false);
104 for (
unsigned iDCC(0); iDCC <
nDCC; ++iDCC ) {
105 if ( sBXSRP.getBinContent(iDCC + 1) > 50. || sBXTCC.getBinContent(iDCC + 1) > 50. )
106 hasMismatchDCC[iDCC] =
true;
112 DetId id(qItr->getId());
113 unsigned iDCC(
dccId(
id) - 1);
117 int integrity(sIntegrity ? sIntegrity->getBinContent(
id) :
kUnknown);
118 int presample(sPresample ? sPresample->getBinContent(
id) :
kUnknown);
119 int hotcell(sHotCell ? sHotCell->getBinContent(
id) :
kUnknown);
120 int timing(sTiming ? sTiming->getBinContent(
id) :
kUnknown);
121 int trigprim(sTriggerPrimitives ? sTriggerPrimitives->getBinContent(
id) :
kUnknown);
122 int rawdata(sRawData.getBinContent(
id));
126 if(integrity ==
kBad && integrityByLumi[iDCC] == 0. && !hasMismatchDCC[iDCC]) integrity =
kGood;
139 qItr->setBinContent(integrity);
142 qItr->setBinContent(status);
153 dccChannels[iDCC] += 1.;
159 dccGoodRaw[iDCC] += 1.;
230 for(
unsigned iDCC(0); iDCC <
nDCC; ++iDCC){
231 if(dccChannels[iDCC] < 1.)
continue;
234 float frac(dccGood[iDCC] / dccChannels[iDCC]);
235 float fracRaw(dccGoodRaw[iDCC] / dccChannels[iDCC]);
236 meReportSummaryMap.setBinContent(dccid, frac);
238 meReportSummaryContents.fill(dccid, fracLS);
243 float totalGoodLS(
onlineMode_ ? totalGood : totalGoodRaw);
244 if(totalChannels > 0.) meReportSummary.fill(totalGoodLS / totalChannels);
247 if(totalChannels > 0.)
MEs_.at(
"GlobalSummary").setBinContent(1, totalGood / totalChannels);
248 MEs_.at(
"NBadFEDs").setBinContent(1, nBad);
void setParams(edm::ParameterSet const &) override
T getUntrackedParameter(std::string const &, T const &) const
#define DEFINE_ECALDQM_WORKER(TYPE)
bool using_(std::string const &_name, ProcessType _type=kJob) const
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
const_iterator & toNextChannel()
uint32_t rawId() const
get the raw id
void producePlots(ProcessType) override
std::set< std::string > qualitySummaries_
virtual double getBinContent(DetId const &, int=0) const
EcalTrigTowerDetId tower() const
get the HCAL/trigger iphi of this crystal
unsigned dccId(DetId const &)