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 18 of file JetTracksAssociationDRVertexAssigned.cc.

References reco::deltaR2(), reco::TrackBase::dz(), eta(), reco::LeafCandidate::eta(), reco::TrackBase::eta(), spr::find(), i, 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(), and edmStreamStallGrapher::t.

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 24 of file JetTracksAssociationDRVertexAssigned.h.

Referenced by produce().