GsfExtrapolator.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
 */

 
#ifndef TrkGsfExtrapolator_H
#define TrkGsfExtrapolator_H
 
#include "TrkGaussianSumFilter/IMultiStateExtrapolator.h"
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "Gaudi/Accumulators.h"
#include "GaudiKernel/ToolHandle.h"
 
#include "TrkExInterfaces/INavigator.h"
#include "TrkExInterfaces/IPropagator.h"
#include "TrkGaussianSumFilter/IMaterialMixtureConvolution.h"
 
#include "GeoPrimitives/GeoPrimitives.h"
#include "TrkGeometry/MagneticFieldProperties.h"
#include "TrkMaterialOnTrack/MaterialEffectsOnTrack.h"
#include "TrkParameters/TrackParameters.h"
 
#include <memory>
#include <string>
#include <vector>
 
namespace Trk {
 
class Layer;
class Surface;
class TrackingVolume;
class TrackingGeometry;
class TrackStateOnSurface;
class MaterialProperties;
class GsfExtrapolator final
  : public AthAlgTool
  , virtual public IMultiStateExtrapolator
{
 
public:
  GsfExtrapolator(const std::string&, const std::string&, const IInterface*);
  virtual ~GsfExtrapolator() override final;
  virtual StatusCode initialize() override final;

  virtual MultiComponentState extrapolate(
    const EventContext& ctx,
    Cache&,
    const MultiComponentState&,
    const Surface&,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const override final;

  virtual MultiComponentState extrapolateDirectly(
    const EventContext& ctx,
    const MultiComponentState&,
    const Surface&,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const override final;

  virtual Trk::ParticleHypothesis particleHypothesis() const override final{
    return m_materialUpdator->particleHypothesis();
  }
 
 
private:
  MultiComponentState extrapolateImpl(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState&,
    const Surface&,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const;

  MultiComponentState extrapolateDirectlyImpl(
    const EventContext& ctx,
    const MultiComponentState&,
    const Surface&,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const;

  void extrapolateToVolumeBoundary(const EventContext& ctx,
                                   Cache& cache,
                                   const MultiComponentState&,
                                   const Layer*,
                                   const TrackingVolume&,
                                   PropDirection direction) const;
 
  MultiComponentState extrapolateInsideVolume(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState&,
    const Surface&,
    const Layer*,
    const TrackingVolume&,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const;

  MultiComponentState extrapolateFromLayerToLayer(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState&,
    const TrackingVolume&,
    const Layer* startLayer,
    const Layer* destinationLayer,
    PropDirection direction) const;

  MultiComponentState extrapolateToIntermediateLayer(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState&,
    const Layer&,
    const TrackingVolume&,
    PropDirection direction) const;

  MultiComponentState extrapolateToDestinationLayer(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState&,
    const Surface&,
    const Layer&,
    const Layer*,
    PropDirection direction,
    const BoundaryCheck& boundaryCheck) const;

  std::unique_ptr<Trk::TrackParameters> initialiseNavigation(
    const EventContext& ctx,
    Cache& cache,
    const MultiComponentState& initialState,
    const Surface& surface,
    const Layer*& associatedLayer,
    const TrackingVolume*& currentVolume,
    const TrackingVolume*& destinationVolume,
    PropDirection& direction) const;
 
  ToolHandle<IPropagator> m_propagator{this, "Propagator", "", ""};
  ToolHandle<INavigator> m_navigator{this, "Navigator",
                                     "Trk::Navigator/Navigator", ""};
  ToolHandle<IMaterialMixtureConvolution> m_materialUpdator{
      this, "GsfMaterialConvolution", "", "Gsf Material effects"};
  BooleanProperty m_fastField{this, "UseFastField", false};
 
  Trk::MagneticFieldProperties m_fieldProperties = Trk::FullField;
};
 
} // end namespace Trk
 
#endif