HBHEDataFrame.h

Go to the documentation of this file.

#ifndef DIGIHCAL_HBHEDATAFRAME_H
#define DIGIHCAL_HBHEDATAFRAME_H
 
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/HcalDetId/interface/HcalElectronicsId.h"
#include "DataFormats/HcalDigi/interface/HcalQIESample.h"
#include <ostream>
 
class HBHEDataFrame {
public:
  typedef HcalDetId key_type;  
 
  constexpr HBHEDataFrame() : id_(0), size_(0), hcalPresamples_(0) {}
  constexpr explicit HBHEDataFrame(const HcalDetId& id) : id_(id), size_(0), hcalPresamples_(0) {
    // TODO : test id for HcalBarrel or HcalEndcap
  }
 
  constexpr const HcalDetId& id() const { return id_; }
  constexpr const HcalElectronicsId& elecId() const { return electronicsId_; }
 
  constexpr int size() const { return size_ & 0xF; }
  constexpr int presamples() const { return hcalPresamples_ & 0xF; }
  constexpr bool zsMarkAndPass() const { return (hcalPresamples_ & 0x10); }
  constexpr bool zsUnsuppressed() const { return (hcalPresamples_ & 0x20); }
  constexpr uint32_t zsCrossingMask() const { return (hcalPresamples_ & 0x3FF000) >> 12; }
 
  constexpr HcalQIESample const& operator[](int i) const { return data_[i]; }
  constexpr HcalQIESample const& sample(int i) const { return data_[i]; }
 
  constexpr int fiberIdleOffset() const {
    int val = (hcalPresamples_ & 0xF00) >> 8;
    return (val == 0) ? (-1000) : (((val & 0x8) == 0) ? (-(val & 0x7)) : (val & 0x7));
  }
 
  constexpr bool validate(int firstSample = 0, int nSamples = 100) const {
    int capid = -1;
    bool ok = true;
    for (int i = 0; ok && i < nSamples && i + firstSample < size_; i++) {
      if (data_[i + firstSample].er() || !data_[i + firstSample].dv())
        ok = false;
      if (i == 0)
        capid = data_[i + firstSample].capid();
      if (capid != data_[i + firstSample].capid())
        ok = false;
      capid = (capid + 1) % 4;
    }
    return ok;
  }
 
  constexpr void setSize(int size) {
    if (size > MAXSAMPLES)
      size_ = MAXSAMPLES;
    else if (size <= 0)
      size_ = 0;
    else
      size_ = size;
  }
  constexpr void setPresamples(int ps) { hcalPresamples_ |= ps & 0xF; }
  constexpr void setZSInfo(bool unsuppressed, bool markAndPass, uint32_t crossingMask = 0) {
    hcalPresamples_ &= 0x7FC00F0F;  // preserve actual presamples and fiber idle offset
    if (markAndPass)
      hcalPresamples_ |= 0x10;
    if (unsuppressed)
      hcalPresamples_ |= 0x20;
    hcalPresamples_ |= (crossingMask & 0x3FF) << 12;
  }
  constexpr void setSample(int i, const HcalQIESample& sam) { data_[i] = sam; }
  constexpr void setReadoutIds(const HcalElectronicsId& eid) { electronicsId_ = eid; }
  constexpr void setFiberIdleOffset(int offset) {
    hcalPresamples_ &= 0x7FFFF0FF;
    if (offset >= 7)
      hcalPresamples_ |= 0xF00;
    else if (offset >= 0)
      hcalPresamples_ |= (0x800) | (offset << 8);
    else if (offset >= -7)
      hcalPresamples_ |= ((-offset) << 8);
    else
      hcalPresamples_ |= 0x700;
  }
 
  static const int MAXSAMPLES = 10;
 
private:
  HcalDetId id_;
  HcalElectronicsId electronicsId_;
  int size_;
  int hcalPresamples_;
  HcalQIESample data_[MAXSAMPLES];
};
 
std::ostream& operator<<(std::ostream&, const HBHEDataFrame&);
 
#endif