RKPropagatorInR Class Reference

#include <RKPropagatorInR.h>

Inheritance diagram for RKPropagatorInR:

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
	RKPropagatorInR (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 RKPropagatorInR.h.

Constructor & Destructor Documentation

RKPropagatorInR()

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

inline

Definition at line 9 of file RKPropagatorInR.h.

Referenced by clone().

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

Member Function Documentation

clone()

Propagator * RKPropagatorInR::clone ( void  ) const

overridevirtual

Implements Propagator.

Definition at line 66 of file RKPropagatorInR.cc.

References RKPropagatorInR().

{ return new RKPropagatorInR(*this); }

magneticField()

const MagneticField* RKPropagatorInR::magneticField ( ) const

inlineoverridevirtual

Implements Propagator.

Definition at line 23 of file RKPropagatorInR.h.

References theVolume.

{ return theVolume; }

myPropagate() [1/2]

TrajectoryStateOnSurface RKPropagatorInR::myPropagate	(	const FreeTrajectoryState &	,
		const Plane &
	)		const

Definition at line 52 of file RKPropagatorInR.cc.

                                                                                                    {
  return TrajectoryStateOnSurface();
}

myPropagate() [2/2]

TrajectoryStateOnSurface RKPropagatorInR::myPropagate	(	const FreeTrajectoryState &	ts,
		const Cylinder &	cyl
	)		const

Definition at line 14 of file RKPropagatorInR.cc.

References FreeTrajectoryState::charge(), MillePedeFileConverter_cfg::e, HLT_2022v14_cff::eps, FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), Cylinder::radius(), command_line::start, theVolume, Surface::toGlobal(), and GloballyPositioned< T >::toLocal().

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

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

  LocalPoint startpos = cyl.toLocal(pos);
  LocalVector startmom = cyl.toLocal(mom);

  CylindricalState startState(startpos, startmom, ts.charge());
  const RKVector& start = startState.parameters();

  RKLocalFieldProvider localField(*theVolume, cyl);
  CylindricalLorentzForce<double, 5> deriv(localField);
  RKCylindricalDistance<double, 5> dist;
  double eps = 1.e-5;
  Solver solver;
  try {
    Scalar step = cyl.radius() - startState.rho();
    RKVector rkresult = solver(startState.rho(), start, step, deriv, dist, eps);
    CylindricalState endState(cyl.radius(), rkresult, startState.prSign());
    return TrajectoryStateOnSurface(GlobalTrajectoryParameters(cyl.toGlobal(endState.position()),
                                                               cyl.toGlobal(endState.momentum()),
                                                               TrackCharge(endState.charge()),
                                                               theVolume),
                                    cyl);
  } catch (CylindricalLorentzForceException& 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();
  }
}

propagateWithPath() [1/2]

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

overridevirtual

Implements Propagator.

Definition at line 56 of file RKPropagatorInR.cc.

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

propagateWithPath() [2/2]

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

overridevirtual

Implements Propagator.

Definition at line 61 of file RKPropagatorInR.cc.

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

Member Data Documentation

theVolume

const MagVolume* RKPropagatorInR::theVolume

private

Definition at line 26 of file RKPropagatorInR.h.

Referenced by magneticField(), and myPropagate().