SiPixelPhase1RawData.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     SiPixelPhase1RawData
// Class:       SiPixelPhase1RawData
//
 
// Original Author: Marcel Schneider
 
#include "DQM/SiPixelPhase1Common/interface/SiPixelPhase1Base.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiPixelRawData/interface/SiPixelRawDataError.h"
#include "FWCore/Framework/interface/MakerMacros.h"
 
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
 
namespace {
 
  class SiPixelPhase1RawData final : public SiPixelPhase1Base {
    enum { NERRORS, FIFOFULL, TBMMESSAGE, TBMTYPE, TYPE_NERRORS };
 
  public:
    explicit SiPixelPhase1RawData(const edm::ParameterSet& conf);
    void analyze(const edm::Event&, const edm::EventSetup&) override;
 
  private:
    edm::EDGetTokenT<edm::DetSetVector<SiPixelRawDataError>> srcToken_;
  };
 
  SiPixelPhase1RawData::SiPixelPhase1RawData(const edm::ParameterSet& iConfig) : SiPixelPhase1Base(iConfig) {
    srcToken_ = consumes<edm::DetSetVector<SiPixelRawDataError>>(iConfig.getParameter<edm::InputTag>("src"));
  }
 
  void SiPixelPhase1RawData::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
    if (!checktrigger(iEvent, iSetup, DCS))
      return;
 
    edm::Handle<edm::DetSetVector<SiPixelRawDataError>> input;
    iEvent.getByToken(srcToken_, input);
    if (!input.isValid())
      return;
 
    for (auto it = input->begin(); it != input->end(); ++it) {
      for (auto& siPixelRawDataError : *it) {
        int fed = siPixelRawDataError.getFedId();
        int type = siPixelRawDataError.getType();
        DetId id = it->detId();
 
        // encoding of the channel number within the FED error word
        const uint32_t LINK_bits = 6;
        const uint32_t LINK_shift = 26;
        const uint64_t LINK_mask = (1 << LINK_bits) - 1;
 
        uint64_t errorWord = 0;
        // use 64bit word for some error types
        // invalid header, invalid trailer, size mismatch
        if (type == 32 || type == 33 || type == 34) {
          errorWord = siPixelRawDataError.getWord64();
        } else {
          errorWord = siPixelRawDataError.getWord32();
        }
 
        int32_t chanNmbr = (errorWord >> LINK_shift) & LINK_mask;
        // timeout
        if (type == 29)
          chanNmbr = -1;  // TODO: different formula needed.
 
        uint32_t error_data = errorWord & 0xFF;
 
        if (type == 28) {  // overflow.
          for (uint32_t i = 0; i < 8; i++) {
            if (error_data & (1 << i))
              histo[FIFOFULL].fill(i, id, &iEvent, fed, chanNmbr);
          }
        }
 
        if (type == 30) {                                      // TBM stuff.
          uint32_t statemachine_state = errorWord >> 8 & 0xF;  // next 4 bits after data
          const uint32_t tbm_types[16] = {0, 1, 2, 4, 2, 4, 2, 4, 3, 1, 4, 4, 4, 4, 4, 4};
 
          histo[TBMTYPE].fill(tbm_types[statemachine_state], id, &iEvent, fed, chanNmbr);
 
          for (uint32_t i = 0; i < 8; i++) {
            if (error_data & (1 << i))
              histo[TBMMESSAGE].fill(i, id, &iEvent, fed, chanNmbr);
          }
          continue;  // we don't really consider these as errors.
        }
 
        // note that a DetId of 0xFFFFFFFF can mean 'no DetId'.
        // We hijack column and row for FED and chan in this case,
        // the GeometryInterface does understand that.
 
        histo[NERRORS].fill(id, &iEvent, fed, chanNmbr);
        histo[TYPE_NERRORS].fill(type, id, &iEvent, fed, chanNmbr);
      }
    }
 
    histo[NERRORS].executePerEventHarvesting(&iEvent);
  }
 
}  //namespace
 
DEFINE_FWK_MODULE(SiPixelPhase1RawData);