ElectronCombinedMaterialEffects.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TrkElectronCombinedMaterialEffects_H
#define TrkElectronCombinedMaterialEffects_H
 
#include "TrkGaussianSumFilterUtils/GsfConstants.h"
#include "TrkGaussianSumFilterUtils/GsfMaterial.h"
//
#include "TrkEventPrimitives/PropDirection.h"
#include "TrkGeometry/MaterialProperties.h"
#include "TrkParameters/ComponentParameters.h"
 
namespace Trk {
 
class ElectronCombinedMaterialEffects
{
public:
  struct ComponentValues
  {
    double weight;
    double mean;
    double variance;
  };
  using MixtureParameters =
    std::array<ComponentValues, GSFConstants::maxNumberofMatComponents>;
 
  using Polynomial = std::array<double, GSFConstants::polynomialCoefficients>;
 
  // ctor with arguments
  ElectronCombinedMaterialEffects(
    const std::string& parameterisationFileName,
    const std::string& parameterisationFileNameHighX0);
 
  // ctor with arguments
  ElectronCombinedMaterialEffects() = default;
  ElectronCombinedMaterialEffects(const ElectronCombinedMaterialEffects&) =
    default;
  ElectronCombinedMaterialEffects(ElectronCombinedMaterialEffects&&) = default;
  ElectronCombinedMaterialEffects& operator=(
    const ElectronCombinedMaterialEffects&) = default;
  ElectronCombinedMaterialEffects& operator=(
    ElectronCombinedMaterialEffects&&) = default;
  ~ElectronCombinedMaterialEffects() = default;
 
  void compute(GsfMaterial::Combined&,
               const Trk::ComponentParameters&,
               const Trk::MaterialProperties&,
               double,
               Trk::PropDirection = anyDirection) const;
 
private:
  // Electron enrrgy loss due to Bremsstrahlung
  void BetheHeitler(GsfMaterial::EnergyLoss& cache,
                    const ComponentParameters& componentParameters,
                    const MaterialProperties& materialProperties,
                    double pathLenght,
                    PropDirection direction = anyDirection) const;
 
  int m_BHnumberOfComponents{};
  int m_BHtransformationCode{};
  int m_BHnumberOfComponentsHighX0{};
  int m_BHtransformationCodeHighX0{};
 
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialWeights{};
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialMeans{};
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialVariances{};
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialWeightsHighX0{};
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialMeansHighX0{};
  std::array<Polynomial, GSFConstants::maxNumberofMatComponents>
    m_BHpolynomialVariancesHighX0{};
};
 
} // end Trk namespace
 
#endif