PfBlockBasedIsolation Class Reference

#include <PfBlockBasedIsolation.h>

Public Member Functions
std::vector< reco::PFCandidateRef >	calculate (math::XYZTLorentzVectorD p4, const reco::PFCandidateRef pfEGCand, const edm::Handle< reco::PFCandidateCollection > pfCandidateHandle)
	PfBlockBasedIsolation ()
void	setup (const edm::ParameterSet &conf)
	~PfBlockBasedIsolation ()

Private Attributes
double	coneSize_

Detailed Description

Definition at line 34 of file PfBlockBasedIsolation.h.

Constructor & Destructor Documentation

PfBlockBasedIsolation::PfBlockBasedIsolation

(

)

Definition at line 22 of file PfBlockBasedIsolation.cc.

                                             {
  // Default Constructor.
}

PfBlockBasedIsolation::~PfBlockBasedIsolation

(

)

Definition at line 30 of file PfBlockBasedIsolation.cc.

{

}

Member Function Documentation

std::vector< reco::PFCandidateRef > PfBlockBasedIsolation::calculate	(	math::XYZTLorentzVectorD	p4,
		const reco::PFCandidateRef	pfEGCand,
		const edm::Handle< reco::PFCandidateCollection >	pfCandidateHandle
	)

Definition at line 43 of file PfBlockBasedIsolation.cc.

References coneSize_, deltaR(), and PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR.

Referenced by ParticleBasedIsoProducer::produce().

                                                                                                                                                                                          {
  
  std::vector<reco::PFCandidateRef> myVec;
  
  math::XYZVector candidateMomentum(p4.px(),p4.py(),p4.pz());
  math::XYZVector candidateDirection=candidateMomentum.Unit();

  const reco::PFCandidate::ElementsInBlocks& theElementsInpfEGcand = (*pfEGCand).elementsInBlocks();
  reco::PFCandidate::ElementsInBlocks::const_iterator ieg = theElementsInpfEGcand.begin();
  const reco::PFBlockRef egblock = ieg->first;


  unsigned nObj = pfCandidateHandle->size();
  for(unsigned int lCand=0; lCand < nObj; lCand++) {

    reco::PFCandidateRef pfCandRef(reco::PFCandidateRef(pfCandidateHandle,lCand));

    float dR = 0.0;
    if( coneSize_ < 10.0 ) {
      dR = deltaR(candidateDirection.Eta(), candidateDirection.Phi(),  pfCandRef->eta(),   pfCandRef->phi());         
      if ( dR> coneSize_ ) continue;
    }

    const reco::PFCandidate::ElementsInBlocks& theElementsInPFcand = pfCandRef->elementsInBlocks();

    bool elementFound=false;
    for (reco::PFCandidate::ElementsInBlocks::const_iterator ipf = theElementsInPFcand.begin(); ipf<theElementsInPFcand.end(); ++ipf) {
 
     if ( ipf->first == egblock && !elementFound ) {

      for (ieg = theElementsInpfEGcand.begin(); ieg<theElementsInpfEGcand.end(); ++ieg) {
        if ( ipf->second == ieg->second && !elementFound  ) {
          elementFound=true;
          myVec.push_back(pfCandRef);    
           
        }
      }
    
    
    
      }
    }

    

  }
  


  return myVec;



 }

void PfBlockBasedIsolation::setup

(

const edm::ParameterSet &

conf

)

Definition at line 36 of file PfBlockBasedIsolation.cc.

References coneSize_, and edm::ParameterSet::getParameter().

Referenced by ParticleBasedIsoProducer::beginRun().

                                                                 {

  coneSize_             = conf.getParameter<double>("coneSize");  

}

Member Data Documentation

double PfBlockBasedIsolation::coneSize_

private

Definition at line 59 of file PfBlockBasedIsolation.h.

Referenced by calculate(), and setup().