/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/CondFormats/L1TObjects/src/L1CaloEtScale.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:     L1TObjects
// Class  :     L1CaloEtScale
// 
// Implementation:
//     <Notes on implementation>
//
// Author:      
// Created:     Wed Sep 27 17:18:27 CEST 2006
// $Id: 

#include "CondFormats/L1TObjects/interface/L1CaloEtScale.h"

#include "FWCore/Utilities/interface/Exception.h"
#include <stdexcept>

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

// default constructor (testing only!)
L1CaloEtScale::L1CaloEtScale() :
  m_linScaleMax(0x3ff),
  m_rankScaleMax(0x3f),
  m_linearLsb(1.0)
{
  for (unsigned i=0; i<m_rankScaleMax; i++) {
    m_thresholds.push_back(m_linearLsb * i);
  }
}

// ctor that provides backwards compatibility with fixed max scale values
// OK to use this with e/gamma and rank scales
L1CaloEtScale::L1CaloEtScale(const double linearLsbInGeV, const vector<double> thresholdsInGeV) :
  m_linScaleMax(0x3ff),
  m_rankScaleMax(0x3f),
  m_linearLsb(linearLsbInGeV),
  m_thresholds(thresholdsInGeV) {

  // protect against too many thresholds!
  //  while ( m_threshold.size() > (L1GctJetScale::maxRank+1) ) {
  //    m_thresholds.pop_back();
  //  }

}


// ctor that sets scale max values
L1CaloEtScale::L1CaloEtScale(const unsigned linScaleMax, const unsigned rankScaleMax, const double linearLsbInGeV, const vector<double> thresholdsInGeV) :
  m_linScaleMax(linScaleMax),
  m_rankScaleMax(rankScaleMax),
  m_linearLsb(linearLsbInGeV),
  m_thresholds(thresholdsInGeV) {

}


L1CaloEtScale::~L1CaloEtScale() {

}

// convert from linear Et to rank
uint16_t L1CaloEtScale::rank(const uint16_t linear) const {

  return rank( (linear & m_linScaleMax) * m_linearLsb);

}

uint16_t L1CaloEtScale::rank(const double EtInGeV) const {

  uint16_t out = 0;

  for (unsigned i=0; i<m_thresholds.size() && i<(unsigned)(m_rankScaleMax+1); i++) {
    if ( EtInGeV >= m_thresholds.at(i) ) { out = i; }
  }

  return out & m_rankScaleMax;
}

// convert from rank to Et/GeV
double L1CaloEtScale::et(const uint16_t rank) const {

  // return bin centre, except for highest bin
//   if (rank < m_thresholds.size()-1) {
//     return (m_thresholds[rank+1]+m_thresholds[rank]) / 2;
//   }
//   else {
//     return m_thresholds.back();
//   }

// return bin lower edge
  try {
    return m_thresholds.at(rank);
  }
  catch(std::out_of_range) {
    throw cms::Exception("OutOfRange") << "Index out of range in L1CaloEtScale::et(rank)" << std::endl;
  }

}

void L1CaloEtScale::print(ostream& s) const {
  s << "L1CaloEtScale :" << endl;
  s << "  Input scale max = " << m_linScaleMax << endl;
  s << "  Input LSB       = " << m_linearLsb << " GeV" << endl;
  s << "  Rank scale max  = " << m_rankScaleMax << endl;
  for (unsigned i=0; i<m_thresholds.size(); i++) {
    s << "  Threshold " << i << " = " << m_thresholds[i] << " GeV" << endl;
  }
}

std::ostream& operator << (std::ostream& os, const L1CaloEtScale obj) {
  os << "L1CaloEtScale :" << endl;
  os << "  Input scale max = " << obj.linScaleMax() << endl;
  os << "  Input LSB       = " << obj.linearLsb() << " GeV" << endl;
  os << "  Rank scale max  = " << obj.rankScaleMax() << endl;
  for (unsigned i=0; i<obj.getThresholds().size(); i++) {
    os << "  Threshold " << i << " = " << obj.getThresholds().at(i) << " GeV" << endl;
  }
  return os ;
}