PNIntegrityClient.cc

Go to the documentation of this file.

#include "DQM/EcalMonitorClient/interface/PNIntegrityClient.h"

#include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"

#include "DataFormats/EcalDetId/interface/EcalPnDiodeDetId.h"

#include "CondFormats/EcalObjects/interface/EcalDQMStatusHelper.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

namespace ecaldqm
{
  PNIntegrityClient::PNIntegrityClient() :
    DQWorkerClient(),
    errFractionThreshold_(0.)
  {
    qualitySummaries_.insert("QualitySummary");
  }

  void
  PNIntegrityClient::setParams(edm::ParameterSet const& _params)
  {
    errFractionThreshold_ = _params.getUntrackedParameter<double>("errFractionThreshold");
  }

  void
  PNIntegrityClient::producePlots(ProcessType)
  {
    uint32_t mask(0x1 << EcalDQMStatusHelper::TT_SIZE_ERROR |
                  0x1 << EcalDQMStatusHelper::TT_ID_ERROR |
                  0x1 << EcalDQMStatusHelper::CH_ID_ERROR |
                  0x1 << EcalDQMStatusHelper::CH_GAIN_ZERO_ERROR);

    MESet& meQualitySummary(MEs_.at("QualitySummary"));

    MESet const& sOccupancy(sources_.at("Occupancy"));
    MESet const& sMEMChId(sources_.at("MEMChId"));
    MESet const& sMEMGain(sources_.at("MEMGain"));
    MESet const& sMEMBlockSize(sources_.at("MEMBlockSize"));
    MESet const& sMEMTowerId(sources_.at("MEMTowerId"));

    for(unsigned iDCC(0); iDCC < nDCC; ++iDCC){
      if(memDCCIndex(iDCC + 1) == unsigned(-1)) continue;
      for(unsigned iPN(0); iPN < 10; ++iPN){
        int subdet(0);
        if(iDCC >= kEBmLow && iDCC <= kEBpHigh) subdet = EcalBarrel;
        else subdet = EcalEndcap;

        EcalPnDiodeDetId id(subdet, iDCC + 1, iPN + 1);

        bool doMask(meQualitySummary.maskMatches(id, mask, statusManager_));

        float entries(sOccupancy.getBinContent(id));

        float chid(sMEMChId.getBinContent(id));
        float gain(sMEMGain.getBinContent(id));

        float blocksize(sMEMBlockSize.getBinContent(id));
        float towerid(sMEMTowerId.getBinContent(id));

        if(entries + gain + chid + blocksize + towerid < 1.){
          meQualitySummary.setBinContent(id, doMask ? kMUnknown : kUnknown);
          continue;
        }

        float chErr((gain + chid + blocksize + towerid) / (entries + gain + chid + blocksize + towerid));

        if(chErr > errFractionThreshold_)
          meQualitySummary.setBinContent(id, doMask ? kMBad : kBad);
        else
          meQualitySummary.setBinContent(id, doMask ? kMGood : kGood);
      }
    }
  }

  DEFINE_ECALDQM_WORKER(PNIntegrityClient);
}