Reco::TrackToVertex Class Reference

#include <TrackToVertex.h>

Inheritance diagram for Reco::TrackToVertex:

Collaboration diagram for Reco::TrackToVertex:

Public Member Functions
	TrackToVertex (const std::string &, const std::string &, const IInterface *)
	AlgTool like constructor. More...
virtual	~TrackToVertex ()=default
	Virtual destructor. More...
virtual StatusCode	initialize () override final
	AlgTool initailize method. More...
virtual std::unique_ptr< Trk::StraightLineSurface >	GetBeamLine (const InDet::BeamSpotData *) const override final
	Use this for MT Coding. More...
virtual std::unique_ptr< Trk::Perigee >	perigeeAtVertex (const EventContext &ctx, const xAOD::TrackParticle &tp, const Amg::Vector3D &gp) const override final
	Interface method for use with xAOD::Trackparticle and given vertex position. More...
virtual std::unique_ptr< Trk::Perigee >	perigeeAtVertex (const EventContext &ctx, const xAOD::TrackParticle &tp) const override final
	Interface method for use with TrackParticle and default primary vertex from TrackParticle - xAOD. More...
virtual std::unique_ptr< Trk::Perigee >	perigeeAtVertex (const EventContext &ctx, const Trk::Track &trk, const Amg::Vector3D &gp) const override final
	Interface method for use with Track and given vertex position - ESD. More...
virtual std::unique_ptr< Trk::Perigee >	perigeeAtBeamline (const EventContext &ctx, const Trk::Track &trk, const InDet::BeamSpotData *) const override final
	Interface method for use with Track and the beamline. More...
virtual std::unique_ptr< Trk::TrackParameters >	trackAtBeamline (const EventContext &ctx, const xAOD::TrackParticle &tp, const InDet::BeamSpotData *) const override final
	Interface method for use with TrackParticle and the beamline from the BeamSpotSvc - xAOD. More...
virtual std::unique_ptr< Trk::TrackParameters >	trackAtBeamline (const EventContext &ctx, const Trk::Track &trk, const Trk::StraightLineSurface *beamline) const override final
	Interface method for use with Track and the beamline from the BeamSpotSvc - ESD. More...
virtual std::unique_ptr< Trk::TrackParameters >	trackAtBeamline (const EventContext &ctx, const Trk::TrackParameters &tpars, const Trk::StraightLineSurface *beamline) const override final
	Interface method for use with Track and the beamline from the BeamSpotSvc - TrackParameters More...

Private Member Functions
bool	startAtOriginalPerigee (const Trk::Track &track) const

Private Attributes
ToolHandle< Trk::IExtrapolator >	m_extrapolator
	ToolHandle for Extrapolator. More...
Gaudi::Property< bool >	m_startTRTSAAtPerigee

Static Private Attributes
static const Amg::Vector3D	s_origin
	static origin More...

Detailed Description

Definition at line 39 of file TrackToVertex.h.

Constructor & Destructor Documentation

TrackToVertex()

TrackToVertex::TrackToVertex	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

AlgTool like constructor.

Definition at line 23 of file TrackToVertex.cxx.

                                                                                            :
  base_class(t,n,p)
{
}

~TrackToVertex()

virtual Reco::TrackToVertex::~TrackToVertex ( )

virtualdefault

Virtual destructor.

Member Function Documentation

GetBeamLine()

std::unique_ptr< Trk::StraightLineSurface > TrackToVertex::GetBeamLine ( const InDet::BeamSpotData *  beamSpotHandle ) const

finaloverridevirtual

Use this for MT Coding.

Definition at line 35 of file TrackToVertex.cxx.

                                                   {
   // get the transform
    Amg::Transform3D beamTransform = Amg::Transform3D(Amg::AngleAxis3D(beamSpotHandle->beamTilt(0),Amg::Vector3D(0.,1.,0.)));
    beamTransform *= Amg::AngleAxis3D(beamSpotHandle->beamTilt(1),Amg::Vector3D(1.,0.,0.));
    beamTransform.pretranslate(beamSpotHandle->beamPos());
    // create the new beam line
    return std::make_unique< Trk::StraightLineSurface >(beamTransform);
}

initialize()

StatusCode TrackToVertex::initialize ( )

finaloverridevirtual

AlgTool initailize method.

Definition at line 29 of file TrackToVertex.cxx.

{
    ATH_CHECK(m_extrapolator.retrieve());
    return StatusCode::SUCCESS;
}

perigeeAtBeamline()

std::unique_ptr< Trk::Perigee > TrackToVertex::perigeeAtBeamline ( const EventContext &  ctx, const Trk::Track &  trk, const InDet::BeamSpotData *  beamspotptr  ) const

finaloverridevirtual

Interface method for use with Track and the beamline.

Definition at line 104 of file TrackToVertex.cxx.

{
 
  Amg::Vector3D beamspot(s_origin);
  float tiltx = 0.0;
  float tilty = 0.0;
  if (beamspotptr) {
    beamspot = Amg::Vector3D(beamspotptr->beamVtx().position());
    tiltx =  beamspotptr->beamTilt(0);
    tilty =  beamspotptr->beamTilt(1);
  }
  Amg::Translation3D amgtranslation(beamspot);
  Amg::Transform3D pAmgTransf =
    amgtranslation * Amg::RotationMatrix3D::Identity();
  pAmgTransf *= Amg::AngleAxis3D(tilty, Amg::Vector3D(0., 1., 0.));
  pAmgTransf *= Amg::AngleAxis3D(tiltx, Amg::Vector3D(1., 0., 0.));
  // preparation
  Trk::PerigeeSurface persf(pAmgTransf);
 
  std::unique_ptr<Trk::Perigee> vertexPerigee;
  if(track.perigeeParameters() && std::abs(track.perigeeParameters()->position().z())>5000.){
    ATH_MSG_DEBUG("Perigee well outside of tracking detector, skipping extrapolation");
  }
  else{
    std::unique_ptr<Trk::TrackParameters> extrapResult =
      !startAtOriginalPerigee(track)
      ? m_extrapolator->extrapolateTrack(ctx,track, persf)
      : m_extrapolator->extrapolate(ctx,*(track.perigeeParameters()), persf);
    if (extrapResult && extrapResult->surfaceType() == Trk::SurfaceType::Perigee) {
      vertexPerigee.reset(static_cast<Trk::Perigee*>(extrapResult.release()));
    }
  }
  if (!vertexPerigee) {
    // workaround.
    // try again using the first track parameter set, since the current extrapolator will
    // use "the closest" track parameter set which is not necessarily the most useful one to
    // start the extrapolation with.
    const DataVector<const Trk::TrackParameters> *track_parameter_list= track.trackParameters();
    if (track_parameter_list) {
      for(const Trk::TrackParameters *trk_params: *track_parameter_list) {
        if (!trk_params || std::abs(trk_params->position().z())>5000.) {
          continue;
        }
        std::unique_ptr<Trk::TrackParameters> extrapResult = m_extrapolator->extrapolate(ctx,*trk_params, persf);
        if (extrapResult &&
            extrapResult->surfaceType() == Trk::SurfaceType::Perigee) {
          vertexPerigee.reset(static_cast<Trk::Perigee*>(extrapResult.release()));
        }
        break;
      }
    }
  }
  if (!vertexPerigee) {
    const Trk::Perigee* trackPerigee = track.perigeeParameters();
    if (trackPerigee && trackPerigee->associatedSurface() == persf) {
      ATH_MSG_DEBUG("Perigee of Track is already expressed to given vertex, a "
                    "copy is returned.");
      vertexPerigee.reset(trackPerigee->clone());
    } else {
      ATH_MSG_DEBUG(
        "Extrapolation to Beamline Perigee failed, NULL pointer is returned.");
    }
  }
  return (vertexPerigee);
}

perigeeAtVertex() [1/3]

std::unique_ptr< Trk::Perigee > TrackToVertex::perigeeAtVertex ( const EventContext &  ctx, const Trk::Track &  trk, const Amg::Vector3D &  gp  ) const

finaloverridevirtual

Interface method for use with Track and given vertex position - ESD.

Definition at line 76 of file TrackToVertex.cxx.

                                 {
 
  Trk::PerigeeSurface persf(gp);
  std::unique_ptr<Trk::Perigee> vertexPerigee;
  std::unique_ptr<Trk::TrackParameters> extrapResult =
     !startAtOriginalPerigee(track)
     ? m_extrapolator->extrapolateTrack(ctx,track, persf)
     : m_extrapolator->extrapolate(ctx, *(track.perigeeParameters()), persf);
  if (extrapResult && extrapResult->surfaceType() == Trk::SurfaceType::Perigee) {
    vertexPerigee.reset( static_cast<Trk::Perigee*>(extrapResult.release()));
  }
  if (!vertexPerigee) {
    const Trk::Perigee* trackPerigee = track.perigeeParameters();
    if (trackPerigee && trackPerigee->associatedSurface() == persf) {
      ATH_MSG_DEBUG("Perigee of Track is already expressed to given vertex, a "
                    "copy is returned.");
      vertexPerigee.reset(trackPerigee->clone());
    } else{
      ATH_MSG_DEBUG(
        "Extrapolation to Perigee failed, NULL pointer is returned.");
    }
  }
  return (vertexPerigee);
}

perigeeAtVertex() [2/3]

std::unique_ptr< Trk::Perigee > TrackToVertex::perigeeAtVertex ( const EventContext &  ctx, const xAOD::TrackParticle &  tp  ) const

finaloverridevirtual

Interface method for use with TrackParticle and default primary vertex from TrackParticle - xAOD.

Definition at line 45 of file TrackToVertex.cxx.

                                                                {
  return perigeeAtVertex(ctx, tp, Amg::Vector3D(tp.vx(), tp.vy(), tp.vz()));
}

perigeeAtVertex() [3/3]

std::unique_ptr< Trk::Perigee > TrackToVertex::perigeeAtVertex ( const EventContext &  ctx, const xAOD::TrackParticle &  tp, const Amg::Vector3D &  gp  ) const

finaloverridevirtual

Interface method for use with xAOD::Trackparticle and given vertex position.

Definition at line 50 of file TrackToVertex.cxx.

                                 {
 
  // preparation
  Trk::PerigeeSurface persf(gp);
  std::unique_ptr<Trk::Perigee> vertexPerigee;
  // retrieve the Perigee from the track particle
  const Trk::Perigee& trackparPerigee = tp.perigeeParameters();
  if (trackparPerigee.associatedSurface() == persf) {
    ATH_MSG_DEBUG("Perigee of TrackParticle is already expressed to given "
                  "vertex, a copy is returned.");
    return std::unique_ptr<Trk::Perigee>(trackparPerigee.clone());
  } else {
    std::unique_ptr<Trk::TrackParameters> extrapResult =
      m_extrapolator->extrapolateDirectly(ctx,trackparPerigee, persf);
    if (extrapResult && extrapResult->surfaceType() == Trk::SurfaceType::Perigee) {
      vertexPerigee.reset(static_cast<Trk::Perigee*>(extrapResult.release()));
    }
  }
  if (!vertexPerigee)
    ATH_MSG_DEBUG(
      "Extrapolation to Perigee failed, a NULL pointer is returned.");
  return vertexPerigee;
}

startAtOriginalPerigee()

bool Reco::TrackToVertex::startAtOriginalPerigee ( const Trk::Track &  track ) const

inlineprivate

Definition at line 104 of file TrackToVertex.h.

                                                                     {
       return m_startTRTSAAtPerigee.value() &&
              track.info().patternRecognition().test(
                  Trk::TrackInfo::TRTStandalone) &&
              track.perigeeParameters();
     }

trackAtBeamline() [1/3]

std::unique_ptr< Trk::TrackParameters > TrackToVertex::trackAtBeamline ( const EventContext &  ctx, const Trk::Track &  trk, const Trk::StraightLineSurface *  beamline  ) const

finaloverridevirtual

Interface method for use with Track and the beamline from the BeamSpotSvc - ESD.

Definition at line 206 of file TrackToVertex.cxx.

                                                  {
 
  return !startAtOriginalPerigee(trk)
     ? m_extrapolator->extrapolateTrack(ctx, trk, *beamline)
     : m_extrapolator->extrapolate(ctx,*(trk.perigeeParameters()), *beamline);
}

trackAtBeamline() [2/3]

std::unique_ptr< Trk::TrackParameters > TrackToVertex::trackAtBeamline ( const EventContext &  ctx, const Trk::TrackParameters &  tpars, const Trk::StraightLineSurface *  beamline  ) const

finaloverridevirtual

Interface method for use with Track and the beamline from the BeamSpotSvc - TrackParameters

Definition at line 215 of file TrackToVertex.cxx.

                                                  {
  return m_extrapolator->extrapolate(ctx, tpars, *beamline);
}

trackAtBeamline() [3/3]

std::unique_ptr< Trk::TrackParameters > TrackToVertex::trackAtBeamline ( const EventContext &  ctx, const xAOD::TrackParticle &  tp, const InDet::BeamSpotData *  beamspotptr  ) const

finaloverridevirtual

Interface method for use with TrackParticle and the beamline from the BeamSpotSvc - xAOD.

Definition at line 173 of file TrackToVertex.cxx.

                                                {
 
  Amg::Vector3D beamspot(s_origin);
  float tiltx = 0.0;
  float tilty = 0.0;
  if (beamspotptr) {
    beamspot = Amg::Vector3D(beamspotptr->beamVtx().position());
    tiltx =  beamspotptr->beamTilt(0);
    tilty =  beamspotptr->beamTilt(1);
  }
  Amg::Transform3D amgTransf;
  Amg::Translation3D amgtranslation(beamspot);
  amgTransf = amgtranslation * Amg::RotationMatrix3D::Identity();
  amgTransf *= Amg::AngleAxis3D(tilty, Amg::Vector3D(0.,1.,0.));
  amgTransf *= Amg::AngleAxis3D(tiltx, Amg::Vector3D(1.,0.,0.));
 // preparation
  Trk::PerigeeSurface persf(amgTransf);
  std::unique_ptr<Trk::TrackParameters> vertexPerigee;
  // retrieve the Perigee from the track particle
  const Trk::Perigee& trackparPerigee = tp.perigeeParameters();
  if ( trackparPerigee.associatedSurface() == persf) {
       ATH_MSG_DEBUG("Perigee of TrackParticle is already expressed to given vertex, a copy is returned.");
       return std::unique_ptr<Trk::TrackParameters>(trackparPerigee.clone());
  } else
      vertexPerigee = m_extrapolator->extrapolateDirectly(ctx,trackparPerigee, persf);
  if (!vertexPerigee){
     ATH_MSG_DEBUG("Extrapolation to Beam Line failed, a NULL pointer is returned.");
  }
  return vertexPerigee;
}

Member Data Documentation

m_extrapolator

ToolHandle<Trk::IExtrapolator> Reco::TrackToVertex::m_extrapolator

private

Initial value:

{
         this, "Extrapolator",
         "Trk::Extrapolator/AtlasExtrapolator"}

ToolHandle for Extrapolator.

Definition at line 111 of file TrackToVertex.h.

m_startTRTSAAtPerigee

Gaudi::Property<bool> Reco::TrackToVertex::m_startTRTSAAtPerigee

private

Initial value:

{
         this, "StartTRTStandaloneTracksAtOriginalPerigee", false,
         "When extrapolating TRT standalone start at their original perigee "
         "which may have a more realistic covariance"}

Definition at line 115 of file TrackToVertex.h.

s_origin

const Amg::Vector3D TrackToVertex::s_origin

staticprivate

static origin

Definition at line 120 of file TrackToVertex.h.

---

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

• TrackToVertex.h
• TrackToVertex.cxx