IsoPhotonEBSelector Class Reference

Public Member Functions
const float	getEffectiveArea (float eta) const
	IsoPhotonEBSelector (const edm::ParameterSet &, edm::ConsumesCollector &iC)
void	newEvent (const edm::Event &, const edm::EventSetup &)
bool	operator() (const reco::Photon &) const
void	printEffectiveAreas () const

Public Attributes
const std::vector< double >	absEtaMax_
const std::vector< double >	absEtaMin_
const std::vector< double >	effectiveAreaValues_
const vector< double >	hOverECut_
const edm::ParameterSet	phIDWP_
const vector< double >	relCombIso_
edm::Handle< double >	rhoHandle_
const vector< double >	sigmaIEtaIEtaCut_
const edm::EDGetTokenT< double >	theRhoToken_

Detailed Description

Definition at line 42 of file IsoPhotonEBSelector.cc.

Constructor & Destructor Documentation

IsoPhotonEBSelector()

IsoPhotonEBSelector::IsoPhotonEBSelector	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &	iC
	)

Definition at line 84 of file IsoPhotonEBSelector.cc.

    : theRhoToken_(iC.consumes<double>(cfg.getParameter<edm::InputTag>("rho"))),
      absEtaMin_(cfg.getParameter<std::vector<double> >("absEtaMin")),
      absEtaMax_(cfg.getParameter<std::vector<double> >("absEtaMax")),
      effectiveAreaValues_(cfg.getParameter<std::vector<double> >("effectiveAreaValues")),
      phIDWP_(cfg.getParameter<edm::ParameterSet>("phID")),
      sigmaIEtaIEtaCut_(phIDWP_.getParameter<std::vector<double> >("full5x5_sigmaIEtaIEtaCut")),
      hOverECut_(phIDWP_.getParameter<std::vector<double> >("hOverECut")),
      relCombIso_(phIDWP_.getParameter<std::vector<double> >("relCombIsolationWithEACut")) {
  //printEffectiveAreas();
}

Member Function Documentation

getEffectiveArea()

const float IsoPhotonEBSelector::getEffectiveArea

(

float

eta

)

const

Definition at line 71 of file IsoPhotonEBSelector.cc.

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

Referenced by operator()().

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

newEvent()

void IsoPhotonEBSelector::newEvent	(	const edm::Event &	ev,
		const edm::EventSetup &
	)

Definition at line 96 of file IsoPhotonEBSelector.cc.

References makeMEIFBenchmarkPlots::ev, rhoHandle_, and theRhoToken_.

                                                                           {
  ev.getByToken(theRhoToken_, rhoHandle_);
}

operator()()

bool IsoPhotonEBSelector::operator()

(

const reco::Photon &

ph

)

const

Definition at line 100 of file IsoPhotonEBSelector.cc.

References dqmMemoryStats::float, getEffectiveArea(), hOverECut_, edm::HandleBase::isValid(), SiStripPI::max, position, edm::Handle< T >::product(), relCombIso_, rho, rhoHandle_, and sigmaIEtaIEtaCut_.

                                                               {
  float pt_e = ph.pt();
  unsigned int ind = 0;
  float abseta = fabs((ph.superCluster().get())->position().eta());

  if (ph.isEB()) {
    if (abseta > 1.479)
      return false;  // check if it is really needed
  }
  if (ph.isEE()) {
    ind = 1;
    if (abseta <= 1.479 || abseta >= 2.5)
      return false;  // check if it is really needed
  }

  if (ph.full5x5_sigmaIetaIeta() > sigmaIEtaIEtaCut_[ind])
    return false;
  if (ph.hadronicOverEm() > hOverECut_[ind])
    return false;
  const float eA = getEffectiveArea(abseta);
  const float rho = rhoHandle_.isValid() ? (float)(*rhoHandle_.product()) : 0;
  if ((ph.getPflowIsolationVariables().chargedHadronIso +
       std::max(float(0.0),
                ph.getPflowIsolationVariables().neutralHadronIso + ph.getPflowIsolationVariables().photonIso -
                    eA * rho)) > relCombIso_[ind] * pt_e)
    return false;

  return true;
}

printEffectiveAreas()

void IsoPhotonEBSelector::printEffectiveAreas

(

)

const

Definition at line 64 of file IsoPhotonEBSelector.cc.

References l1tTowerCalibrationProducer_cfi::iEta, L1TMuonDQMOffline_cfi::nEtaBins, and parallelization::uint.

                                                    {
  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_

const std::vector<double> IsoPhotonEBSelector::absEtaMax_

Definition at line 54 of file IsoPhotonEBSelector.cc.

absEtaMin_

const std::vector<double> IsoPhotonEBSelector::absEtaMin_

Definition at line 53 of file IsoPhotonEBSelector.cc.

effectiveAreaValues_

const std::vector<double> IsoPhotonEBSelector::effectiveAreaValues_

Definition at line 55 of file IsoPhotonEBSelector.cc.

hOverECut_

const vector<double> IsoPhotonEBSelector::hOverECut_

Definition at line 60 of file IsoPhotonEBSelector.cc.

Referenced by operator()().

phIDWP_

const edm::ParameterSet IsoPhotonEBSelector::phIDWP_

Definition at line 57 of file IsoPhotonEBSelector.cc.

relCombIso_

const vector<double> IsoPhotonEBSelector::relCombIso_

Definition at line 61 of file IsoPhotonEBSelector.cc.

Referenced by operator()().

rhoHandle_

edm::Handle<double> IsoPhotonEBSelector::rhoHandle_

Definition at line 51 of file IsoPhotonEBSelector.cc.

Referenced by newEvent(), and operator()().

sigmaIEtaIEtaCut_

const vector<double> IsoPhotonEBSelector::sigmaIEtaIEtaCut_

Definition at line 59 of file IsoPhotonEBSelector.cc.

Referenced by operator()().

theRhoToken_

const edm::EDGetTokenT<double> IsoPhotonEBSelector::theRhoToken_

Definition at line 50 of file IsoPhotonEBSelector.cc.

Referenced by newEvent().