JetTracksAssociationDR.cc

Go to the documentation of this file.

// Associate jets with tracks by simple "dR" criteria
// Fedor Ratnikov (UMd), Aug. 28, 2007

#include "RecoJets/JetAssociationAlgorithms/interface/JetTracksAssociationDR.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

// -----------------------------------------------------------------------------
//
JetTracksAssociationDR::JetTracksAssociationDR(double fDr) : mDeltaR2Threshold(fDr * fDr) { ; }

// -----------------------------------------------------------------------------
//
JetTracksAssociationDR::~JetTracksAssociationDR() { ; }

// -----------------------------------------------------------------------------
//
void JetTracksAssociationDR::associateTracksToJets(Association* fAssociation,
                                                   const JetRefs& fJets,
                                                   const TrackRefs& fTracks) {
  JetRefs::const_iterator ii = fJets.begin();
  JetRefs::const_iterator jj = fJets.end();
  for (; ii != jj; ++ii) {
    reco::TrackRefVector associated;
    associateTracksToJet(associated, **ii, fTracks);
    reco::JetTracksAssociation::setValue(fAssociation, *ii, associated);
  }
}

// -----------------------------------------------------------------------------
//
void JetTracksAssociationDR::createJetRefs(JetRefs& output, const Jets& input) {
  output.clear();
  output.reserve(input->size());
  for (unsigned ii = 0; ii < input->size(); ++ii) {
    output.push_back(input->refAt(ii));
  }
}

// -----------------------------------------------------------------------------
//
void JetTracksAssociationDR::createTrackRefs(TrackRefs& output, const Tracks& input, const TrackQuality& quality) {
  if (quality == reco::TrackBase::undefQuality) {
    edm::LogError("JetTracksAssociationDR") << " Unknown TrackQuality value: " << static_cast<int>(quality)
                                            << ". See possible values in 'reco::TrackBase::TrackQuality'";
  }

  output.clear();
  output.reserve(input->size());
  for (unsigned ii = 0; ii < input->size(); ++ii) {
    if ((*input)[ii].quality(quality)) {
      output.push_back(reco::TrackRef(input, ii));
    }
  }
}