13 #include "cetlib_except/exception.h" 73 double sinphi = std::sin(
fPhi);
74 double cosphi = std::cos(
fPhi);
77 uvw[0] = xyz[0] -
fX0;
80 uvw[1] = (xyz[1] -
fY0) * cosphi + (xyz[2] -
fZ0) * sinphi;
83 uvw[2] = -(xyz[1] -
fY0) * sinphi + (xyz[2] -
fZ0) * cosphi;
95 double sinphi = std::sin(
fPhi);
96 double cosphi = std::cos(
fPhi);
99 xyz[0] =
fX0 + uvw[0];
102 xyz[1] =
fY0 + uvw[1] * cosphi - uvw[2] * sinphi;
105 xyz[2] =
fZ0 + uvw[1] * sinphi + uvw[2] * cosphi;
120 while(
result(2) <= -TMath::Pi())
121 result(2) += TMath::TwoPi();
122 while(
result(2) > TMath::Pi())
123 result(2) -= TMath::TwoPi();
139 double sinphi = std::sin(phi);
140 double cosphi = std::cos(phi);
143 uvw[0] = -vec(0) * sinphi;
145 uvw[2] = vec(0) * cosphi;
175 double sinphi = std::sin(phi);
176 double cosphi = std::cos(phi);
177 double sh = 1./std::cosh(eta);
178 double th = std::tanh(eta);
182 double pu = p * cosphi * sh;
184 double pw = p * sinphi * sh;
188 double sinfphi = std::sin(
fPhi);
189 double cosfphi = std::cos(
fPhi);
192 mom[1] = pv * cosfphi - pw * sinfphi;
193 mom[2] = pv * sinfphi + pw * cosfphi;
220 double delta_phi = TVector2::Phi_mpi_pi(
fPhi - psurf->
phi());
243 throw cet::exception(
"SurfYZLine") <<
"Attempt to find distance to non-parallel surface.\n";
248 double otheruvw[3] = {0., 0., 0.};
281 double delta_phi = TVector2::Phi_mpi_pi(
fPhi - psurf->
phi());
282 double dx =
fX0 - psurf->
x0();
283 double dy =
fY0 - psurf->
y0();
284 double dz =
fZ0 - psurf->
z0();
297 out <<
"SurfYZLine{ x0=" <<
fX0 <<
", y0=" <<
fY0 <<
", z0=" <<
fZ0 <<
", phi=" <<
fPhi <<
"}";
TrackDirection
Track direction enum.
virtual void toGlobal(const double uvw[3], double xyz[3]) const
Transform local to global coordinates.
virtual bool isEqual(const Surface &surf) const
Test two surfaces for equality, within tolerance.
virtual void toLocal(const double xyz[3], double uvw[3]) const
Transform global to local coordinates.
virtual bool isTrackValid(const TrackVector &vec) const
Surface-specific tests of validity of track parameters.
std::enable_if_t< std::is_arithmetic_v< T >, T > hypot(T x, T y)
virtual void getMomentum(const TrackVector &vec, double mom[3], TrackDirection dir=UNKNOWN) const
Get momentum vector of track.
virtual bool isParallel(const Surface &surf) const
Test whether two surfaces are parallel, within tolerance.
virtual double distanceTo(const Surface &surf) const
Find perpendicular distance to a parallel surface.
virtual ~SurfYZLine()
Destructor.
static double fPhiTolerance
Phi tolerance for parallel.
virtual std::ostream & Print(std::ostream &out) const
Printout.
KVector< 5 >::type TrackVector
Track state vector, dimension 5.
static int max(int a, int b)
double z0() const
Z origin.
static double fSepTolerance
Separation tolerance for equal.
double x0() const
X origin.
virtual TrackVector getDiff(const TrackVector &vec1, const TrackVector &vec2) const
Calculate difference of two track parameter vectors.
double phi() const
Rotation angle about x-axis.
double y0() const
Y origin.
Line surface perpendicular to x-axis.
double fPhi
Rotation angle about x-axis.
virtual Surface * clone() const
Clone method.
virtual void getPosition(const TrackVector &vec, double xyz[3]) const
Get position of track.
virtual void toGlobal(const double uvw[3], double xyz[3]) const =0
Transform local to global coordinates.
cet::coded_exception< error, detail::translate > exception
SurfYZLine()
Default constructor.