SiPixelErrorsDigisToCalibDigis.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiPixelErrorsDigisToCalibDigis
// Class:      SiPixelErrorsDigisToCalibDigis
//
//
// Original Author:  Ricardo Vasquez Sierra
//         Created:  Wed Apr  9 12:43:02 CEST 2008
//
//

// system include files
#include <memory>

#include "SiPixelErrorsDigisToCalibDigis.h"

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

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

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
SiPixelErrorsDigisToCalibDigis::SiPixelErrorsDigisToCalibDigis(const edm::ParameterSet& iConfig)

{
  siPixelProducerLabel_ = iConfig.getParameter<edm::InputTag>("SiPixelProducerLabelTag");
  createOutputFile_ = iConfig.getUntrackedParameter<bool>("saveFile", false);
  outputFilename_ = iConfig.getParameter<std::string>("outputFilename");
  daqBE_ = &*edm::Service<DQMStore>();
  folderMaker_ = new SiPixelFolderOrganizer();

  tPixelCalibDigiError = consumes<edm::DetSetVector<SiPixelCalibDigiError> >(siPixelProducerLabel_);

  //  std::cout<<"siPixelProducerLabel_ = "<<siPixelProducerLabel_<<std::endl;
  //  std::cout<<"createOutputFile_= "<< createOutputFile_<<std::endl;
  //  std::cout<<"outpuFilename_= "<< outputFilename_<< std::endl;
}

SiPixelErrorsDigisToCalibDigis::~SiPixelErrorsDigisToCalibDigis() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

//
// member functions
//

// ------------ method called to for each event  ------------
void SiPixelErrorsDigisToCalibDigis::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  using namespace edm;

  static int first(1);
  if (1 == first) {
    first = 0;
    iSetup.get<TrackerDigiGeometryRecord>().get(geom_);
    theHistogramIdWorker_ = new SiPixelHistogramId(siPixelProducerLabel_.label());
  }

  Handle<DetSetVector<SiPixelCalibDigiError> > thePlaquettes;
  iEvent.getByToken(tPixelCalibDigiError, thePlaquettes);
  // iEvent.getByLabel("siPixelCalibDigis", thePlaquettes);

  DetSetVector<SiPixelCalibDigiError>::const_iterator digiIter;

  for (digiIter = thePlaquettes->begin(); digiIter != thePlaquettes->end(); digiIter++) {
    uint32_t detId = digiIter->id;

    DetSet<SiPixelCalibDigiError>::const_iterator ipix;
    //loop over pixel errors pulsed in the current plaquette

    MonitorElement* temp_;

    std::map<uint32_t, MonitorElement*>::iterator mapIterator =
        SiPixelErrorsDigisToCalibDigis_2DErrorInformation_.find(detId);

    if (digiIter->begin() != digiIter->end()) {
      if (mapIterator == SiPixelErrorsDigisToCalibDigis_2DErrorInformation_.end()) {
        //      std::cout << "This is the beginning of an error 2d histo booking: "<<std::endl;
        setDQMDirectory(detId);
        temp_ = bookDQMHistoPlaquetteSummary2D(detId, "SiPixelErrorsCalibDigis", "SiPixelErrorsDigisToCalibDigis");
        SiPixelErrorsDigisToCalibDigis_2DErrorInformation_.insert(std::make_pair(detId, temp_));
      } else {
        //      std::cout << "This one was already booked."<<std::endl;
        temp_ = (*mapIterator).second;
      }

      for (ipix = digiIter->begin(); ipix != digiIter->end(); ++ipix) {
        temp_->Fill(ipix->getCol(), ipix->getRow());
        //      std::cout << "detId: " << detId << " " << ipix->getRow() << " " << ipix->getCol() << std::endl;
      }

    }  // end of the if statement asking if the plaquette in question has any errors in it

  }  // end of the for loop that goes through all plaquettes
}

// ------------ method called once each job just before starting event loop  ------------
void SiPixelErrorsDigisToCalibDigis::beginJob() {}

// ------------ method called once each job just after ending the event loop  ------------
void SiPixelErrorsDigisToCalibDigis::endJob() {
  if (!outputFilename_.empty() && createOutputFile_) {
    edm::LogInfo("SiPixelErrorCalibDigis") << "Writing ROOT file to: " << outputFilename_ << std::endl;
    if (edm::Service<DQMStore>().isAvailable())
      edm::Service<DQMStore>()->save(outputFilename_);
  }
}

// ------------ helper functions ---------------------------------------------------------

SiPixelErrorsDigisToCalibDigis::MonitorElement* SiPixelErrorsDigisToCalibDigis::bookDQMHistogram2D(uint32_t detid,
                                                                                                   std::string name,
                                                                                                   std::string title,
                                                                                                   int nchX,
                                                                                                   double lowX,
                                                                                                   double highX,
                                                                                                   int nchY,
                                                                                                   double lowY,
                                                                                                   double highY) {
  std::string hid = theHistogramIdWorker_->setHistoId(name, detid);
  return daqBE_->book2D(hid, title, nchX, lowX, highX, nchY, lowY, highY);
}

SiPixelErrorsDigisToCalibDigis::MonitorElement* SiPixelErrorsDigisToCalibDigis::bookDQMHistoPlaquetteSummary2D(
    uint32_t detid, std::string name, std::string title) {
  DetId detId(detid);
  const TrackerGeometry& theTracker(*geom_);
  const PixelGeomDetUnit* theGeomDet = dynamic_cast<const PixelGeomDetUnit*>(theTracker.idToDet(detId));
  int maxcol = theGeomDet->specificTopology().ncolumns();
  int maxrow = theGeomDet->specificTopology().nrows();

  std::string hid = theHistogramIdWorker_->setHistoId(name, detid);
  return daqBE_->book2D(hid, title, maxcol, 0, maxcol, maxrow, 0, maxrow);
}

bool SiPixelErrorsDigisToCalibDigis::setDQMDirectory(std::string dirName) {
  daqBE_->setCurrentFolder(dirName);
  return daqBE_->dirExists(dirName);
}

bool SiPixelErrorsDigisToCalibDigis::setDQMDirectory(uint32_t detID) { return folderMaker_->setModuleFolder(detID, 0); }

// -- define this as a plug-in
DEFINE_FWK_MODULE(SiPixelErrorsDigisToCalibDigis);