trkf::InteractGeneral Class Reference

#include "InteractGeneral.h"

Inheritance diagram for trkf::InteractGeneral:

Public Member Functions
	InteractGeneral (double tcut)
	Constructor. More...
virtual	~InteractGeneral ()
	Destructor. More...
Interactor *	clone () const
	Clone method. More...
virtual bool	noise (const KTrack &trk, double s, TrackError &noise_matrix) const
	Calculate noise matrix. More...
double	getTcut () const

Private Attributes
InteractPlane	fInteract
	Plane interactor. More...
PropAny	fProp
	Propagator. More...

Detailed Description

Definition at line 24 of file InteractGeneral.h.

Constructor & Destructor Documentation

trkf::InteractGeneral::InteractGeneral

(

double

tcut

)

Constructor.

Constructor.

Arguments:

tcut - Maximum delta ray energy.

Definition at line 24 of file InteractGeneral.cxx.

Referenced by clone().

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

trkf::InteractGeneral::~InteractGeneral ( )

virtual

Destructor.

Definition at line 31 of file InteractGeneral.cxx.

  {}

Member Function Documentation

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

inlinevirtual

Clone method.

Implements trkf::Interactor.

Definition at line 37 of file InteractGeneral.h.

References InteractGeneral(), noise(), and s.

{return new InteractGeneral(*this);}

double trkf::Interactor::getTcut ( ) const

inlineinherited

Definition at line 39 of file Interactor.h.

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

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

{return fTcut;}

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

virtual

Calculate noise matrix.

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 53 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;
    boost::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

Plane interactor.

Definition at line 47 of file InteractGeneral.h.

Referenced by noise().

PropAny trkf::InteractGeneral::fProp

private

Propagator.

Definition at line 50 of file InteractGeneral.h.

Referenced by noise().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/lardata/v07_01_01/source/lardata/RecoObjects/InteractGeneral.h
• /cvmfs/larsoft.opensciencegrid.org/products/lardata/v07_01_01/source/lardata/RecoObjects/InteractGeneral.cxx