EffectiveAreas Class Reference

#include <EffectiveAreas.h>

Public Member Functions
void	checkConsistency () const
	EffectiveAreas (const std::string &filename)
const float	getEffectiveArea (float eta) const
void	printEffectiveAreas () const

Private Attributes
std::vector< float >	absEtaMax_
std::vector< float >	absEtaMin_
std::vector< float >	effectiveAreaValues_
const std::string	filename_

Detailed Description

Definition at line 7 of file EffectiveAreas.h.

Constructor & Destructor Documentation

EffectiveAreas()

EffectiveAreas::EffectiveAreas

(

const std::string &

filename

)

Definition at line 9 of file EffectiveAreas.cc.

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

                                                        : filename_(filename) {
  // 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;
  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);
  }

  // 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 67 of file EffectiveAreas.cc.

References absEtaMax_, absEtaMin_, effectiveAreaValues_, Exception, filename_, L1TowerCalibrationProducer_cfi::iEta, L1TMuonDQMOffline_cfi::nEtaBins, and parallelization::uint.

Referenced by EffectiveAreas().

                                            {
  // There should be at least one eta range with one constant
  if (effectiveAreaValues_.empty())
    throw cms::Exception("EffectiveAreas config failure")
        << "found no effective area constans 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;
  }
}

getEffectiveArea()

const float EffectiveAreas::getEffectiveArea

(

float

eta

)

const

Definition at line 44 of file EffectiveAreas.cc.

References funct::abs(), absEtaMax_, absEtaMin_, effectiveAreaValues_, PVValHelper::eta, L1TowerCalibrationProducer_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;
}

printEffectiveAreas()

void EffectiveAreas::printEffectiveAreas

(

)

const

Definition at line 57 of file EffectiveAreas.cc.

References absEtaMax_, absEtaMin_, effectiveAreaValues_, filename_, L1TowerCalibrationProducer_cfi::iEta, L1TMuonDQMOffline_cfi::nEtaBins, and parallelization::uint.

                                               {
  printf("EffectiveAreas: source file %s\n", filename_.c_str());
  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]);
  }
}

Member Data Documentation

absEtaMax_

std::vector<float> EffectiveAreas::absEtaMax_

private

Definition at line 23 of file EffectiveAreas.h.

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

absEtaMin_

std::vector<float> EffectiveAreas::absEtaMin_

private

Definition at line 22 of file EffectiveAreas.h.

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

effectiveAreaValues_

std::vector<float> EffectiveAreas::effectiveAreaValues_

private

Definition at line 24 of file EffectiveAreas.h.

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

filename_

const std::string EffectiveAreas::filename_

private

Definition at line 21 of file EffectiveAreas.h.

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