91 LogVerbatim(
"DTDQM|DTMonitorClient|DTRunConditionVarClient") <<
"DTRunConditionVarClient: end job";
96 ibooker.
book2D(
"VDriftGlbSummary",
"# of MBs with good mean and good sigma of vDrift", 12, 1, 13, 5, -2, 3);
97 glbT0Summary = ibooker.
book2D(
"T0GlbSummary",
"# of MBs with good mean and good sigma of t0", 12, 1, 13, 5, -2, 3);
102 ibooker.
book2D(
"MeanVDriftGlbSummary",
"mean VDrift average per sector", 12, 1., 13., 5, -2., 3.);
104 ibooker.
book1D(
"VDriftMeanAllWheels",
"mean VDrift for all chambers", 60, 0.0048, 0.006);
109 ibooker.
book2D(
"SigmaVDriftGlbSummary",
"# of Chambers with good sigma VDrift", 12, 1., 13., 5, -2., 3.);
111 ibooker.
book1D(
"VDriftSigmaAllWheels",
"sigma VDrift for all chambers", 30, 0., 0.0006);
116 ibooker.
book2D(
"MeanT0GlbSummary",
"mean T0 average per sector", 12, 1., 13., 5, -2., 3.);
117 allwheelHistos[
"allMeanT0"] = ibooker.
book1D(
"T0MeanAllWheels",
"mean T0 for all chambers", 100, -25., 25.);
122 ibooker.
book2D(
"SigmaT0GlbSummary",
"# of Chambers with good sigma T0", 12, 1., 13., 5, -2., 3.);
123 allwheelHistos[
"allSigmaT0"] = ibooker.
book1D(
"T0SigmaAllWheels",
"sigma T0 for all chambers", 50, 0, 25);
125 for (
int wh = -2; wh <= 2; wh++) {
126 bookWheelHistos(ibooker,
"MeanVDrift",
"02-MeanVDrift", wh, 60, 0.0048, 0.006,
true);
127 bookWheelHistos(ibooker,
"SigmaVDrift",
"02-SigmaVDrift", wh, 30, 0., 0.0006);
128 bookWheelHistos(ibooker,
"MeanT0",
"03-MeanT0", wh, 100, -25., 25.,
false,
true);
129 bookWheelHistos(ibooker,
"SigmaT0",
"03-SigmaT0", wh, 50, 0, 25,
false,
true);
132 for (
int wheel = -2; wheel <= 2; wheel++) {
133 for (
int sec = 1; sec <= 14; sec++) {
135 if ((sec == 13 || sec == 14) &&
stat != 4)
142 if (!VDriftME || !T0ME) {
143 edm::LogWarning(
"DTRunConditionVarClient") <<
"ME not available" << std::endl;
148 float vDriftMean = VDriftME->
getMean();
150 float t0Mean = T0ME->
getMean();
153 float vDriftSigma = VDriftME->
getRMS();
154 float t0Sigma = T0ME->
getRMS();
174 float vDriftDev(0.), errvDriftDev(0.);
175 percDevVDrift(indexCh, vDriftMean, vDriftSigma, vDriftDev, errvDriftDev);
178 if (sec == 13 || sec == 14)
179 sec_ = (sec == 13) ? 4 : 10;
190 if (sec == 13 || sec == 14) {
191 float binVDriftDev = (
wheelHistos[wheel])[
"MeanVDriftSummary"]->getBinContent(sec_,
stat);
192 binVDriftDev = (fabs(binVDriftDev) > fabs(fillvDriftDev)) ? binVDriftDev : fillvDriftDev;
193 (
wheelHistos[wheel])[
"MeanVDriftSummary"]->setBinContent(sec_,
stat, binVDriftDev);
195 float binT0MeanVal = (
wheelHistos[wheel])[
"MeanT0Summary"]->getBinContent(sec_,
stat);
196 binT0MeanVal = (fabs(binT0MeanVal) > fabs(fillT0Mean)) ? binT0MeanVal : fillT0Mean;
197 (
wheelHistos[wheel])[
"MeanT0Summary"]->setBinContent(sec_,
stat, binT0MeanVal);
199 float binVDriftSigmVal = (
wheelHistos[wheel])[
"SigmaVDriftSummary"]->getBinContent(sec_,
stat);
200 binVDriftSigmVal = (binVDriftSigmVal > 0. && binVDriftSigmVal < vDriftSigmQ) ? binVDriftSigmVal : vDriftSigmQ;
201 (
wheelHistos[wheel])[
"SigmaVDriftSummary"]->setBinContent(sec_,
stat, binVDriftSigmVal);
203 float binT0SigmVal = (
wheelHistos[wheel])[
"SigmaT0Summary"]->getBinContent(sec_,
stat);
204 binT0SigmVal = (binT0SigmVal > 0. && binT0SigmVal < t0SigmQ) ? binT0SigmVal : t0SigmQ;
205 (
wheelHistos[wheel])[
"SigmaT0Summary"]->setBinContent(sec_,
stat, binT0SigmVal);
208 (
wheelHistos[wheel])[
"MeanVDriftSummary"]->setBinContent(sec_,
stat, fillvDriftDev);
209 (
wheelHistos[wheel])[
"MeanT0Summary"]->setBinContent(sec_,
stat, fillT0Mean);
210 (
wheelHistos[wheel])[
"SigmaVDriftSummary"]->setBinContent(sec_,
stat, vDriftSigmQ);
211 (
wheelHistos[wheel])[
"SigmaT0Summary"]->setBinContent(sec_,
stat, t0SigmQ);
215 if ((sec_ == 4 || sec_ == 10) &&
stat == 4)
218 if (vDriftDevQ > 0.85 && vDriftSigmQ > 0.85) {
220 summaryHistos[
"MeanVDriftGlbSummary"]->Fill(sec_, wheel, weight);
221 summaryHistos[
"SigmaVDriftGlbSummary"]->Fill(sec_, wheel, weight);
224 if (vDriftDevQ > 0.85 && vDriftSigmQ < 0.85) {
225 summaryHistos[
"MeanVDriftGlbSummary"]->Fill(sec_, wheel, weight);
227 if (vDriftDevQ < 0.85 && vDriftSigmQ > 0.85) {
228 summaryHistos[
"SigmaVDriftGlbSummary"]->Fill(sec_, wheel, weight);
232 if (t0MeanQ > 0.85 && t0SigmQ > 0.85) {
234 summaryHistos[
"MeanT0GlbSummary"]->Fill(sec_, wheel, weight);
235 summaryHistos[
"SigmaT0GlbSummary"]->Fill(sec_, wheel, weight);
237 if (t0MeanQ > 0.85 && t0SigmQ < 0.85) {
238 summaryHistos[
"MeanT0GlbSummary"]->Fill(sec_, wheel, weight);
240 if (t0MeanQ < 0.85 && t0SigmQ > 0.85) {
241 summaryHistos[
"SigmaT0GlbSummary"]->Fill(sec_, wheel, weight);
Log< level::Info, true > LogVerbatim
std::map< std::string, MonitorElement * > summaryHistos
virtual void setCurrentFolder(std::string const &fullpath)
MonitorElement * getChamberHistos(DQMStore::IGetter &, const DTChamberId &, std::string)
std::map< std::string, MonitorElement * > allwheelHistos
float varQuality(float var, float maxGood, float minBad)
MonitorElement * glbVDriftSummary
std::map< int, std::map< int, std::map< std::string, MonitorElement * > > > wheelRingHistos
virtual void setAxisRange(double xmin, double xmax, int axis=1)
set x-, y- or z-axis range (axis=1, 2, 3 respectively)
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
virtual double getEntries() const
get # of entries
void percDevVDrift(DTChamberId indexCh, float meanVD, float sigmaVD, float &devVD, float &errdevVD)
MonitorElement * glbT0Summary
virtual double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
void bookWheelHistos(DQMStore::IBooker &, std::string histoType, std::string subfolder, int wh, int nbins, float min, float max, bool isVDCorr=false, bool makeRings=false)
book the report summary
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
virtual double getRMS(int axis=1) const
get RMS of histogram along x, y or z axis (axis=1, 2, 3 respectively)
std::map< int, std::map< std::string, MonitorElement * > > wheelHistos
Log< level::Warning, false > LogWarning
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())