CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
JetTracksAssociationDRVertexAssigned.cc
Go to the documentation of this file.
1 // Associate jets with tracks by simple "dR" criteria
2 // Fedor Ratnikov (UMd), Aug. 28, 2007
3 
5 
8 
13 
15 
18  const std::vector<reco::TrackRef>& fTracks,
19  const reco::VertexCollection& vertices) const {
20  // cache tracks kinematics
21  std::vector<math::RhoEtaPhiVector> trackP3s;
22  std::map<int, double> trackvert;
23 
24  // std::cout<<" Number of vertices "<<vertices.size()<<std::endl;
25 
26  trackP3s.reserve(fTracks.size());
27  for (unsigned i = 0; i < fTracks.size(); ++i) {
28  const reco::Track* track = &*(fTracks[i]);
29  trackP3s.push_back(math::RhoEtaPhiVector(track->p(), track->eta(), track->phi()));
30 
31  // OK: Look for the tracks not associated with vertices
32 
33  const reco::TrackBaseRef ttr1(fTracks[i]);
34 
35  int trackhasvert = -1;
36  for (reco::VertexCollection::const_iterator iv = vertices.begin(); iv != vertices.end(); iv++) {
37  std::vector<reco::TrackBaseRef>::const_iterator rr = find((*iv).tracks_begin(), (*iv).tracks_end(), ttr1);
38  if (rr != (*iv).tracks_end()) {
39  trackhasvert = 1;
40  trackvert[i] = (*iv).position().z();
41  // std::cout<<" Z "<<i<<" "<<trackhasvert<<" "<<(*iv).position().z()<<std::endl;
42  break;
43  }
44  } // all vertices
45  if (trackhasvert < 0) {
46  // Take impact parameter of the track as vertex position
47  math::XYZPoint ppt(0., 0., 0.);
48  trackvert[i] = track->dz(ppt);
49  // std::cout<<" Z "<<i<<" "<<trackhasvert<<" "<<track->dz(ppt)<<" "<<track->vz()<<" "<<track->vx()<<" "<<
50  // track->vy()<<" "<<track->pz()<<std::endl;
51  }
52  // OK
53  } // tracks
54 
55  for (unsigned j = 0; j < fJets.size(); ++j) {
56  reco::TrackRefVector assoTracks;
57  const reco::Jet* jet = &*(fJets[j]);
58  double jetEta = jet->eta();
59  double jetPhi = jet->phi();
60  double neweta = 0;
61  for (unsigned t = 0; t < fTracks.size(); ++t) {
62  std::map<int, double>::iterator cur = trackvert.find(t);
63  if (cur != trackvert.end()) {
64  neweta = jet->physicsEta((*cur).second, jetEta);
65  } else {
66  neweta = jetEta;
67  //std::cout<<" Lost track - not in map "<<std::endl;
68  }
69 
70  //std::cout<<" Old eta-new eta "<<(*cur).second<<" "<<jetEta<<" "<<neweta<<" Track "<<t<<" "<<trackP3s[t].eta()<<std::endl;
71 
72  double dR2 = deltaR2(neweta, jetPhi, trackP3s[t].eta(), trackP3s[t].phi());
73  if (dR2 < mDeltaR2Threshold)
74  assoTracks.push_back(fTracks[t]);
75  }
76 
77  reco::JetTracksAssociation::setValue(fAssociation, fJets[j], assoTracks);
78  }
79 }
double p() const
momentum vector magnitude
Definition: TrackBase.h:631
int32_t *__restrict__ iv
Base class for all types of Jets.
Definition: Jet.h:20
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:649
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
double mDeltaR2Threshold
fidutial dR between track in the vertex and jet&#39;s reference direction
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
RhoEtaPhiVectorD RhoEtaPhiVector
spatial vector with cylindrical internal representation using pseudorapidity
Definition: Vector3D.h:33
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
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
bool setValue(Container &, const reco::JetBaseRef &, reco::TrackRefVector)
associate jet with value. Returns false and associate nothing if jet is already associated ...
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:622
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
void push_back(value_type const &ref)
Add a Ref&lt;C, T&gt; to the RefVector.
Definition: RefVector.h:67
double phi() const final
momentum azimuthal angle
static float physicsEta(float fZVertex, float fDetectorEta)
static function to convert detector eta to physics eta
double eta() const final
momentum pseudorapidity