SiPixelPhase1Digis.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiPixelPhase1Digis
// Class:      SiPixelPhase1Digis
//

// Original Author: Marcel Schneider

// C++ stuff
#include <iostream>

// CMSSW stuff
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

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

// Input data stuff
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiPixelDigi/interface/PixelDigi.h"

// PixelDQM Framework
#include "DQM/SiPixelPhase1Common/interface/SiPixelPhase1Base.h"

namespace {

  class SiPixelPhase1Digis final : public SiPixelPhase1Base {
    // List of quantities to be plotted.
    enum {
      ADC,              // digi ADC readouts
      NDIGIS,           // number of digis per event and module
      NDIGISINCLUSIVE,  //Total number of digis in BPix and FPix
      NDIGIS_FED,       // number of digis per event and FED
      NDIGIS_FEDtrend,  // number of digis per event and FED
      EVENT,            // event frequency
      MAP,              // digi hitmap per module
      OCCUPANCY,        // like map but coarser

      MAX_HIST  // a sentinel that gives the number of quantities (not a plot).
    };

  public:
    explicit SiPixelPhase1Digis(const edm::ParameterSet& conf);

    void analyze(const edm::Event&, const edm::EventSetup&) override;

  private:
    edm::EDGetTokenT<edm::DetSetVector<PixelDigi>> srcToken_;
  };

  SiPixelPhase1Digis::SiPixelPhase1Digis(const edm::ParameterSet& iConfig) : SiPixelPhase1Base(iConfig) {
    srcToken_ = consumes<edm::DetSetVector<PixelDigi>>(iConfig.getParameter<edm::InputTag>("src"));
  }

  void SiPixelPhase1Digis::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
    if (!checktrigger(iEvent, iSetup, DCS))
      return;

    edm::Handle<edm::DetSetVector<PixelDigi>> input;
    iEvent.getByToken(srcToken_, input);
    if (!input.isValid())
      return;
    bool hasDigis = false;

    edm::DetSetVector<PixelDigi>::const_iterator it;
    for (it = input->begin(); it != input->end(); ++it) {
      for (PixelDigi const& digi : *it) {
        hasDigis = true;
        histo[ADC].fill((double)digi.adc(), DetId(it->detId()), &iEvent, digi.column(), digi.row());
        histo[MAP].fill(DetId(it->detId()), &iEvent, digi.column(), digi.row());
        histo[OCCUPANCY].fill(DetId(it->detId()), &iEvent, digi.column(), digi.row());
        histo[NDIGIS].fill(DetId(it->detId()), &iEvent);           // count
        histo[NDIGISINCLUSIVE].fill(DetId(it->detId()), &iEvent);  // count
        histo[NDIGIS_FED].fill(DetId(it->detId()), &iEvent);
        histo[NDIGIS_FEDtrend].fill(DetId(it->detId()), &iEvent);
      }
    }
    if (hasDigis)
      histo[EVENT].fill(DetId(0), &iEvent);
    histo[NDIGIS].executePerEventHarvesting(&iEvent);
    histo[NDIGISINCLUSIVE].executePerEventHarvesting(&iEvent);
    histo[NDIGIS_FED].executePerEventHarvesting(&iEvent);
    histo[NDIGIS_FEDtrend].executePerEventHarvesting(&iEvent);
  }

}  //namespace

DEFINE_FWK_MODULE(SiPixelPhase1Digis);