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

IsoDeposits

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

Definition at line 33 of file PFIsoDepositAlgo.h.

Particle

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

Definition at line 29 of file PFIsoDepositAlgo.h.

ParticleCollection

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

Definition at line 30 of file PFIsoDepositAlgo.h.

Constructor & Destructor Documentation

PFIsoDepositAlgo()

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

explicit

Definition at line 15 of file PFIsoDepositAlgo.cc.

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

~PFIsoDepositAlgo()

PFIsoDepositAlgo::~PFIsoDepositAlgo

(

)

Definition at line 20 of file PFIsoDepositAlgo.cc.

{}

Member Function Documentation

buildIsoDeposit()

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

private

build the IsoDeposit for "particle"

Definition at line 43 of file PFIsoDepositAlgo.cc.

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

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

Referenced by produce().

produce()

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 22 of file PFIsoDepositAlgo.cc.

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

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

sameParticle()

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

private

checks if the 2 particles are in fact the same

Definition at line 81 of file PFIsoDepositAlgo.cc.

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

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

Referenced by buildIsoDeposit().

Member Data Documentation

isoDeposits_

IsoDeposits pf2pat::PFIsoDepositAlgo::isoDeposits_

private

IsoDeposits computed in the produce function.

Definition at line 53 of file PFIsoDepositAlgo.h.

Referenced by produce().

verbose_

bool pf2pat::PFIsoDepositAlgo::verbose_

private

Definition at line 55 of file PFIsoDepositAlgo.h.

Referenced by buildIsoDeposit(), and produce().