ElectronMaterialMixtureConvolution.h

Go to the documentation of this file.

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

 
#ifndef TrkElectronMaterialMixtureConvolution_H
#define TrkElectronMaterialMixtureConvolution_H
 
#include "TrkGaussianSumFilter/IMaterialMixtureConvolution.h"
//
#include "TrkGaussianSumFilterUtils/ElectronCombinedMaterialEffects.h"
#include "TrkGaussianSumFilterUtils/GsfMaterial.h"
#include "TrkParameters/ComponentParameters.h"
//
#include "AthenaBaseComps/AthAlgTool.h"
#include "GaudiKernel/ToolHandle.h"
 
namespace Trk {
class Layer;
class MaterialProperties;
 
class ElectronMaterialMixtureConvolution final
  : public AthAlgTool
  , virtual public IMaterialMixtureConvolution
{
 
public:
  enum MaterialUpdateType
  {
    Normal = 0,
    Preupdate = 1,
    Postupdate = 2
  };

  ElectronMaterialMixtureConvolution(const std::string&,
                                     const std::string&,
                                     const IInterface*);

  virtual ~ElectronMaterialMixtureConvolution();

  virtual StatusCode initialize() override final;

  virtual MultiComponentState update(std::vector<GsfMaterial::Combined>&,
                                     const MultiComponentState&,
                                     const Layer&,
                                     PropDirection direction = anyDirection) const override final;

  virtual MultiComponentState preUpdate(std::vector<GsfMaterial::Combined>&,
                                        const MultiComponentState&,
                                        const Layer&,
                                        PropDirection direction = anyDirection) const override final;

  virtual MultiComponentState postUpdate(std::vector<GsfMaterial::Combined>&,
                                         const MultiComponentState&,
                                         const Layer&,
                                         PropDirection direction = anyDirection) const override final;
 
  virtual Trk::ParticleHypothesis particleHypothesis() const override final {
    return Trk::electron;
  }
 
 private:
  Trk::MultiComponentState update(std::vector<GsfMaterial::Combined>&,
                                  const Trk::MultiComponentState& inputState,
                                  const Trk::Layer& layer,
                                  Trk::PropDirection direction,
                                  double updateFactor) const;
 
  //pimpl
  std::unique_ptr<ElectronCombinedMaterialEffects> m_materialEffects;
 
  Gaudi::Property<unsigned int> m_maximumNumberOfComponents{
    this,
    "MaximumNumberOfComponents",
    12,
    "Maximum number of components"
  };
 
  Gaudi::Property<std::string> m_parameterisationFileName{
    this,
    "BetheHeitlerParameterisationFileName",
    "GeantSim_LT01_cdf_nC6_O5.par",
    "Parametrization of Bethe Heitler material effects"
  };
 
  Gaudi::Property<std::string> m_parameterisationFileNameHighX0{
    this,
    "BetheHeitlerParameterisationFileNameHighX0",
    "GeantSim_GT01_cdf_nC6_O5.par",
    "Parametrization of Bethe Heitler material effects for high X0"
  };
};
 
} // end Trk namespace
 
#endif