JetTracksAssociationDRVertexAssigned Class Reference

#include <JetTracksAssociationDRVertexAssigned.h>

Public Member Functions
	JetTracksAssociationDRVertexAssigned (double fDr)
void	produce (reco::JetTracksAssociation::Container *fAssociation, const std::vector< edm::RefToBase< reco::Jet > > &fJets, const std::vector< reco::TrackRef > &fTracks, const reco::VertexCollection &vertices) const
	~JetTracksAssociationDRVertexAssigned ()

Private Attributes
double	mDeltaR2Threshold
	fidutial dR between track in the vertex and jet's reference direction More...

Detailed Description

Definition at line 13 of file JetTracksAssociationDRVertexAssigned.h.

Constructor & Destructor Documentation

JetTracksAssociationDRVertexAssigned::JetTracksAssociationDRVertexAssigned

(

double

fDr

)

Definition at line 14 of file JetTracksAssociationDRVertexAssigned.cc.

: mDeltaR2Threshold(fDr * fDr) {}

JetTracksAssociationDRVertexAssigned::~JetTracksAssociationDRVertexAssigned ( )

inline

Definition at line 16 of file JetTracksAssociationDRVertexAssigned.h.

{}

Member Function Documentation

void JetTracksAssociationDRVertexAssigned::produce	(	reco::JetTracksAssociation::Container *	fAssociation,
		const std::vector< edm::RefToBase< reco::Jet > > &	fJets,
		const std::vector< reco::TrackRef > &	fTracks,
		const reco::VertexCollection &	vertices
	)		const

Definition at line 16 of file JetTracksAssociationDRVertexAssigned.cc.

References reco::deltaR2(), reco::TrackBase::dz(), PVValHelper::eta, reco::LeafCandidate::eta(), reco::TrackBase::eta(), spr::find(), mps_fire::i, gpuVertexFinder::iv, dqmiolumiharvest::j, metsig::jet, reco::btau::jetEta, reco::btau::jetPhi, mDeltaR2Threshold, reco::TrackBase::p(), phi, reco::LeafCandidate::phi(), reco::TrackBase::phi(), reco::Jet::physicsEta(), edm::RefVector< C, T, F >::push_back(), findQualityFiles::rr, reco::JetTracksAssociation::setValue(), submitPVValidationJobs::t, and HLT_FULL_cff::track.

Referenced by JetTracksAssociatorAtVertex::produce().

                                                                                               {
  // cache tracks kinematics
  std::vector<math::RhoEtaPhiVector> trackP3s;
  std::map<int, double> trackvert;

  // std::cout<<" Number of vertices "<<vertices.size()<<std::endl;

  trackP3s.reserve(fTracks.size());
  for (unsigned i = 0; i < fTracks.size(); ++i) {
    const reco::Track* track = &*(fTracks[i]);
    trackP3s.push_back(math::RhoEtaPhiVector(track->p(), track->eta(), track->phi()));

    // OK: Look for the tracks not associated with vertices

    const reco::TrackBaseRef ttr1(fTracks[i]);

    int trackhasvert = -1;
    for (reco::VertexCollection::const_iterator iv = vertices.begin(); iv != vertices.end(); iv++) {
      std::vector<reco::TrackBaseRef>::const_iterator rr = find((*iv).tracks_begin(), (*iv).tracks_end(), ttr1);
      if (rr != (*iv).tracks_end()) {
        trackhasvert = 1;
        trackvert[i] = (*iv).position().z();
        // std::cout<<" Z "<<i<<" "<<trackhasvert<<" "<<(*iv).position().z()<<std::endl;
        break;
      }
    }  // all vertices
    if (trackhasvert < 0) {
      // Take impact parameter of the track as vertex position
      math::XYZPoint ppt(0., 0., 0.);
      trackvert[i] = track->dz(ppt);
      //    std::cout<<" Z "<<i<<" "<<trackhasvert<<" "<<track->dz(ppt)<<" "<<track->vz()<<" "<<track->vx()<<" "<<
      //    track->vy()<<" "<<track->pz()<<std::endl;
    }
    // OK
  }  // tracks

  for (unsigned j = 0; j < fJets.size(); ++j) {
    reco::TrackRefVector assoTracks;
    const reco::Jet* jet = &*(fJets[j]);
    double jetEta = jet->eta();
    double jetPhi = jet->phi();
    double neweta = 0;
    for (unsigned t = 0; t < fTracks.size(); ++t) {
      std::map<int, double>::iterator cur = trackvert.find(t);
      if (cur != trackvert.end()) {
        neweta = jet->physicsEta((*cur).second, jetEta);
      } else {
        neweta = jetEta;
        //std::cout<<" Lost track - not in map "<<std::endl;
      }

      //std::cout<<" Old eta-new eta "<<(*cur).second<<" "<<jetEta<<" "<<neweta<<" Track "<<t<<" "<<trackP3s[t].eta()<<std::endl;

      double dR2 = deltaR2(neweta, jetPhi, trackP3s[t].eta(), trackP3s[t].phi());
      if (dR2 < mDeltaR2Threshold)
        assoTracks.push_back(fTracks[t]);
    }

    reco::JetTracksAssociation::setValue(fAssociation, fJets[j], assoTracks);
  }
}

Member Data Documentation

double JetTracksAssociationDRVertexAssigned::mDeltaR2Threshold

private

fidutial dR between track in the vertex and jet's reference direction

Definition at line 25 of file JetTracksAssociationDRVertexAssigned.h.

Referenced by produce().