HFRecalibration.cc

Go to the documentation of this file.

// File: HFRecalibration.cc
// Description: simple helper class containing parameterized 
//              function for HF damade recovery for Upgrade studies  
//              evaluated using SimG4CMS/Calo/ HFDarkening   

#include "DataFormats/HcalCalibObjects/interface/HFRecalibration.h"

// CMSSW Headers
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace edm;

HFRecalibration::HFRecalibration(const edm::ParameterSet& pset)
{
  //HFParsAB[Depth=0/1][A=0/B=1]
  HFParsAB[0][0] = pset.getParameter<vecOfDoubles>("HFdepthOneParameterA");
  HFParsAB[0][1] = pset.getParameter<vecOfDoubles>("HFdepthOneParameterB");
  HFParsAB[1][0] = pset.getParameter<vecOfDoubles>("HFdepthTwoParameterA");
  HFParsAB[1][1] = pset.getParameter<vecOfDoubles>("HFdepthTwoParameterB");
}

HFRecalibration::~HFRecalibration()
{ }

double HFRecalibration::getCorr(int ieta, int depth, double lumi)
{

  // parameterizations provided by James Wetzel
 
  ieta = (abs(ieta) - loweriEtaBin);
  
  if (ieta < 0 || ieta > 11 || depth < 1 || depth > 2)
  {
    return 1.0;
  }

  switch (depth)
  {
    case 1:
      reCalFactor = (1 + HFParsAB[0][0][ieta] * sqrt(lumi) + HFParsAB[0][1][ieta] * lumi);
      
    case 2:
      reCalFactor = (1 + HFParsAB[1][0][ieta] * sqrt(lumi) + HFParsAB[1][1][ieta] * lumi);
  }

  return reCalFactor;
}