TowerStatusTask.cc

Go to the documentation of this file.

#include "../interface/TowerStatusTask.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "CondFormats/EcalObjects/interface/EcalDAQTowerStatus.h"
#include "CondFormats/DataRecord/interface/EcalDAQTowerStatusRcd.h"
#include "CondFormats/EcalObjects/interface/EcalDCSTowerStatus.h"
#include "CondFormats/DataRecord/interface/EcalDCSTowerStatusRcd.h"

#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "DataFormats/EcalDetId/interface/EcalScDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"

namespace ecaldqm {

  TowerStatusTask::TowerStatusTask() :
    DQWorkerClient(),
    doDAQInfo_(false),
    doDCSInfo_(false)
  {
    std::fill_n(daqStatus_, nDCC, 0.);
    std::fill_n(dcsStatus_, nDCC, 0.);
  }

  void
  TowerStatusTask::setParams(edm::ParameterSet const& _params)
  {
    doDAQInfo_ = _params.getUntrackedParameter<bool>("doDAQInfo");
    doDCSInfo_ = _params.getUntrackedParameter<bool>("doDCSInfo");

    if(doDAQInfo_ && doDCSInfo_) return;
    if(doDAQInfo_){
      MEs_.erase(std::string("DCSSummary"));
      MEs_.erase(std::string("DCSSummaryMap"));
      MEs_.erase(std::string("DCSContents"));
    }
    else if(doDCSInfo_){
      MEs_.erase(std::string("DAQSummary"));
      MEs_.erase(std::string("DAQSummaryMap"));
      MEs_.erase(std::string("DAQContents"));
    }
    else
      throw cms::Exception("InvalidConfiguration") << "Nothing to do in TowerStatusTask";
  }

  void
  TowerStatusTask::endLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const& _es)
  {
    if(doDAQInfo_){
      std::fill_n(daqStatus_, nDCC, 1.);

      edm::ESHandle<EcalDAQTowerStatus> daqHndl;
      _es.get<EcalDAQTowerStatusRcd>().get(daqHndl);
      if(daqHndl.isValid()){
        for(unsigned id(0); id < EcalTrigTowerDetId::kEBTotalTowers; id++){
          if(daqHndl->barrel(id).getStatusCode() != 0){
            EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(id));
            daqStatus_[dccId(ttid) - 1] -= 25. / 1700.;
          }
        }
        for(unsigned id(0); id < EcalScDetId::kSizeForDenseIndexing; id++){
          if(daqHndl->endcap(id).getStatusCode() != 0){
            EcalScDetId scid(EcalScDetId::unhashIndex(id));
            unsigned dccid(dccId(scid));
            daqStatus_[dccid - 1] -= double(scConstituents(scid).size()) / nCrystals(dccid);
          }
        }
      }
      else
        edm::LogWarning("EventSetup") << "EcalDAQTowerStatus record not valid";
    }

    if(doDCSInfo_){
      std::fill_n(dcsStatus_, nDCC, 1.);

      edm::ESHandle<EcalDCSTowerStatus> dcsHndl;
      _es.get<EcalDCSTowerStatusRcd>().get(dcsHndl);
      if(dcsHndl.isValid()){
        for(unsigned id(0); id < EcalTrigTowerDetId::kEBTotalTowers; id++){
          if(dcsHndl->barrel(id).getStatusCode() != 0){
            EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(id));
            dcsStatus_[dccId(ttid) - 1] -= 25. / 1700.;
          }
        }
        for(unsigned id(0); id < EcalScDetId::kSizeForDenseIndexing; id++){
          if(dcsHndl->endcap(id).getStatusCode() != 0){
            EcalScDetId scid(EcalScDetId::unhashIndex(id));
            unsigned dccid(dccId(scid));
            dcsStatus_[dccid - 1] -= double(scConstituents(scid).size()) / nCrystals(dccid);
          }
        }
      }
      else
        edm::LogWarning("EventSetup") << "EcalDCSTowerStatus record not valid";
    }
  }

  void
  TowerStatusTask::producePlots(ProcessType)
  {
    if(doDAQInfo_) producePlotsTask_(daqStatus_, "DAQ");
    if(doDCSInfo_) producePlotsTask_(dcsStatus_, "DCS");
  }

  void
  TowerStatusTask::producePlotsTask_(float const* _status, std::string const& _type)
  {
    MESet* meSummary(0);
    MESet* meSummaryMap(0);
    MESet* meContents(0);
    meSummary = &MEs_.at(_type + "Summary");
    meSummaryMap = &MEs_.at(_type + "SummaryMap");
    meContents = &MEs_.at(_type + "Contents");

    meSummary->reset(-1.);
    meSummaryMap->resetAll(-1.);
    meSummaryMap->reset();
    meContents->reset(-1.);

    float totalFraction(0.);
    for(int iDCC(0); iDCC < nDCC; iDCC++){
      meSummaryMap->setBinContent(iDCC + 1, _status[iDCC]);
      meContents->fill(iDCC + 1, _status[iDCC]);
      totalFraction += _status[iDCC] / nCrystals(iDCC + 1);
    }

    meSummary->fill(totalFraction);
  }

  DEFINE_ECALDQM_WORKER(TowerStatusTask);
}