L1GctInternJetData.cc

Go to the documentation of this file.

#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctInternJetData.h"

L1GctInternJetData::L1GctInternJetData() : capBlock_(0), capIndex_(0), bx_(0), type_(null), data_(0) {}

L1GctInternJetData::L1GctInternJetData(L1CaloRegionDetId rgn,
                                       uint16_t capBlock,
                                       uint16_t capIndex,
                                       int16_t bx,
                                       uint8_t sgnEta,
                                       uint8_t oflow,
                                       uint16_t et,
                                       uint8_t eta,
                                       uint8_t phi,
                                       uint8_t tauVeto,
                                       uint8_t rank)
    : regionId_(rgn), capBlock_(capBlock), capIndex_(capIndex), bx_(bx), type_(null), data_(0) {
  setData(sgnEta, oflow, et, eta, phi, tauVeto, rank);
}

// 'named' constructors to avoid confusion

// emulator calibrated jet ctor
L1GctInternJetData L1GctInternJetData::fromEmulator(
    L1CaloRegionDetId rgn, int16_t bx, uint16_t et, bool oflow, bool tauVeto, uint8_t eta, uint8_t phi, uint16_t rank) {
  L1GctInternJetData d;

  d.setType(emulator);
  d.setRegionId(rgn);
  d.setData(0, (oflow ? 1 : 0), et, eta, phi, (tauVeto ? 1 : 0), rank);
  d.setBx(bx);

  return d;
}

L1GctInternJetData L1GctInternJetData::fromJetCluster(
    L1CaloRegionDetId rgn, uint16_t capBlock, uint16_t capIndex, int16_t bx, uint32_t data) {
  L1GctInternJetData d;

  d.setRegionId(rgn);
  d.setCapBlock(capBlock);
  d.setCapIndex(capIndex);
  d.setBx(bx);
  d.setType(jet_cluster);
  d.setRawData(data & 0x1fffffff);

  return d;
}

L1GctInternJetData L1GctInternJetData::fromJetPreCluster(
    L1CaloRegionDetId rgn, uint16_t capBlock, uint16_t capIndex, int16_t bx, uint32_t data) {
  L1GctInternJetData d;

  d.setRegionId(rgn);
  d.setCapBlock(capBlock);
  d.setCapIndex(capIndex);
  d.setBx(bx);
  d.setType(jet_precluster);
  d.setData(0, (data >> 10) & 0x1, data & 0x3ff, (data >> 12) & 0xf, 0, (data >> 11) & 0x1, 0);

  return d;
}

L1GctInternJetData L1GctInternJetData::fromJetClusterMinimal(
    L1CaloRegionDetId rgn, uint16_t capBlock, uint16_t capIndex, int16_t bx, uint32_t data) {
  L1GctInternJetData d;

  d.setRegionId(rgn);
  d.setCapBlock(capBlock);
  d.setCapIndex(capIndex);
  d.setBx(bx);
  d.setType(jet_cluster_minimal);
  d.setData(0, (data >> 10) & 0x1, 0, (data >> 6) & 0xf, (data >> 12) & 0x1, (data >> 11) & 0x1, data & 0x3f);

  return d;
}

L1GctInternJetData L1GctInternJetData::fromGctTrigObject(
    L1CaloRegionDetId rgn, uint16_t capBlock, uint16_t capIndex, int16_t bx, uint32_t data) {
  L1GctInternJetData d;

  d.setRegionId(rgn);
  d.setCapBlock(capBlock);
  d.setCapIndex(capIndex);
  d.setBx(bx);
  d.setType(gct_trig_object);
  d.setData((data >> 10) & 0x1, 0, 0, (data >> 6) & 0xf, (data >> 11) & 0x1f, 0, data & 0x3f);

  return d;
}

L1GctInternJetData::~L1GctInternJetData() {}

void L1GctInternJetData::setData(
    uint8_t sgnEta, uint8_t oflow, uint16_t et, uint8_t eta, uint8_t phi, uint8_t tauVeto, uint8_t rank) {
  data_ = 0;
  data_ = rank & 0x3f;
  data_ |= (tauVeto & 0x1) << 6;
  data_ |= (phi & 0x1f) << 7;
  data_ |= (eta & 0xf) << 12;
  data_ |= (et & 0xfff) << 16;
  data_ |= (oflow & 0x1) << 28;
  data_ |= (sgnEta & 0x1) << 29;
}


bool L1GctInternJetData::operator==(const L1GctInternJetData& c) const {
  return (type_ == c.type() && data_ == c.raw() && regionId_ == c.regionId() && bx_ == c.bx());
}

std::ostream& operator<<(std::ostream& s, const L1GctInternJetData& c) {
  s << "L1GctInternJetData :";
  if (c.empty()) {
    s << " empty!";
  }
  if (c.type() == L1GctInternJetData::jet_cluster) {
    s << " type=jet_cluster";
    s << " oflow=" << c.oflow();
    s << " et=" << c.et();
    s << " eta=" << c.eta();
    s << " phi=" << c.phi();
    s << " tauVeto=" << c.tauVeto();
    s << " rank=" << c.rank();
  } else if (c.type() == L1GctInternJetData::jet_precluster) {
    s << " type=jet_precluster";
    s << " oflow=" << c.oflow();
    s << " et=" << c.et();
    s << " eta=" << c.eta();
    s << " tauVeto=" << c.tauVeto();
  } else if (c.type() == L1GctInternJetData::jet_cluster_minimal) {
    s << " type=jet_cluster_minimal";
    s << " oflow=" << c.oflow();
    s << " eta=" << c.eta();
    s << " phi=" << c.phi();
    s << " tauVeto=" << c.tauVeto();
    s << " rank=" << c.rank();
  } else if (c.type() == L1GctInternJetData::gct_trig_object) {
    s << " type=gct_trig_object";
    s << " eta=" << c.eta();
    s << " phi=" << c.phi();
    s << " rank=" << c.rank();
  }
  s << " cap block=" << std::hex << c.capBlock();
  s << " index=" << std::dec << c.capIndex();
  s << " BX=" << c.bx();

  return s;
}