GXFMaterialEffects.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRKGXFMATERIALEFFECTS_H
#define TRKGXFMATERIALEFFECTS_H
 
#include "TrkMaterialOnTrack/EnergyLoss.h"
 
#include <memory>
 
namespace Trk {
  class MaterialEffectsOnTrack;
  class MaterialEffectsBase;
  class EnergyLoss;
  class Surface;
  class MaterialProperties;

  class GXFMaterialEffects {
 
  public:
   GXFMaterialEffects(const MaterialEffectsOnTrack &);
   GXFMaterialEffects() = default;
   GXFMaterialEffects(const GXFMaterialEffects &rhs);
   GXFMaterialEffects &operator=(const GXFMaterialEffects &rhs);
   GXFMaterialEffects(GXFMaterialEffects &&rhs) = default;
   GXFMaterialEffects &operator=(GXFMaterialEffects &&rhs) = default;
   ~GXFMaterialEffects() = default;
 
   void setScatteringAngles(double, double);
   void setScatteringSigmas(double, double);
   double x0() const;
   void setX0(double);
   double deltaPhi() const;
   double measuredDeltaPhi() const;
   void setMeasuredDeltaPhi(double);
   double deltaTheta() const;
   double sigmaDeltaPhi() const;
   double sigmaDeltaTheta() const;
   double deltaE() const;

   void setEloss(std::unique_ptr<EnergyLoss>);
   double sigmaDeltaE() const;
   double sigmaDeltaEPos() const;
   double sigmaDeltaENeg() const;
   double sigmaDeltaEAve() const;
   void setDeltaE(double);
   void setSigmaDeltaE(double);
   void setSigmaDeltaEPos(double);
   void setSigmaDeltaENeg(double);
   double delta_p() const;
   void setdelta_p(double);
   void setKink(bool);
   bool isKink() const;
   void setMeasuredEloss(bool);
   bool isMeasuredEloss() const;
   const Surface &associatedSurface() const;

   void setSurface(const Surface *);
   std::unique_ptr<MaterialEffectsBase> makeMEOT() const;
   const MaterialProperties *materialProperties() const;

   void setMaterialProperties(const MaterialProperties *);
 
  private:
   double m_scatphi = 0;
   double m_scattheta = 0;
   double m_sigmascatphi = 0;
   double m_sigmascattheta = 0;
   double m_x0 = 0;
   double m_deltap = 0;
   double m_deltae = 0;
   double m_sigmadeltae = 0;
   double m_sigmadeltaepos = 0;
   double m_sigmadeltaeneg = 0;
   std::unique_ptr<const EnergyLoss> m_eloss;
   //not owning ptrs
   const Surface *m_surf = nullptr;
   const MaterialProperties *m_matprop = nullptr;
   bool m_iskink = false;
   bool m_ismeasuredeloss = false;
   double m_measscatphi = 0;  // fudge to stabilize fit in muon system
   double m_sintheta = 1;
  };
}
#endif