HighETPhotonsFilterAlgo Class Reference

#include <HighETPhotonsFilterAlgo.h>

Public Member Functions
bool	filter (const edm::Event &iEvent)
	HighETPhotonsFilterAlgo (const edm::ParameterSet &)
	~HighETPhotonsFilterAlgo ()

Private Attributes
float	FILTER_ETA_MAX_
float	isoConeSize_
float	nonPhotETMax_
float	seedETThreshold_
float	sumETThreshold_

Detailed Description

HighETPhotonsFilterAlgo a gen-level filter that selects events that will reconstruct an isolated photon only tested with high ET thresholds (aiming for 100 GeV photons)

Author

J Lamb, UCSB

Definition at line 28 of file HighETPhotonsFilterAlgo.h.

Constructor & Destructor Documentation

HighETPhotonsFilterAlgo::HighETPhotonsFilterAlgo

(

const edm::ParameterSet &

iConfig

)

Definition at line 8 of file HighETPhotonsFilterAlgo.cc.

References edm::ParameterSet::getParameter().

                                                                               { 

  //set constants
  FILTER_ETA_MAX_=2.5;
  //parameters
  sumETThreshold_=(float)iConfig.getParameter<double>("sumETThreshold");
  seedETThreshold_=(float)iConfig.getParameter<double>("seedETThreshold");
  nonPhotETMax_=(float)iConfig.getParameter<double>("nonPhotETMax");
  isoConeSize_=(float)iConfig.getParameter<double>("isoConeSize");


}

HighETPhotonsFilterAlgo::~HighETPhotonsFilterAlgo

(

)

Definition at line 21 of file HighETPhotonsFilterAlgo.cc.

                                                  {
}

Member Function Documentation

bool HighETPhotonsFilterAlgo::filter

(

const edm::Event &

iEvent

)

Definition at line 25 of file HighETPhotonsFilterAlgo.cc.

References deltaR(), reco::LeafCandidate::et(), edm::Event::getByLabel(), reco::LeafCandidate::pdgId(), and reco::LeafCandidate::status().

                                                            {

  Handle<reco::GenParticleCollection> genParsHandle;
  iEvent.getByLabel("genParticles",genParsHandle);
  vector<reco::GenParticle> genPars=*genParsHandle;
  
  for (uint32_t ig=0;ig<genPars.size();ig++) {
    reco::GenParticle gp=genPars.at(ig);
    //look for seeds
    if (fabs(gp.pdgId())!=22) continue;
    if (gp.status()!=1) continue;
    if (gp.et()<seedETThreshold_) continue;
    //found a seed
    float photetsum=gp.et();
    float nonphotetsum=0;
    for (uint32_t jg=0;jg<genPars.size();jg++) {
      if (ig==jg) continue;
      reco::GenParticle opar=genPars.at(jg);
      if (opar.status()!=1) continue;
      float dr=deltaR(gp,opar);
      if (dr>=isoConeSize_) continue;
      if (fabs(opar.pdgId())==22) {
    photetsum+=opar.et();
      } else {
    nonphotetsum+=opar.et();
      }
    }
    //check that photon energy exceeds threshold, and isolation
    //energy is below threshold
    if (photetsum>sumETThreshold_ && nonphotetsum<nonPhotETMax_) return true;
  }
  return false;

}

Member Data Documentation

float HighETPhotonsFilterAlgo::FILTER_ETA_MAX_

private

Definition at line 42 of file HighETPhotonsFilterAlgo.h.

float HighETPhotonsFilterAlgo::isoConeSize_

private

Definition at line 47 of file HighETPhotonsFilterAlgo.h.

float HighETPhotonsFilterAlgo::nonPhotETMax_

private

Definition at line 46 of file HighETPhotonsFilterAlgo.h.

float HighETPhotonsFilterAlgo::seedETThreshold_

private

Definition at line 45 of file HighETPhotonsFilterAlgo.h.

float HighETPhotonsFilterAlgo::sumETThreshold_

private

Definition at line 44 of file HighETPhotonsFilterAlgo.h.