Skip to content

Commit

Permalink
---
Browse files Browse the repository at this point in the history
yaml
---
r: 32744
b: "refs/heads/CMSSW_7_1_X"
c: 845cee8
h: "refs/heads/CMSSW_7_1_X"
v: v3
  • Loading branch information
Andrea Bocci committed Oct 24, 2007
1 parent 875bb31 commit 3bf830b
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 7 deletions.
2 changes: 1 addition & 1 deletion [refs]
Original file line number Diff line number Diff line change
@@ -1,3 +1,3 @@
---
refs/heads/gh-pages: 09c786f70121f131b3715aaf3464996502bbeb7e
"refs/heads/CMSSW_7_1_X": 0fe471bd94cb1310bf1296bafff618cff7ec3532
"refs/heads/CMSSW_7_1_X": 845cee8e93ae68ce0119971274b359909b4a5a8f
12 changes: 6 additions & 6 deletions trunk/RecoBTag/ImpactParameter/plugins/TrackIPProducer.cc
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,7 @@
//
// Original Author: Andrea Rizzi
// Created: Thu Apr 6 09:56:23 CEST 2006
// $Id: TrackIPProducer.cc,v 1.8 2007/10/08 18:45:22 arizzi Exp $
// $Id: TrackIPProducer.cc,v 1.9 2007/10/12 09:31:37 arizzi Exp $
//
//

Expand Down Expand Up @@ -162,7 +162,7 @@ TrackIPProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
int ind =0;
for (TrackRefVector::const_iterator itTrack = tracks.begin(); itTrack != tracks.end(); ++itTrack) {
const Track & track = **itTrack;
const TransientTrack & transientTrack = builder->build(&(**itTrack));
const TransientTrack & transientTrack = builder->build(&(**itTrack));
/* cout << " pt " << track.pt() <<
" d0 " << fabs(track.d0()) <<
" #hit " << track.hitPattern().numberOfValidHits()<<
Expand All @@ -188,7 +188,7 @@ TrackIPProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
trackIp.ip3d=IPTools::signedImpactParameter3D(transientTrack,direction,*pv).second;
trackIp.ip2d=IPTools::signedTransverseImpactParameter(transientTrack,direction,*pv).second;
TrajectoryStateOnSurface closest = IPTools::closestApproachToJet(transientTrack.impactPointState(), *pv, direction,transientTrack.field());
if(closes.isValid()) trackIp.closestToJetAxis=closest.globalPosition();
if (closest.isValid()) trackIp.closestToJetAxis = closest.globalPosition();
//TODO:cross check if it is the same using other methods
trackIp.distanceToJetAxis=IPTools::jetTrackDistance(transientTrack,direction,*pv).second.value();

Expand Down Expand Up @@ -221,17 +221,17 @@ TrackIPProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
int second=last->second;

for(int n=0;n< ipData.size();n++)
{
{
int use;
if(n==first) use = second; else use = first;
TrajectoryStateOnSurface trackState1 = builder->build(selectedTracks[n]).impactPointState();
TrajectoryStateOnSurface trackState2 = builder->build(selectedTracks[use]).impactPointState();
std::pair<GlobalPoint,GlobalPoint> points = minDist.points(trackState1,trackState2);
float distance = ( points.first - points.second ).mag();
ipData[n].closestToFirstTrack=points.first;
ipData[n].closestToFirstTrack=points.first;
ipData[n].distanceToFirstTrack=distance;

}
}
}
TrackIPTagInfo tagInfo(ipData,prob2D,prob3D,selectedTracks,
edm::Ref<JetTracksAssociationCollection>(jetTracksAssociation,i),
Expand Down

0 comments on commit 3bf830b

Please sign in to comment.