pf2pat::PFIsoDepositAlgo Class Reference

#include <PFIsoDepositAlgo.h>

Public Types
typedef std::vector< reco::IsoDeposit >	IsoDeposits
typedef reco::PFCandidate	Particle
typedef std::vector< Particle >	ParticleCollection

Public Member Functions
	PFIsoDepositAlgo (const edm::ParameterSet &)
const IsoDeposits &	produce (const ParticleCollection &toBeIsolated, const ParticleCollection &forIsolation)
	~PFIsoDepositAlgo ()

Private Member Functions
reco::IsoDeposit	buildIsoDeposit (const Particle &particle, const ParticleCollection &forIsolation) const
	build the IsoDeposit for "particle" More...
bool	sameParticle (const Particle &particle1, const Particle &particle2) const
	checks if the 2 particles are in fact the same More...

Private Attributes
IsoDeposits	isoDeposits_
	IsoDeposits computed in the produce function. More...
bool	verbose_

Detailed Description

Definition at line 26 of file PFIsoDepositAlgo.h.

Member Typedef Documentation

typedef std::vector< reco::IsoDeposit > pf2pat::PFIsoDepositAlgo::IsoDeposits

Definition at line 34 of file PFIsoDepositAlgo.h.

typedef reco::PFCandidate pf2pat::PFIsoDepositAlgo::Particle

Definition at line 30 of file PFIsoDepositAlgo.h.

typedef std::vector< Particle > pf2pat::PFIsoDepositAlgo::ParticleCollection

Definition at line 31 of file PFIsoDepositAlgo.h.

Constructor & Destructor Documentation

PFIsoDepositAlgo::PFIsoDepositAlgo ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 16 of file PFIsoDepositAlgo.cc.

                                                                : 
  verbose_ ( iConfig.getUntrackedParameter<bool>("verbose",false) )

{}

PFIsoDepositAlgo::~PFIsoDepositAlgo

(

)

Definition at line 23 of file PFIsoDepositAlgo.cc.

{ }

Member Function Documentation

IsoDeposit PFIsoDepositAlgo::buildIsoDeposit ( const Particle &  particle, const ParticleCollection &  forIsolation  ) const

private

build the IsoDeposit for "particle"

Definition at line 52 of file PFIsoDepositAlgo.cc.

References gather_cfg::cout, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::LeafCandidate::eta(), mps_fire::i, reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), EnergyCorrector::pt, sameParticle(), and verbose_.

Referenced by produce().

                                                                         {
  

  reco::isodeposit::Direction pfDir(particle.eta(), 
                    particle.phi());
//   reco::IsoDeposit::Veto veto;
//   veto.vetoDir = pfDir;
//   veto.dR = 0.05; 

  IsoDeposit isoDep( pfDir );
  
  for( unsigned i=0; i<forIsolation.size(); i++ ) {
    
    const reco::PFCandidate& pfc = forIsolation[i];

    // need to remove "particle"!

    if( sameParticle( particle, pfc ) ) continue; 


    XYZTLorentzVector pvi(pfc.p4());
    reco::isodeposit::Direction dirPfc(pfc.eta(), pfc.phi());
    double dR = pfDir.deltaR(dirPfc);

    //COLIN make a parameter
    double maxDeltaRForIsoDep_ = 1;
    if(dR > maxDeltaRForIsoDep_) {
      //      if( verbose_ ) cout<<"OUT OF CONE"<<endl;
      continue;
    }
    //    else if(verbose_) cout<<endl;

    if(verbose_ ) 
      cout<<"\t"<<pfc<<endl;

    double pt = pvi.Pt();
    isoDep.addDeposit(dirPfc, pt); 
  }

  return isoDep;
}

const PFIsoDepositAlgo::IsoDeposits & PFIsoDepositAlgo::produce	(	const ParticleCollection &	toBeIsolated,
		const ParticleCollection &	forIsolation
	)

all the filtering is done before could gain in performance by having a single loop, and by producing all isodeposits at the same time? however, would not gain in ease of maintenance, and in flexibility.

Definition at line 27 of file PFIsoDepositAlgo.cc.

References buildIsoDeposit(), gather_cfg::cout, mps_fire::i, isoDeposits_, and verbose_.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                       {
  

  isoDeposits_.clear();
  isoDeposits_.reserve( toBeIsolated.size() );

  for( unsigned i=0; i<toBeIsolated.size(); i++ ) {
    const reco::PFCandidate& toBeIso = toBeIsolated[i];

    if(verbose_ ) 
      cout<<"to be isolated: "<<toBeIso<<endl;

    isoDeposits_.push_back( buildIsoDeposit( toBeIso, forIsolation ) ); 
  }
  

  if(verbose_) {
    cout<<"PFIsoDepositAlgo "<<endl;
  }

  return isoDeposits_;
}

bool PFIsoDepositAlgo::sameParticle ( const Particle &  particle1, const Particle &  particle2  ) const

private

checks if the 2 particles are in fact the same

Definition at line 96 of file PFIsoDepositAlgo.cc.

References MillePedeFileConverter_cfg::e, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), and reco::PFCandidate::particleId().

Referenced by buildIsoDeposit().

                                                       {

  double smallNumber = 1e-15;
  
  if( particle1.particleId() != particle2.particleId() ) return false;
  else if( fabs( particle1.energy() - particle2.energy() ) > smallNumber ) return false;
  else if( fabs( particle1.eta() - particle2.eta() ) > smallNumber ) return false;
  else if( fabs( particle1.eta() - particle2.eta() ) > smallNumber ) return false;
  else return true; 
  
}

Member Data Documentation

IsoDeposits pf2pat::PFIsoDepositAlgo::isoDeposits_

private

IsoDeposits computed in the produce function.

Definition at line 59 of file PFIsoDepositAlgo.h.

Referenced by produce().

bool pf2pat::PFIsoDepositAlgo::verbose_

private

Definition at line 61 of file PFIsoDepositAlgo.h.

Referenced by buildIsoDeposit(), and produce().