PropAny.cxx

Go to the documentation of this file.

#include "lardata/RecoObjects/PropAny.h"
#include "cetlib_except/exception.h"
#include "lardata/RecoObjects/InteractPlane.h"
#include "lardata/RecoObjects/SurfXYZPlane.h"
#include "lardata/RecoObjects/SurfYZLine.h"
#include "lardata/RecoObjects/SurfYZPlane.h"

namespace trkf {

  PropAny::PropAny(detinfo::DetectorPropertiesData const& detProp, double tcut, bool doDedx)
    : Propagator(detProp,
                 tcut,
                 doDedx,
                 (tcut >= 0 ? std::make_shared<InteractPlane const>(detProp, tcut) :
                              std::shared_ptr<Interactor const>{}))
    , fPropYZLine(detProp, tcut, doDedx)
    , fPropYZPlane(detProp, tcut, doDedx)
    , fPropXYZPlane(detProp, tcut, doDedx)
  {}

  std::optional<double> PropAny::short_vec_prop(KTrack& trk,
                                                const std::shared_ptr<const Surface>& psurf,
                                                Propagator::PropDirection dir,
                                                bool doDedx,
                                                TrackMatrix* prop_matrix,
                                                TrackError* noise_matrix) const
  {
    // Default result.

    std::optional<double> result{std::nullopt};

    // Test the type of the destination surface.

    if (dynamic_cast<const SurfYZLine*>(&*psurf))
      result = fPropYZLine.short_vec_prop(trk, psurf, dir, doDedx, prop_matrix, noise_matrix);
    else if (dynamic_cast<const SurfYZPlane*>(&*psurf))
      result = fPropYZPlane.short_vec_prop(trk, psurf, dir, doDedx, prop_matrix, noise_matrix);
    else if (dynamic_cast<const SurfXYZPlane*>(&*psurf))
      result = fPropXYZPlane.short_vec_prop(trk, psurf, dir, doDedx, prop_matrix, noise_matrix);
    else
      throw cet::exception("PropAny") << "Destination surface has unknown type.\n";

    // Done.

    return result;
  }

  std::optional<double> PropAny::origin_vec_prop(KTrack& trk,
                                                 const std::shared_ptr<const Surface>& porient,
                                                 TrackMatrix* prop_matrix) const
  {
    // Default result.

    std::optional<double> result{std::nullopt};

    // Test the type of the destination surface.

    if (dynamic_cast<const SurfYZLine*>(&*porient))
      result = fPropYZLine.origin_vec_prop(trk, porient, prop_matrix);
    else if (dynamic_cast<const SurfYZPlane*>(&*porient))
      result = fPropYZPlane.origin_vec_prop(trk, porient, prop_matrix);
    else if (dynamic_cast<const SurfXYZPlane*>(&*porient))
      result = fPropXYZPlane.origin_vec_prop(trk, porient, prop_matrix);
    else
      throw cet::exception("PropAny") << "Destination surface has unknown type.\n";

    // Done.

    return result;
  }

} // end namespace trkf