HFValueStruct.cc

Go to the documentation of this file.

#include "HFValueStruct.h"
#include <cstdio>
#include <cstdlib>
//version -1 will take information from DB (NOT DONE YET)
//version 0 has energy corrections on, everything else off
//version 1 has energy correction, and pile up slope and interceept on
// version 2+ will use defaults of 'do nothing'
//
//
//

reco::HFValueStruct::HFValueStruct(const int& version, const std::vector<double>& vect) : v_(version), hfvv_(vect) {
  //if(v_==-1) hfvv_=SetHfvvFromDB_();
  //v==99 will always give defaults

  //version control, add in versions as they appear!!
  if (v_ == 0 || v_ == 1)
    doEnCor_ = true;
  else
    doEnCor_ = false;

  if (v_ == 1)
    doPU_ = true;
  else
    doPU_ = false;
}

int reco::HFValueStruct::indexByIeta(int& ieta) const { return (ieta > 0) ? (abs(ieta) - 29 + 13) : (41 - abs(ieta)); }
int reco::HFValueStruct::ietaByIndex(int& indx) const { return (indx > 13) ? (indx + 29 - 13) : (indx - 41); }
//version 0
// EnCor=energy corrections,default 1.0, 26 slots

//version 1
//PUSlope= slope for pile up corrections,default 0.0, 26 slots, 0,1,12,13,24,25 are all defaults always

//PUIntercept= intercept  slope for pile up corrections,default 1.0, 26 slots, 0,1,12,13,24,25 are all defaults always

// returns single value by index

double reco::HFValueStruct::EnCor(int ieta) const {
  int indx = indexByIeta(ieta);
  if (doEnCor_)
    return hfvv_[indx];
  else
    return 1.0;
}
double reco::HFValueStruct::PUSlope(int ieta) const {
  int indx = indexByIeta(ieta) + 26;
  if (doPU_)
    return hfvv_[indx];
  else
    return 0.0;
}
double reco::HFValueStruct::PUIntercept(int ieta) const {
  int indx = indexByIeta(ieta) + 52;
  if (doPU_)
    return hfvv_[indx];
  else
    return 1.0;
}

// sets single value by index
void reco::HFValueStruct::setEnCor(int ieta, double val) {
  int indx = indexByIeta(ieta);
  hfvv_[indx] = val;
}
void reco::HFValueStruct::setPUSlope(int ieta, double val) {
  int indx = indexByIeta(ieta) + 26;
  hfvv_[indx] = val;
}
void reco::HFValueStruct::setPUIntercept(int ieta, double val) {
  int indx = indexByIeta(ieta) + 52;
  hfvv_[indx] = val;
}

// returns whole vector
std::vector<double> reco::HFValueStruct::EnCor() const {
  std::vector<double> vct;
  if (doEnCor_) {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(hfvv_[ii]);
  } else {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(1.0);
  }
  return vct;
}

std::vector<double> reco::HFValueStruct::PUSlope() const {
  std::vector<double> vct;
  if (doPU_) {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(hfvv_[ii + 26]);
  } else {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(0.0);
  }
  return vct;
}

std::vector<double> reco::HFValueStruct::PUIntercept() const {
  std::vector<double> vct;
  if (doPU_) {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(hfvv_[ii + 52]);
  } else {
    for (int ii = 0; ii < 13; ii++)
      vct.push_back(1.0);
  }
  return vct;
}

// set whole vector
void reco::HFValueStruct::setEnCor(const std::vector<double>& val) {
  for (int ii = 0; ii < 13; ii++)
    hfvv_[ii] = val[ii];
}
void reco::HFValueStruct::setPUSlope(const std::vector<double>& val) {
  for (int ii = 0; ii < 13; ii++)
    hfvv_[ii + 26] = val[ii];
}
void reco::HFValueStruct::setPUIntercept(const std::vector<double>& val) {
  for (int ii = 0; ii < 13; ii++)
    hfvv_[ii + 52] = val[ii];
}