RKPropagatorInZ Class Reference

#include <RKPropagatorInZ.h>

Inheritance diagram for RKPropagatorInZ:

Public Member Functions
Propagator *	clone () const override
const MagneticField *	magneticField () const override
TrajectoryStateOnSurface	myPropagate (const FreeTrajectoryState &, const Plane &) const
TrajectoryStateOnSurface	myPropagate (const FreeTrajectoryState &, const Cylinder &) const
std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Plane &) const override
std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Cylinder &) const override
	RKPropagatorInZ (const MagVolume &vol, PropagationDirection dir=alongMomentum)
Public Member Functions inherited from Propagator
template<typename STA , typename SUR >
TrajectoryStateOnSurface	propagate (STA const &state, SUR const &surface) const
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest) const final
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest1, const GlobalPoint &pDest2) const final
virtual FreeTrajectoryState	propagate (const FreeTrajectoryState &ftsStart, const reco::BeamSpot &beamSpot) const final
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const FreeTrajectoryState &, const Surface &) const final
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Surface &sur) const final
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Plane &sur) const
virtual std::pair < TrajectoryStateOnSurface, double >	propagateWithPath (const TrajectoryStateOnSurface &tsos, const Cylinder &sur) const
virtual std::pair < FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest) const
virtual std::pair < FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const GlobalPoint &pDest1, const GlobalPoint &pDest2) const
	Propagate to PCA to a line (given by 2 points) given a starting point. More...
virtual std::pair < FreeTrajectoryState, double >	propagateWithPath (const FreeTrajectoryState &ftsStart, const reco::BeamSpot &beamSpot) const
	Propagate to PCA to a line (given by beamSpot position and slope) given a starting point. More...
virtual PropagationDirection	propagationDirection () const final
	Propagator (PropagationDirection dir=alongMomentum)
virtual bool	setMaxDirectionChange (float phiMax)
virtual void	setPropagationDirection (PropagationDirection dir)
virtual	~Propagator ()

Private Attributes
const MagVolume *	theVolume

Detailed Description

Definition at line 7 of file RKPropagatorInZ.h.

Constructor & Destructor Documentation

RKPropagatorInZ::RKPropagatorInZ ( const MagVolume &  vol, PropagationDirection  dir = alongMomentum  )

inline

Definition at line 9 of file RKPropagatorInZ.h.

Referenced by clone().

: Propagator(dir), theVolume(&vol) {}

Member Function Documentation

Propagator * RKPropagatorInZ::clone ( void  ) const

overridevirtual

Implements Propagator.

Definition at line 73 of file RKPropagatorInZ.cc.

References RKPropagatorInZ().

{ return new RKPropagatorInZ(*this); }

const MagneticField* RKPropagatorInZ::magneticField ( ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 23 of file RKPropagatorInZ.h.

References theVolume.

{ return theVolume; }

TrajectoryStateOnSurface RKPropagatorInZ::myPropagate	(	const FreeTrajectoryState &	ts,
		const Plane &	plane
	)		const

Definition at line 12 of file RKPropagatorInZ.cc.

References FreeTrajectoryState::charge(), alignCSCRings::e, gpuVertexFinder::eps, PV3DBase< T, PVType, FrameType >::mag(), FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), command_line::start, theVolume, GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                             {
  //typedef RK4PreciseSolver<double,5>           Solver;
  //typedef RKAdaptiveSolver<double,RKOne4OrderStep, 5>   Solver;
  typedef RKAdaptiveSolver<double, RKOneCashKarpStep, 5> Solver;
  typedef Solver::Vector RKVector;

  GlobalPoint pos(ts.position());
  GlobalVector mom(ts.momentum());

  // cout << "RKPropagatorInZ: starting from FTS " << ts << endl;

  LocalPoint startpos = plane.toLocal(pos);
  LocalVector startmom = plane.toLocal(mom);
  double pzSign = startmom.z() > 0 ? 1.0 : -1.0;

  // cout << "In local plane coordinates: " << startpos << ", momentum " << startmom << endl;

  RKVector start;
  start(0) = startpos.x();
  start(1) = startpos.y();
  start(2) = startmom.x() / startmom.z();
  start(3) = startmom.y() / startmom.z();
  start(4) = pzSign * ts.charge() / startmom.mag();

  // cout << "RKPropagatorInZ: Solving with par " <<  startpos.z() << " and state " << start << endl;

  RKLocalFieldProvider localField(*theVolume, plane);

  CurvilinearLorentzForce<double, 5> deriv(localField);
  RKCurvilinearDistance<double, 5> dist;
  double eps = 1.e-5;
  Solver solver;
  try {
    RKVector rkresult = solver(startpos.z(), start, -startpos.z(), deriv, dist, eps);

    return TrajectoryStateOnSurface(
        LocalTrajectoryParameters(rkresult(4), rkresult(2), rkresult(3), rkresult(0), rkresult(1), pzSign),
        plane,
        theVolume);
  } catch (CurvilinearLorentzForceException& e) {
    // the propagation failed due to momentum almost parallel to the plane.
    // This does not mean the propagation is impossible, but it should be done
    // in a different parametrization (e.g. s)
    return TrajectoryStateOnSurface();
  }
}

TrajectoryStateOnSurface RKPropagatorInZ::myPropagate	(	const FreeTrajectoryState &	,
		const Cylinder &
	)		const

Definition at line 59 of file RKPropagatorInZ.cc.

                                                                                                       {
  return TrajectoryStateOnSurface();
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInZ::propagateWithPath ( const FreeTrajectoryState &  , const Plane &    ) const

overridevirtual

Implements Propagator.

Definition at line 63 of file RKPropagatorInZ.cc.

                                                                                                   {
  return std::pair<TrajectoryStateOnSurface, double>();
}

std::pair< TrajectoryStateOnSurface, double > RKPropagatorInZ::propagateWithPath ( const FreeTrajectoryState &  , const Cylinder &    ) const

overridevirtual

Implements Propagator.

Definition at line 68 of file RKPropagatorInZ.cc.

                                                                                                      {
  return std::pair<TrajectoryStateOnSurface, double>();
}

Member Data Documentation

const MagVolume* RKPropagatorInZ::theVolume

private

Definition at line 26 of file RKPropagatorInZ.h.

Referenced by magneticField(), and myPropagate().