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.

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().

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

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.

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

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

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.

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_

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().