HGCal_helpers::SimpleTrackPropagator Class Reference

Public Member Functions
bool	propagate (const double px, const double py, const double pz, const double x, const double y, const double z, const float charge, Coordinates &coords) const
bool	propagate (const math::XYZTLorentzVectorD &momentum, const math::XYZTLorentzVectorD &position, const float charge, Coordinates &coords) const
void	setPropagationTargetZ (const float &z)
	SimpleTrackPropagator (const MagneticField *f)

Private Member Functions
const RKPropagatorInS &	RKProp () const
	SimpleTrackPropagator ()

Private Attributes
float	absz_target_
CurvilinearTrajectoryError	err_
const MagneticField *	field_
defaultRKPropagator::Product	prod_
Plane::PlanePointer	targetPlaneBackward_
Plane::PlanePointer	targetPlaneForward_

Detailed Description

Definition at line 34 of file HGCalTriggerNtupleGen.cc.

Constructor & Destructor Documentation

HGCal_helpers::SimpleTrackPropagator::SimpleTrackPropagator ( const MagneticField *  f )

inline

Definition at line 36 of file HGCalTriggerNtupleGen.cc.

References gen::C, ALCARECOTkAlJpsiMuMu_cff::charge, triggerObjects_cff::id, position, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, HGCal_helpers::Coordinates::x, HGCal_helpers::Coordinates::y, and HGCal_helpers::Coordinates::z.

                                                  : field_(f), prod_(field_, alongMomentum), absz_target_(0) {
      ROOT::Math::SMatrixIdentity id;
      AlgebraicSymMatrix55 C(id);
      const float uncert = 0.001;
      C *= uncert;
      err_ = CurvilinearTrajectoryError(C);
    }

HGCal_helpers::SimpleTrackPropagator::SimpleTrackPropagator ( )

inlineprivate

Definition at line 60 of file HGCalTriggerNtupleGen.cc.

: field_(nullptr), prod_(field_, alongMomentum), absz_target_(0) {}

Member Function Documentation

bool HGCal_helpers::SimpleTrackPropagator::propagate	(	const double	px,
		const double	py,
		const double	pz,
		const double	x,
		const double	y,
		const double	z,
		const float	charge,
		Coordinates &	coords
	)		const

Definition at line 74 of file HGCalTriggerNtupleGen.cc.

References Plane::build(), HGCal_helpers::Coordinates::Coordinates(), HGCal_helpers::Coordinates::eta, PV3DBase< T, PVType, FrameType >::eta(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), HGCal_helpers::Coordinates::phi, PV3DBase< T, PVType, FrameType >::phi(), HGCal_helpers::Coordinates::x, PV3DBase< T, PVType, FrameType >::x(), HGCal_helpers::Coordinates::y, PV3DBase< T, PVType, FrameType >::y(), HGCal_helpers::Coordinates::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by HGCalTriggerNtupleGen::fill().

                                                                   {
    output = Coordinates();

    typedef TrajectoryStateOnSurface TSOS;
    GlobalPoint startingPosition(x, y, z);
    GlobalVector startingMomentum(px, py, pz);
    Plane::PlanePointer startingPlane = Plane::build(Plane::PositionType(x, y, z), Plane::RotationType());
    TSOS startingStateP(
        GlobalTrajectoryParameters(startingPosition, startingMomentum, charge, field_), err_, *startingPlane);

    TSOS trackStateP;
    if (pz > 0) {
      trackStateP = RKProp().propagate(startingStateP, *targetPlaneForward_);
    } else {
      trackStateP = RKProp().propagate(startingStateP, *targetPlaneBackward_);
    }
    if (trackStateP.isValid()) {
      output.x = trackStateP.globalPosition().x();
      output.y = trackStateP.globalPosition().y();
      output.z = trackStateP.globalPosition().z();
      output.phi = trackStateP.globalPosition().phi();
      output.eta = trackStateP.globalPosition().eta();
      return true;
    }
    return false;
  }

bool HGCal_helpers::SimpleTrackPropagator::propagate	(	const math::XYZTLorentzVectorD &	momentum,
		const math::XYZTLorentzVectorD &	position,
		const float	charge,
		Coordinates &	coords
	)		const

Definition at line 108 of file HGCalTriggerNtupleGen.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, and convertSQLitetoXML_cfg::output.

                                                                   {
    return propagate(
        momentum.px(), momentum.py(), momentum.pz(), position.x(), position.y(), position.z(), charge, output);
  }

const RKPropagatorInS& HGCal_helpers::SimpleTrackPropagator::RKProp ( ) const

inlineprivate

Definition at line 61 of file HGCalTriggerNtupleGen.cc.

{ return prod_.propagator; }

void HGCal_helpers::SimpleTrackPropagator::setPropagationTargetZ

(

const float &

z

)

Definition at line 69 of file HGCalTriggerNtupleGen.cc.

References funct::abs(), and Plane::build().

Referenced by HGCalTriggerNtupleGen::fill().

                                                                  {
    targetPlaneForward_ = Plane::build(Plane::PositionType(0, 0, std::abs(z)), Plane::RotationType());
    targetPlaneBackward_ = Plane::build(Plane::PositionType(0, 0, -std::abs(z)), Plane::RotationType());
    absz_target_ = std::abs(z);
  }

Member Data Documentation

float HGCal_helpers::SimpleTrackPropagator::absz_target_

private

Definition at line 66 of file HGCalTriggerNtupleGen.cc.

CurvilinearTrajectoryError HGCal_helpers::SimpleTrackPropagator::err_

private

Definition at line 64 of file HGCalTriggerNtupleGen.cc.

const MagneticField* HGCal_helpers::SimpleTrackPropagator::field_

private

Definition at line 63 of file HGCalTriggerNtupleGen.cc.

defaultRKPropagator::Product HGCal_helpers::SimpleTrackPropagator::prod_

private

Definition at line 65 of file HGCalTriggerNtupleGen.cc.

Plane::PlanePointer HGCal_helpers::SimpleTrackPropagator::targetPlaneBackward_

private

Definition at line 62 of file HGCalTriggerNtupleGen.cc.

Plane::PlanePointer HGCal_helpers::SimpleTrackPropagator::targetPlaneForward_

private

Definition at line 62 of file HGCalTriggerNtupleGen.cc.