22 std::vector<reco::TransientTrack>
result;
23 TwoTrackMinimumDistance dist;
27 float pvDistance = ipSeed.second.value();
29 float densityFactor = 2./
sqrt(20.*80);
30 for(std::vector<reco::TransientTrack>::const_iterator
tt =
tracks.begin();
tt!=
tracks.end(); ++
tt ) {
32 if(*
tt==seed)
continue;
38 GlobalError ttPointErr =
tt->impactPointState().cartesianError().position();
41 Measurement1D m = distanceComputer.
distance(VertexState(seedPosition,seedPositionErr), VertexState(ttPoint, ttPointErr));
45 float distanceFromPV = (dist.points().second-pv).
mag();
46 float distance = dist.distance();
47 GlobalVector trackDir2D(
tt->impactPointState().globalDirection().x(),
tt->impactPointState().globalDirection().y(),0.);
51 float dotprodTrack = (dist.points().first-pv).
unit().dot(
tt->impactPointState().globalDirection().unit());
54 float w = distanceFromPV*distanceFromPV/(pvDistance*distance);
68 distance*
clusterScale <<
" < " << densityFactor*distanceFromPV <<
" crossingtoPV: " << distanceFromPV <<
" dis*scal " << distance*
clusterScale <<
" < " << densityFactor*distanceFromPV <<
" dist: " << distance <<
" < " <<
clusterMaxDistance << std::endl;
72 result.push_back(*
tt);
73 seedingPoint =
GlobalPoint(
cp.x()*w+seedingPoint.
x(),
cp.y()*w+seedingPoint.
y(),
cp.z()*w+seedingPoint.
z());
79 seedingPoint =
GlobalPoint(seedingPoint.
x()/sumWeights,seedingPoint.
y()/sumWeights,seedingPoint.
z()/sumWeights);
80 return std::pair<std::vector<reco::TransientTrack>,
GlobalPoint>(
result,seedingPoint);
double clusterMaxSignificance
virtual Measurement1D distance(const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Global3DPoint GlobalPoint
const CartesianTrajectoryError cartesianError() const
const Point & position() const
position
double clusterMaxDistance
double clusterMinAngleCosine
Vector3DBase unit() const
double significance() const
const GlobalError position() const
Position error submatrix.
TrajectoryStateOnSurface impactPointState() const
GlobalVector globalDirection() const