Skip to content

Commit

Permalink
Switch from global to inner track for muon IP evaluation
Browse files Browse the repository at this point in the history
  • Loading branch information
Ivan Marchesini committed Aug 7, 2015
1 parent 6bb2ee3 commit 977c94f
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions plugins/BTagAnalyzer.cc
Original file line number Diff line number Diff line change
Expand Up @@ -2148,18 +2148,18 @@ void BTagAnalyzerT<IPTI,VTX>::processJets(const edm::Handle<PatJetCollection>& j
//NOTE: temporary workaround for an issue in https://github.com/cms-btv-pog/cmssw/blob/CMSSW_7_5_X/DataFormats/BTauReco/interface/TemplatedSoftLeptonTagInfo.h#L166
//where the signed IP is actually filled with the IP significance. To be changed when this is fixed in 76X and the IP branches can be filled directly from the taginfos

reco::TrackRef bestTrackmuon = muonPtr->muonBestTrack();
reco::TransientTrack tmuon = trackBuilder->build(bestTrackmuon);

reco::TrackRef innerTrackmuon( muonPtr->innerTrack() );
reco::TransientTrack tmuon = trackBuilder->build(innerTrackmuon);
GlobalVector directionformuon(pjet->px(), pjet->py(), pjet->pz());

Measurement1D ip2dmuon = IPTools::signedTransverseImpactParameter(tmuon, directionformuon, *pv).second;
Measurement1D ip3dmuon = IPTools::signedImpactParameter3D(tmuon, directionformuon, *pv).second;

JetInfo[iJetColl].PFMuon_IP[JetInfo[iJetColl].nPFMuon] = (ip3dmuon.value());
JetInfo[iJetColl].PFMuon_IP2D[JetInfo[iJetColl].nPFMuon] = (ip2dmuon.value());
JetInfo[iJetColl].PFMuon_IPsig[JetInfo[iJetColl].nPFMuon] = (ip3dmuon.value())/(ip3dmuon.error());
JetInfo[iJetColl].PFMuon_IP2Dsig[JetInfo[iJetColl].nPFMuon] = (ip2dmuon.value())/(ip2dmuon.error());
JetInfo[iJetColl].PFMuon_IPsig[JetInfo[iJetColl].nPFMuon] = (ip3dmuon.significance());
JetInfo[iJetColl].PFMuon_IP2Dsig[JetInfo[iJetColl].nPFMuon] = (ip2dmuon.significance());

JetInfo[iJetColl].PFMuon_nMuHit[JetInfo[iJetColl].nPFMuon] = muonPtr->outerTrack()->hitPattern().numberOfValidMuonHits();
JetInfo[iJetColl].PFMuon_nTkHit[JetInfo[iJetColl].nPFMuon] = muonPtr->innerTrack()->hitPattern().numberOfValidHits();
Expand Down

0 comments on commit 977c94f

Please sign in to comment.