EnergyScaleCorrection_class Class Reference

Read and get energy scale and smearings from .dat files. More...

#include <Calibration/ZFitter/src/EnergyScaleCorrection_class.cc>

Public Types
enum	fileFormat_t { UNKNOWN = 0, GLOBE, ECALELF_TOY, ECALELF }
enum	paramSmear_t { kNone = 0, kRho, kPhi, kNParamSmear }

Public Member Functions
	EnergyScaleCorrection_class (std::string correctionFileName, unsigned int genSeed=0)
	EnergyScaleCorrection_class ()
float	getSmearingRho (int runNumber, bool isEBEle, float R9Ele, float etaSCEle, float EtEle) const
	public for sigmaE estimate More...
float	getSmearingSigma (int runNumber, bool isEBEle, float R9Ele, float etaSCEle, float EtEle, paramSmear_t par, float nSigma=0.) const
float	getSmearingSigma (int runNumber, bool isEBEle, float R9Ele, float etaSCEle, float EtEle, float nSigma_rho, float nSigma_phi) const
float	ScaleCorrection (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
	method to get energy scale corrections More...
float	ScaleCorrectionUncertainty (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
	method to get scale correction uncertainties: it's stat+syst in eta x R9 categories More...
void	SetSmearingType (fileFormat_t value)
	~EnergyScaleCorrection_class (void)
	dummy constructor needed in ElectronEnergyCalibratorRun2 More...

Public Attributes
bool	doScale
bool	doSmearings

Private Member Functions
void	AddScale (TString category_, int runMin_, int runMax_, double deltaP_, double err_deltaP_, double err_syst_deltaP)
void	AddSmearing (TString category_, int runMin_, int runMax_, double rho, double err_rho, double phi, double err_phi, double Emean, double err_Emean)
correctionValue_class	getScaleCorrection (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
	returns the correction value class More...
float	getScaleOffset (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
float	getScaleStatUncertainty (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
float	getScaleSystUncertainty (unsigned int runNumber, bool isEBEle, double R9Ele, double etaSCEle, double EtEle) const
void	ReadFromFile (TString filename)
	category "runNumber" runMin runMax deltaP err_deltaP_per_bin err_deltaP_stat err_deltaP_syst More...
void	ReadSmearingFromFile (TString filename)
	File structure: category constTerm alpha;. More...

Private Attributes
correction_map_t	scales
correction_map_t	scales_not_defined
correction_map_t	smearings
correction_map_t	smearings_not_defined
fileFormat_t	smearingType_

Detailed Description

Read and get energy scale and smearings from .dat files.

Author

Shervin Nourbakhsh

Definition at line 118 of file EnergyScaleCorrection_class.h.

Member Enumeration Documentation

fileFormat_t

enum EnergyScaleCorrection_class::fileFormat_t

Enumerator
UNKNOWN
GLOBE
ECALELF_TOY
ECALELF

Definition at line 120 of file EnergyScaleCorrection_class.h.

{ UNKNOWN = 0, GLOBE, ECALELF_TOY, ECALELF };

paramSmear_t

enum EnergyScaleCorrection_class::paramSmear_t

Enumerator
kNone
kRho
kPhi
kNParamSmear

Definition at line 122 of file EnergyScaleCorrection_class.h.

{ kNone = 0, kRho, kPhi, kNParamSmear };

Constructor & Destructor Documentation

EnergyScaleCorrection_class() [1/2]

EnergyScaleCorrection_class::EnergyScaleCorrection_class	(	std::string	correctionFileName,
		unsigned int	genSeed = 0
	)

Definition at line 19 of file EnergyScaleCorrection_class.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, beamvalidation::exit(), corrVsCorr::filename, ReadFromFile(), ReadSmearingFromFile(), scales, smearings, and AlCaHLTBitMon_QueryRunRegistry::string.

    : doScale(false), doSmearings(false), smearingType_(ECALELF) {
  if (!correctionFileName.empty()) {
    std::string filename = correctionFileName + "_scales.dat";
    ReadFromFile(filename);
    if (scales.empty()) {
      std::cerr << "[ERROR] scale correction map empty" << std::endl;
      exit(1);
    }
  }

  if (!correctionFileName.empty()) {
    std::string filename = correctionFileName + "_smearings.dat";
    ReadSmearingFromFile(filename);
    if (smearings.empty()) {
      std::cerr << "[ERROR] smearing correction map empty" << std::endl;
      exit(1);
    }
  }

  return;
}

EnergyScaleCorrection_class() [2/2]

EnergyScaleCorrection_class::EnergyScaleCorrection_class ( )

inline

Definition at line 128 of file EnergyScaleCorrection_class.h.

{};  

~EnergyScaleCorrection_class()

EnergyScaleCorrection_class::~EnergyScaleCorrection_class

(

void

)

dummy constructor needed in ElectronEnergyCalibratorRun2

Definition at line 42 of file EnergyScaleCorrection_class.cc.

{ return; }

Member Function Documentation

AddScale()

void EnergyScaleCorrection_class::AddScale ( TString  category_, int  runMin_, int  runMax_, double  deltaP_, double  err_deltaP_, double  err_syst_deltaP  )

private

Definition at line 175 of file EnergyScaleCorrection_class.cc.

References eostools::cat(), EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, alignCSCRings::corr, gather_cfg::cout, beamvalidation::exit(), and scales.

Referenced by ReadFromFile().

                                                                                                             {
  correctionCategory_class cat(category_);  // build the category from the string
  cat.runmin = runMin_;
  cat.runmax = runMax_;

  // the following check is not needed, can be removed
  if (scales.count(cat) != 0) {
    std::cerr << "[ERROR] Category already defined!" << std::endl;
    std::cerr << "        Adding category:  " << cat << std::endl;
    std::cerr << "        Defined category: " << scales[cat] << std::endl;
    exit(1);
  }

  correctionValue_class corr;  // define the correction values
  corr.scale = deltaP_;
  corr.scale_err = err_deltaP_;
  corr.scale_err_syst = err_syst_deltaP;
  scales[cat] = corr;

#ifdef PEDANTIC_OUTPUT
  std::cout << "[INFO:scale correction] " << cat << corr << std::endl;
#endif
  return;
}

AddSmearing()

void EnergyScaleCorrection_class::AddSmearing ( TString  category_, int  runMin_, int  runMax_, double  rho, double  err_rho, double  phi, double  err_phi, double  Emean, double  err_Emean  )

private

Definition at line 202 of file EnergyScaleCorrection_class.cc.

References eostools::cat(), EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, alignCSCRings::corr, gather_cfg::cout, beamvalidation::exit(), phi, rho, and smearings.

Referenced by ReadSmearingFromFile().

                                                                {
  correctionCategory_class cat(category_);
  cat.runmin = (runMin_ < 0) ? 0 : runMin_;
  cat.runmax = runMax_;

  if (smearings.count(cat) != 0) {
    std::cerr << "[ERROR] Smearing category already defined!" << std::endl;
    std::cerr << "        Adding category:  " << cat << std::endl;
    std::cerr << "        Defined category: " << smearings[cat] << std::endl;
    exit(1);
  }

  correctionValue_class corr;
  corr.rho = rho;
  corr.rho_err = err_rho;
  corr.phi = phi;
  corr.phi_err = err_phi;
  corr.Emean = Emean;
  corr.Emean_err = err_Emean;
  smearings[cat] = corr;

#ifdef PEDANTIC_OUTPUT
#ifndef CMSSW
  std::cout << "[INFO:smearings] " << cat << corr << std::endl;
#else
  edm::LogInfo("[INFO:smearings] ") << cat << corr;
#endif
#endif

  return;
}

getScaleCorrection()

correctionValue_class EnergyScaleCorrection_class::getScaleCorrection ( unsigned int  runNumber, bool  isEBEle, double  R9Ele, double  etaSCEle, double  EtEle  ) const

private

returns the correction value class

Definition at line 71 of file EnergyScaleCorrection_class.cc.

References taus_updatedMVAIds_cff::category, gather_cfg::cout, convertSQLiteXML::runNumber, and scales.

Referenced by getScaleStatUncertainty(), and getScaleSystUncertainty().

                                                                                             {
  // buld the category based on the values of the object
  correctionCategory_class category(runNumber, etaSCEle, R9Ele, EtEle);
  correction_map_t::const_iterator corr_itr =
      scales.find(category);  // find the correction value in the map that associates the category to the correction

  if (corr_itr == scales.end()) {  // if not in the standard classes, add it in the list of not defined classes

    // this part is commented because it makes the method not constant
    // if(scales_not_defined.count(category) == 0) {
    //  correctionValue_class corr;
    //  scales_not_defined[category] = corr;
    // }
    // corr_itr = scales_not_defined.find(category);
    std::cout << "[ERROR] Scale category not found: " << std::endl;
    std::cout << category << std::endl;
    std::cout << "Returning uncorrected value." << std::endl;
    //     exit(1);
    correctionValue_class nocorr;
    std::cout << nocorr << std::endl;
    return nocorr;
  }

#ifdef DEBUG
  std::cout << "[DEBUG] Checking scale correction for category: " << category << std::endl;
  std::cout << "[DEBUG] Correction is: " << corr_itr->second << std::endl
            << "        given for category " << corr_itr->first << std::endl;
  ;
#endif
  return corr_itr->second;
}

getScaleOffset()

float EnergyScaleCorrection_class::getScaleOffset ( unsigned int  runNumber, bool  isEBEle, double  R9Ele, double  etaSCEle, double  EtEle  ) const

private

Definition at line 105 of file EnergyScaleCorrection_class.cc.

References taus_updatedMVAIds_cff::category, gather_cfg::cout, convertSQLiteXML::runNumber, correctionValue_class::scale, and scales.

Referenced by ScaleCorrection().

                                                                                             {
  // buld the category based on the values of the object
  correctionCategory_class category(runNumber, etaSCEle, R9Ele, EtEle);
  correction_map_t::const_iterator corr_itr =
      scales.find(category);  // find the correction value in the map that associates the category to the correction

  if (corr_itr == scales.end()) {  // if not in the standard classes, add it in the list of not defined classes

    // this part is commented because it makes the method not constant
    // if(scales_not_defined.count(category) == 0) {
    //  correctionValue_class corr;
    //  scales_not_defined[category] = corr;
    // }
    // corr_itr = scales_not_defined.find(category);
    std::cout << "[ERROR] Scale offset category not found: " << std::endl;
    std::cout << category << std::endl;
    std::cout << "Returning uncorrected value." << std::endl;
    //     exit(1);
    correctionValue_class nocorr;
    return nocorr.scale;
  }

#ifdef DEBUG
  std::cout << "[DEBUG] Checking scale offset correction for category: " << category << std::endl;
  std::cout << "[DEBUG] Correction is: " << corr_itr->second << std::endl
            << "        given for category " << corr_itr->first << std::endl;
  ;
#endif

  return corr_itr->second.scale;
}

getScaleStatUncertainty()

float EnergyScaleCorrection_class::getScaleStatUncertainty ( unsigned int  runNumber, bool  isEBEle, double  R9Ele, double  etaSCEle, double  EtEle  ) const

private

Definition at line 61 of file EnergyScaleCorrection_class.cc.

References getScaleCorrection(), convertSQLiteXML::runNumber, and correctionValue_class::scale_err.

Referenced by ScaleCorrectionUncertainty().

                                                                                             {
  return getScaleCorrection(runNumber, isEBEle, R9Ele, etaSCEle, EtEle).scale_err;
}

getScaleSystUncertainty()

float EnergyScaleCorrection_class::getScaleSystUncertainty ( unsigned int  runNumber, bool  isEBEle, double  R9Ele, double  etaSCEle, double  EtEle  ) const

private

Definition at line 66 of file EnergyScaleCorrection_class.cc.

References getScaleCorrection(), convertSQLiteXML::runNumber, and correctionValue_class::scale_err_syst.

Referenced by ScaleCorrectionUncertainty().

                                                                                             {
  return getScaleCorrection(runNumber, isEBEle, R9Ele, etaSCEle, EtEle).scale_err_syst;
}

getSmearingRho()

float EnergyScaleCorrection_class::getSmearingRho	(	int	runNumber,
		bool	isEBEle,
		float	R9Ele,
		float	etaSCEle,
		float	EtEle
	)		const

public for sigmaE estimate

Definition at line 376 of file EnergyScaleCorrection_class.cc.

References taus_updatedMVAIds_cff::category, convertSQLiteXML::runNumber, smearings, and smearings_not_defined.

                                                                                 {
  correctionCategory_class category(runNumber, etaSCEle, R9Ele, EtEle);
  correction_map_t::const_iterator corr_itr = smearings.find(category);
  if (corr_itr == smearings.end()) {  // if not in the standard classes, add it in the list of not defined classes
    // if(smearings_not_defined.count(category) == 0) {
    //  correctionValue_class corr;
    //  smearings_not_defined[category] = corr;
    // }
    corr_itr = smearings_not_defined.find(category);
  }

  return corr_itr->second.rho;
}

getSmearingSigma() [1/2]

float EnergyScaleCorrection_class::getSmearingSigma	(	int	runNumber,
		bool	isEBEle,
		float	R9Ele,
		float	etaSCEle,
		float	EtEle,
		paramSmear_t	par,
		float	nSigma = 0.
	)		const

Definition at line 336 of file EnergyScaleCorrection_class.cc.

References kPhi, kRho, HLTSiStripMonitoring_cff::nSigma, and convertSQLiteXML::runNumber.

Referenced by ElectronEnergyCalibratorRun2::calibrate(), and PhotonEnergyCalibratorRun2::calibrate().

                                                                                                                 {
  if (par == kRho)
    return getSmearingSigma(runNumber, isEBEle, R9Ele, etaSCEle, EtEle, nSigma, 0.);
  if (par == kPhi)
    return getSmearingSigma(runNumber, isEBEle, R9Ele, etaSCEle, EtEle, 0., nSigma);
  return getSmearingSigma(runNumber, isEBEle, R9Ele, etaSCEle, EtEle, 0., 0.);
}

getSmearingSigma() [2/2]

float EnergyScaleCorrection_class::getSmearingSigma	(	int	runNumber,
		bool	isEBEle,
		float	R9Ele,
		float	etaSCEle,
		float	EtEle,
		float	nSigma_rho,
		float	nSigma_phi
	)		const

Definition at line 345 of file EnergyScaleCorrection_class.cc.

References alpha, taus_updatedMVAIds_cff::category, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, funct::cos(), gather_cfg::cout, phi, rho, convertSQLiteXML::runNumber, funct::sin(), smearings, smearings_not_defined, and mathSSE::sqrt().

                                                                                                                     {
  correctionCategory_class category(runNumber, etaSCEle, R9Ele, EtEle);
  correction_map_t::const_iterator corr_itr = smearings.find(category);
  if (corr_itr == smearings.end()) {  // if not in the standard classes, add it in the list of not defined classes
    // the following commented part makes the method non const
    // if(smearings_not_defined.count(category) == 0) {
    //  correctionValue_class corr;
    //  smearings_not_defined[category] = corr;
    // }
    corr_itr = smearings_not_defined.find(category);
    std::cerr << "[WARNING] Smearing category not found: " << std::endl;
    std::cerr << category << std::endl;
    //     exit(1);
  }

#ifdef DEBUG
  std::cout << "[DEBUG] Checking smearing correction for category: " << category << std::endl;
  std::cout << "[DEBUG] Correction is: " << corr_itr->second << std::endl
            << "        given for category " << corr_itr->first;
#endif

  double rho = corr_itr->second.rho + corr_itr->second.rho_err * nSigma_rho;
  double phi = corr_itr->second.phi + corr_itr->second.phi_err * nSigma_phi;

  double constTerm = rho * sin(phi);
  double alpha = rho * corr_itr->second.Emean * cos(phi);

  return sqrt(constTerm * constTerm + alpha * alpha / EtEle);
}

ReadFromFile()

void EnergyScaleCorrection_class::ReadFromFile ( TString  filename )

private

category "runNumber" runMin runMax deltaP err_deltaP_per_bin err_deltaP_stat err_deltaP_syst

Input file structure: category "runNumber" runMin runMax deltaP err_deltaP(stat on single bins) err_deltaP_stat(to be used) err_deltaP_syst(to be used)

Definition at line 144 of file EnergyScaleCorrection_class.cc.

References AddScale(), taus_updatedMVAIds_cff::category, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, gather_cfg::cout, beamvalidation::exit(), corrVsCorr::filename, contentValuesFiles::fullPath, test_db_connect::runMax, and test_db_connect::runMin.

Referenced by EnergyScaleCorrection_class().

                                                               {
#ifdef PEDANTIC_OUTPUT
  std::cout << "[STATUS] Reading energy scale correction  values from file: " << filename << std::endl;
#endif
  std::cout << "[STATUS] Reading energy scale correction  values from file: " << filename << std::endl;

  //std::ifstream Ccufile(edm::FileInPath(Ccufilename).fullPath().c_str(),std::ios::in);
  std::ifstream f_in(edm::FileInPath(filename).fullPath().c_str());

  if (!f_in.good()) {
    std::cerr << "[ERROR] file " << filename << " not readable" << std::endl;
    exit(1);
    return;
  }

  int runMin, runMax;
  TString category, region2;
  double deltaP, err_deltaP, err_deltaP_stat, err_deltaP_syst;

  for (f_in >> category; f_in.good(); f_in >> category) {
    f_in >> region2 >> runMin >> runMax >> deltaP >> err_deltaP >> err_deltaP_stat >> err_deltaP_syst;

    AddScale(category, runMin, runMax, deltaP, err_deltaP_stat, err_deltaP_syst);
  }

  f_in.close();

  return;
}

ReadSmearingFromFile()

void EnergyScaleCorrection_class::ReadSmearingFromFile ( TString  filename )

private

File structure: category constTerm alpha;.

File structure: EBlowEtaBad8TeV 0 0.0 1.0 -999. 0.94 -999999 999999 6.73 0. 7.7e-3 6.32e-4 0.00 0.16 EBlowEtaGold8TeV 0 0.0 1.0 0.94 999. -999999 999999 6.60 0. 7.4e-3 6.50e-4 0.00 0.16 EBhighEtaBad8TeV 0 1.0 1.5 -999. 0.94 -999999 999999 6.73 0. 1.26e-2 1.03e-3 0.00 0.07 EBhighEtaGold8TeV 0 1.0 1.5 0.94 999. -999999 999999 6.52 0. 1.12e-2 1.32e-3 0.00 0.22

EElowEtaBad8TeV 0 1.5 2.0 -999. 0.94 -999999 999999 0. 0. 1.98e-2 3.03e-3 0. 0. EElowEtaGold8TeV 0 1.5 2.0 0.94 999. -999999 999999 0. 0. 1.63e-2 1.22e-3 0. 0. EEhighEtaBad8TeV 0 2.0 3.0 -999. 0.94 -999999 999999 0. 0. 1.92e-2 9.22e-4 0. 0. EEhighEtaGold8TeV 0 2.0 3.0 0.94 999. -999999 999999 0. 0. 1.86e-2 7.81e-4 0. 0.

Definition at line 257 of file EnergyScaleCorrection_class.cc.

References AddSmearing(), cms::cuda::assert(), taus_updatedMVAIds_cff::category, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, gather_cfg::cout, ECALELF, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, beamvalidation::exit(), corrVsCorr::filename, contentValuesFiles::fullPath, GLOBE, M_PI_2, phi, cosmicPhotonAnalyzer_cfi::r9Max, photonAnalyzer_cfi::r9Min, rho, test_db_connect::runMax, test_db_connect::runMin, smearingType_, AlCaHLTBitMon_QueryRunRegistry::string, and UNKNOWN.

Referenced by EnergyScaleCorrection_class().

                                                                       {
#ifdef PEDANTIC_OUTPUT
  std::cout << "[STATUS] Reading smearing values from file: " << filename << std::endl;
#endif
  //edm::FileInPath(Ccufilename).fullPath().c_str(),std::ios::in); .fullPath().c_str()
  std::ifstream f_in(edm::FileInPath(filename).fullPath().c_str());
  if (!f_in.good()) {
    std::cerr << "[ERROR] file " << filename << " not readable" << std::endl;
    exit(1);
  }

  int runMin = 0, runMax = 900000;
  int unused = 0;
  TString category, region2;
  //double smearing, err_smearing;
  double rho, phi, Emean, err_rho, err_phi, err_Emean;
  double etaMin, etaMax, r9Min, r9Max;
  std::string phi_string, err_phi_string;

  while (f_in.peek() != EOF && f_in.good()) {
    if (f_in.peek() == 10) {  // 10 = \n
      f_in.get();
      continue;
    }

    if (f_in.peek() == 35) {  // 35 = #
      f_in.ignore(1000, 10);  // ignore the rest of the line until \n
      continue;
    }

    if (smearingType_ == UNKNOWN) {  // trying to guess: not recommended
      std::cerr << "[ERROR] Not implemented" << std::endl;
      assert(false);

    } else if (smearingType_ == GLOBE) {
      f_in >> category >> unused >> etaMin >> etaMax >> r9Min >> r9Max >> runMin >> runMax >> Emean >> err_Emean >>
          rho >> err_rho >> phi >> err_phi;

      AddSmearing(category, runMin, runMax, rho, err_rho, phi, err_phi, Emean, err_Emean);

    } else if (smearingType_ == ECALELF) {
      f_in >> category >> Emean >> err_Emean >> rho >> err_rho >> phi_string >> err_phi_string;
#ifdef DEBUG
      std::cout
          << category
          //<< "\t" << etaMin << "\t" << etaMax << "\t" << r9Min << "\t" << r9Max << "\t" << runMin << "\t" << runMax
          << "\tEmean=" << Emean << "\t" << rho << "\t" << err_rho << "\tphi_string=" << phi_string
          << "#\terr_phi_string=" << err_phi_string << std::endl;
#endif

      if (phi_string == "M_PI_2")
        phi = M_PI_2;
      else
        phi = std::stod(phi_string);

      if (err_phi_string == "M_PI_2")
        err_phi = M_PI_2;
      else
        err_phi = std::stod(err_phi_string);

      AddSmearing(category, runMin, runMax, rho, err_rho, phi, err_phi, Emean, err_Emean);

    } else {
      f_in >> category >> rho >> phi;
      err_rho = err_phi = Emean = err_Emean = 0;
      AddSmearing(category, runMin, runMax, rho, err_rho, phi, err_phi, Emean, err_Emean);
    }
#ifdef DEBUG
    std::cout << category << "\t" << etaMin << "\t" << etaMax << "\t" << r9Min << "\t" << r9Max << "\t" << runMin
              << "\t" << runMax << "\tEmean=" << Emean << "\t" << rho << "\t" << phi << std::endl;
#endif
  }

  f_in.close();
  //  runCorrection_itr=runMin_map.begin();

  return;
}

ScaleCorrection()

float EnergyScaleCorrection_class::ScaleCorrection	(	unsigned int	runNumber,
		bool	isEBEle,
		double	R9Ele,
		double	etaSCEle,
		double	EtEle
	)		const

method to get energy scale corrections

Definition at line 44 of file EnergyScaleCorrection_class.cc.

References pfMETCorrectionType0_cfi::correction, doScale, getScaleOffset(), and convertSQLiteXML::runNumber.

Referenced by ElectronEnergyCalibratorRun2::calibrate(), and PhotonEnergyCalibratorRun2::calibrate().

                                                                                             {
  float correction = 1;
  if (doScale)
    correction *= getScaleOffset(runNumber, isEBEle, R9Ele, etaSCEle, EtEle);

  return correction;
}

ScaleCorrectionUncertainty()

float EnergyScaleCorrection_class::ScaleCorrectionUncertainty	(	unsigned int	runNumber,
		bool	isEBEle,
		double	R9Ele,
		double	etaSCEle,
		double	EtEle
	)		const

method to get scale correction uncertainties: it's stat+syst in eta x R9 categories

Definition at line 53 of file EnergyScaleCorrection_class.cc.

References getScaleStatUncertainty(), getScaleSystUncertainty(), convertSQLiteXML::runNumber, and mathSSE::sqrt().

                                                                                             {
  float statUncert = getScaleStatUncertainty(runNumber, isEBEle, R9Ele, etaSCEle, EtEle);
  float systUncert = getScaleSystUncertainty(runNumber, isEBEle, R9Ele, etaSCEle, EtEle);

  return sqrt(statUncert * statUncert + systUncert * systUncert);
}

SetSmearingType()

void EnergyScaleCorrection_class::SetSmearingType ( fileFormat_t  value )

inline

Definition at line 194 of file EnergyScaleCorrection_class.h.

References smearingType_, UNKNOWN, and relativeConstraints::value.

                                                  {
    if (value <= 1) {
      smearingType_ = value;
    } else {
      smearingType_ = UNKNOWN;
    }
  };

Member Data Documentation

doScale

bool EnergyScaleCorrection_class::doScale

Definition at line 124 of file EnergyScaleCorrection_class.h.

Referenced by ElectronEnergyCalibratorRun2::ElectronEnergyCalibratorRun2(), PhotonEnergyCalibratorRun2::PhotonEnergyCalibratorRun2(), and ScaleCorrection().

doSmearings

bool EnergyScaleCorrection_class::doSmearings

Definition at line 124 of file EnergyScaleCorrection_class.h.

Referenced by ElectronEnergyCalibratorRun2::ElectronEnergyCalibratorRun2(), and PhotonEnergyCalibratorRun2::PhotonEnergyCalibratorRun2().

scales

correction_map_t EnergyScaleCorrection_class::scales

private

Definition at line 180 of file EnergyScaleCorrection_class.h.

Referenced by AddScale(), EnergyScaleCorrection_class(), getScaleCorrection(), and getScaleOffset().

scales_not_defined

correction_map_t EnergyScaleCorrection_class::scales_not_defined

private

Definition at line 180 of file EnergyScaleCorrection_class.h.

smearings

correction_map_t EnergyScaleCorrection_class::smearings

private

Definition at line 181 of file EnergyScaleCorrection_class.h.

Referenced by AddSmearing(), EnergyScaleCorrection_class(), getSmearingRho(), and getSmearingSigma().

smearings_not_defined

correction_map_t EnergyScaleCorrection_class::smearings_not_defined

private

Definition at line 181 of file EnergyScaleCorrection_class.h.

Referenced by getSmearingRho(), and getSmearingSigma().

smearingType_

fileFormat_t EnergyScaleCorrection_class::smearingType_

private

Definition at line 178 of file EnergyScaleCorrection_class.h.

Referenced by ReadSmearingFromFile(), and SetSmearingType().