CalibrationSummaryClient.cc

Go to the documentation of this file.

#include "../interface/CalibrationSummaryClient.h"

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

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

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

#include <algorithm>

namespace ecaldqm
{
  CalibrationSummaryClient::CalibrationSummaryClient() :
    DQWorkerClient(),
    laserWlToME_(),
    ledWlToME_(),
    tpGainToME_(),
    tpPNGainToME_(),
    pedGainToME_(),
    pedPNGainToME_()
  {
  }

  void
  CalibrationSummaryClient::setParams(edm::ParameterSet const& _params)
  {
    std::vector<std::string> sourceList(_params.getUntrackedParameter<std::vector<std::string> >("activeSources"));
    if(std::find(sourceList.begin(), sourceList.end(), "Laser") == sourceList.end()){
      sources_.erase(std::string("Laser"));
      sources_.erase(std::string("LaserPN"));
    }
    if(std::find(sourceList.begin(), sourceList.end(), "Led") == sourceList.end()){
      sources_.erase(std::string("Led"));
      sources_.erase(std::string("LedPN"));
    }
    if(std::find(sourceList.begin(), sourceList.end(), "TestPulse") == sourceList.end()){
      sources_.erase(std::string("TestPulse"));
      sources_.erase(std::string("TestPulsePN"));
    }
    if(std::find(sourceList.begin(), sourceList.end(), "Pedestal") == sourceList.end()){
      sources_.erase(std::string("Pedestal"));
      sources_.erase(std::string("PedestalPN"));
    }

    MESet::PathReplacements repl;

    if(using_("Laser")){
      repl.clear();
      std::vector<int> laserWavelengths(_params.getUntrackedParameter<std::vector<int> >("laserWavelengths"));
      MESetMulti const& laser(static_cast<MESetMulti const&>(sources_.at("Laser")));
      unsigned nWL(laserWavelengths.size());
      for(unsigned iWL(0); iWL != nWL; ++iWL){
        int wl(laserWavelengths[iWL]);
        if(wl <= 0 || wl >= 5) throw cms::Exception("InvalidConfiguration") << "Laser Wavelength";
        repl["wl"] = std::to_string(wl);
        laserWlToME_[wl] = laser.getIndex(repl);
      }
    }

    if(using_("Led")){
      repl.clear();
      std::vector<int> ledWavelengths(_params.getUntrackedParameter<std::vector<int> >("ledWavelengths"));
      MESetMulti const& led(static_cast<MESetMulti const&>(sources_.at("Led")));
      unsigned nWL(ledWavelengths.size());
      for(unsigned iWL(0); iWL != nWL; ++iWL){
        int wl(ledWavelengths[iWL]);
        if(wl <= 0 || wl >= 5) throw cms::Exception("InvalidConfiguration") << "Led Wavelength";
        repl["wl"] = std::to_string(wl);
        ledWlToME_[wl] = led.getIndex(repl);
      }
    }

    if(using_("TestPulse")){
      repl.clear();
      std::vector<int> tpMGPAGains(_params.getUntrackedParameter<std::vector<int> >("testPulseMGPAGains"));
      MESetMulti const& tp(static_cast<MESetMulti const&>(sources_.at("TestPulse")));
      unsigned nG(tpMGPAGains.size());
      for(unsigned iG(0); iG != nG; ++iG){
        int gain(tpMGPAGains[iG]);
        if(gain != 1 && gain != 6 && gain != 12) throw cms::Exception("InvalidConfiguration") << "MGPA gain";
        repl["gain"] = std::to_string(gain);
        tpGainToME_[gain] = tp.getIndex(repl);
      }

      repl.clear();
      std::vector<int> tpMGPAGainsPN(_params.getUntrackedParameter<std::vector<int> >("testPulseMGPAGainsPN"));
      MESetMulti const& tppn(static_cast<MESetMulti const&>(sources_.at("TestPulsePN")));
      unsigned nGPN(tpMGPAGainsPN.size());
      for(unsigned iG(0); iG != nGPN; ++iG){
        int gain(tpMGPAGainsPN[iG]);
        if(gain != 1 && gain != 16) throw cms::Exception("InvalidConfiguration") << "PN MGPA gain";
        repl["pngain"] = std::to_string(gain);
        tpPNGainToME_[gain] = tppn.getIndex(repl);
      }
    }

    if(using_("Pedestal")){
      repl.clear();
      std::vector<int> pedMGPAGains(_params.getUntrackedParameter<std::vector<int> >("pedestalMGPAGains"));
      MESetMulti const& ped(static_cast<MESetMulti const&>(sources_.at("Pedestal")));
      unsigned nG(pedMGPAGains.size());
      for(unsigned iG(0); iG != nG; ++iG){
        int gain(pedMGPAGains[iG]);
        if(gain != 1 && gain != 6 && gain != 12) throw cms::Exception("InvalidConfiguration") << "MGPA gain";
        repl["gain"] = std::to_string(gain);
        pedGainToME_[gain] = ped.getIndex(repl);
      }

      repl.clear();
      std::vector<int> pedMGPAGainsPN(_params.getUntrackedParameter<std::vector<int> >("pedestalMGPAGainsPN"));
      MESetMulti const& pedpn(static_cast<MESetMulti const&>(sources_.at("PedestalPN")));
      unsigned nGPN(pedMGPAGainsPN.size());
      for(unsigned iG(0); iG != nGPN; ++iG){
        int gain(pedMGPAGainsPN[iG]);
        if(gain != 1 && gain != 16) throw cms::Exception("InvalidConfiguration") << "PN MGPA gain";
        repl["pngain"] = std::to_string(gain);
        pedPNGainToME_[gain] = pedpn.getIndex(repl);
      }
    }

    qualitySummaries_.insert("QualitySummary");
    qualitySummaries_.insert("PNQualitySummary");
  }

  void
  CalibrationSummaryClient::producePlots(ProcessType)
  {
    using namespace std;

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

    MESetMulti const* sLaser(using_("Laser") ? static_cast<MESetMulti const*>(&sources_.at("Laser")) : nullptr);
    MESetMulti const* sLaserPN(using_("LaserPN") ? static_cast<MESetMulti const*>(&sources_.at("LaserPN")) : nullptr);
    MESetMulti const* sLed(using_("Led") ? static_cast<MESetMulti const*>(&sources_.at("Led")) : nullptr);
    MESetMulti const* sLedPN(using_("LedPN") ? static_cast<MESetMulti const*>(&sources_.at("LedPN")) : nullptr);
    MESetMulti const* sTestPulse(using_("TestPulse") ? static_cast<MESetMulti const*>(&sources_.at("TestPulse")) : nullptr);
    MESetMulti const* sTestPulsePN(using_("TestPulsePN") ? static_cast<MESetMulti const*>(&sources_.at("TestPulsePN")) : nullptr);
    MESetMulti const* sPedestal(using_("Pedestal") ? static_cast<MESetMulti const*>(&sources_.at("Pedestal")) : nullptr);
    MESetMulti const* sPedestalPN(using_("PedestalPN") ? static_cast<MESetMulti const*>(&sources_.at("PedestalPN")) : nullptr);
    MESet const& sPNIntegrity(sources_.at("PNIntegrity"));

    MESet::iterator qEnd(meQualitySummary.end());
    for(MESet::iterator qItr(meQualitySummary.beginChannel()); qItr != qEnd; qItr.toNextChannel()){
      DetId id(qItr->getId());

      int status(kGood);

      if(status == kGood && sLaser){
        for(map<int, unsigned>::iterator wlItr(laserWlToME_.begin()); wlItr != laserWlToME_.end(); ++wlItr){
          sLaser->use(wlItr->second);
          if(sLaser->getBinContent(id) == kBad){
            status = kBad;
            break;
          }
        }
      }

      if(status == kGood && sLed){
        DetId id(qItr->getId());
        if(id.subdetId() == EcalEndcap){
          for(map<int, unsigned>::iterator wlItr(ledWlToME_.begin()); wlItr != ledWlToME_.end(); ++wlItr){
            sLed->use(wlItr->second);
            if(sLed->getBinContent(id) == kBad){
              status = kBad;
              break;
            }
          }
        }
      }

      if(status == kGood && sTestPulse){
        for(map<int, unsigned>::iterator gainItr(tpGainToME_.begin()); gainItr != tpGainToME_.end(); ++gainItr){
          sTestPulse->use(gainItr->second);
          if(sTestPulse->getBinContent(id) == kBad){
            status = kBad;
            break;
          }
        }
      }

      if(status == kGood && sPedestal){
        for(map<int, unsigned>::iterator gainItr(pedGainToME_.begin()); gainItr != pedGainToME_.end(); ++gainItr){
          sPedestal->use(gainItr->second);
          if(sPedestal->getBinContent(id) == kBad){
            status = kBad;
            break;
          }
        }
      }

      qItr->setBinContent(status);
    }

    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);

        int status(kGood);

        if(sPNIntegrity.getBinContent(id) == kBad) status = kBad;

        if(status == kGood && sLaserPN){
          for(map<int, unsigned>::iterator wlItr(laserWlToME_.begin()); wlItr != laserWlToME_.end(); ++wlItr){
            sLaserPN->use(wlItr->second);
            if(sLaserPN->getBinContent(id) == kBad){
              status = kBad;
              break;
            }
          }
        }

        if(status == kGood && sLedPN){
          for(map<int, unsigned>::iterator wlItr(ledWlToME_.begin()); wlItr != ledWlToME_.end(); ++wlItr){
            sLedPN->use(wlItr->second);
            if(sLedPN->getBinContent(id) == kBad){
              status = kBad;
              break;
            }
          }
        }

        if(status == kGood && sTestPulsePN){
          for(map<int, unsigned>::iterator gainItr(tpPNGainToME_.begin()); gainItr != tpPNGainToME_.end(); ++gainItr){
            sTestPulsePN->use(gainItr->second);
            if(sTestPulsePN->getBinContent(id) == kBad){
              status = kBad;
              break;
            }
          }
        }

        if(status == kGood && sPedestalPN){
          for(map<int, unsigned>::iterator gainItr(pedPNGainToME_.begin()); gainItr != pedPNGainToME_.end(); ++gainItr){
            sPedestalPN->use(gainItr->second);
            if(sPedestalPN->getBinContent(id) == kBad){
              status = kBad;
              break;
            }
          }
        }

        mePNQualitySummary.setBinContent(id, status);
      }
    }
  }

  DEFINE_ECALDQM_WORKER(CalibrationSummaryClient);
}