CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
Classes | Functions
GhostTrackVertexFinder.cc File Reference
#include <algorithm>
#include <iterator>
#include <memory>
#include <map>
#include <set>
#include <vector>
#include <Math/SVector.h>
#include <Math/SMatrix.h>
#include <Math/MatrixFunctions.h>
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/GlobalError.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "TrackingTools/TrajectoryParametrization/interface/GlobalTrajectoryParameters.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/GeomPropagators/interface/AnalyticalImpactPointExtrapolator.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "TrackingTools/TransientTrack/interface/TransientTrackFromFTSFactory.h"
#include "RecoVertex/VertexPrimitives/interface/TransientVertex.h"
#include "RecoVertex/VertexPrimitives/interface/VertexFitter.h"
#include "RecoVertex/VertexPrimitives/interface/VertexState.h"
#include "RecoVertex/VertexPrimitives/interface/CachingVertex.h"
#include "RecoVertex/VertexPrimitives/interface/LinearizedTrackState.h"
#include "RecoVertex/VertexPrimitives/interface/ConvertError.h"
#include "RecoVertex/VertexPrimitives/interface/ConvertToFromReco.h"
#include "RecoVertex/VertexTools/interface/LinearizedTrackStateFactory.h"
#include "RecoVertex/VertexTools/interface/VertexTrackFactory.h"
#include "RecoVertex/VertexTools/interface/VertexDistance3D.h"
#include "RecoVertex/VertexTools/interface/GeometricAnnealing.h"
#include "RecoVertex/KalmanVertexFit/interface/KalmanVertexFitter.h"
#include "RecoVertex/AdaptiveVertexFit/interface/AdaptiveVertexFitter.h"
#include "RecoVertex/GhostTrackFitter/interface/GhostTrack.h"
#include "RecoVertex/GhostTrackFitter/interface/GhostTrackState.h"
#include "RecoVertex/GhostTrackFitter/interface/GhostTrackFitter.h"
#include "RecoVertex/GhostTrackFitter/interface/GhostTrackPrediction.h"
#include "RecoVertex/GhostTrackFitter/interface/SequentialGhostTrackFitter.h"
#include "RecoVertex/GhostTrackFitter/interface/KalmanGhostTrackUpdater.h"
#include "RecoVertex/GhostTrackFitter/interface/GhostTrackVertexFinder.h"

Go to the source code of this file.

Classes

struct  reco::GhostTrackVertexFinder::FinderInfo
 

Functions

static bool covarianceUpdate (Matrix3S &cov, const Vector3 &residual, const Matrix3S &error, double &chi2, double theta, double phi)
 
static GhostTrackPrediction dummyPrediction (const Vertex &primaryVertex, const Track &ghostTrack)
 
static double fitChi2 (const CachingVertex< 5 > &vtx)
 
static void mergeTrackHelper (const std::vector< RefCountedVertexTrack > &tracks, std::vector< RefCountedVertexTrack > &newTracks, const VertexState &state, const VtxTrackIs &ghostTrackFinder, RefCountedVertexTrack &ghostTrack, const VertexTrackFactory< 5 > &factory)
 
static RefCountedVertexTrack relinearizeTrack (const RefCountedVertexTrack &track, const VertexState &state, const VertexTrackFactory< 5 > factory)
 
static std::vector
< RefCountedVertexTrack > 		relinearizeTracks (const std::vector< RefCountedVertexTrack > &tracks, const VertexState &state)
 
static VertexState stateMean (const VertexState &v1, const VertexState &v2)
 
static double trackVertexCompat (const CachingVertex< 5 > &vtx, const RefCountedVertexTrack &vertexTrack)
 
static TransientTrack transientGhostTrack (const GhostTrackPrediction &pred, const MagneticField *field)
 
static CachingVertex< 5 > vertexAtState (const TransientTrack &ghostTrack, const GhostTrackPrediction &pred, const GhostTrackState &state)
 
static double vtxErrorLong (const GlobalError &error, const GlobalVector &dir)
 
static GlobalPoint vtxMean (const GlobalPoint &p1, const GlobalError &e1, const GlobalPoint &p2, const GlobalError &e2)
 

Function Documentation

static bool covarianceUpdate ( Matrix3S &  cov,
const Vector3 &  residual,
const Matrix3S &  error,
double &  chi2,
double  theta,
double  phi 
)
static

Definition at line 138 of file GhostTrackVertexFinder.cc.

References funct::cos(), and funct::sin().

Referenced by vertexAtState().

139  {
140  using namespace ROOT::Math;
141 
142  Matrix23 jacobian;
143  jacobian(0, 0) = std::cos(phi) * std::cos(theta);
144  jacobian(0, 1) = std::sin(phi) * std::cos(theta);
145  jacobian(0, 2) = -std::sin(theta);
146  jacobian(1, 0) = -std::sin(phi);
147  jacobian(1, 1) = std::cos(phi);
148 
149  Matrix2S measErr = Similarity(jacobian, error);
150  Matrix2S combErr = Similarity(jacobian, cov) + measErr;
151  if (!measErr.Invert() || !combErr.Invert())
152  return false;
153 
154  cov -= SimilarityT(jacobian * cov, combErr);
155  chi2 += Similarity(jacobian * residual, measErr);
156 
157  return true;
158 }
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:150
HLT_FULL_cff.chi2
tuple chi2
Definition: HLT_FULL_cff.py:10024
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
static GhostTrackPrediction dummyPrediction ( const Vertex &  primaryVertex,
const Track &  ghostTrack 
)
static

Definition at line 535 of file GhostTrackVertexFinder.cc.

References RecoVertex::convertError(), RecoVertex::convertPos(), reco::Vertex::error(), reco::Vertex::position(), reco::TrackBase::px(), reco::TrackBase::py(), and reco::TrackBase::pz().

Referenced by reco::GhostTrackVertexFinder::vertices().

535  {
536  return GhostTrackPrediction(RecoVertex::convertPos(primaryVertex.position()),
537  RecoVertex::convertError(primaryVertex.error()),
538  GlobalVector(ghostTrack.px(), ghostTrack.py(), ghostTrack.pz()),
539  0.05);
540 }
RecoVertex::convertPos
reco::Vertex::Point convertPos(const GlobalPoint &p)
Definition: ConvertToFromReco.h:7
RecoVertex::convertError
reco::Vertex::Error convertError(const GlobalError &ge)
Definition: ConvertError.h:8
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:640
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:127
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:646
reco::Vertex::error
Error error() const
return SMatrix
Definition: Vertex.h:163
reco::TrackBase::py
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:643
GlobalVector
Global3DVector GlobalVector
Definition: GlobalVector.h:10
static double fitChi2 ( const CachingVertex< 5 > &  vtx)
static

Definition at line 638 of file GhostTrackVertexFinder.cc.

References CachingVertex< N >::degreesOfFreedom(), and CachingVertex< N >::totalChiSquared().

Referenced by PFGsfHelper::PFGsfHelper(), TtEvtBuilder< C >::produce(), and reco::GhostTrackVertexFinder::recursiveMerge().

638 { return vtx.totalChiSquared() / vtx.degreesOfFreedom(); }
CachingVertex::totalChiSquared
float totalChiSquared() const
Definition: CachingVertex.h:159
CachingVertex::degreesOfFreedom
float degreesOfFreedom() const
Definition: CachingVertex.cc:348
static void mergeTrackHelper ( const std::vector< RefCountedVertexTrack > &  tracks,
std::vector< RefCountedVertexTrack > &  newTracks,
const VertexState &  state,
const VtxTrackIs &  ghostTrackFinder,
RefCountedVertexTrack &  ghostTrack,
const VertexTrackFactory< 5 > &  factory 
)
static

Definition at line 661 of file GhostTrackVertexFinder.cc.

References submitPVValidationJobs::gt, and relinearizeTrack().

Referenced by reco::GhostTrackVertexFinder::mergeVertices().

666  {
667  for (std::vector<RefCountedVertexTrack>::const_iterator iter = tracks.begin(); iter != tracks.end(); ++iter) {
668  bool gt = ghostTrackFinder(*iter);
669  if (gt && ghostTrack)
670  continue;
671 
672  RefCountedVertexTrack track = relinearizeTrack(*iter, state, factory);
673 
674  if (gt)
675  ghostTrack = *iter;
676  else
677  newTracks.push_back(*iter);
678  }
679 }
relinearizeTrack
static RefCountedVertexTrack relinearizeTrack(const RefCountedVertexTrack &track, const VertexState &state, const VertexTrackFactory< 5 > factory)
Definition: GhostTrackVertexFinder.cc:211
tracks
auto const & tracks
cannot be loose
Definition: CAHitNtupletGeneratorKernelsImpl.h:146
HLT_FULL_cff.track
tuple track
Definition: HLT_FULL_cff.py:11953
static RefCountedVertexTrack relinearizeTrack ( const RefCountedVertexTrack &  track,
const VertexState &  state,
const VertexTrackFactory< 5 >  factory 
)
static

Definition at line 211 of file GhostTrackVertexFinder.cc.

References VertexState::position(), and VertexTrackFactory< N >::vertexTrack().

Referenced by mergeTrackHelper(), reco::GhostTrackVertexFinder::reassignTracks(), and relinearizeTracks().

213  {
214  RefCountedLinearizedTrackState linTrack = track->linearizedTrack();
215  linTrack = linTrack->stateWithNewLinearizationPoint(state.position());
216  return factory.vertexTrack(linTrack, state);
217 }
VertexTrackFactory::vertexTrack
RefCountedVertexTrack vertexTrack(const RefCountedLinearizedTrackState lt, const VertexState vs, float weight=1.0) const
Definition: VertexTrackFactory.h:27
HLT_FULL_cff.track
tuple track
Definition: HLT_FULL_cff.py:11953
VertexState::position
GlobalPoint position() const
Definition: VertexState.h:62
static std::vector<RefCountedVertexTrack> relinearizeTracks ( const std::vector< RefCountedVertexTrack > &  tracks,
const VertexState &  state 
)
static

Definition at line 219 of file GhostTrackVertexFinder.cc.

References relinearizeTrack().

Referenced by reco::GhostTrackVertexFinder::reassignTracks().

220  {
221  VertexTrackFactory<5> vertexTrackFactory;
222 
223  std::vector<RefCountedVertexTrack> finalTracks;
224  finalTracks.reserve(tracks.size());
225 
226  for (std::vector<RefCountedVertexTrack>::const_iterator iter = tracks.begin(); iter != tracks.end(); ++iter)
227  finalTracks.push_back(relinearizeTrack(*iter, state, vertexTrackFactory));
228 
229  return finalTracks;
230 }
relinearizeTrack
static RefCountedVertexTrack relinearizeTrack(const RefCountedVertexTrack &track, const VertexState &state, const VertexTrackFactory< 5 > factory)
Definition: GhostTrackVertexFinder.cc:211
tracks
auto const & tracks
cannot be loose
Definition: CAHitNtupletGeneratorKernelsImpl.h:146
static VertexState stateMean ( const VertexState &  v1,
const VertexState &  v2 
)
static

Definition at line 134 of file GhostTrackVertexFinder.cc.

References VertexState::error(), VertexState::position(), and vtxMean().

Referenced by reco::GhostTrackVertexFinder::mergeVertices().

134  {
135  return VertexState(vtxMean(v1.position(), v1.error(), v2.position(), v2.error()), v1.error() + v2.error());
136 }
VertexState::position
GlobalPoint position() const
Definition: VertexState.h:62
vtxMean
static GlobalPoint vtxMean(const GlobalPoint &p1, const GlobalError &e1, const GlobalPoint &p2, const GlobalError &e2)
Definition: GhostTrackVertexFinder.cc:123
VertexState::error
GlobalError error() const
Definition: VertexState.h:64
static double trackVertexCompat ( const CachingVertex< 5 > &  vtx,
const RefCountedVertexTrack &  vertexTrack 
)
static

Definition at line 240 of file GhostTrackVertexFinder.cc.

References conv, DeadROC_duringRun::dir, relativeConstraints::error, CachingVertex< N >::error(), GlobalErrorBase< T, ErrorWeightType >::matrix(), point, CachingVertex< N >::position(), and HLT_FULL_cff::track.

Referenced by reco::GhostTrackVertexFinder::reassignTracks().

240  {
241  using namespace ROOT::Math;
242 
243  TransientTrack track = vertexTrack->linearizedTrack()->track();
244  GlobalPoint point = vtx.position();
245  AnalyticalImpactPointExtrapolator extrap(track.field());
246  TrajectoryStateOnSurface tsos = extrap.extrapolate(track.impactPointState(), point);
247 
248  if (!tsos.isValid())
249  return 1.0e6;
250 
251  GlobalPoint point1 = vtx.position();
252  GlobalPoint point2 = tsos.globalPosition();
253  Vector3 dir = conv(point2 - point1);
254  Matrix3S error = vtx.error().matrix() + tsos.cartesianError().matrix().Sub<Matrix3S>(0, 0);
255  if (!error.Invert())
256  return 1.0e6;
257 
258  return ROOT::Math::Similarity(error, dir);
259 }
GlobalErrorBase::matrix
const AlgebraicSymMatrix33 matrix() const
Definition: GlobalErrorBase.h:121
HLT_FULL_cff.track
tuple track
Definition: HLT_FULL_cff.py:11953
DeadROC_duringRun.dir
string dir
Definition: DeadROC_duringRun.py:23
conv
EPOS::IO_EPOS conv
Definition: ReggeGribovPartonMCHadronizer.cc:42
CachingVertex::position
GlobalPoint position() const
Definition: CachingVertex.cc:308
CachingVertex::error
GlobalError error() const
Definition: CachingVertex.cc:318
point
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
static TransientTrack transientGhostTrack ( const GhostTrackPrediction &  pred,
const MagneticField *  field 
)
static

Definition at line 115 of file GhostTrackVertexFinder.cc.

References TransientTrackFromFTSFactory::build(), and reco::GhostTrackPrediction::fts().

Referenced by reco::GhostTrackVertexFinder::refitGhostTrack(), and reco::GhostTrackVertexFinder::vertices().

115  {
116  return TransientTrackFromFTSFactory().build(pred.fts(field));
117 }
reco::GhostTrackPrediction::fts
FreeTrajectoryState fts(const MagneticField *fieldProvider) const
Definition: GhostTrackPrediction.cc:185
TransientTrackFromFTSFactory::build
reco::TransientTrack build(const FreeTrajectoryState &fts) const
Definition: TransientTrackFromFTSFactory.cc:7
static CachingVertex<5> vertexAtState ( const TransientTrack &  ghostTrack,
const GhostTrackPrediction &  pred,
const GhostTrackState &  state 
)
static

Definition at line 160 of file GhostTrackVertexFinder.cc.

References reco::GhostTrackState::cartesianError(), HLT_FULL_cff::chi2, conv, covarianceUpdate(), relativeConstraints::error, reco::GhostTrackState::globalPosition(), reco::GhostTrackState::isValid(), sistrip::SpyUtilities::isValid(), reco::GhostTrackState::lambda(), LinearizedTrackStateFactory::linearizedTrackState(), GlobalErrorBase< T, ErrorWeightType >::matrix(), point, reco::GhostTrackPrediction::position(), reco::GhostTrackPrediction::positionError(), reco::GhostTrackState::track(), VertexTrackFactory< N >::vertexTrack(), and vtxMean().

Referenced by reco::GhostTrackVertexFinder::initialVertices(), reco::GhostTrackVertexFinder::reassignTracks(), and reco::GhostTrackVertexFinder::refitGhostTrack().

162  {
163  LinearizedTrackStateFactory linTrackFactory;
164  VertexTrackFactory<5> vertexTrackFactory;
165 
166  if (!state.isValid())
167  return CachingVertex<5>();
168 
169  GlobalPoint pca1 = pred.position(state.lambda());
170  GlobalError err1 = pred.positionError(state.lambda());
171 
172  GlobalPoint pca2 = state.globalPosition();
173  GlobalError err2 = state.cartesianError();
174 
175  GlobalPoint point = vtxMean(pca1, err1, pca2, err2);
176 
177  const TransientTrack &recTrack = state.track();
178 
179  RefCountedLinearizedTrackState linState[2] = {linTrackFactory.linearizedTrackState(point, ghostTrack),
180  linTrackFactory.linearizedTrackState(point, recTrack)};
181  if (!linState[0]->isValid() || !linState[1]->isValid())
182  return CachingVertex<5>();
183 
184  Matrix3S cov = SMatrixIdentity();
185  cov *= 10000;
186  double chi2 = 0.;
187  if (!covarianceUpdate(cov,
188  conv(pca1 - point),
189  err1.matrix(),
190  chi2,
191  linState[0]->predictedStateParameters()[1],
192  linState[0]->predictedStateParameters()[2]) ||
193  !covarianceUpdate(cov,
194  conv(pca2 - point),
195  err2.matrix(),
196  chi2,
197  linState[1]->predictedStateParameters()[1],
198  linState[1]->predictedStateParameters()[2]))
199  return CachingVertex<5>();
200 
201  GlobalError error(cov);
202  VertexState vtxState(point, error);
203 
204  std::vector<RefCountedVertexTrack> linTracks(2);
205  linTracks[0] = vertexTrackFactory.vertexTrack(linState[0], vtxState);
206  linTracks[1] = vertexTrackFactory.vertexTrack(linState[1], vtxState);
207 
208  return CachingVertex<5>(point, error, linTracks, chi2);
209 }
LinearizedTrackStateFactory::linearizedTrackState
RefCountedLinearizedTrackState linearizedTrackState(const GlobalPoint &linP, const reco::TransientTrack &track) const override
Definition: LinearizedTrackStateFactory.cc:9
sistrip::SpyUtilities::isValid
const bool isValid(const Frame &aFrame, const FrameQuality &aQuality, const uint16_t aExpectedPos)
Definition: SiStripSpyUtilities.cc:124
covarianceUpdate
static bool covarianceUpdate(Matrix3S &cov, const Vector3 &residual, const Matrix3S &error, double &chi2, double theta, double phi)
Definition: GhostTrackVertexFinder.cc:138
reco::GhostTrackState::track
const TransientTrack & track() const
Definition: GhostTrackState.cc:56
VertexTrackFactory::vertexTrack
RefCountedVertexTrack vertexTrack(const RefCountedLinearizedTrackState lt, const VertexState vs, float weight=1.0) const
Definition: VertexTrackFactory.h:27
GlobalErrorBase::matrix
const AlgebraicSymMatrix33 matrix() const
Definition: GlobalErrorBase.h:121
reco::GhostTrackState::cartesianError
GlobalError cartesianError() const
Definition: GhostTrackState.h:37
HLT_FULL_cff.chi2
tuple chi2
Definition: HLT_FULL_cff.py:10024
vtxMean
static GlobalPoint vtxMean(const GlobalPoint &p1, const GlobalError &e1, const GlobalPoint &p2, const GlobalError &e2)
Definition: GhostTrackVertexFinder.cc:123
reco::GhostTrackState::lambda
double lambda() const
Definition: GhostTrackState.h:40
conv
EPOS::IO_EPOS conv
Definition: ReggeGribovPartonMCHadronizer.cc:42
reco::GhostTrackPrediction::position
GlobalPoint position(double lambda=0.) const
Definition: GhostTrackPrediction.h:67
reco::GhostTrackState::globalPosition
GlobalPoint globalPosition() const
Definition: GhostTrackState.h:36
point
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
Definition: invegas.h:5
reco::GhostTrackPrediction::positionError
GlobalError positionError(double lambda=0.) const
Definition: GhostTrackPrediction.cc:131
static double vtxErrorLong ( const GlobalError &  error,
const GlobalVector &  dir 
)
static

Definition at line 119 of file GhostTrackVertexFinder.cc.

References conv, and GlobalErrorBase< T, ErrorWeightType >::matrix().

Referenced by vtxMean().

119  {
120  return ROOT::Math::Similarity(conv(dir), error.matrix());
121 }
GlobalErrorBase::matrix
const AlgebraicSymMatrix33 matrix() const
Definition: GlobalErrorBase.h:121
conv
EPOS::IO_EPOS conv
Definition: ReggeGribovPartonMCHadronizer.cc:42
static GlobalPoint vtxMean ( const GlobalPoint &  p1,
const GlobalError &  e1,
const GlobalPoint &  p2,
const GlobalError &  e2 
)
static

Definition at line 123 of file GhostTrackVertexFinder.cc.

References change_name::diff, fireworks::p1, vtxErrorLong(), and histoStyle::weight.

Referenced by GlobalTrackingRegion::checkRZ(), RectangularEtaPhiTrackingRegion::checkRZOld(), stateMean(), and vertexAtState().

123  {
124  GlobalVector diff = p2 - p1;
125 
126  double err1 = vtxErrorLong(e1, diff);
127  double err2 = vtxErrorLong(e2, diff);
128 
129  double weight = err1 / (err1 + err2);
130 
131  return p1 + weight * diff;
132 }
change_name.diff
list diff
Definition: change_name.py:13
fireworks::p1
const TString p1
Definition: fwPaths.cc:12
vtxErrorLong
static double vtxErrorLong(const GlobalError &error, const GlobalVector &dir)
Definition: GhostTrackVertexFinder.cc:119
histoStyle.weight
int weight
Definition: histoStyle.py:51
Generated for CMSSW Reference Manual by doxygen 1.8.5