EffectiveAreas Class Reference

#include <EffectiveAreas.h>

Public Member Functions
void	checkConsistency () const
	EffectiveAreas (const std::string &filename, const bool quadraticEAflag=false)
const float	getEffectiveArea (float eta) const
const float	getLinearEA (float eta) const
const float	getQuadraticEA (float eta) const
void	printEffectiveAreas () const

Private Attributes
std::vector< float >	absEtaMax_
std::vector< float >	absEtaMin_
std::vector< float >	effectiveAreaValues_
const std::string	filename_
std::vector< float >	linearEffectiveAreaValues_
const bool	quadraticEAflag_
std::vector< float >	quadraticEffectiveAreaValues_

Detailed Description

Definition at line 7 of file EffectiveAreas.h.

Constructor & Destructor Documentation

EffectiveAreas()

EffectiveAreas::EffectiveAreas	(	const std::string &	filename,
		const bool	quadraticEAflag = false
	)

Definition at line 9 of file EffectiveAreas.cc.

References absEtaMax_, absEtaMin_, checkConsistency(), effectiveAreaValues_, ALCARECOTkAlBeamHalo_cff::etaMax, ALCARECOTkAlBeamHalo_cff::etaMin, Exception, filename_, makeListRunsInFiles::inputFile, mps_splice::line, linearEffectiveAreaValues_, quadraticEAflag_, quadraticEffectiveAreaValues_, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, and trackerHitRTTI::undef.

    : filename_(filename), quadraticEAflag_(quadraticEAflag) {
  // Open the file with the effective area constants
  std::ifstream inputFile;
  inputFile.open(filename_.c_str());
  if (!inputFile.is_open())
    throw cms::Exception("EffectiveAreas config failure") << "failed to open the file " << filename_ << std::endl;

  // Read file line by line
  std::string line;
  const float undef = -999;

  if (!quadraticEAflag_) {
    while (getline(inputFile, line)) {
      if (line[0] == '#')
        continue;  // skip the comments lines
      float etaMin = undef, etaMax = undef, effArea = undef;
      std::stringstream ss(line);
      ss >> etaMin >> etaMax >> effArea;
      // In case if the format is messed up, there are letters
      // instead of numbers, or not exactly three numbers in the line,
      // it is likely that one or more of these vars never changed
      // the original "undef" value:
      if (etaMin == undef || etaMax == undef || effArea == undef)
        throw cms::Exception("EffectiveAreas config failure")
            << "wrong file format, file name " << filename_ << std::endl;

      absEtaMin_.push_back(etaMin);
      absEtaMax_.push_back(etaMax);
      effectiveAreaValues_.push_back(effArea);
    }

  }

  else {
    while (getline(inputFile, line)) {
      if (line[0] == '#')
        continue;  // skip the comments lines

      float etaMin = undef, etaMax = undef;
      float linEffArea = undef, quadEffArea = undef;

      std::stringstream ss(line);

      ss >> etaMin >> etaMax >> linEffArea >> quadEffArea;

      // In case if the format is messed up, there are letters
      // instead of numbers, or not exactly three numbers in the line,
      // it is likely that one or more of these vars never changed
      // the original "undef" value:
      if (etaMin == undef || etaMax == undef || linEffArea == undef || quadEffArea == undef)
        throw cms::Exception("EffectiveAreas config failure")
            << "wrong file format, file name " << filename_ << std::endl;

      absEtaMin_.push_back(etaMin);
      absEtaMax_.push_back(etaMax);
      linearEffectiveAreaValues_.push_back(linEffArea);
      quadraticEffectiveAreaValues_.push_back(quadEffArea);
    }
  }

  // Extra consistency checks are in the function below.
  // If any of them fail, an exception is thrown.
  checkConsistency();
}

Member Function Documentation

checkConsistency()

void EffectiveAreas::checkConsistency

(

void

)

const

Definition at line 145 of file EffectiveAreas.cc.

References absEtaMax_, absEtaMin_, effectiveAreaValues_, Exception, filename_, l1tTowerCalibrationProducer_cfi::iEta, linearEffectiveAreaValues_, L1TMuonDQMOffline_cfi::nEtaBins, quadraticEAflag_, quadraticEffectiveAreaValues_, and parallelization::uint.

Referenced by EffectiveAreas().

                                            {
  // There should be at least one eta range with one constant

  if (!quadraticEAflag_) {
    if (effectiveAreaValues_.empty())
      throw cms::Exception("EffectiveAreas config failure")
          << "found no effective area constants in the file " << filename_ << std::endl;

    uint nEtaBins = absEtaMin_.size();

    for (uint iEta = 0; iEta < nEtaBins; iEta++) {
      // The low limit should be lower than the upper limit
      if (!(absEtaMin_[iEta] < absEtaMax_[iEta]))
        throw cms::Exception("EffectiveAreas config failure")
            << "eta ranges improperly defined (min>max) in the file" << filename_ << std::endl;

      // The low limit of the next range should be (near) equal to the
      // upper limit of the previous range
      if (iEta != nEtaBins - 1)  // don't do the check for the last bin
        if (!(absEtaMin_[iEta + 1] - absEtaMax_[iEta] < 0.0001))
          throw cms::Exception("EffectiveAreas config failure")
              << "eta ranges improperly defined (disjointed) in the file " << filename_ << std::endl;

      // The effective area should be non-negative number,
      // and should be less than the whole calorimeter area
      // eta range -2.5 to 2.5, phi 0 to 2pi => Amax = 5*2*pi ~= 31.4
      if (!(effectiveAreaValues_[iEta] >= 0 && effectiveAreaValues_[iEta] < 31.4))
        throw cms::Exception("EffectiveAreas config failure")
            << "effective area values are too large or negative in the file" << filename_ << std::endl;
    }
  }

  else {
    if (linearEffectiveAreaValues_.empty() || quadraticEffectiveAreaValues_.empty())
      throw cms::Exception("EffectiveAreas config failure")
          << "found no effective area constants (linear and/or quadratic) in the file " << filename_ << std::endl;

    uint nEtaBins = absEtaMin_.size();

    for (uint iEta = 0; iEta < nEtaBins; iEta++) {
      // The low limit should be lower than the upper limit
      if (!(absEtaMin_[iEta] < absEtaMax_[iEta]))
        throw cms::Exception("EffectiveAreas config failure")
            << "eta ranges improperly defined (min>max) in the file" << filename_ << std::endl;

      // The low limit of the next range should be (near) equal to the
      // upper limit of the previous range
      if (iEta != nEtaBins - 1)  // don't do the check for the last bin
        if (!(absEtaMin_[iEta + 1] - absEtaMax_[iEta] < 0.0001))
          throw cms::Exception("EffectiveAreas config failure")
              << "eta ranges improperly defined (disjointed) in the file " << filename_ << std::endl;

      // The linear effective area should be non-negative number,
      // and should be less than the whole calorimeter area
      // eta range -2.5 to 2.5, phi 0 to 2pi => Amax = 5*2*pi ~= 31.4
      if (!(linearEffectiveAreaValues_[iEta] >= 0 && linearEffectiveAreaValues_[iEta] < 31.4))
        throw cms::Exception("EffectiveAreas config failure")
            << "effective area values are too large or negative in the file" << filename_ << std::endl;
    }
  }
}

getEffectiveArea()

const float EffectiveAreas::getEffectiveArea

(

float

eta

)

const

Definition at line 76 of file EffectiveAreas.cc.

References funct::abs(), absEtaMax_, absEtaMin_, effectiveAreaValues_, PVValHelper::eta, l1tTowerCalibrationProducer_cfi::iEta, L1TMuonDQMOffline_cfi::nEtaBins, and parallelization::uint.

Referenced by ElectronIdentifier::isolation(), GsfEleEffAreaPFIsoCut::operator()(), PhoAnyPFIsoWithEAAndExpoScalingEBCut::operator()(), PhoAnyPFIsoWithEAAndQuadScalingCut::operator()(), PhoAnyPFIsoWithEAAndExpoScalingCut::operator()(), PhoAnyPFIsoWithEACut::operator()(), GsfEleCalPFClusterIsoCut::operator()(), PhoGenericRhoPtScaledCut::operator()(), PhoAnyPFIsoWithEAAndExpoScalingCut::value(), GsfEleEffAreaPFIsoCut::value(), PhoAnyPFIsoWithEAAndExpoScalingEBCut::value(), PhoAnyPFIsoWithEACut::value(), PhoAnyPFIsoWithEAAndQuadScalingCut::value(), GsfEleRelPFIsoScaledCut::value(), and GsfEleCalPFClusterIsoCut::value().

                                                            {
  float effArea = 0;
  uint nEtaBins = absEtaMin_.size();
  for (uint iEta = 0; iEta < nEtaBins; iEta++) {
    if (std::abs(eta) >= absEtaMin_[iEta] && std::abs(eta) < absEtaMax_[iEta]) {
      effArea = effectiveAreaValues_[iEta];
      break;
    }
  }

  return effArea;
}

getLinearEA()

const float EffectiveAreas::getLinearEA

(

float

eta

)

const

Definition at line 90 of file EffectiveAreas.cc.

References funct::abs(), absEtaMax_, absEtaMin_, PVValHelper::eta, l1tTowerCalibrationProducer_cfi::iEta, linearEffectiveAreaValues_, L1TMuonDQMOffline_cfi::nEtaBins, and parallelization::uint.

Referenced by PhoGenericQuadraticRhoPtScaledCut::operator()().

                                                       {
  float linEffArea = 0;
  uint nEtaBins = absEtaMin_.size();

  for (uint iEta = 0; iEta < nEtaBins; iEta++) {
    if (std::abs(eta) >= absEtaMin_[iEta] && std::abs(eta) < absEtaMax_[iEta]) {
      linEffArea = linearEffectiveAreaValues_[iEta];
      break;
    }
  }

  return linEffArea;
}

getQuadraticEA()

const float EffectiveAreas::getQuadraticEA

(

float

eta

)

const

Definition at line 105 of file EffectiveAreas.cc.

References funct::abs(), absEtaMax_, absEtaMin_, PVValHelper::eta, l1tTowerCalibrationProducer_cfi::iEta, L1TMuonDQMOffline_cfi::nEtaBins, quadraticEffectiveAreaValues_, and parallelization::uint.

Referenced by PhoGenericQuadraticRhoPtScaledCut::operator()().

                                                          {
  float quadEffArea = 0;
  uint nEtaBins = absEtaMin_.size();

  for (uint iEta = 0; iEta < nEtaBins; iEta++) {
    if (std::abs(eta) >= absEtaMin_[iEta] && std::abs(eta) < absEtaMax_[iEta]) {
      quadEffArea = quadraticEffectiveAreaValues_[iEta];
      break;
    }
  }

  return quadEffArea;
}

printEffectiveAreas()

void EffectiveAreas::printEffectiveAreas

(

)

const

Definition at line 119 of file EffectiveAreas.cc.

References absEtaMax_, absEtaMin_, effectiveAreaValues_, filename_, l1tTowerCalibrationProducer_cfi::iEta, linearEffectiveAreaValues_, L1TMuonDQMOffline_cfi::nEtaBins, quadraticEAflag_, quadraticEffectiveAreaValues_, and parallelization::uint.

                                               {
  printf("EffectiveAreas: source file %s\n", filename_.c_str());

  if (!quadraticEAflag_) {
    printf("  eta_min   eta_max    effective area\n");
    uint nEtaBins = absEtaMin_.size();

    for (uint iEta = 0; iEta < nEtaBins; iEta++) {
      printf("  %8.4f    %8.4f   %8.5f\n", absEtaMin_[iEta], absEtaMax_[iEta], effectiveAreaValues_[iEta]);
    }
  }

  else {
    printf("  eta_min   eta_max    EA linear term    EA quadratic term\n");
    uint nEtaBins = absEtaMin_.size();
    for (uint iEta = 0; iEta < nEtaBins; iEta++) {
      printf("  %8.4f    %8.4f   %8.5f    %8.5f\n",
             absEtaMin_[iEta],
             absEtaMax_[iEta],
             linearEffectiveAreaValues_[iEta],
             quadraticEffectiveAreaValues_[iEta]);
    }
  }
}

Member Data Documentation

absEtaMax_

std::vector<float> EffectiveAreas::absEtaMax_

private

Definition at line 25 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), getEffectiveArea(), getLinearEA(), getQuadraticEA(), and printEffectiveAreas().

absEtaMin_

std::vector<float> EffectiveAreas::absEtaMin_

private

Definition at line 24 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), getEffectiveArea(), getLinearEA(), getQuadraticEA(), and printEffectiveAreas().

effectiveAreaValues_

std::vector<float> EffectiveAreas::effectiveAreaValues_

private

Definition at line 26 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), getEffectiveArea(), and printEffectiveAreas().

filename_

const std::string EffectiveAreas::filename_

private

Definition at line 23 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), and printEffectiveAreas().

linearEffectiveAreaValues_

std::vector<float> EffectiveAreas::linearEffectiveAreaValues_

private

Definition at line 30 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), getLinearEA(), and printEffectiveAreas().

quadraticEAflag_

const bool EffectiveAreas::quadraticEAflag_

private

Definition at line 29 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), and printEffectiveAreas().

quadraticEffectiveAreaValues_

std::vector<float> EffectiveAreas::quadraticEffectiveAreaValues_

private

Definition at line 31 of file EffectiveAreas.h.

Referenced by checkConsistency(), EffectiveAreas(), getQuadraticEA(), and printEffectiveAreas().