/data/refman/pasoursint/CMSSW_5_3_3/src/DQM/EcalBarrelMonitorClient/src/IntegrityClient.cc

Go to the documentation of this file.

#include "../interface/IntegrityClient.h"

#include "DQM/EcalBarrelMonitorTasks/interface/OccupancyTask.h"
#include "DQM/EcalBarrelMonitorTasks/interface/IntegrityTask.h"

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

namespace ecaldqm {

  IntegrityClient::IntegrityClient(const edm::ParameterSet& _params, const edm::ParameterSet& _paths) :
    DQWorkerClient(_params, _paths, "IntegrityClient"),
    errFractionThreshold_(0.)
  {
    edm::ParameterSet const& taskParams(_params.getUntrackedParameterSet(name_));
    errFractionThreshold_ = taskParams.getUntrackedParameter<double>("errFractionThreshold");

    edm::ParameterSet const& sources(_params.getUntrackedParameterSet("sources"));
    source_(sOccupancy, "OccupancyTask", OccupancyTask::kDigi, sources);
    source_(sGain, "IntegrityTask", IntegrityTask::kGain, sources);
    source_(sChId, "IntegrityTask", IntegrityTask::kChId, sources);
    source_(sGainSwitch, "IntegrityTask", IntegrityTask::kGainSwitch, sources);
    source_(sTowerId, "IntegrityTask", IntegrityTask::kTowerId, sources);
    source_(sBlockSize, "IntegrityTask", IntegrityTask::kBlockSize, sources);
  }

  void
  IntegrityClient::bookMEs()
  {
    DQWorker::bookMEs();

    MEs_[kQuality]->resetAll(-1.);
    MEs_[kQualitySummary]->resetAll(-1.);
  }

  void
  IntegrityClient::producePlots()
  {
    uint32_t mask(1 << EcalDQMStatusHelper::CH_ID_ERROR |
                  1 << EcalDQMStatusHelper::CH_GAIN_ZERO_ERROR |
                  1 << EcalDQMStatusHelper::CH_GAIN_SWITCH_ERROR |
                  1 << EcalDQMStatusHelper::TT_ID_ERROR |
                  1 << EcalDQMStatusHelper::TT_SIZE_ERROR);

    for(unsigned dccid(1); dccid <= 54; dccid++){
      for(unsigned tower(1); tower <= getNSuperCrystals(dccid); tower++){
        std::vector<DetId> ids(getElectronicsMap()->dccTowerConstituents(dccid, tower));

        if(ids.size() == 0) continue;

        float towerEntries(0.);
        bool towerGood(true);

        for(std::vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
          float entries(sources_[sOccupancy]->getBinContent(*idItr));
          towerEntries += entries;

          float gain(sources_[sGain]->getBinContent(*idItr));
          float chid(sources_[sChId]->getBinContent(*idItr));
          float gainswitch(sources_[sGainSwitch]->getBinContent(*idItr));

          if(entries + gain + chid + gainswitch < 1.){
            fillQuality_(kQuality, *idItr, mask, 2.);
            continue;
          }

          float chErr((gain + chid + gainswitch) / (entries + gain + chid + gainswitch));

          if(chErr > errFractionThreshold_){
            fillQuality_(kQuality, *idItr, mask, 0.);
            towerGood = false;
          }
          else
            fillQuality_(kQuality, *idItr, mask, 1.);
        }

        float towerid(sources_[sTowerId]->getBinContent(ids[0]));
        float blocksize(sources_[sBlockSize]->getBinContent(ids[0]));

        float quality(-1.);

        if(towerEntries + towerid + blocksize > 1.){
          float towerErr((towerid + blocksize) / (towerEntries + towerid + blocksize));
          if(towerErr > errFractionThreshold_) towerGood = false;

          quality = towerGood ? 1. : 0.;
        }
        else{
          quality = 2.;
        }

        if(dccid <= 9 || dccid >= 46){
          std::vector<EcalScDetId> scs(getElectronicsMap()->getEcalScDetId(dccid, tower));
          for(std::vector<EcalScDetId>::iterator scItr(scs.begin()); scItr != scs.end(); ++scItr)
            fillQuality_(kQualitySummary, *scItr, mask, quality);
        }
        else
          fillQuality_(kQualitySummary, ids[0], mask, quality);
      }
    }
  }

  /*static*/
  void
  IntegrityClient::setMEData(std::vector<MEData>& _data)
  {
    _data[kQuality] = MEData("Quality", BinService::kSM, BinService::kCrystal, MonitorElement::DQM_KIND_TH2F);
    _data[kQualitySummary] = MEData("QualitySummary", BinService::kEcal2P, BinService::kSuperCrystal, MonitorElement::DQM_KIND_TH2F);
  }

  DEFINE_ECALDQM_WORKER(IntegrityClient);
}