/data/refman/pasoursint/CMSSW_5_3_4/src/DQM/EcalCommon/src/DQWorkerClient.cc

Go to the documentation of this file.

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

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

#include "FWCore/Utilities/interface/Exception.h"

namespace ecaldqm {

  EcalDQMChannelStatus const* DQWorkerClient::channelStatus(0);
  EcalDQMTowerStatus const* DQWorkerClient::towerStatus(0);

  DQWorkerClient::DQWorkerClient(const edm::ParameterSet &_params, const edm::ParameterSet& _paths, std::string const& _name) :
    DQWorker(_params, _paths, _name),
    sources_(0)
  {
  }

  void
  DQWorkerClient::endLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
  {
    for(std::vector<MESet const*>::iterator sItr(sources_.begin()); sItr != sources_.end(); ++sItr){
      MESetChannel const* meset(dynamic_cast<MESetChannel const*>(*sItr));
      if(meset) meset->checkDirectory();
    }
  }

  void
  DQWorkerClient::reset()
  {
    DQWorker::reset();
    for(std::vector<MESet const*>::iterator sItr(sources_.begin()); sItr != sources_.end(); ++sItr)
      (*sItr)->clear();
  }

  void
  DQWorkerClient::initialize()
  {
    initialized_ = true;
    for(std::vector<MESet const*>::iterator sItr(sources_.begin()); sItr != sources_.end(); ++sItr)
      initialized_ &= (*sItr)->retrieve();
  }

  void
  DQWorkerClient::source_(unsigned _iS, std::string const& _worker, unsigned _iW, edm::ParameterSet const& _sources)
  {
    if(_iS >= sources_.size()) sources_.resize(_iS + 1, 0);

    std::map<std::string, std::vector<MEData> >::const_iterator dataItr(meData.find(_worker));
    if(dataItr == meData.end())
      throw cms::Exception("InvalidCall") << "DQWorker " << _worker << " is not defined";

    MEData const& data(dataItr->second.at(_iW));

    edm::ParameterSet const& workerPaths(_sources.getUntrackedParameterSet(_worker));

    std::string fullpath(workerPaths.getUntrackedParameter<std::string>(data.pathName));

    sources_.at(_iS) = createMESet_(fullpath, data, true);
  }

  void
  DQWorkerClient::fillQuality_(unsigned _iME, DetId const& _id, uint32_t _mask, float _quality)
  {
    using namespace std;

    bool doMask(false);

    // turn off masking for good channel for the time being
    // update the RP then enable again
    if(_quality != 1. && channelStatus && towerStatus){
      BinService::BinningType btype(MEs_[_iME]->getBinType());
      bool searchTower(btype == BinService::kTriggerTower || btype == BinService::kSuperCrystal);

      switch(_id.subdetId()){
      case EcalBarrel:
        if(searchTower){
          EcalTrigTowerDetId ttid(EBDetId(_id).tower());
          int tccid(getElectronicsMap()->TCCid(ttid));
          int itt(getElectronicsMap()->iTT(ttid));
          vector<DetId> ids(getElectronicsMap()->ttConstituents(tccid, itt));
          for(vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
            if(doMask) break;
            EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(idItr->rawId()));
            if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;
          }
        }
        else{
          EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(_id.rawId()));
          if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;
        }

        if(!doMask){
          EcalDQMTowerStatus::const_iterator tItr(towerStatus->find(EBDetId(_id).tower().rawId()));
          if(tItr != towerStatus->end()) doMask |= tItr->getStatusCode() & _mask;
        }

        break;

      case EcalEndcap:
        if(isEcalScDetId(_id)){
          EcalScDetId scid(_id);
          for(int ix(1); ix <= 5; ix++){
            for(int iy(1); iy <= 5; iy++){
              if(doMask) break;
              int iix((scid.ix() - 1) * 5 + ix);
              int iiy((scid.iy() - 1) * 5 + iy);
              if(!EEDetId::validDetId(iix, iiy, scid.zside())) continue;
              EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(EEDetId(iix, iiy, scid.zside()).rawId()));
              if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;
            }
          }

          if(!doMask){
            EcalDQMTowerStatus::const_iterator tItr(towerStatus->find(_id.rawId()));
            if(tItr != towerStatus->end()) doMask |= tItr->getStatusCode() & _mask;
          }
        }
        else{
          if(searchTower){
            EcalScDetId scid(EEDetId(_id).sc());
            for(int ix(1); ix <= 5; ix++){
              for(int iy(1); iy <= 5; iy++){
                if(doMask) break;
                int iix((scid.ix() - 1) * 5 + ix);
                int iiy((scid.iy() - 1) * 5 + iy);
                if(!EEDetId::validDetId(iix, iiy, scid.zside())) continue;
                EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(EEDetId(iix, iiy, scid.zside()).rawId()));
                if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;
              }
            }
          }
          else{
            EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(_id.rawId()));
            if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;
          }

          if(!doMask){
            EcalDQMTowerStatus::const_iterator tItr(towerStatus->find(EEDetId(_id).sc().rawId()));
            if(tItr != towerStatus->end()) doMask |= tItr->getStatusCode() & _mask;
          }
        }

        break;

      case EcalTriggerTower:
        {
          EcalTrigTowerDetId ttid(_id);
          vector<DetId> ids(getTrigTowerMap()->constituentsOf(ttid));
          for(vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
            if(doMask) break;
            EcalDQMChannelStatus::const_iterator cItr(channelStatus->find(idItr->rawId()));
            if(cItr != channelStatus->end()) doMask |= cItr->getStatusCode() & _mask;

            if(doMask) break;
            if(idItr->subdetId() == EcalBarrel){
              if(idItr != ids.begin()) continue;
              EcalDQMTowerStatus::const_iterator tItr(towerStatus->find(EBDetId(*idItr).tower().rawId()));
              if(tItr != towerStatus->end()) doMask |= tItr->getStatusCode() & _mask;
            }
            else{
              EcalDQMTowerStatus::const_iterator tItr(towerStatus->find(EEDetId(*idItr).sc().rawId()));
              if(tItr != towerStatus->end()) doMask |= tItr->getStatusCode() & _mask;
            }
          }
        }

        break;

      default:
        break;
      }
    }

    float quality(doMask ? _quality + 3. : _quality);

    MEs_[_iME]->setBinContent(_id, quality);
  }

}