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;
158 float vDriftMean = VDriftME->
getMean();
159 float t0Mean = T0ME->
getMean();
162 float vDriftSigma = VDriftME->
getRMS();
163 float t0Sigma = T0ME->
getRMS();
189 float vDriftDev(0.), errvDriftDev(0.);
190 percDevVDrift(indexCh, vDriftMean, vDriftSigma, vDriftDev, errvDriftDev);
193 if( sec == 13 || sec == 14 ) sec_ = ( sec == 13 ) ? 4 : 10;
204 if( sec == 13 || sec == 14 ) {
206 float binVDriftDev = (
wheelHistos[wheel])[
"MeanVDriftSummary"]->getBinContent(sec_,stat);
207 binVDriftDev = (fabs(binVDriftDev) > fabs(fillvDriftDev)) ? binVDriftDev : fillvDriftDev;
208 (
wheelHistos[wheel])[
"MeanVDriftSummary"] -> setBinContent(sec_,stat,binVDriftDev);
210 float binT0MeanVal = (
wheelHistos[wheel])[
"MeanT0Summary"] -> getBinContent(sec_,stat);
211 binT0MeanVal = (fabs(binT0MeanVal) > fabs(fillT0Mean)) ? binT0MeanVal : fillT0Mean;
212 (
wheelHistos[wheel])[
"MeanT0Summary"] -> setBinContent(sec_,stat,binT0MeanVal);
214 float binVDriftSigmVal = (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> getBinContent(sec_,stat);
215 binVDriftSigmVal = (binVDriftSigmVal > 0. && binVDriftSigmVal < vDriftSigmQ) ? binVDriftSigmVal : vDriftSigmQ;
216 (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> setBinContent(sec_,stat,binVDriftSigmVal);
218 float binT0SigmVal = (
wheelHistos[wheel])[
"SigmaT0Summary"] -> getBinContent(sec_,stat);
219 binT0SigmVal = (binT0SigmVal > 0. && binT0SigmVal < t0SigmQ) ? binT0SigmVal : t0SigmQ;
220 (
wheelHistos[wheel])[
"SigmaT0Summary"] -> setBinContent(sec_,stat,binT0SigmVal);
224 (
wheelHistos[wheel])[
"MeanVDriftSummary"] -> setBinContent(sec_,stat,fillvDriftDev);
225 (
wheelHistos[wheel])[
"MeanT0Summary"] -> setBinContent(sec_,stat,fillT0Mean);
226 (
wheelHistos[wheel])[
"SigmaVDriftSummary"] -> setBinContent(sec_,stat,vDriftSigmQ);
227 (
wheelHistos[wheel])[
"SigmaT0Summary"] -> setBinContent(sec_,stat,t0SigmQ);
232 if(( sec_ == 4 || sec_ == 10) && stat == 4) weight = 1/8.;
234 if( vDriftDevQ > 0.85 && vDriftSigmQ > 0.85 ) {
240 if( vDriftDevQ > 0.85 && vDriftSigmQ < 0.85 ) {
243 if( vDriftDevQ < 0.85 && vDriftSigmQ > 0.85 ) {
248 if( t0MeanQ > 0.85 && t0SigmQ > 0.85 ) {
253 if( t0MeanQ > 0.85 && t0SigmQ < 0.85 ) {
257 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_
const T & max(const T &a, const T &b)
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)