#include <TrackParamsEstimationTool.h>

Inheritance diagram for ActsTrk::TrackParamsEstimationTool:

Collaboration diagram for ActsTrk::TrackParamsEstimationTool:

Public Member Functions
	TrackParamsEstimationTool (const std::string &type, const std::string &name, const IInterface *parent)

virtual	~TrackParamsEstimationTool ()=default

virtual StatusCode	initialize () override

virtual std::optional< Acts::BoundTrackParameters >	estimateTrackParameters (const EventContext &ctx, const ActsTrk::Seed &seed, const Acts::GeometryContext &geoContext, const Acts::MagneticFieldContext &magFieldContext, std::function< const Acts::Surface &(const ActsTrk::Seed &)> retrieveSurface) const override

virtual std::optional< Acts::BoundTrackParameters >	estimateTrackParameters (const EventContext &ctx, const ActsTrk::Seed &seed, const Acts::GeometryContext &geoContext, const Acts::Surface &surface, const Acts::Vector3 &bField) const override

Private Member Functions
const Acts::Logger &	logger () const
	Private access to the logger. More...

Private Attributes
Gaudi::Property< double >	m_sigmaLoc0

Gaudi::Property< double >	m_sigmaLoc1

Gaudi::Property< double >	m_sigmaPhi

Gaudi::Property< double >	m_sigmaTheta

Gaudi::Property< double >	m_sigmaQOverP

Gaudi::Property< double >	m_sigmaT0

Gaudi::Property< double >	m_initialSigmaPtRel

Gaudi::Property< std::vector< double > >	m_initialVarInflation

Gaudi::Property< bool >	m_useTopSp

std::unique_ptr< const Acts::Logger >	m_logger {nullptr}
	logging instance More...

Detailed Description

Definition at line 15 of file TrackParamsEstimationTool.h.

Constructor & Destructor Documentation

◆ TrackParamsEstimationTool()

ActsTrk::TrackParamsEstimationTool::TrackParamsEstimationTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 11 of file TrackParamsEstimationTool.cxx.

14 : base_class(type, name, parent)

15 {}

◆ ~TrackParamsEstimationTool()

virtual ActsTrk::TrackParamsEstimationTool::~TrackParamsEstimationTool ( )

virtualdefault

Member Function Documentation

◆ estimateTrackParameters() [1/2]

std::optional< Acts::BoundTrackParameters > ActsTrk::TrackParamsEstimationTool::estimateTrackParameters	(	const EventContext &	ctx,
		const ActsTrk::Seed &	seed,
		const Acts::GeometryContext &	geoContext,
		const Acts::MagneticFieldContext &	magFieldContext,
		std::function< const Acts::Surface &(const ActsTrk::Seed &)>	retrieveSurface
	)		const

overridevirtual

Definition at line 37 of file TrackParamsEstimationTool.cxx.

   {
     const auto& sp_collection = seed.sp();
     if ( sp_collection.size() < 3 ) return std::nullopt;
     const auto& bottom_sp = m_useTopSp ? sp_collection.back() : sp_collection.front();
  
     // Magnetic Field
     ATLASMagneticFieldWrapper magneticField;
     Acts::MagneticFieldProvider::Cache magFieldCache = magneticField.makeCache( magFieldContext );
     Acts::Vector3 bField = *magneticField.getField( Acts::Vector3(bottom_sp->x(), bottom_sp->y(), bottom_sp->z()),
                                                     magFieldCache );
  
     // Get the surface
     const Acts::Surface& surface = retrieveSurface(seed);
  
     return estimateTrackParameters(ctx,
                    seed,
                    geoContext,
                    surface,
                    bField);
   }

◆ estimateTrackParameters() [2/2]

std::optional< Acts::BoundTrackParameters > ActsTrk::TrackParamsEstimationTool::estimateTrackParameters	(	const EventContext &	ctx,
		const ActsTrk::Seed &	seed,
		const Acts::GeometryContext &	geoContext,
		const Acts::Surface &	surface,
		const Acts::Vector3 &	bField
	)		const

overridevirtual

Definition at line 64 of file TrackParamsEstimationTool.cxx.

   {
     // Get SPs
     const auto& sp_collection = seed.sp();
     if ( sp_collection.size() < 3 ) return std::nullopt;
     
     // Compute Bound parameters at surface
     auto params_result = m_useTopSp ?
       Acts::estimateTrackParamsFromSeed(geoContext,
                                         std::ranges::views::reverse(sp_collection),
                                         surface,
                                         bField) :
       Acts::estimateTrackParamsFromSeed(geoContext,
                                         sp_collection,
                                         surface,
                                         bField);
     
     if ( not params_result.ok() ) {
       return std::nullopt;
     }
  
     auto& params = params_result.value();
  
     if (m_useTopSp) {
       // reverse direction so momentum vector pointing outwards
       params = Acts::reflectBoundParameters(params);
     }
  
     Acts::EstimateTrackParamCovarianceConfig covarianceEstimationConfig = {
       .initialSigmas = {m_sigmaLoc0, m_sigmaLoc1, m_sigmaPhi, m_sigmaTheta, m_sigmaQOverP, m_sigmaT0},
       .initialSigmaPtRel = m_initialSigmaPtRel,
       .initialVarInflation = Eigen::Map<const Acts::BoundVector>(m_initialVarInflation.value().data()),
       .noTimeVarInflation = 1.0,
     };
     Acts::BoundMatrix covariance = Acts::estimateTrackParamCovariance(covarianceEstimationConfig, params, false);
  
     // Create BoundTrackParameters
     return Acts::BoundTrackParameters(surface.getSharedPtr(),
                                       params,
                                       covariance,
                                       Acts::ParticleHypothesis::pion());
   }

◆ initialize()

StatusCode ActsTrk::TrackParamsEstimationTool::initialize ( )

overridevirtual

Definition at line 17 of file TrackParamsEstimationTool.cxx.

   {
     ATH_MSG_INFO( "Initializing " << name() << "..." );
  
     ATH_MSG_DEBUG( "Properties Summary:" );
     ATH_MSG_DEBUG( "   " << m_sigmaLoc0 );
     ATH_MSG_DEBUG( "   " << m_sigmaLoc1 );
     ATH_MSG_DEBUG( "   " << m_sigmaPhi );
     ATH_MSG_DEBUG( "   " << m_sigmaTheta );
     ATH_MSG_DEBUG( "   " << m_sigmaQOverP );
     ATH_MSG_DEBUG( "   " << m_sigmaT0 );
     ATH_MSG_DEBUG( "   " << m_initialVarInflation );
     ATH_MSG_DEBUG( "   " << m_useTopSp );
  
     // Make the logger And Propagate to ACTS routines
     m_logger = makeActsAthenaLogger(this, "Acts");
     return StatusCode::SUCCESS;
   }