trkf::InteractGeneral Class Reference

#include "InteractGeneral.h"

Inheritance diagram for trkf::InteractGeneral:

Public Member Functions
	InteractGeneral (detinfo::DetectorPropertiesData const &detProp, double tcut)
Interactor *	clone () const override
	Clone method. More...
bool	noise (const KTrack &trk, double s, TrackError &noise_matrix) const override
double	getTcut () const

Private Attributes
InteractPlane	fInteract
PropAny	fProp

Detailed Description

Definition at line 24 of file InteractGeneral.h.

Constructor & Destructor Documentation

trkf::InteractGeneral::InteractGeneral ( detinfo::DetectorPropertiesData const &  detProp, double  tcut  )

explicit

Constructor.

Arguments:

tcut - Maximum delta ray energy.

Definition at line 23 of file InteractGeneral.cxx.

Referenced by clone().

    : Interactor(tcut), fInteract(detProp, tcut), fProp(detProp, -1., false)
  {}

Member Function Documentation

Interactor* trkf::InteractGeneral::clone ( ) const

inlineoverridevirtual

Clone method.

Implements trkf::Interactor.

Definition at line 28 of file InteractGeneral.h.

References InteractGeneral(), and noise().

{ return new InteractGeneral(*this); }

double trkf::Interactor::getTcut ( ) const

inlineinherited

Definition at line 32 of file Interactor.h.

References trkf::Interactor::clone(), trkf::Interactor::fTcut, and trkf::Interactor::noise().

Referenced by trkf::InteractPlane::noise().

{ return fTcut; }

bool trkf::InteractGeneral::noise ( const KTrack &  trk, double  s, TrackError &  noise_matrix  ) const

overridevirtual

Calculate noise matrix.

Arguments:

trk - Original track. s - Path distance. noise_matrix - Resultant noise matrix.

Returns: True if success.

Currently calculate noise from multiple scattering only.

Note about multiple scattering calculation:

We make a zero distance propagation to a plane surface (SurfXYZPlane) that is normal to the track direction. Then calculate the noise matrix on that surface and transform back to the original surface.

Implements trkf::Interactor.

Definition at line 46 of file InteractGeneral.cxx.

References fInteract, fProp, trkf::KTrack::getMomentum(), trkf::KTrack::getPosition(), trkf::invert(), trkf::InteractPlane::noise(), trkf::PropAny::short_vec_prop(), and trkf::Propagator::UNKNOWN.

Referenced by clone().

  {
    // Get track position and direction.

    double xyz[3];
    double mom[3];
    trk.getPosition(xyz);
    trk.getMomentum(mom);

    // Generate a SurfXYZPlane with origin at current track position, and
    // normal to current track direction.

    std::shared_ptr<Surface> psurf(
      new SurfXYZPlane(xyz[0], xyz[1], xyz[2], mom[0], mom[1], mom[2]));

    // Propagate track to newly created surface (zero-distance propagation).

    TrackMatrix prop_matrix;
    KTrack temp_trk = trk;
    std::optional<double> result =
      fProp.short_vec_prop(temp_trk, psurf, Propagator::UNKNOWN, false, &prop_matrix);

    // Return failure if propagation did not succeed.

    if (!result) return false;

    // Calculate noise on plane surface.

    TrackError plane_noise(5);
    fInteract.noise(temp_trk, s, plane_noise);

    // Transform noise matrix to original surface using inverse of propagation matrix.

    invert(prop_matrix);
    TrackMatrix temp = prod(plane_noise, trans(prop_matrix));
    TrackMatrix temp2 = prod(prop_matrix, temp);
    noise_matrix = ublas::symmetric_adaptor<TrackMatrix>(temp2);

    // Done (success).

    return true;
  }

Member Data Documentation

InteractPlane trkf::InteractGeneral::fInteract

private

Definition at line 32 of file InteractGeneral.h.

Referenced by noise().

PropAny trkf::InteractGeneral::fProp

private

Definition at line 33 of file InteractGeneral.h.

Referenced by noise().

---

The documentation for this class was generated from the following files:

• lardata/v09_16_03/source/lardata/RecoObjects/InteractGeneral.h
• lardata/v09_16_03/source/lardata/RecoObjects/InteractGeneral.cxx