SiStripMonitorRawData.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripMonitorRawData
// Class:      SiStripMonitorRawData
//
//
// Original Author:  Suchandra Dutta
//         Created:  Fri June  1 17:00:00 CET 2007
//
//

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"

#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiStripDigi/interface/SiStripRawDigi.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <FWCore/Framework/interface/EventSetup.h>

#include "DQM/SiStripMonitorPedestals/interface/SiStripMonitorRawData.h"

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

SiStripMonitorRawData::SiStripMonitorRawData(edm::ParameterSet const &iConfig)
    : BadFedNumber(nullptr),
      dqmStore_(edm::Service<DQMStore>().operator->()),
      conf_(iConfig),
      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("SiStripMonitorRawData") << "SiStripMonitorRawData  "
                                        << " Constructing....... ";
}

SiStripMonitorRawData::~SiStripMonitorRawData() {
  edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData  "
                                        << " Destructing....... ";
}
//
// -- BeginRun
//

void SiStripMonitorRawData::bookHistograms(DQMStore::IBooker &ibooker,
                                           const edm::Run &run,
                                           const edm::EventSetup &eSetup) {
  if (BadFedNumber)
    BadFedNumber->Reset();

  if (detCablingWatcher_.check(eSetup)) {
    SelectedDetIds.clear();
    const auto &detcabling = eSetup.getData(detCablingToken_);
    detcabling.addActiveDetectorsRawIds(SelectedDetIds);

    edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::bookHistograms: "
                                          << " Creating MEs for new Cabling ";
    ibooker.setCurrentFolder("Track/GlobalParameter");
    if (!BadFedNumber) {
      BadFedNumber =
          ibooker.book1D("FaultyFedNumberAndChannel", "Faulty Fed Id and Channel and Numbers", 60000, 0.5, 600.5);
      BadFedNumber->setAxisTitle("Fed Id and Channel numbers", 1);
    }
  }
}

// ------------ method called to produce the data  ------------
void SiStripMonitorRawData::analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) {
  edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::analyze: Run " << iEvent.id().run() << " Event "
                                        << iEvent.id().event();

  const auto &detcabling = iSetup.getData(detCablingToken_);

  // get DigiCollection object from Event
  edm::Handle<edm::DetSetVector<SiStripRawDigi>> digi_collection;
  iEvent.getByToken(digiToken_, digi_collection);

  for (std::vector<uint32_t>::const_iterator idetid = SelectedDetIds.begin(), iEnd = SelectedDetIds.end();
       idetid != iEnd;
       ++idetid) {
    std::vector<edm::DetSet<SiStripRawDigi>>::const_iterator digis = digi_collection->find((*idetid));
    if (digis == digi_collection->end() || digis->data.empty() || digis->data.size() > 768) {
      std::vector<const FedChannelConnection *> fed_conns = detcabling.getConnections((*idetid));
      for (unsigned int k = 0; k < fed_conns.size(); k++) {
        if (fed_conns[k] && fed_conns[k]->isConnected()) {
          float fed_id = fed_conns[k]->fedId() + 0.01 * fed_conns[k]->fedCh();
          BadFedNumber->Fill(fed_id);
        }
      }
      continue;
    }
  }
}
//
// -- End Run
//
void SiStripMonitorRawData::dqmEndRun(edm::Run const &run, edm::EventSetup const &eSetup) {
  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
  std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
  if (outputMEsInRootFile) {
    dqmStore_->save(outputFileName);
  }
}
//
// -- End Job
//
void SiStripMonitorRawData::endJob(void) {
  edm::LogInfo("SiStripMonitorRawData") << "SiStripMonitorRawData::EndJob: "
                                        << " Finishing!! ";
}
DEFINE_FWK_MODULE(SiStripMonitorRawData);