SiStripMonitorPedestals.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripMonitorPedestals
// Class:      SiStripMonitorPedestals
//
//
// Original Author:  Simone Gennai and Suchandra Dutta
//         Created:  Sat Feb  4 20:49:10 CET 2006
//
//

#include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
#include "CalibTracker/Records/interface/SiStripQualityRcd.h"
#include "CalibTracker/SiStripAPVAnalysis/interface/ApvAnalysisFactory.h"
#include "CondFormats/DataRecord/interface/SiStripNoisesRcd.h"
#include "CondFormats/DataRecord/interface/SiStripPedestalsRcd.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
#include "CondFormats/SiStripObjects/interface/SiStripPedestals.h"
#include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
#include "DQM/SiStripCommon/interface/SiStripHistoId.h"
#include "DQM/SiStripMonitorPedestals/interface/SiStripMonitorPedestals.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"

// ROOT profile
#include "TProfile.h"

// std includes
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <numeric>

const std::string SiStripMonitorPedestals::RunMode1 = "ConDBPlotsOnly";
const std::string SiStripMonitorPedestals::RunMode2 = "CalculatedPlotsOnly";
const std::string SiStripMonitorPedestals::RunMode3 = "AllPlots";

SiStripMonitorPedestals::SiStripMonitorPedestals(edm::ParameterSet const &iConfig)
    : dqmStore_(edm::Service<DQMStore>().operator->()),
      conf_(iConfig),
      pedsPSet_(iConfig.getParameter<edm::ParameterSet>("PedestalsPSet")),
      analyzed(false),
      firstEvent(true),
      signalCutPeds_(4),
      nEvTot_(0),
      nIteration_(0),
      apvFactory_(nullptr),
      tTopoToken_(esConsumes<edm::Transition::BeginRun>()),
      detCablingToken_(esConsumes<edm::Transition::BeginRun>()) {
  // retrieve producer name of input StripDigiCollection
  std::string digiProducer = conf_.getParameter<std::string>("DigiProducer");
  std::string digiType = "VirginRaw";
  digiToken_ = consumes<edm::DetSetVector<SiStripRawDigi>>(edm::InputTag(digiProducer, digiType));

  edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals  "
                                          << " Constructing....... ";

  theEventInitNumber_ = pedsPSet_.getParameter<int>("NumberOfEventsForInit");
  theEventIterNumber_ = pedsPSet_.getParameter<int>("NumberOfEventsForIteration");
  NumCMstripsInGroup_ = pedsPSet_.getParameter<int>("NumCMstripsInGroup");
  runTypeFlag_ = conf_.getParameter<std::string>("RunTypeFlag");

  if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
    pedestalToken_ = esConsumes<edm::Transition::BeginRun>();
    noiseToken_ = esConsumes<edm::Transition::BeginRun>();
    qualityToken_ = esConsumes<edm::Transition::BeginRun>(
        edm::ESInputTag{"", iConfig.getParameter<std::string>("StripQualityLabel")});
  }
}
//
// -- Destructor
//
SiStripMonitorPedestals::~SiStripMonitorPedestals() {
  edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals  "
                                          << " Destructing...... ";
  if (apvFactory_) {
    delete apvFactory_;
  }
}
//
// -- BeginRun
//

void SiStripMonitorPedestals::bookHistograms(DQMStore::IBooker &ibooker,
                                             const edm::Run &run,
                                             const edm::EventSetup &eSetup) {
  if (detCablingWatcher_.check(eSetup)) {
    detcabling = &eSetup.getData(detCablingToken_);
    edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::bookHistograms: "
                                            << " Creating MEs for new Cabling ";
    createMEs(ibooker, eSetup);
  }

  if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3)
    fillCondDBMEs(eSetup);
}

//
// -- Create Monitor Elements
//
void SiStripMonitorPedestals::createMEs(DQMStore::IBooker &ibooker, const edm::EventSetup &es) {
  const auto tTopo = &es.getData(tTopoToken_);

  std::vector<uint32_t> SelectedDetIds;

  // ApvAnalysisFactory
  if (apvFactory_) {
    delete apvFactory_;
  }
  apvFactory_ = new ApvAnalysisFactory(pedsPSet_);

  detcabling->addActiveDetectorsRawIds(SelectedDetIds);

  // use SistripHistoId for producing histogram id (and title)
  SiStripHistoId hidmanager;
  // create SiStripFolderOrganizer
  SiStripFolderOrganizer folder_organizer;

  edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: "
                                          << "Number of Detector Present in cabling " << SelectedDetIds.size();

  for (std::vector<uint32_t>::const_iterator idetid = SelectedDetIds.begin(), iEnd = SelectedDetIds.end();
       idetid != iEnd;
       ++idetid) {
    uint32_t detid = *idetid;

    // Check consistency in DetId
    if (detid == 0 || detid == 0xFFFFFFFF) {
      edm::LogError("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: "
                                               << "Wrong DetId !!!!!! " << detid << " Neglecting !!!!!! ";
      continue;
    }

    unsigned int apv_pairs = detcabling->nApvPairs(detid);

    // Check consistency in Apv numbers
    if (apv_pairs < 1 || apv_pairs > 3) {
      edm::LogError("SiStripMonitorPedestals") << "SiStripMonitorPedestals::createMEs: Wrong APV Pairs  => detId "
                                               << detid << " APV pairs " << apv_pairs << " Neglecting !!!!!! ";
      continue;
    }
    unsigned int napvs = apv_pairs * 2;
    unsigned int nStrip = napvs * 128;

    bool newDetId = apvFactory_->instantiateApvs(detid, napvs);

    if (newDetId) {
      ModMEs local_modmes;
      local_modmes.PedsPerStrip = nullptr;
      local_modmes.PedsDistribution = nullptr;
      local_modmes.PedsEvolution = nullptr;
      local_modmes.CMSubNoisePerStrip = nullptr;
      local_modmes.RawNoisePerStrip = nullptr;
      local_modmes.CMSubNoiseProfile = nullptr;
      local_modmes.RawNoiseProfile = nullptr;
      local_modmes.NoisyStrips = nullptr;
      local_modmes.NoisyStripDistribution = nullptr;
      local_modmes.CMDistribution = nullptr;
      local_modmes.CMSlopeDistribution = nullptr;
      local_modmes.PedsPerStripDB = nullptr;
      local_modmes.CMSubNoisePerStripDB = nullptr;
      local_modmes.BadStripsDB = nullptr;

      std::string hid;
      // set appropriate folder using SiStripFolderOrganizer
      folder_organizer.setDetectorFolder(detid, tTopo);  // pass the detid to this method

      // if the deid already exists in the map, then reset MEs otherwise create
      // them
      resetMEs(detid);

      if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
        // Pedestals histos
        hid = hidmanager.createHistoId("PedestalFromCondDB", "det", detid);
        local_modmes.PedsPerStripDB =
            ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);  // to modify the size binning
        (local_modmes.PedsPerStripDB)->setAxisTitle("Pedestal from CondDB(ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("NoiseFromCondDB", "det", detid);
        local_modmes.CMSubNoisePerStripDB = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
        (local_modmes.CMSubNoisePerStripDB)->setAxisTitle("CMSubNoise from CondDB(ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("BadStripFlagCondDB", "det", detid);
        local_modmes.BadStripsDB = ibooker.book2D(hid, hid, nStrip, 0.5, nStrip + 0.5, 6, -0.5, 5.5);
        (local_modmes.BadStripsDB)->setAxisTitle("Strip Flag from CondDB(ADC) vs Strip Number", 1);
      }
      if (runTypeFlag_ == RunMode2 || runTypeFlag_ == RunMode3) {
        // Pedestals histos
        hid = hidmanager.createHistoId("PedsPerStrip", "det", detid);
        local_modmes.PedsPerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);  // to modify the size binning
        (local_modmes.PedsPerStrip)->setAxisTitle("Pedestal (ADC)  vs Strip Number ", 1);

        hid = hidmanager.createHistoId("PedsDistribution", "det", detid);
        local_modmes.PedsDistribution = ibooker.book2D(hid,
                                                       hid,
                                                       napvs,
                                                       -0.5,
                                                       napvs - 0.5,
                                                       300,
                                                       200,
                                                       500);  // to modify the size binning
        (local_modmes.PedsDistribution)->setAxisTitle("Apv Number", 1);
        (local_modmes.PedsDistribution)->setAxisTitle("Mean Pedestal Value (ADC)", 2);

        hid = hidmanager.createHistoId("PedsEvolution", "det", detid);
        local_modmes.PedsEvolution = ibooker.book2D(hid,
                                                    hid,
                                                    napvs,
                                                    -0.5,
                                                    napvs - 0.5,
                                                    50,
                                                    0.,
                                                    50.);  // to modify the size binning
        (local_modmes.PedsEvolution)->setAxisTitle("Apv Number", 1);
        (local_modmes.PedsEvolution)->setAxisTitle("Iteration Number", 2);

        // Noise histos
        hid = hidmanager.createHistoId("CMSubNoisePerStrip", "det", detid);
        local_modmes.CMSubNoisePerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
        (local_modmes.CMSubNoisePerStrip)->setAxisTitle("CMSubNoise (ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("RawNoisePerStrip", "det", detid);
        local_modmes.RawNoisePerStrip = ibooker.book1D(hid, hid, nStrip, 0.5, nStrip + 0.5);
        (local_modmes.RawNoisePerStrip)->setAxisTitle("RawNoise(ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("CMSubNoiseProfile", "det", detid);
        local_modmes.CMSubNoiseProfile = ibooker.bookProfile(hid, hid, nStrip, 0.5, nStrip + 0.5, 100, 0., 100.);
        (local_modmes.CMSubNoiseProfile)->setAxisTitle("Mean of CMSubNoise (ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("RawNoiseProfile", "det", detid);
        local_modmes.RawNoiseProfile = ibooker.bookProfile(hid, hid, nStrip, 0.5, nStrip + 0.5, 100, 0., 100.);
        (local_modmes.RawNoiseProfile)->setAxisTitle("Mean of RawNoise (ADC) vs Strip Number", 1);

        hid = hidmanager.createHistoId("NoisyStrips", "det", detid);
        local_modmes.NoisyStrips = ibooker.book2D(hid, hid, nStrip, 0.5, nStrip + 0.5, 6, -0.5, 5.5);
        (local_modmes.NoisyStrips)->setAxisTitle("Strip Number", 1);
        (local_modmes.NoisyStrips)->setAxisTitle("Flag Value", 2);

        hid = hidmanager.createHistoId("NoisyStripDistribution", "det", detid);
        local_modmes.NoisyStripDistribution = ibooker.book1D(hid, hid, 11, -0.5, 10.5);
        (local_modmes.NoisyStripDistribution)->setAxisTitle("Flag Value", 1);

        // Common Mode histos
        hid = hidmanager.createHistoId("CMDistribution", "det", detid);
        local_modmes.CMDistribution = ibooker.book2D(hid, hid, napvs, -0.5, napvs - 0.5, 150, -15., 15.);
        (local_modmes.CMDistribution)->setAxisTitle("Common Mode (ADC) vs APV Number", 1);

        hid = hidmanager.createHistoId("CMSlopeDistribution", "det", detid);
        local_modmes.CMSlopeDistribution = ibooker.book2D(hid, hid, napvs, -0.5, napvs - 0.5, 100, -0.05, 0.05);
        (local_modmes.CMSlopeDistribution)->setAxisTitle("Common Mode Slope vs APV Number", 1);
      }
      // data from CondDB
      // append to PedMEs
      PedMEs.insert(std::make_pair(detid, local_modmes));
    }  // newDetId
  }
  edm::LogInfo("SiStripMonitorPedestals")
      << "SiStripMonitorPedestals::createMEs: Number of DETS used " << PedMEs.size();
}
// ------------ method called to produce the data  ------------
void SiStripMonitorPedestals::analyze(const edm::Event &iEvent, const edm::EventSetup &eSetup) {
  edm::LogInfo("SiStripMonitorPedestals")
      << "SiStripMonitorPedestals::analyze: Run " << iEvent.id().run() << " Event " << iEvent.id().event();

  detcabling = &eSetup.getData(detCablingToken_);

  if (runTypeFlag_ != RunMode2 && runTypeFlag_ != RunMode3)
    return;

  // Increment # of Events
  nEvTot_++;

  // get DigiCollection object from Event
  edm::Handle<edm::DetSetVector<SiStripRawDigi>> digi_collection;
  // you have a collection as there are all the digis for the event for every
  // detector
  iEvent.getByToken(digiToken_, digi_collection);

  // Increase the number of iterations ...
  if ((nEvTot_ - theEventInitNumber_) % theEventIterNumber_ == 1)
    nIteration_++;

  // loop over all MEs
  for (std::map<uint32_t, ModMEs>::const_iterator i = PedMEs.begin(); i != PedMEs.end(); i++) {
    uint32_t detid = i->first;
    ModMEs local_modmes = i->second;
    // get iterators for digis belonging to one DetId, it is an iterator, i.e.
    // one element of the vector
    std::vector<edm::DetSet<SiStripRawDigi>>::const_iterator digis = digi_collection->find(detid);
    if (digis == digi_collection->end() || digis->data.empty() || digis->data.size() > 768) {
      if (digis == digi_collection->end()) {
        edm::LogError("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: Event " << nEvTot_ << " DetId "
                                                 << detid << " at the end of Digi Collection!!!";
      } else {
        edm::LogError("SiStripMonitorPedestals") << " [SiStripMonitorPedestals::analyze: Event " << nEvTot_ << " DetId "
                                                 << detid << " # of Digis " << digis->data.size();
      }
      std::vector<const FedChannelConnection *> fed_conns = detcabling->getConnections(detid);
      bool firstchannel(true);
      for (unsigned int k = 0; k < fed_conns.size(); k++) {
        if (fed_conns[k] && fed_conns[k]->isConnected()) {
          if (firstchannel) {
            edm::LogError("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: Fed Id "
                                                     << fed_conns[k]->fedId() << " Channel " << fed_conns[k]->fedCh();
            firstchannel = false;
          } else
            edm::LogError("SiStripMonitorPedestals")
                << "  SiStripMonitorPedestals::analyze: Channel " << fed_conns[k]->fedCh();
        }
      }
      std::cout << std::endl;
      continue;
    }

    if (digis->data.empty()) {
      edm::LogError("MonitorDigi_tmp") << "[SiStripRawDigiToRaw::createFedBuffers] Zero digis found!";
    }
    uint32_t id = detid;
    //    cout <<"Data size "<<digis->data.size()<<endl;
    apvFactory_->update(id, (*digis));

    if (nEvTot_ > theEventInitNumber_) {
      if (local_modmes.CMDistribution != nullptr) {
        std::vector<float> tmp;
        tmp.clear();
        apvFactory_->getCommonMode(id, tmp);
        // unpacking the info looking at the right topology
        int numberCMBlocks = int(128. / NumCMstripsInGroup_);
        int ibin = 0;
        for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
          int iapv = int(ibin / numberCMBlocks);
          (local_modmes.CMDistribution)->Fill(iapv, static_cast<float>(*iped));
          ibin++;
        }
      }
      if (local_modmes.CMSlopeDistribution != nullptr) {
        std::vector<float> tmp;
        tmp.clear();
        int iapv = 0;
        apvFactory_->getCommonModeSlope(id, tmp);
        for (std::vector<float>::const_iterator it = tmp.begin(); it != tmp.end(); it++) {
          (local_modmes.CMSlopeDistribution)->Fill(iapv, static_cast<float>(*it));
          iapv++;
        }
      }
    }

    // asking for the status
    if ((nEvTot_ - theEventInitNumber_ - theEventIterNumber_) % theEventIterNumber_ == 1) {
      std::vector<float> tmp;
      tmp.clear();
      apvFactory_->getPedestal(id, tmp);
      if (local_modmes.PedsPerStrip != nullptr) {
        int numberOfApvs = int(tmp.size() / 128.);
        for (int i = 0; i < numberOfApvs; i++) {
          std::vector<float> myPedPerApv;
          apvFactory_->getPedestal(id, i, myPedPerApv);
          float avarage = 0;
          avarage = std::accumulate(myPedPerApv.begin(), myPedPerApv.end(), avarage);
          avarage = avarage / 128.;
          (local_modmes.PedsEvolution)->setBinContent(i + 1, nIteration_, avarage);
        }
        int ibin = 0;

        for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
          int napv = int(ibin / 128.);
          ibin++;
          float last_value = (local_modmes.PedsPerStrip)->getBinContent(ibin);
          if (last_value != 0.) {
            (local_modmes.PedsPerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
          } else {
            (local_modmes.PedsPerStrip)->setBinContent(ibin, static_cast<float>(*iped));
          }
          (local_modmes.PedsDistribution)->Fill(napv, static_cast<float>(*iped));
        }
      }

      if (local_modmes.CMSubNoisePerStrip != nullptr && local_modmes.CMSubNoiseProfile != nullptr) {
        tmp.clear();
        apvFactory_->getNoise(id, tmp);
        int ibin = 0;
        for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
          ibin++;
          (local_modmes.CMSubNoiseProfile)->Fill(static_cast<double>(ibin * 1.), static_cast<float>(*iped));

          float last_value = (local_modmes.CMSubNoisePerStrip)->getBinContent(ibin);
          if (last_value != 0.) {
            (local_modmes.CMSubNoisePerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
          } else {
            (local_modmes.CMSubNoisePerStrip)->setBinContent(ibin, static_cast<float>(*iped));
          }
        }
      }

      if (local_modmes.RawNoisePerStrip != nullptr && local_modmes.RawNoiseProfile != nullptr) {
        tmp.clear();
        apvFactory_->getRawNoise(id, tmp);
        int ibin = 0;
        for (std::vector<float>::const_iterator iped = tmp.begin(); iped != tmp.end(); iped++) {
          ibin++;
          (local_modmes.RawNoiseProfile)->Fill(static_cast<double>(ibin * 1.), static_cast<float>(*iped));
          float last_value = (local_modmes.RawNoisePerStrip)->getBinContent(ibin);
          if (last_value != 0.) {
            (local_modmes.RawNoisePerStrip)->setBinContent(ibin, (static_cast<float>(*iped) + last_value) / 2.);
          } else {
            (local_modmes.RawNoisePerStrip)->setBinContent(ibin, static_cast<float>(*iped));
          }
        }
      }

      if (local_modmes.NoisyStrips != nullptr) {
        TkApvMask::MaskType temp;
        apvFactory_->getMask(id, temp);
        int ibin = 0;
        for (TkApvMask::MaskType::const_iterator iped = temp.begin(); iped != temp.end(); iped++) {
          ibin++;

          if (nIteration_ < 2) {
            if (*iped == 1)
              (local_modmes.NoisyStrips)->Fill(ibin, 3.);
            if (*iped == 2)
              (local_modmes.NoisyStrips)->Fill(ibin, 4.);
            if (*iped == 0)
              (local_modmes.NoisyStrips)->Fill(ibin, 0.);
          } else {
            (local_modmes.NoisyStrips)->Fill(ibin, static_cast<float>(*iped));
            (local_modmes.NoisyStripDistribution)->Fill(static_cast<float>(*iped));
          }
        }
      }
    }
  }
  if (firstEvent)
    firstEvent = false;
}
//
// -- End Run
//
void SiStripMonitorPedestals::dqmEndRun(edm::Run const &run, edm::EventSetup const &eSetup) {
  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
  if (outputMEsInRootFile) {
    std::string outPutFileName = conf_.getParameter<std::string>("OutPutFileName");
    dqmStore_->save(outPutFileName);
  }
}
//
// -- End Job
//
void SiStripMonitorPedestals::endJob(void) {
  edm::LogInfo("SiStripMonitorPedestals") << "SiStripMonitorPedestals::EndJob: "
                                          << " Finishing!! ";
}
//
// -- Reset Monitor Elements corresponding to a detetor
//
void SiStripMonitorPedestals::resetMEs(uint32_t idet) {
  std::map<uint32_t, ModMEs>::iterator pos = PedMEs.find(idet);
  if (pos != PedMEs.end()) {
    ModMEs mod_me = pos->second;

    if (runTypeFlag_ == RunMode1 || runTypeFlag_ == RunMode3) {
      mod_me.PedsPerStripDB->Reset();
      mod_me.CMSubNoisePerStripDB->Reset();
      mod_me.BadStripsDB->Reset();
    }
    if (runTypeFlag_ == RunMode2 || runTypeFlag_ == RunMode3) {
      mod_me.PedsPerStrip->Reset();
      mod_me.PedsDistribution->Reset();
      mod_me.PedsEvolution->Reset();
      mod_me.CMSubNoisePerStrip->Reset();
      mod_me.RawNoisePerStrip->Reset();
      mod_me.CMSubNoiseProfile->Reset();
      mod_me.RawNoiseProfile->Reset();
      mod_me.NoisyStrips->Reset();
      mod_me.CMDistribution->Reset();
      mod_me.CMSlopeDistribution->Reset();
    }
  }
}
//
// -- Fill CondDB Monitoring Elements
//
void SiStripMonitorPedestals::fillCondDBMEs(edm::EventSetup const &eSetup) {
  const auto &pedestals = eSetup.getData(pedestalToken_);
  const auto &noises = eSetup.getData(noiseToken_);
  const auto &quality = eSetup.getData(qualityToken_);

  for (std::map<uint32_t, ModMEs>::const_iterator i = PedMEs.begin(); i != PedMEs.end(); i++) {
    uint32_t detid = i->first;
    ModMEs local_modmes = i->second;
    edm::LogInfo("SiStripMonitorPedestals") << " SiStripMonitorPedestals::analyze: "
                                            << " Get Ped/Noise/Bad Strips from CondDb for DetId " << detid;
    int nStrip = detcabling->nApvPairs(detid) * 256;
    // Get range of pedestal and noise for the detid
    SiStripNoises::Range noiseRange = noises.getRange(detid);
    SiStripPedestals::Range pedRange = pedestals.getRange(detid);
    SiStripQuality::Range qualityRange = quality.getRange(detid);

    for (int istrip = 0; istrip < nStrip; ++istrip) {
      try {
        // Fill Pedestals
        (local_modmes.PedsPerStripDB)->Fill(istrip + 1, pedestals.getPed(istrip, pedRange));
      } catch (cms::Exception &e) {
        edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                    "SiStripPedestalsService_.getPedestal("
                                                 << detid << "," << istrip << ") :  "
                                                 << " " << e.what();
      }
      try {
        // Fill Noises
        (local_modmes.CMSubNoisePerStripDB)->Fill(istrip + 1, noises.getNoise(istrip, noiseRange));

      } catch (cms::Exception &e) {
        edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                    "SiStripNoiseService_.getNoise("
                                                 << detid << "," << istrip << ") :  "
                                                 << " " << e.what();
      }
      try {
        // Fill BadStripsNoise
        (local_modmes.BadStripsDB)->Fill(istrip + 1, quality.IsStripBad(qualityRange, istrip) ? 1. : 0.);

      } catch (cms::Exception &e) {
        edm::LogError("SiStripMonitorPedestals") << "[SiStripMonitorPedestals::analyze]  cms::Exception accessing "
                                                    "SiStripNoiseService_.getDisable("
                                                 << detid << "," << istrip << ") :  "
                                                 << " " << e.what();
      }
    }  // close istrip loop
  }
}
DEFINE_FWK_MODULE(SiStripMonitorPedestals);