ExtrapolationEngine.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// ExtrapolationEngine.cxx, (c) ATLAS Detector software
 
// STL
#include <sstream>
// Trk include
#include "TrkExEngine/ExtrapolationEngine.h"
#include "TrkDetDescrUtils/GeometrySignature.h"
 
// constructor
Trk::ExtrapolationEngine::ExtrapolationEngine(const std::string& t, const std::string& n, const IInterface* p)
: AthCheckedComponent<AthAlgTool>(t,n,p),   
  m_forceSearchInit(false)
{
    declareInterface<Trk::IExtrapolationEngine>(this);
    // steering of the screen outoput (SOP)
    declareProperty("OutputPrefix"                          , m_sopPrefix);
    declareProperty("OutputPostfix"                         , m_sopPostfix);
    // the properties to be given 
    declareProperty("ForceSearchAtInit"                     , m_forceSearchInit);
}
 
// destructor
Trk::ExtrapolationEngine::~ExtrapolationEngine()
= default;
 
 
// the interface method initialize
StatusCode Trk::ExtrapolationEngine::initialize()
{
            
    EX_MSG_DEBUG( "", "initialize", "", "starting initialize()" );
    // get the TrackingGeometrySvc
    ATH_CHECK( m_trackingGeometryReadKey.initialize(!m_trackingGeometryReadKey.key().empty()) );
 
    // retrieve the ExtrapolationEngines
    if (m_extrapolationEngines.retrieve().isFailure()){
        
    } else {
        EX_MSG_DEBUG( "", "initialize", "", "Successfully retrieved " << m_extrapolationEngines.size() << " ExtrapolationEngines. Ordering them now." );
        m_eeAccessor = std::vector<const Trk::IExtrapolationEngine*>(int(Trk::NumberOfGeometryTypes), (const Trk::IExtrapolationEngine*)nullptr);
        for (auto& ee : m_extrapolationEngines){
            EX_MSG_DEBUG( "", "initialize", "", "Registering " << ee->name() << " - for GeometryType : "  << ee->geometryType() );
            m_eeAccessor[ee->geometryType()] = (&*ee);
        }
    }
 
    if (m_propagationEngine.retrieve().isFailure()){
      EX_MSG_FATAL("", "initialize", "", "failed to retrieve propagation engine '"<< m_propagationEngine << "'. Aborting." );
      return StatusCode::FAILURE;
    } else 
      EX_MSG_DEBUG("", "initialize", "", "successfully  propagation engine '" << m_propagationEngine << "'." );
    
    if (m_navigationEngine.retrieve().isFailure()){
        EX_MSG_FATAL("", "initialize", "", "failed to retrieve navigation engine '"<< m_navigationEngine << "'. Aborting." );
        return StatusCode::FAILURE;
    } else 
        EX_MSG_DEBUG("", "initialize", "", "successfully retrieved '" << m_navigationEngine << "'." );
    
    return StatusCode::SUCCESS;
}    
 
// the interface method finalize
StatusCode Trk::ExtrapolationEngine::finalize()
{    
    EX_MSG_DEBUG( "", "finalize", "", "successful" );    
    return StatusCode::SUCCESS;
}
 
 
Trk::ExtrapolationCode Trk::ExtrapolationEngine::extrapolate(ExCellCharged& ecCharged,
                                                        const Surface* sf,
                                                        const BoundaryCheck& bcheck) const
{ return extrapolateT<TrackParameters>(ecCharged,sf,ecCharged.propDirection,bcheck); }
 
 
Trk::ExtrapolationCode Trk::ExtrapolationEngine::extrapolate(ExCellNeutral& ecNeutral,
                                                        const Surface* sf,
                                                        const BoundaryCheck& bcheck) const
{ return extrapolateT<NeutralParameters>(ecNeutral,sf,ecNeutral.propDirection,bcheck); }
 
void Trk::ExtrapolationEngine::throwFailedToGetTrackingGeomtry() const {
   std::stringstream msg;
   msg << "Failed to get conditions data " << m_trackingGeometryReadKey.key() << ".";
   throw std::runtime_error(msg.str());
}