21 minCrossingRadius = pset.
getParameter<
double>(
"minCrossingRadius");
22 maxCrossingRadius = pset.
getParameter<
double>(
"maxCrossingRadius");
23 maxImpactMother = pset.
getParameter<
double>(
"maxImpactMother");
28 theMagField = magField.product();
49 track.
charge(),theMagField);
84 pair<GlobalPoint, GlobalPoint> points = theMinimum.
points();
87 pair<GlobalVector,GlobalVector> momenta;
104 crossing.
perp() > minCrossingRadius &&
105 crossing.
perp() < maxCrossingRadius)
111 if(vertices->size() > 0)
114 for(reco::VertexCollection::const_iterator
115 vertex = vertices->begin(); vertex!= vertices->end(); vertex++)
119 crossing.
z() - vertex->position().z());
125 if(im < impact || vertex == vertices->
begin())
144 if(impact < maxImpactMother)
148 (momenta.first.cross(momenta.second)).
mag()/momentum.
mag();
150 (momenta.first.mag2() - momenta.second.mag2())/momentum.
mag2();
virtual float distance() const
T getParameter(std::string const &) const
math::GlobalPoint crossingPoint
virtual bool calculate(const TrajectoryStateOnSurface &sta, const TrajectoryStateOnSurface &stb)
const Vector & momentum() const
track momentum vector
virtual GlobalPoint crossingPoint() const
GlobalTrajectoryParameters getGlobalTrajectoryParameters(const reco::Track &track)
Sin< T >::type sin(const T &t)
std::vector< Vertex > VertexCollection
collection of Vertex objects
std::pair< GlobalVector, GlobalVector > momenta
static int position[TOTALCHAMBERS][3]
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
GlobalVector momentum() const
Cos< T >::type cos(const T &t)
double firstAngle() const
double secondAngle() const
const Point & vertex() const
reference point on the track. This method is DEPRECATED, please use referencePoint() instead ...
GlobalVector rotate(const GlobalVector &p, double a)
VZeroFinder(const edm::EventSetup &es, const edm::ParameterSet &pset)
int charge() const
track electric charge
bool checkTrackPair(const reco::Track &posTrack, const reco::Track &negTrack, const reco::VertexCollection *vertices, reco::VZeroData &data)
virtual std::pair< GlobalPoint, GlobalPoint > points() const
Global3DVector GlobalVector