PFIsoDepositAlgo.cc

Go to the documentation of this file.

#include "CommonTools/ParticleFlow/interface/PFIsoDepositAlgo.h"

#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"

#include "DataFormats/Math/interface/LorentzVector.h"

#include "FWCore/Framework/interface/EventSetup.h"

using namespace std;
using namespace edm;
using namespace reco;
using namespace math;
using namespace pf2pat;

PFIsoDepositAlgo::PFIsoDepositAlgo(const edm::ParameterSet& iConfig)
    : verbose_(iConfig.getUntrackedParameter<bool>("verbose", false))

{}

PFIsoDepositAlgo::~PFIsoDepositAlgo() {}

const PFIsoDepositAlgo::IsoDeposits& PFIsoDepositAlgo::produce(const ParticleCollection& toBeIsolated,
                                                               const ParticleCollection& forIsolation) {
  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_;
}

IsoDeposit PFIsoDepositAlgo::buildIsoDeposit(const Particle& particle, const ParticleCollection& forIsolation) const {
  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;
}

bool PFIsoDepositAlgo::sameParticle(const Particle& particle1, const Particle& particle2) const {
  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;
}