/data/refman/pasoursint/CMSSW_5_3_3/src/DataFormats/L1GlobalCaloTrigger/src/L1GctJetCounts.cc

Go to the documentation of this file.

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

using std::vector;
using std::ostream;
using std::endl;

const unsigned L1GctJetCounts::MAX_TOTAL_COUNTS=12;
const unsigned L1GctJetCounts::MAX_TRUE_COUNTS=6;

// default constructor
L1GctJetCounts::L1GctJetCounts() :
  m_data0(0),
  m_data1(0),
  m_bx(0)
{

}

// constructor for unpacking
L1GctJetCounts::L1GctJetCounts(uint32_t data0, uint32_t data1) :
  m_data0(data0 & 0x7fff7fff), // Mask off bits 15 and 31 for better compression and consistency
  m_data1(data1 & 0x7fff7fff),  // with emulator constructor - these bits are not jet count data!
  m_bx(0)
{
}

// constructor for unpacking
L1GctJetCounts::L1GctJetCounts(uint32_t data0, uint32_t data1, int16_t bx) :
  m_data0(data0 & 0x7fff7fff), // Mask off bits 15 and 31 for better compression and consistency
  m_data1(data1 & 0x7fff7fff),  // with emulator constructor - these bits are not jet count data!
  m_bx(bx)
{
}

// constructor for emulator
L1GctJetCounts::L1GctJetCounts(vector<unsigned> counts) :
  m_data0(0), 
  m_data1(0),
  m_bx(0)
{
  // Assumes all required output data has been packed 
  // into 12 5-bit fields
  if (counts.size() != MAX_TOTAL_COUNTS) { }
  else {
    for (unsigned int i=0; i<3; ++i) {
      m_data0 += (counts[i]   << (5*i));
      m_data0 += (counts[i+3] << (5*i + 16));
      m_data1 += (counts[i+6] << (5*i));
      m_data1 += (counts[i+9] << (5*i + 16));
    }
  }
}

// constructor for emulator
L1GctJetCounts::L1GctJetCounts(vector<unsigned> counts, int16_t bx) :
  m_data0(0), 
  m_data1(0),
  m_bx(bx)
{
  if (counts.size() != MAX_TOTAL_COUNTS) { }
  else {
    for (unsigned int i=0; i<3; ++i) {
      m_data0 += (counts[i]   << (5*i));
      m_data0 += (counts[i+3] << (5*i + 16));
      m_data1 += (counts[i+6] << (5*i));
      m_data1 += (counts[i+9] << (5*i + 16));
    }
  }
}

// destructor
L1GctJetCounts::~L1GctJetCounts()
{

}

// return counts by index
unsigned L1GctJetCounts::count(unsigned i) const
{
  if (i<6){ return ((m_data0 >> (i<3 ? (5*i) : ((5*i)+1))) & 0x1f); }
  else if (i < MAX_TOTAL_COUNTS) { return ((m_data1 >> (i<9 ? ((5*i)-30) : ((5*i)-29))) & 0x1f); }    
  else { return 0; }
}

// pretty print
ostream& operator<<(ostream& s, const L1GctJetCounts& c) {
  s << "L1GctJetCounts : ";
  for (unsigned int i=0 ; i<12 ; ++i) {
    s << "\n     count " << i<< "=" << c.count(i);
  }
  return s;
}