#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.

                                                         : 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();
 }

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.

Member Function Documentation

◆ checkConsistency()

void EffectiveAreas::checkConsistency ( void ) const

Definition at line 67 of file EffectiveAreas.cc.

                                             {
   // 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;
   }
 }

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

Referenced by EffectiveAreas().

◆ getEffectiveArea()

const float EffectiveAreas::getEffectiveArea ( float eta ) const

Definition at line 44 of file EffectiveAreas.cc.

                                                             {
   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;
 }

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().

◆ printEffectiveAreas()

void EffectiveAreas::printEffectiveAreas ( ) const

Definition at line 57 of file EffectiveAreas.cc.

                                                {
   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]);
   }
 }