MaterialStep.h

Go to the documentation of this file.

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

// MaterialStep.h, (c) ATLAS Detector software
 
#ifndef TRKGEOMETRY_MATERIALSTEP_H
#define TRKGEOMETRY_MATERIALSTEP_H
 
#include <iostream>
#include <vector>
 
#include "TrkGeometry/Material.h"
 
class MsgStream;
 
class MaterialStepCnv_p1;
 
namespace Trk {

class MaterialStep {
 public:
  MaterialStep() = default;

  MaterialStep(float x, float y, float z, float t, float matX0, float matL0,
               float matA, float matZ, float matRho, float matdEdX = 0.);

  MaterialStep(float x, float y, float z, float t, float matX0, float matL0,
               float matA, float matZ, float matRho,
               const std::vector<unsigned char>& elements,
               const std::vector<unsigned char>& fractions, float matdEdX = 0.);

  MaterialStep(const MaterialStep& mstep) = default;
  MaterialStep(MaterialStep&& mstep) = default;

  MaterialStep& operator=(const MaterialStep& mstep) = default;
  MaterialStep& operator=(MaterialStep&& mstep) = default;

  ~MaterialStep() = default;

  MsgStream& dump(MsgStream& sl) const;

  std::ostream& dump(std::ostream& sl) const;

  double steplength() const;

  double steplengthInX0() const;

  double steplengthInL0() const;

  double hitX() const;
  double hitY() const;
  double hitZ() const;
  double hitR() const;

  double x0() const;
  double l0() const;
  double A() const;
  double Z() const;
  double rho() const;
 
  const Material& fullMaterial() const;
 
 protected:
  friend class ::MaterialStepCnv_p1;
  float m_steplength = 0;
  float m_hitX = 0;
  float m_hitY = 0;
  float m_hitZ = 0;
  float m_hitR = 0;
  Material m_material{};
};
 
inline double MaterialStep::steplength() const { return m_steplength; }
 
inline double MaterialStep::steplengthInX0() const {
  return (m_steplength / m_material.X0);
}
 
inline double MaterialStep::steplengthInL0() const {
  return (m_steplength / m_material.L0);
}
 
inline double MaterialStep::hitX() const { return m_hitX; }
 
inline double MaterialStep::hitY() const { return m_hitY; }
 
inline double MaterialStep::hitZ() const { return m_hitZ; }
 
inline double MaterialStep::hitR() const { return m_hitR; }
 
inline double MaterialStep::x0() const { return m_material.X0; }
 
inline double MaterialStep::l0() const { return m_material.L0; }
 
inline double MaterialStep::A() const { return m_material.A; }
 
inline double MaterialStep::Z() const { return m_material.Z; }
 
inline double MaterialStep::rho() const { return m_material.rho; }
 
inline const Material& MaterialStep::fullMaterial() const { return m_material; }

MsgStream& operator<<(MsgStream& sl, const MaterialStep& mstep);
std::ostream& operator<<(std::ostream& sl, const MaterialStep& mstep);
 
}  // namespace Trk
 
#endif