26 #include "cetlib_except/exception.h" 38 throw cet::exception(
"TrackProjectedLength") <<
"cannot provide projected length for " 44 double angleToVert = 0.;
46 if (plane.View() == view) {
47 angleToVert = plane.Wire(0).ThetaZ(
false) - 0.5 * ::util::pi<>();
59 std::sqrt(std::pow(pos_cur.x() - pos_prev.x(), 2) + std::pow(pos_cur.y() - pos_prev.y(), 2) +
60 std::pow(pos_cur.z() - pos_prev.z(), 2));
66 std::abs(std::sin(angleToVert) * dir_prev.Y() + std::cos(angleToVert) * dir_prev.Z());
69 length += dist / cosgamma;
78 size_t trajectory_point )
90 <<
"ERROR: Asking for trajectory point #" << trajectory_point
97 auto const& geom = *(lar::providerFrom<geo::Geometry>());
111 throw cet::exception(
"Track") <<
"track at point #" << trajectory_point
112 <<
" is almost parallel to the wires in view " 115 << point.
direction() <<
", its projection on plane " << plane.
ID()
116 <<
" is " <<
proj <<
").\n";
A point in the trajectory, with position and momentum.
Utilities related to art service access.
static std::string ViewName(View_t view)
Returns the name of the specified view.
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
Provides simple real number checks.
constexpr bool zero(Value_t value) const
Returns whether the value is no farther from 0 than the threshold.
Point_t const & LocationAtPoint(size_t i) const
Access to track position at different points.
double TrackProjectedLength(recob::Track const &track, geo::View_t view)
Returns the length of the projection of a track on a view.
size_t NumberTrajectoryPoints() const
Various functions related to the presence and the number of (valid) points.
constexpr auto abs(T v)
Returns the absolute value of the argument.
cout<< "Opened file "<< fin<< " ixs= "<< ixs<< endl;if(ixs==0) hhh=(TH1F *) fff-> Get("h1")
Class for approximate comparisons.
Access the description of the physical detector geometry.
Definitions of geometry vector data types.
Vector_t const & GetWireDirection() const
Returns the direction of the wires.
Geometry information for a single wire plane.The plane is represented in the geometry by a solid whic...
TrajectoryPoint_t TrajectoryPoint(size_t i) const
Access to i-th TrajectoryPoint or its Flags.
Vector_t direction() const
Returns the direction of the trajectory (unit vector of the momentum).
Encapsulate the construction of a single detector plane .
double dist(const TReal *x, const TReal *y, const unsigned int dimension)
Provides recob::Track data product.
double TrackPitchInView(recob::Track const &track, geo::View_t view, size_t trajectory_point=0U)
Returns the projected length of track on a wire pitch step [cm].
Vector_t DirectionAtPoint(size_t i) const
Access to track direction at different points.
WireCoordProjection_t Projection(Point_t const &point) const
Returns the projection of the specified point on the plane.
Collection of Physical constants used in LArSoft.
PlaneID const & ID() const
Returns the identifier of this plane.
Utility functions to extract information from recob::Track
art framework interface to geometry description
Track from a non-cascading particle.A recob::Track consists of a recob::TrackTrajectory, plus additional members relevant for a "fitted" track:
cet::coded_exception< error, detail::translate > exception
double WirePitch() const
Return the wire pitch (in centimeters). It is assumed constant.
Encapsulate the construction of a single detector plane .
The data type to uniquely identify a cryostat.
Point_t position
position in the trajectory [cm].