d3/d97/GsfExtrapolator_8h_source.html

/*

   Copyright (C) 2002-2022 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,

    ParticleHypothesis particleHypothesis =

      nonInteracting) const override final;


  virtual MultiComponentState extrapolateDirectly(

    const EventContext& ctx,

    const MultiComponentState&,

    const Surface&,

    PropDirection direction,

    const BoundaryCheck& boundaryCheck,

    ParticleHypothesis particleHypothesis =

      nonInteracting) const override final;


private:

  MultiComponentState extrapolateImpl(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState&,

    const Surface&,

    PropDirection direction,

    const BoundaryCheck& boundaryCheck,

    ParticleHypothesis particleHypothesis) const;


  MultiComponentState extrapolateDirectlyImpl(

    const EventContext& ctx,

    const MultiComponentState&,

    const Surface&,

    PropDirection direction = anyDirection,

    const BoundaryCheck& boundaryCheck = true,

    ParticleHypothesis particleHypothesis = nonInteracting) const;


  void extrapolateToVolumeBoundary(const EventContext& ctx,

                                   Cache& cache,

                                   const MultiComponentState&,

                                   const Layer*,

                                   const TrackingVolume&,

                                   PropDirection direction,

                                   ParticleHypothesis particleHypothesis) const;


  MultiComponentState extrapolateInsideVolume(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState&,

    const Surface&,

    const Layer*,

    const TrackingVolume&,

    PropDirection direction,

    const BoundaryCheck& boundaryCheck,

    ParticleHypothesis particleHypothesis) const;


  MultiComponentState extrapolateFromLayerToLayer(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState&,

    const TrackingVolume&,

    const Layer* startLayer,

    const Layer* destinationLayer,

    PropDirection direction,

    ParticleHypothesis particleHypothesis) const;


  MultiComponentState extrapolateToIntermediateLayer(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState&,

    const Layer&,

    const TrackingVolume&,

    PropDirection direction,

    ParticleHypothesis particleHypothesis,

    bool perpendicularCheck = true) const;


  MultiComponentState extrapolateToDestinationLayer(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState&,

    const Surface&,

    const Layer&,

    const Layer*,

    PropDirection direction,

    const BoundaryCheck& boundaryCheck,

    ParticleHypothesis particleHypothesis) const;


  void initialiseNavigation(

    const EventContext& ctx,

    Cache& cache,

    const MultiComponentState& initialState,

    const Surface& surface,

    const Layer*& associatedLayer,

    const TrackingVolume*& currentVolume,

    const TrackingVolume*& destinationVolume,

    std::unique_ptr<TrackParameters>& referenceParameters,

    PropDirection direction) const;


  ToolHandle<IPropagator> m_propagator{ this, "Propagator", "", "" };

  ToolHandle<INavigator> m_navigator{ this,

                                      "Navigator",

                                      "Trk::Navigator/Navigator",

                                      "" };

  ToolHandle<IMaterialMixtureConvolution> m_materialUpdator{

    this,

    "GsfMaterialConvolution",

    "Trk::GsfMaterialMixtureConvolution/GsfMaterialMixtureConvolution",

    ""

  };


  bool m_fastField;

  Trk::MagneticFieldProperties m_fieldProperties;

};


} // end namespace Trk


#endif