140 LogVerbatim (
"DTDQM|DTMonitorClient|DTRunConditionVarClient")
141 <<
"DTRunConditionVarClient: endRun";
147 for(
int wheel=-2;wheel<=2;wheel++){
148 for(
int sec=1; sec<=14; sec++) {
149 for(
int stat=1; stat<=4; stat++) {
151 if( (sec == 13 || sec == 14) && stat != 4 )
continue;
157 if (!VDriftME || !T0ME) {
158 edm::LogWarning(
"DTRunConditionVarClient") <<
"ME not available" << std::endl;
164 float vDriftMean = VDriftME->
getMean();
165 float t0Mean = T0ME->
getMean();
168 float vDriftSigma = VDriftME->
getRMS();
169 float t0Sigma = T0ME->
getRMS();
195 float vDriftDev(0.), errvDriftDev(0.);
196 percDevVDrift(indexCh, vDriftMean, vDriftSigma, vDriftDev, errvDriftDev);
199 if( sec == 13 || sec == 14 ) sec_ = ( sec == 13 ) ? 4 : 10;
210 if( sec == 13 || sec == 14 ) {
212 float binVDriftDev = (
wheelHistos[wheel])[
"MeanVDriftSummary"]->getBinContent(sec_,stat);
213 binVDriftDev = (fabs(binVDriftDev) > fabs(fillvDriftDev)) ? binVDriftDev : fillvDriftDev;
214 (
wheelHistos[wheel])[
"MeanVDriftSummary"] -> setBinContent(sec_,stat,binVDriftDev);
216 float binT0MeanVal = (
wheelHistos[wheel])[
"MeanT0Summary"] -> getBinContent(sec_,stat);
217 binT0MeanVal = (fabs(binT0MeanVal) > fabs(fillT0Mean)) ? binT0MeanVal : fillT0Mean;
218 (
wheelHistos[wheel])[
"MeanT0Summary"] -> setBinContent(sec_,stat,binT0MeanVal);
220 float binVDriftSigmVal = (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> getBinContent(sec_,stat);
221 binVDriftSigmVal = (binVDriftSigmVal > 0. && binVDriftSigmVal < vDriftSigmQ) ? binVDriftSigmVal : vDriftSigmQ;
222 (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> setBinContent(sec_,stat,binVDriftSigmVal);
224 float binT0SigmVal = (
wheelHistos[wheel])[
"SigmaT0Summary"] -> getBinContent(sec_,stat);
225 binT0SigmVal = (binT0SigmVal > 0. && binT0SigmVal < t0SigmQ) ? binT0SigmVal : t0SigmQ;
226 (
wheelHistos[wheel])[
"SigmaT0Summary"] -> setBinContent(sec_,stat,binT0SigmVal);
230 (
wheelHistos[wheel])[
"MeanVDriftSummary"] -> setBinContent(sec_,stat,fillvDriftDev);
231 (
wheelHistos[wheel])[
"MeanT0Summary"] -> setBinContent(sec_,stat,fillT0Mean);
232 (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> setBinContent(sec_,stat,vDriftSigmQ);
233 (
wheelHistos[wheel])[
"SigmaT0Summary"] -> setBinContent(sec_,stat,t0SigmQ);
238 if(( sec_ == 4 || sec_ == 10) && stat == 4) weight = 1/8.;
240 if( vDriftDevQ > 0.85 && vDriftSigmQ > 0.85 ) {
246 if( vDriftDevQ > 0.85 && vDriftSigmQ < 0.85 ) {
249 if( vDriftDevQ < 0.85 && vDriftSigmQ > 0.85 ) {
254 if( t0MeanQ > 0.85 && t0SigmQ > 0.85 ) {
259 if( t0MeanQ > 0.85 && t0SigmQ < 0.85 ) {
263 if( t0MeanQ < 0.85 && t0SigmQ > 0.85 ) {
std::map< std::string, MonitorElement * > summaryHistos
double getEntries(void) const
get # of entries
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
float varQuality(float var, float maxGood, float minBad)
MonitorElement * glbVDriftSummary
MonitorElement * getChamberHistos(const DTChamberId &, std::string)
const DTMtime * mTimeMap_
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
void percDevVDrift(DTChamberId indexCh, float meanVD, float sigmaVD, float &devVD, float &errdevVD)
edm::ESHandle< DTMtime > mTime
MonitorElement * glbT0Summary
std::map< std::string, MonitorElement * > allwheelHistos
std::map< int, std::map< std::string, MonitorElement * > > wheelHistos
double getRMS(int axis=1) const
get RMS of histogram along x, y or z axis (axis=1, 2, 3 respectively)