/data/refman/pasoursint/CMSSW_5_3_10_patch2/src/DQM/EcalBarrelMonitorClient/src/SummaryClient.cc

Go to the documentation of this file.

#include "../interface/SummaryClient.h"

#include "DQM/EcalBarrelMonitorClient/interface/IntegrityClient.h"
#include "DQM/EcalBarrelMonitorClient/interface/PresampleClient.h"
#include "DQM/EcalBarrelMonitorClient/interface/TimingClient.h"
#include "DQM/EcalBarrelMonitorClient/interface/RawDataClient.h"
#include "DQM/EcalBarrelMonitorTasks/interface/OccupancyTask.h"

#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"

namespace ecaldqm {

  SummaryClient::SummaryClient(const edm::ParameterSet& _params, const edm::ParameterSet& _paths) :
    DQWorkerClient(_params, _paths, "SummaryClient")
  {
    edm::ParameterSet const& sources(_params.getUntrackedParameterSet("sources"));
    source_(sIntegrity, "IntegrityClient", IntegrityClient::kQuality, sources);
    source_(sPresample, "PresampleClient", PresampleClient::kQuality, sources);
    source_(sTiming, "TimingClient", TimingClient::kQuality, sources);
    source_(sRawData, "RawDataClient", RawDataClient::kQualitySummary, sources);
    source_(sDigiOccupancy, "OccupancyTask", OccupancyTask::kDigi, sources);
  }

  void
  SummaryClient::bookMEs()
  {
  }

  void
  SummaryClient::beginRun(const edm::Run &, const edm::EventSetup &)
  {
    for(unsigned iME(0); iME < MEs_.size(); ++iME)
      MEs_[iME]->book();

    MEs_[kQualitySummary]->resetAll(-1.);
    MEs_[kQualitySummary]->reset(1.);
    MEs_[kReportSummaryMap]->resetAll(-1.);
    MEs_[kReportSummaryMap]->reset(1.);
    MEs_[kReportSummaryContents]->reset(1.);
    MEs_[kReportSummary]->reset(1.);
  }

  void
  SummaryClient::producePlots()
  {
    MEs_[kQualitySummary]->reset(2.);
    MEs_[kReportSummaryMap]->reset(1.);

    float totalChannels(0.);
    float totalGood(0.);

    for(unsigned dccid(1); dccid <= 54; dccid++){
      MonitorElement const* me(sources_[sDigiOccupancy]->getME(dccid - 1));
      if(me && me->getTH1()->GetEntries() < 1.){
        MEs_[kReportSummaryContents]->fill(dccid, 1.);
        continue;
      }

      float dccChannels(0.);
      float dccGood(0.);

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

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

        float towerChannels(0.);
        float towerGood(0.);

        for(std::vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){

          float integrity(sources_[sIntegrity]->getBinContent(*idItr));

          if(integrity == 2.) continue;

          float presample(sources_[sPresample]->getBinContent(*idItr));
          float timing(sources_[sTiming]->getBinContent(*idItr));
          float rawdata(sources_[sRawData]->getBinContent(*idItr));

          if(integrity > 2.) integrity = 1.;
          if(presample > 2.) presample = 1.;
          if(timing > 2.) timing = 1.;
          if(rawdata > 2.) rawdata = 1.;

          float status(1.);
          if(integrity < 2.) status *= integrity;
          if(presample < 2.) status *= presample;
          if(timing < 2.) status *= timing;
          if(rawdata < 2.) status *= rawdata;

          MEs_[kQualitySummary]->setBinContent(*idItr, status);

          if(status == 1.){
            towerGood += 1.;
            dccGood += 1.;
            totalGood += 1.;
          }
          towerChannels += 1.;
          dccChannels += 1.;
          totalChannels += 1.;
        }

        if(towerChannels < 1.) continue;

        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)
            MEs_[kReportSummaryMap]->setBinContent(*scItr, towerGood / towerChannels);
        }
        else
          MEs_[kReportSummaryMap]->setBinContent(ids[0], towerGood / towerChannels);
      }

      if(dccChannels < 1.) continue;

      MEs_[kReportSummaryContents]->fill(dccid, dccGood / dccChannels);
    }

    if(totalChannels > 0.) MEs_[kReportSummary]->fill(totalGood / totalChannels);

  }

  /*static*/
  void
  SummaryClient::setMEData(std::vector<MEData>& _data)
  {
    _data[kQualitySummary] = MEData("QualitySummary", BinService::kEcal2P, BinService::kCrystal, MonitorElement::DQM_KIND_TH2F);
    _data[kReportSummaryMap] = MEData("ReportSummaryMap", BinService::kEcal, BinService::kSuperCrystal, MonitorElement::DQM_KIND_TH2F);
    _data[kReportSummaryContents] = MEData("ReportSummaryContents", BinService::kSM, BinService::kReport, MonitorElement::DQM_KIND_REAL);
    _data[kReportSummary] = MEData("ReportSummary", BinService::kEcal, BinService::kReport, MonitorElement::DQM_KIND_REAL);
  }

  DEFINE_ECALDQM_WORKER(SummaryClient);
}