A class estimating the transverse impact parameter d0 and its error of a given track wrt to the reconstructed vertex. More...

#include <TrackToVertexIPEstimator.h>

Inheritance diagram for Trk::TrackToVertexIPEstimator:

Collaboration diagram for Trk::TrackToVertexIPEstimator:

Public Member Functions
virtual StatusCode	initialize () override
	Default Athena interface methods. More...

	TrackToVertexIPEstimator (const std::string &t, const std::string &n, const IInterface *p)
	Default Athena interface constructor and destructor. More...

	~TrackToVertexIPEstimator ()

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const xAOD::TrackParticle track, const xAOD::Vertex vtx, bool doRemoval) const override
	Estimate methods returning a d0 and its calculated sigma. More...

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const TrackParameters track, const xAOD::Vertex vtx, bool doRemoval) const override

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const xAOD::TrackParticle track, const xAOD::TrackParticle newtrack, const xAOD::Vertex *vtx, bool doRemoval) const override

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const TrackParameters track, const TrackParameters newtrack, const xAOD::Vertex *vtx, bool doRemoval) const override

virtual xAOD::Vertex *	getUnbiasedVertex (const xAOD::TrackParticle track, const xAOD::Vertex vtx) const override

virtual xAOD::Vertex *	getUnbiasedVertex (const TrackParameters track, const xAOD::Vertex vtx) const override

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const xAOD::TrackParticle track, const xAOD::Vertex vtx) const override

virtual std::unique_ptr< ImpactParametersAndSigma >	estimate (const TrackParameters track, const xAOD::Vertex vtx) const override

virtual double	get3DLifetimeSignOfTrack (const TrackParameters &track, const Amg::Vector3D &jetDirection, const xAOD::Vertex &primaryVertex) const override

virtual double	get3DLifetimeSignOfTrack (const TrackParameters &track, const CLHEP::Hep3Vector &jetDirection, const xAOD::Vertex &primaryVertex) const override

virtual double	get2DLifetimeSignOfTrack (const TrackParameters &track, const Amg::Vector3D &jetDirection, const xAOD::Vertex &primaryVertex) const override

virtual double	get2DLifetimeSignOfTrack (const TrackParameters &track, const CLHEP::Hep3Vector &jetDirection, const xAOD::Vertex &primaryVertex) const override

virtual double	getZLifetimeSignOfTrack (const TrackParameters &track, const Amg::Vector3D &jetDirection, const xAOD::Vertex &primaryVertex) const override

virtual double	getZLifetimeSignOfTrack (const TrackParameters &track, const CLHEP::Hep3Vector &jetDirection, const xAOD::Vertex &primaryVertex) const override

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
std::unique_ptr< ImpactParametersAndSigma >	calculate (const TrackParameters *track, const xAOD::Vertex &vtx) const
	A method calculating the do and its error. More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
ToolHandle< Trk::IExtrapolator >	m_extrapolator

ToolHandle< Trk::IVertexUpdator >	m_Updator

ToolHandle< Trk::IVertexLinearizedTrackFactory >	m_linFactory

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

A class estimating the transverse impact parameter d0 and its error of a given track wrt to the reconstructed vertex.

Algorithm:

(a) Provided track or track particle are checked, whether they do already belong to the vertex. (b) In the positive case, the vertex is reduced by removing the track or track particle in question. (c) The trajectory state is extrapolated to the perigee with the vertex in question, the d0 and its error is estimated. The error on the vertex position is included in this calculation.

Kiril.nosp@m.l.Pr.nosp@m.okofi.nosp@m.ev@c.nosp@m.ern.c.nosp@m.h February 2008

Changes:

Giaci.nosp@m.nto..nosp@m.Piacq.nosp@m.uadi.nosp@m.o@phy.nosp@m.sik..nosp@m.uni-f.nosp@m.reib.nosp@m.urg.d.nosp@m.e (update of 11.September 2008)

In addition the tool is now able to provide the impact parameter z0 and its error.

At the same time an option can be passed not to remove the track from the vertex (non-default)

In addition, the lifetime sign of the Impact Parameter given a lifetime direction of flight is also provided with 3 methods: 2D, Z and full 3D.

Johanna Fleckner (03.12.2009) allow for estimation of IP parameter even if vertex doesn't contain any tracks as can happen in BeamSpot or early running when PV is not reconstructed degrade some INFO/WARNINGS to DEBUG

David Shope david.nosp@m..ric.nosp@m.hard..nosp@m.shop.nosp@m.e@cer.nosp@m.n.ch (2016-04-19) EDM Migration to xAOD - Remove methods involving Trk::VxCandidate and Trk::RecVertex

Definition at line 71 of file TrackToVertexIPEstimator.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ TrackToVertexIPEstimator()

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

Default Athena interface constructor and destructor.

Definition at line 19 of file TrackToVertexIPEstimator.cxx.

                                                                                                                :
   AthAlgTool(t,n,p)
  {
   declareInterface<ITrackToVertexIPEstimator>(this);
  }

◆ ~TrackToVertexIPEstimator()

Trk::TrackToVertexIPEstimator::~TrackToVertexIPEstimator ( )

default

Member Function Documentation

◆ calculate()

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::calculate	(	const TrackParameters *	track,
		const xAOD::Vertex &	vtx
	)		const

private

A method calculating the do and its error.

At this stage the trajectory is removed from the vertex.

Definition at line 122 of file TrackToVertexIPEstimator.cxx.

  {
   //estimating the d0 and its significance by propagating the trajectory state towards
   //the vertex position. By this time the vertex should NOT contain this trajectory anymore
  
   //estrapolating to the  perigee of the reconstructed vertex
   const Amg::Vector3D & lp = vtx.position();
   PerigeeSurface perigeeSurface(lp);
   const Trk::Perigee* extrapolatedParameters =
     dynamic_cast<const Trk::Perigee*>(
       m_extrapolator
         ->extrapolate(Gaudi::Hive::currentContext(), *track, perigeeSurface)
         .release());
   if (extrapolatedParameters && extrapolatedParameters->covariance()) {
  
     //actual calculation of d0 and sigma.
     const AmgVector(5) & par = extrapolatedParameters->parameters();
     const double d0  = par[Trk::d0];
     const double z0  = par[Trk::z0];
     const double phi = par[Trk::phi];
     const double theta = par[Trk::theta];
  
     AmgSymMatrix(2) vrtXYCov = vtx.covariancePosition().block<2,2>(0,0);
  
     const AmgSymMatrix(5) & perigeeCov = *(extrapolatedParameters->covariance());
     //   std::cout<<"Perigee covariance: "<<perigeeCov<<std::endl;
  
     //d0phi->cartesian Jacobian
     Amg::Vector2D d0JacXY(-sin(phi), cos(phi));
  
  
     auto newIPandSigma=std::make_unique<ImpactParametersAndSigma>();
     newIPandSigma->IPd0=d0;
     double d0_PVcontrib=d0JacXY.transpose()*(vrtXYCov*d0JacXY);
     if (d0_PVcontrib>=0)
     {
       newIPandSigma->sigmad0=sqrt(d0_PVcontrib+perigeeCov(Trk::d0,Trk::d0));
       newIPandSigma->PVsigmad0=sqrt(d0_PVcontrib);
     }
     else
     {
       msg(MSG::WARNING) << " The contribution to d0_err: " << d0_PVcontrib << " from PV is negative: critical error in PV error matrix! Removing contribution from PV ... "  << endmsg;
       newIPandSigma->sigmad0=sqrt(perigeeCov(Trk::d0,Trk::d0));
       newIPandSigma->PVsigmad0=0;
     }
  
     AmgSymMatrix(2) covPerigeeZ0Theta;
    covPerigeeZ0Theta(0,0)=perigeeCov(Trk::z0,Trk::z0);
    covPerigeeZ0Theta(0,1)=perigeeCov(Trk::z0,Trk::theta);
    covPerigeeZ0Theta(1,0)=perigeeCov(Trk::theta,Trk::z0);
    covPerigeeZ0Theta(1,1)=perigeeCov(Trk::theta,Trk::theta);
  
    double vrtZZCov = vtx.covariancePosition()(Trk::z,Trk::z);
  
    Amg::Vector2D IPz0JacZ0Theta (sin(theta), z0*cos(theta));
    if (vrtZZCov>=0)
    {
      newIPandSigma->IPz0SinTheta=z0*sin(theta);
      newIPandSigma->sigmaz0SinTheta=
          sqrt(IPz0JacZ0Theta.transpose()*(covPerigeeZ0Theta*IPz0JacZ0Theta)+sin(theta)*vrtZZCov*sin(theta));
      newIPandSigma->PVsigmaz0SinTheta=sqrt(sin(theta)*vrtZZCov*sin(theta));
      newIPandSigma->IPz0 = z0;
      newIPandSigma->sigmaz0 = std::sqrt( vrtZZCov + perigeeCov(Trk::z0,Trk::z0) );
      newIPandSigma->PVsigmaz0 = std::sqrt( vrtZZCov );
    }
    else
    {
      ATH_MSG_WARNING(" The contribution to z0_err: " << vrtZZCov << " from PV is negative: critical error in PV error matrix! Removing contribution from PV ... ");
      newIPandSigma->IPz0SinTheta=z0*sin(theta);
      double temp = (IPz0JacZ0Theta.transpose()*(covPerigeeZ0Theta*IPz0JacZ0Theta));
      newIPandSigma->sigmaz0SinTheta=sqrt(temp);
      newIPandSigma->PVsigmaz0SinTheta=0;
  
      newIPandSigma->IPz0 = z0;
      newIPandSigma->sigmaz0 = std::sqrt( perigeeCov(Trk::z0,Trk::z0) );
      newIPandSigma->PVsigmaz0 = 0;
    }
  
  
    delete extrapolatedParameters;
    return newIPandSigma;
   }
     ATH_MSG_DEBUG ("Cannot extrapolate the trajectory state. Returning null. ");
     return nullptr;
   //end of successfull extrapolation check
  }//end of actual calculation method

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ estimate() [1/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const TrackParameters *	track,
		const TrackParameters *	newtrack,
		const xAOD::Vertex *	vtx,
		bool	doRemoval
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 90 of file TrackToVertexIPEstimator.cxx.

  {
  
    if (vtx==nullptr)
    {
      ATH_MSG_WARNING("Vertex is zero pointer. Will not estimate IP of track.");
      return nullptr;
    }
  
    const xAOD::Vertex *newVertex = vtx;
    if (doRemoval)
    {
      newVertex = getUnbiasedVertex(track,vtx);
      if (newVertex == nullptr) {
        ATH_MSG_WARNING("Unbiasing of vertex failed. Will not estimate IP of track.");
        return nullptr;
      }
    }
  
    std::unique_ptr<ImpactParametersAndSigma>  IPandSigma=calculate(newtrack,*newVertex);
  
    if (doRemoval)
    {
      delete newVertex;
      newVertex=nullptr;
    }
  
    return IPandSigma;
  
  }//end of parameterBase estimate method

◆ estimate() [2/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const TrackParameters *	track,
		const xAOD::Vertex *	vtx
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 441 of file TrackToVertexIPEstimator.cxx.

  {
  
    if(track && vtx ){
      return calculate( track , *vtx);
    }
    return nullptr;
  
  }

◆ estimate() [3/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const TrackParameters *	track,
		const xAOD::Vertex *	vtx,
		bool	doRemoval
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 79 of file TrackToVertexIPEstimator.cxx.

  {
    if(track && vtx){
      return estimate(track,track,vtx,doRemoval);
    }
    ATH_MSG_INFO( "Empty TrackParticle or Vertex pointer passed. Returning zero " );
    return nullptr;
    //end of track particle validity check
  
  }//end of parameterBase estimate method

◆ estimate() [4/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const xAOD::TrackParticle *	track,
		const xAOD::TrackParticle *	newtrack,
		const xAOD::Vertex *	vtx,
		bool	doRemoval
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 65 of file TrackToVertexIPEstimator.cxx.

  {
   if(track && vtx)
   {
     return estimate(&(track->perigeeParameters()),&(newtrack->perigeeParameters()),vtx,doRemoval);
   }
    ATH_MSG_INFO( "Empty TrackParticle or Vertex pointer passed. Returning zero " );
    return nullptr;
   //end of track particle validity check
  }//end of method using track particles

◆ estimate() [5/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const xAOD::TrackParticle *	track,
		const xAOD::Vertex *	vtx
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 430 of file TrackToVertexIPEstimator.cxx.

  {
  
    if(track && vtx ){
      return estimate( &(track->perigeeParameters()), vtx);
    }
    return nullptr;
  
  }

◆ estimate() [6/6]

std::unique_ptr< ImpactParametersAndSigma > Trk::TrackToVertexIPEstimator::estimate	(	const xAOD::TrackParticle *	track,
		const xAOD::Vertex *	vtx,
		bool	doRemoval
	)		const

overridevirtual

Estimate methods returning a d0 and its calculated sigma.

All methods do check whether a track is actually fitted to a vertex

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 54 of file TrackToVertexIPEstimator.cxx.

  {
   if(track && vtx)
   {
    return estimate(&(track->perigeeParameters()),&(track->perigeeParameters()),vtx,doRemoval);
   }
   ATH_MSG_INFO( "Empty TrackParticle or Vertex pointer passed. Returning zero " );
   return nullptr;
   //end of track particle validity check
  }//end of method using track particles

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ get2DLifetimeSignOfTrack() [1/2]

double Trk::TrackToVertexIPEstimator::get2DLifetimeSignOfTrack	(	const TrackParameters &	track,
		const Amg::Vector3D &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 254 of file TrackToVertexIPEstimator.cxx.

   {
     const Amg::Vector3D & lp = primaryVertex.position();
     PerigeeSurface perigeeSurface(lp);
  
     std::unique_ptr<const Trk::TrackParameters> extrapolatedParameters =
       m_extrapolator->extrapolate(
         Gaudi::Hive::currentContext(), track, perigeeSurface);
  
     if (!extrapolatedParameters) return 0.;
  
     double trackD0  = extrapolatedParameters->parameters()[Trk::d0];
     double trackPhi = extrapolatedParameters->parameters()[Trk::phi];
     double vs = sinf( atan2(jetMomentum.y(),jetMomentum.x()) - trackPhi )*trackD0;
  
     return (vs>=0. ? 1. : -1.);
   }

◆ get2DLifetimeSignOfTrack() [2/2]

double Trk::TrackToVertexIPEstimator::get2DLifetimeSignOfTrack	(	const TrackParameters &	track,
		const CLHEP::Hep3Vector &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 246 of file TrackToVertexIPEstimator.cxx.

   {
     Amg::Vector3D eigenJetMomentum(jetMomentum.x(), jetMomentum.y(), jetMomentum.z());
     return get2DLifetimeSignOfTrack(track, eigenJetMomentum, primaryVertex);
   }

◆ get3DLifetimeSignOfTrack() [1/2]

double Trk::TrackToVertexIPEstimator::get3DLifetimeSignOfTrack	(	const TrackParameters &	track,
		const Amg::Vector3D &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 222 of file TrackToVertexIPEstimator.cxx.

   {
     const Amg::Vector3D & lp = primaryVertex.position();
     PerigeeSurface perigeeSurface(lp);
  
     std::unique_ptr<const Trk::TrackParameters> extrapolatedParameters =
       m_extrapolator->extrapolate(
         Gaudi::Hive::currentContext(), track, perigeeSurface);
  
     if (!extrapolatedParameters) return 0.;
  
     const Amg::Vector3D & primaryPos=primaryVertex.position();
     const Amg::Vector3D & trackPos=extrapolatedParameters->position();
     const Amg::Vector3D & trackMom=extrapolatedParameters->momentum();
  
     double sign=(jetMomentum.cross(trackMom)).dot(trackMom.cross(primaryPos-trackPos));
  
     return sign>=0.?1.:-1;
   }

◆ get3DLifetimeSignOfTrack() [2/2]

double Trk::TrackToVertexIPEstimator::get3DLifetimeSignOfTrack	(	const TrackParameters &	track,
		const CLHEP::Hep3Vector &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 214 of file TrackToVertexIPEstimator.cxx.

   {
     Amg::Vector3D eigenJetMomentum(jetMomentum.x(), jetMomentum.y(), jetMomentum.z());
     return get3DLifetimeSignOfTrack(track, eigenJetMomentum, primaryVertex);
   }

◆ getUnbiasedVertex() [1/2]

xAOD::Vertex * Trk::TrackToVertexIPEstimator::getUnbiasedVertex	(	const TrackParameters *	track,
		const xAOD::Vertex *	vtx
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 342 of file TrackToVertexIPEstimator.cxx.

  {
    if (!track) {
      msg(MSG::INFO) << "Empty Trk::TrackParameter pointer passed. Returning zero " << endmsg;
      return nullptr;
    }
    if (!vtx) {
      msg(MSG::INFO) << "Empty xAOD::Vertex pointer passed. Returning zero " << endmsg;
      return nullptr;
    }
  
    if (vtx->nTrackParticles() == 0)
     {
       ATH_MSG_DEBUG("This vertex has no associated tracks. Normal if beam spot is used. Vertex already unbiased");
       return new xAOD::Vertex(*vtx);
     }
  
    //create new vertex for output
    xAOD::Vertex *outputVertex = new xAOD::Vertex(*vtx);
    outputVertex->clearTracks(); //remove all tracks -> will add them back one by one
    if (outputVertex->vxTrackAtVertexAvailable()) outputVertex->vxTrackAtVertex().clear(); //remove all VxTrackAtVertex
  
    bool tmpLinTrack = false; //do we created a new linearised track?
  
    //loop over tracks
    const Amg::Vector3D & pos = track->position();
    const Amg::Vector3D & mom = track->momentum();
    for (unsigned int itrk=0; itrk < vtx->nTrackParticles(); ++itrk) {
      const Perigee& testTP = vtx->trackParticle(itrk)->perigeeParameters();
      if ((testTP.position() == pos) and (testTP.momentum() == mom)) {
        //track found, now unbias the vertex using linearized track
        const LinearizedTrack *linTrack = nullptr;
        double trackWeight(0.0);
        //first check if a VxTrackAtVertex is already available with linearized track
        if (vtx->vxTrackAtVertexAvailable()) {
      if (vtx->vxTrackAtVertex().size() > itrk) {
        linTrack = vtx->vxTrackAtVertex()[itrk].linState();
        trackWeight = vtx->vxTrackAtVertex()[itrk].weight();
      }
        }
        //if linearized track is not available, create it
        if (!linTrack) {
      linTrack = m_linFactory->linearizedTrack(track, vtx->position());
      trackWeight = vtx->trackWeight(itrk);
      if (!linTrack) {
        ATH_MSG_INFO("Failing to linearized track. Returning biased vertex.");
        outputVertex->addTrackAtVertex(vtx->trackParticleLinks()[itrk], vtx->trackWeight(itrk));
        if (vtx->vxTrackAtVertexAvailable()) {
          outputVertex->vxTrackAtVertex().push_back(vtx->vxTrackAtVertex()[itrk]);//will clone everything inside -> output vertex owns all the memory
        }
        continue;
      }
      tmpLinTrack = true;
        }
        //now update vertex position removing the linearized track, and do not add the track back to the output vertex
        const IVertexUpdator::positionUpdateOutcome & reducedVertex = m_Updator->positionUpdate(*vtx, linTrack, trackWeight,IVertexUpdator::removeTrack);
  
        //calculate updated chi2
        double chi2 = vtx->chiSquared();
        double trk_chi = m_Updator->trackParametersChi2( reducedVertex, linTrack );
        chi2 += -1 * (m_Updator->vertexPositionChi2(*vtx, reducedVertex) + trackWeight * trk_chi);
  
        //calculate updated ndf
        double ndf = vtx->numberDoF();
        ndf += -1 * trackWeight * (2.0);
  
        //reducedVertex has the updated position and covariance
        outputVertex->setPosition(reducedVertex.position);
        outputVertex->setCovariancePosition(reducedVertex.covariancePosition);
  
        //chi2 and ndf now store the updated FitQuality
        outputVertex->setFitQuality( chi2, ndf );
  
        if (tmpLinTrack) delete linTrack; //only delete if it was created new, and doesn't belong to a vertex
      } else {
        //track not to be removed. Add back the track to the vertex collection
        outputVertex->addTrackAtVertex(vtx->trackParticleLinks()[itrk], vtx->trackWeight(itrk));
        if (vtx->vxTrackAtVertexAvailable()) {
      outputVertex->vxTrackAtVertex().push_back(vtx->vxTrackAtVertex()[itrk]);//will clone everything inside -> output vertex owns all the memory
        }
      }
    }
    return outputVertex;
  
  }

◆ getUnbiasedVertex() [2/2]

xAOD::Vertex * Trk::TrackToVertexIPEstimator::getUnbiasedVertex	(	const xAOD::TrackParticle *	track,
		const xAOD::Vertex *	vtx
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 331 of file TrackToVertexIPEstimator.cxx.

  {
   if(track)
   {
    return getUnbiasedVertex(&(track->perigeeParameters()),vtx);
   }
    msg(MSG::INFO) << "Empty xAOD::TrackParticle pointer passed. Returning zero " << endmsg;
    return nullptr;
   //end of track particle validity check
  }

◆ getZLifetimeSignOfTrack() [1/2]

double Trk::TrackToVertexIPEstimator::getZLifetimeSignOfTrack	(	const TrackParameters &	track,
		const Amg::Vector3D &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 285 of file TrackToVertexIPEstimator.cxx.

   {
  
     const Amg::Vector3D & lp = primaryVertex.position();
     PerigeeSurface perigeeSurface(lp);
  
     std::unique_ptr<const Trk::TrackParameters> extrapolatedParameters =
       m_extrapolator->extrapolate(
         Gaudi::Hive::currentContext(), track, perigeeSurface);
  
     if (!extrapolatedParameters) return 0.;
  
     double trackTheta  = extrapolatedParameters->parameters()[Trk::theta];
     double trackZ0 = extrapolatedParameters->parameters()[Trk::z0];
     double trackEta = -logf(tanf(trackTheta/2.));
     double jetEta = jetMomentum.eta();
     double zs = (jetEta - trackEta)*trackZ0;
     return (zs>=0. ? 1. : -1.);
   }

◆ getZLifetimeSignOfTrack() [2/2]

double Trk::TrackToVertexIPEstimator::getZLifetimeSignOfTrack	(	const TrackParameters &	track,
		const CLHEP::Hep3Vector &	jetDirection,
		const xAOD::Vertex &	primaryVertex
	)		const

overridevirtual

Implements Trk::ITrackToVertexIPEstimator.

Definition at line 277 of file TrackToVertexIPEstimator.cxx.

   {
     Amg::Vector3D eigenJetMomentum(jetMomentum.x(), jetMomentum.y(), jetMomentum.z());
     return getZLifetimeSignOfTrack(track, eigenJetMomentum, primaryVertex);
   }

◆ initialize()

StatusCode Trk::TrackToVertexIPEstimator::initialize ( )

overridevirtual

Default Athena interface methods.

Definition at line 27 of file TrackToVertexIPEstimator.cxx.

  {
  
   // uploading the corresponding tools
   // extrapolator
   if (m_extrapolator.retrieve().isFailure()) {
     ATH_MSG_FATAL("Failed to retrieve tool " << m_extrapolator);
     return StatusCode::FAILURE;
   }
  
   // updator
   if (m_Updator.retrieve().isFailure()) {
     ATH_MSG_FATAL("Failed to retrieve tool " << m_Updator);
     return StatusCode::FAILURE;
   }
  
   // linearized track factory
   if ( m_linFactory.retrieve().isFailure() )
   {
     ATH_MSG_FATAL("Failed to retrieve tool " << m_linFactory );
     return StatusCode::FAILURE;
   }
  
   return StatusCode::SUCCESS;
  }//end of initialize method

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& Trk::ITrackToVertexIPEstimator::interfaceID ( )

inlinestaticinherited

Definition at line 53 of file ITrackToVertexIPEstimator.h.

                                           {
     return IID_ITrackToVertexCompatibilityEstimator;
   };

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extrapolator

ToolHandle<Trk::IExtrapolator> Trk::TrackToVertexIPEstimator::m_extrapolator

private

Initial value:

{this, "Extrapolator",

"Trk::Extrapolator"}

Definition at line 148 of file TrackToVertexIPEstimator.h.

◆ m_linFactory

ToolHandle<Trk::IVertexLinearizedTrackFactory> Trk::TrackToVertexIPEstimator::m_linFactory

private

Initial value:

{

this, "LinearizedTrackFactory", "Trk::FullLinearizedTrackFactory"}

Definition at line 152 of file TrackToVertexIPEstimator.h.

◆ m_Updator

ToolHandle<Trk::IVertexUpdator> Trk::TrackToVertexIPEstimator::m_Updator

private

Initial value:

{this, "VertexUpdator",

"Trk::KalmanVertexUpdator"}

Definition at line 150 of file TrackToVertexIPEstimator.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ TrackToVertexIPEstimator()

◆ ~TrackToVertexIPEstimator()

Member Function Documentation

◆ calculate()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ estimate() [1/6]

◆ estimate() [2/6]

◆ estimate() [3/6]

◆ estimate() [4/6]

◆ estimate() [5/6]

◆ estimate() [6/6]

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ get2DLifetimeSignOfTrack() [1/2]

◆ get2DLifetimeSignOfTrack() [2/2]

◆ get3DLifetimeSignOfTrack() [1/2]

◆ get3DLifetimeSignOfTrack() [2/2]

◆ getUnbiasedVertex() [1/2]

◆ getUnbiasedVertex() [2/2]

◆ getZLifetimeSignOfTrack() [1/2]

◆ getZLifetimeSignOfTrack() [2/2]

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_evtStore

◆ m_extrapolator

◆ m_linFactory

◆ m_Updator

◆ m_varHandleArraysDeclared

◆ m_vhka