State of the Cerenkov photon in the quartz optics. More...

#include <PhotonState.h>

Public Types
enum	EType { c_Undefined = 0 , c_BarSegment = 1 , c_MirrorSegment = 2 , c_Prism = 3 }
	enumerator for the type of quartz segments More...

Public Member Functions
	PhotonState ()
	Default constructor.

	PhotonState (const ROOT::Math::XYZPoint &position, const ROOT::Math::XYZVector &direction)
	Constructor with position and direction vectors.

	PhotonState (const ROOT::Math::XYZPoint &position, double kx, double ky, double kz)
	Constructor with position vector and direction components.

	PhotonState (const ROOT::Math::XYZPoint &position, const ROOT::Math::XYZVector &trackDir, double thc, double fic)
	Constructor with position, track direction and cerenkov angles.

PhotonState &	flipKy ()
	Changes sign of direction y component.

ROOT::Math::XYZPoint	getPosition () const
	Returns position as 3D vector.

double	getX () const
	Returns position in x.

double	getY () const
	Returns position in y.

double	getZ () const
	Returns position in z.

double	getUnfoldedX (double x) const
	Unfolds the position in x.

double	getUnfoldedY (double y) const
	Unfolds the position in y.

double	getXD () const
	Returns detection position x in unfolded prism.

double	getYD () const
	Returns detection position y in unfolded prism.

double	getZD () const
	Returns detection position z in unfolded prism.

ROOT::Math::XYZVector	getDirection () const
	Returns direction as 3D unit vector.

double	getKx () const
	Returns direction in x.

double	getKy () const
	Returns direction in y.

double	getKz () const
	Returns direction in z.

double	getPropagationLen () const
	Returns total propagation length since initial position.

int	getNx () const
	Returns number of reflections in x at last propagation step.

int	getNy () const
	Returns number of reflections in y at last propagation step.

double	getA () const
	Returns width (dimension in x) of the quartz segment at last propagation step.

double	getB () const
	Returns thickness (dimension in y) of the quartz segment at last propagation step.

EType	getSegmentType () const
	Returns the type of the quartz segment at last propagation.

bool	getPropagationStatus () const
	Returns propagation status.

bool	getTotalReflStatus (double cosTotal) const
	Returns total internal reflection status.

bool	isInside (const RaytracerBase::BarSegment &bar) const
	Checks if photon is inside the bar segment (including surface).

bool	isInside (const RaytracerBase::BarSegment &bar, const RaytracerBase::Mirror &mirror) const
	Checks if photon is inside the mirror segment (including surface).

bool	isInside (const RaytracerBase::Prism &prism) const
	Checks if photon is inside the prism (including surface).

void	propagate (const RaytracerBase::BarSegment &bar)
	Propagate photon to the exit of bar segment.

void	propagateSemiLinear (const RaytracerBase::BarSegment &bar, const RaytracerBase::Mirror &mirror)
	Propagate photon to the mirror and reflect it using semi-linear mirror optics.

void	propagateExact (const RaytracerBase::BarSegment &bar, const RaytracerBase::Mirror &mirror)
	Propagate photon to the mirror and reflect it using exact mirror optics.

void	propagate (const RaytracerBase::Prism &prism)
	Propagate photon in the prism to the detector plane.

Static Public Member Functions
static void	setMaxPropagationLen (double maxLen)
	Sets maximal allowed propagation length.

Private Attributes
double	m_x = 0
	position in x

double	m_y = 0
	position in y

double	m_z = 0
	position in z

double	m_kx = 0
	direction in x

double	m_ky = 0
	direction in y

double	m_kz = 0
	direction in z

double	m_propLen = 0
	propagation length since initial position

int	m_nx = 0
	signed number of reflections in x at last propagation step

int	m_ny = 0
	signed number of reflections in y at last propagation step

double	m_cosx = 0
	maximal cosine of impact angle to surface in x

double	m_cosy = 0
	maximal cosine of impact angle to surface in y

double	m_A = 0
	width of the quartz segment (dimension in x) for unfolding

double	m_B = 0
	thickness of the quartz segment (dimension in y) for unfolding

double	m_y0 = 0
	origin in y for unfolding

double	m_yD = 0
	unfolded prism detection position in y

double	m_zD = 0
	unfolded prism detection position in z

EType	m_type = c_Undefined
	quartz segment type at last propagation step

bool	m_status = false
	propagation status

Static Private Attributes
static double	s_maxLen = 10000
	maximal allowed propagation length

Detailed Description

State of the Cerenkov photon in the quartz optics.

Definition at line 27 of file PhotonState.h.

Member Enumeration Documentation

◆ EType

enum EType

enumerator for the type of quartz segments

Enumerator
c_Undefined	undefined
c_BarSegment	bar segment
c_MirrorSegment	mirror segment
c_Prism	prism

Definition at line 34 of file PhotonState.h.

                 {
        c_Undefined = 0,     
        c_BarSegment = 1,    
        c_MirrorSegment = 2, 
        c_Prism = 3          
      };

Constructor & Destructor Documentation

◆ PhotonState() [1/4]

PhotonState ( )

inline

Default constructor.

Definition at line 44 of file PhotonState.h.

45 {}

◆ PhotonState() [2/4]

PhotonState	(	const ROOT::Math::XYZPoint &	position,
		const ROOT::Math::XYZVector &	direction )

Constructor with position and direction vectors.

Parameters

position	initial photon position (must be inside quartz)
direction	initial direction vector (must be unit vector)

Definition at line 26 of file PhotonState.cc.

                                                                                                :
      m_x(position.X()), m_y(position.Y()), m_z(position.Z()),
      m_kx(direction.X()), m_ky(direction.Y()), m_kz(direction.Z()),
      m_status(true)
    {}

◆ PhotonState() [3/4]

PhotonState	(	const ROOT::Math::XYZPoint &	position,
		double	kx,
		double	ky,
		double	kz )

Constructor with position vector and direction components.

Parameters

position	initial photon position (must be inside quartz)
kx	initial direction x-component (must be unit vector)
ky	initial direction y-component (must be unit vector)
kz	initial direction z-component (must be unit vector)

Definition at line 33 of file PhotonState.cc.

                                                                                             :
      m_x(position.X()), m_y(position.Y()), m_z(position.Z()),
      m_kx(kx), m_ky(ky), m_kz(kz),
      m_status(true)
    {}

◆ PhotonState() [4/4]

PhotonState	(	const ROOT::Math::XYZPoint &	position,
		const ROOT::Math::XYZVector &	trackDir,
		double	thc,
		double	fic )

Constructor with position, track direction and cerenkov angles.

Parameters

position	initial photon position (must be inside quartz)
trackDir	track direction vector (must be unit vector)
thc	Cerenkov (polar) angle
fic	Cerenkov azimuthal angle

Definition at line 39 of file PhotonState.cc.

                                                    :
      m_x(position.X()), m_y(position.Y()), m_z(position.Z()),
      m_status(true)
    {
      ROOT::Math::XYZVector dir(cos(fic) * sin(thc), sin(fic) * sin(thc), cos(thc));
      func::rotateUz(dir, trackDir);
      m_kx = dir.X();
      m_ky = dir.Y();
      m_kz = dir.Z();
    }

Member Function Documentation

◆ flipKy()

PhotonState & flipKy ( )

inline

Changes sign of direction y component.

Definition at line 75 of file PhotonState.h.

75{m_ky = -m_ky; return *this;}

◆ getA()

double getA ( ) const

inline

Returns width (dimension in x) of the quartz segment at last propagation step.

Returns: width

Definition at line 185 of file PhotonState.h.

185{return m_A;}

◆ getB()

double getB ( ) const

inline

Returns thickness (dimension in y) of the quartz segment at last propagation step.

For prism it is the size of exit window.

Returns: thickness

Definition at line 192 of file PhotonState.h.

192{return m_B;}

◆ getDirection()

ROOT::Math::XYZVector getDirection ( ) const

inline

Returns direction as 3D unit vector.

Returns: direction

Definition at line 143 of file PhotonState.h.

143{return ROOT::Math::XYZVector(m_kx, m_ky, m_kz);}

◆ getKx()

double getKx ( ) const

inline

Returns direction in x.

Returns: direction in x

Definition at line 149 of file PhotonState.h.

149{return m_kx;}

◆ getKy()

double getKy ( ) const

inline

Returns direction in y.

Returns: direction in y

Definition at line 155 of file PhotonState.h.

155{return m_ky;}

◆ getKz()

double getKz ( ) const

inline

Returns direction in z.

Returns: direction in z

Definition at line 161 of file PhotonState.h.

161{return m_kz;}

◆ getNx()

int getNx ( ) const

inline

Returns number of reflections in x at last propagation step.

Returns: signed number of reflections

Definition at line 173 of file PhotonState.h.

173{return m_nx;}

◆ getNy()

int getNy ( ) const

inline

Returns number of reflections in y at last propagation step.

Returns: signed number of reflections

Definition at line 179 of file PhotonState.h.

179{return m_ny;}

◆ getPosition()

ROOT::Math::XYZPoint getPosition ( ) const

inline

Returns position as 3D vector.

Returns: position

Definition at line 87 of file PhotonState.h.

87{return ROOT::Math::XYZPoint(m_x, m_y, m_z);}

◆ getPropagationLen()

double getPropagationLen ( ) const

inline

Returns total propagation length since initial position.

Returns: propagation length

Definition at line 167 of file PhotonState.h.

167{return m_propLen;}

◆ getPropagationStatus()

bool getPropagationStatus ( ) const

inline

Returns propagation status.

Returns: true on success

Definition at line 204 of file PhotonState.h.

204{return m_status;}

◆ getSegmentType()

EType getSegmentType ( ) const

inline

Returns the type of the quartz segment at last propagation.

Returns: type quartz segment type

Definition at line 198 of file PhotonState.h.

198{return m_type;}

◆ getTotalReflStatus()

bool getTotalReflStatus ( double cosTotal ) const

inline

Returns total internal reflection status.

Parameters

cosTotal cosine of total reflection angle

Returns: true if totally reflected

Definition at line 211 of file PhotonState.h.

      {
        return ((m_nx == 0 or m_cosx < cosTotal) and (m_ny == 0 or m_cosy < cosTotal));
      }

◆ getUnfoldedX()

double getUnfoldedX ( double x ) const

inline

Unfolds the position in x.

Parameters

x	position to unfold

Returns: unfolded position

Definition at line 112 of file PhotonState.h.

112{return func::unfold(x, m_nx, m_A);}

◆ getUnfoldedY()

double getUnfoldedY ( double y ) const

inline

Unfolds the position in y.

Parameters

y	position to unfold

Returns: unfolded position

Definition at line 119 of file PhotonState.h.

119{return (func::unfold(y - m_y0, m_ny, m_B) + m_y0);}

◆ getX()

double getX ( ) const

inline

Returns position in x.

Returns: position in x

Definition at line 93 of file PhotonState.h.

93{return m_x;}

◆ getXD()

double getXD ( ) const

inline

Returns detection position x in unfolded prism.

Returns: detection position x in unfolded prism

Definition at line 125 of file PhotonState.h.

125{return getUnfoldedX(m_x);}

◆ getY()

double getY ( ) const

inline

Returns position in y.

Returns: position in y

Definition at line 99 of file PhotonState.h.

99{return m_y;}

◆ getYD()

double getYD ( ) const

inline

Returns detection position y in unfolded prism.

Returns: detection position y in unfolded prism

Definition at line 131 of file PhotonState.h.

131{return m_yD;}

◆ getZ()

double getZ ( ) const

inline

Returns position in z.

Returns: position in z

Definition at line 105 of file PhotonState.h.

105{return m_z;}

◆ getZD()

double getZD ( ) const

inline

Returns detection position z in unfolded prism.

Returns: detection position z in unfolded prism

Definition at line 137 of file PhotonState.h.

137{return m_zD;}

◆ isInside() [1/3]

bool isInside ( const RaytracerBase::BarSegment & bar ) const

Checks if photon is inside the bar segment (including surface).

Parameters

bar	bar segment data

Returns: true if inside

Definition at line 51 of file PhotonState.cc.

    {
      if (std::abs(m_x) > bar.A / 2) return false;
      if (std::abs(m_y) > bar.B / 2) return false;
      if (m_z < bar.zL or m_z > bar.zR) return false;
      return true;
    }

◆ isInside() [2/3]

bool isInside	(	const RaytracerBase::BarSegment &	bar,
		const RaytracerBase::Mirror &	mirror ) const

Checks if photon is inside the mirror segment (including surface).

Parameters

bar	bar segment data
mirror	spherical mirror data

Returns: true if inside

Definition at line 60 of file PhotonState.cc.

    {
      if (std::abs(m_x) > bar.A / 2) return false;
      if (std::abs(m_y) > bar.B / 2) return false;
      if (m_z < bar.zL) return false;
      double Rsq = pow(m_x - mirror.xc, 2) + pow(m_y - mirror.yc, 2) + pow(m_z - mirror.zc, 2);
      if (Rsq > pow(mirror.R, 2)) return false;
      return true;
    }

◆ isInside() [3/3]

bool isInside ( const RaytracerBase::Prism & prism ) const

Checks if photon is inside the prism (including surface).

Parameters

prism prism data

Returns: true if inside

Definition at line 71 of file PhotonState.cc.

    {
      if (std::abs(m_x) > prism.A / 2) return false;
      if (m_z < prism.zL or m_z > prism.zR) return false;
      if (m_y > prism.yUp or m_y < prism.yDown) return false;
      double y = prism.yDown + (prism.yDown + prism.B / 2) / (prism.zFlat - prism.zR) * (m_z - prism.zFlat);
      if (m_y < y) return false;
      return true;
    }

◆ propagate() [1/2]

void propagate ( const RaytracerBase::BarSegment & bar )

Propagate photon to the exit of bar segment.

The state is not reliable if propagation fails.

Parameters

bar	bar segment data

Definition at line 82 of file PhotonState.cc.

    {
      if (not m_status) return;
 
      // local copies: stores to the data members would otherwise be assumed to alias bar
      const double A = bar.A;
      const double B = bar.B;
      double x = m_x;
      double y = m_y;
      double kx = m_kx;
      double ky = m_ky;
 
      m_cosx = std::abs(kx);
      m_cosy = std::abs(ky);
      m_A = A;
      m_B = B;
      m_type = c_BarSegment;
 
      double z = bar.zR;
      if (m_kz < 0) z = bar.zL;
      if (z == m_z) return;
 
      m_status = false;
 
      double len = (z - m_z) / m_kz;
      if (len < 0 or len > s_maxLen) return;
 
      func::fold(x + len * kx, A, x, kx, m_nx);
      func::fold(y + len * ky, B, y, ky, m_ny);
 
      m_propLen += len;
      m_x = x;
      m_y = y;
      m_kx = kx;
      m_ky = ky;
      m_z = z;
 
      m_status = true;
    }

◆ propagate() [2/2]

void propagate ( const RaytracerBase::Prism & prism )

Propagate photon in the prism to the detector plane.

The state is not reliable if propagation fails.

Parameters

prism prism data

Definition at line 260 of file PhotonState.cc.

    {
      if (not m_status) return;
 
      m_status = false;
 
      // Work on local copies of the state and of the prism geometry. Writing to the data members
      // inside the loop below would force the compiler to re-load the window data on every
      // iteration (a double stored through 'this' may alias the doubles of the unfolded windows).
 
      double x = m_x;
      double y = m_y;
      double z = m_z;
      double kx = m_kx;
      double ky = m_ky;
      double kz = m_kz;
      double propLen = m_propLen;
      double cosy = m_cosy;
 
      const double prismA = prism.A;
      const double yUp = prism.yUp;
      const double yDown = prism.yDown;
      const double zR = prism.zR;
      const double zFlat = prism.zFlat;
      const double zDet = prism.zD;
 
      m_cosx = std::abs(kx);
      m_A = prismA;
      m_B = yUp - yDown;
      m_y0 = (yUp + yDown) / 2;
      m_type = c_Prism;
 
      if (kz > 0) {
        if (z >= zR or std::abs(ky / kz) < std::abs(prism.slope)) return;
        if (std::abs(y + ky / kz * (zR - z)) < prism.B / 2) return;
      }
      const double ky_in = ky;
      const double kz_in = kz;
 
      double yD = y;
      double zD = z;
 
      if (z > zFlat) {
 
        int step = 1;
        int ii = 0;
        if (ky < 0) {
          step = -1;
          ii = 1;
          y = std::min(y, yUp);
        }
 
        const auto* windows = prism.unfoldedWindows.data();
        const unsigned numWindows = prism.unfoldedWindows.size();
        const int k0 = prism.k0;
 
        bool found = false;
        unsigned k = k0;
        while (k < numWindows) {
          const auto& win = windows[k];
          double s = ky * win.sz - kz * win.sy;
          if (s != 0) {
            double len = ((win.y0 - y) * win.sz - (win.z0 - z) * win.sy) / s;
            yD = y + len * ky;
            zD = z + len * kz;
            double yu = yD - win.y0;
            double zu = zD - win.z0;
            double yw = yu * win.sy + zu * win.sz;
            if (yw >= yDown and yw <= yUp) {
              if (len < 0 or len > s_maxLen) return;
              double kyNew = ky * win.sy + kz * win.sz;
              double kzNew = kz * win.sy - ky * win.sz;
              int ny = k - k0;
              x += len * kx;
              y = yw;
              ky = kyNew;
              m_ny = ny;
              z = zFlat;
              kz = ny % 2 == 0 ? kzNew : -kzNew;
              propLen += len;
              found = true;
              break;
            }
            cosy = std::max(cosy, std::abs(ky_in * win.nsy[ii] + kz_in * win.nsz[ii]));
          }
          k += step;
          ii ^= 1;
        }
 
        if (not found) {
          B2DEBUG(20, "TOP::PhotonState::propagate: unfolded prism window not found"
                  << LogVar("yUp", prism.yUp) << LogVar("yDown", prism.yDown) << LogVar("zR", prism.zR)
                  << LogVar("y", y) << LogVar("z", z)
                  << LogVar("ky", ky_in) << LogVar("kz", kz_in));
          return;
        }
      }
 
      double len = (zDet - z) / kz;
      if (len < 0 or len > s_maxLen) return;
 
      func::fold(x + len * kx, prismA, x, kx, m_nx);
      y += len * ky;
      z = zDet;
      propLen += len;
      yD += len * ky_in;
      zD += len * kz_in;
 
      m_x = x;
      m_y = y;
      m_z = z;
      m_kx = kx;
      m_ky = ky;
      m_kz = kz;
      m_propLen = propLen;
      m_yD = yD;
      m_zD = zD;
      m_cosy = cosy;
 
      m_status = true;
    }

◆ propagateExact()

void propagateExact	(	const RaytracerBase::BarSegment &	bar,
		const RaytracerBase::Mirror &	mirror )

Propagate photon to the mirror and reflect it using exact mirror optics.

The state is not reliable if propagation fails.

Parameters

bar	mirror segment data
mirror	spherical mirror data

Definition at line 189 of file PhotonState.cc.

    {
      if (not m_status) return;
 
      m_cosx = std::abs(m_kx);
      m_cosy = std::abs(m_ky);
      m_A = bar.A;
      m_B = bar.B;
      m_type = c_MirrorSegment;
 
      m_status = false;
 
      if (m_kz < 0) return;
 
      double len = 0;
      if (m_z < mirror.zb) {
        len = (mirror.zb - m_z) / m_kz;
        if (len > s_maxLen) return;
      }
 
      double xm = m_x + len * m_kx;
      int nx = func::lround(xm / bar.A);
      double ym = m_y + len * m_ky;
      int ny = func::lround(ym / bar.B);
      int i = 0;
      while (true) {
        double xc = func::unfold(mirror.xc, nx, bar.A);
        double yc = func::unfold(mirror.yc, ny, bar.B);
        double x = m_x - xc;
        double y = m_y - yc;
        double z = m_z - mirror.zc;
        double rdir = x * m_kx + y * m_ky + z * m_kz;
        double rr = x * x + y * y + z * z;
        double D = rdir * rdir + (mirror.R * mirror.R - rr);
        if (D < 0) return;
        D = sqrt(D);
        len = (D - rdir);
        if (len < 0 or len > s_maxLen) return;
        double xmm = m_x + len * m_kx;
        int nxx = func::lround(xmm / bar.A);
        double ymm = m_y + len * m_ky;
        int nyy = func::lround(ymm / bar.B);
        if (nxx == nx and nyy == ny) break;
        i++;
        if (i == 10) {
          if (std::abs(xmm - xm) < 0.001 and std::abs(ymm - ym) < 0.001) break;
          B2DEBUG(20, "TOP::PhotonState::propagateExact: not converging");
          return;
        }
        nx = nxx;
        ny = nyy;
      }
 
      m_propLen += len;
 
      func::fold(m_x + len * m_kx, bar.A, m_x, m_kx, m_nx);
      func::fold(m_y + len * m_ky, bar.B, m_y, m_ky, m_ny);
      m_z += len * m_kz;
 
      double normX = (m_x - mirror.xc) / mirror.R;
      double normY = (m_y - mirror.yc) / mirror.R;
      double normZ = (m_z - mirror.zc) / mirror.R;
      double s = 2 * (m_kx * normX + m_ky * normY + m_kz * normZ);
      m_kx -= s * normX;
      m_ky -= s * normY;
      m_kz -= s * normZ;
 
      m_status = true;
    }

◆ propagateSemiLinear()

void propagateSemiLinear	(	const RaytracerBase::BarSegment &	bar,
		const RaytracerBase::Mirror &	mirror )

Propagate photon to the mirror and reflect it using semi-linear mirror optics.

The state is not reliable if propagation fails.

Semi-linear: mirror surface approximated with a cylinder along y-axis and using linear optics approximation in y for the photon reflection. Useful to get rid of some of the discontinuities in the TOP image at given wavelength, which are anyway smeared-out by the dispersion and other processes.

Parameters

bar	mirror segment data
mirror	spherical mirror data

Definition at line 123 of file PhotonState.cc.

    {
      if (not m_status) return;
 
      m_cosx = std::abs(m_kx);
      m_cosy = std::abs(m_ky);
      m_A = bar.A;
      m_B = bar.B;
      m_type = c_MirrorSegment;
 
      m_status = false;
 
      if (m_kz < 0) return;
 
      double len = 0;
      if (m_z < mirror.zb) {
        len = (mirror.zb - m_z) / m_kz;
        if (len > s_maxLen) return;
      }
 
      double xm = m_x + len * m_kx;
      int nx = func::lround(xm / bar.A);
      double ss = m_kx * m_kx + m_kz * m_kz;
      if (ss == 0) return;
      int i = 0;
      while (true) {
        double xc = func::unfold(mirror.xc, nx, bar.A);
        double x = m_x - xc;
        double z = m_z - mirror.zc;
        double rdir = x * m_kx + z * m_kz;
        double rr = x * x + z * z;
        double D = rdir * rdir + (mirror.R * mirror.R - rr) * ss;
        if (D < 0) return;
        D = sqrt(D);
        len = (D - rdir) / ss;
        if (len < 0 or len > s_maxLen) return;
        double xmm = m_x + len * m_kx;
        int nxx = func::lround(xmm / bar.A);
        if (nxx == nx) break;
        i++;
        if (i == 10) {
          if (std::abs(xmm - xm) < 0.001) break;
          B2DEBUG(20, "TOP::PhotonState::propagateSemiLinear: not converging");
          return;
        }
        nx = nxx;
        xm = xmm;
      }
 
      m_propLen += len;
 
      func::fold(m_x + len * m_kx, bar.A, m_x, m_kx, m_nx);
      func::fold(m_y + len * m_ky, bar.B, m_y, m_ky, m_ny);
      m_y = mirror.yc;
      m_z += len * m_kz;
 
      double normX = (m_x - mirror.xc) / mirror.R;
      double normZ = (m_z - mirror.zc) / mirror.R;
      double s = 2 * (m_kx * normX + m_kz * normZ);
      m_kx -= s * normX;
      m_kz -= s * normZ;
 
      m_status = true;
    }

◆ setMaxPropagationLen()

static void setMaxPropagationLen ( double maxLen )

inlinestatic

Sets maximal allowed propagation length.

Parameters

maxLen maximal allowed propagation length

Definition at line 81 of file PhotonState.h.

81{s_maxLen = maxLen;}

Member Data Documentation

◆ m_A

double m_A = 0

private

width of the quartz segment (dimension in x) for unfolding

Definition at line 286 of file PhotonState.h.

◆ m_B

double m_B = 0

private

thickness of the quartz segment (dimension in y) for unfolding

Definition at line 287 of file PhotonState.h.

◆ m_cosx

double m_cosx = 0

private

maximal cosine of impact angle to surface in x

Definition at line 284 of file PhotonState.h.

◆ m_cosy

double m_cosy = 0

private

maximal cosine of impact angle to surface in y

Definition at line 285 of file PhotonState.h.

◆ m_kx

double m_kx = 0

private

direction in x

Definition at line 278 of file PhotonState.h.

◆ m_ky

double m_ky = 0

private

direction in y

Definition at line 279 of file PhotonState.h.

◆ m_kz

double m_kz = 0

private

direction in z

Definition at line 280 of file PhotonState.h.

◆ m_nx

int m_nx = 0

private

signed number of reflections in x at last propagation step

Definition at line 282 of file PhotonState.h.

◆ m_ny

int m_ny = 0

private

signed number of reflections in y at last propagation step

Definition at line 283 of file PhotonState.h.

◆ m_propLen

double m_propLen = 0

private

propagation length since initial position

Definition at line 281 of file PhotonState.h.

◆ m_status

bool m_status = false

private

propagation status

Definition at line 292 of file PhotonState.h.

◆ m_type

EType m_type = c_Undefined

private

quartz segment type at last propagation step

Definition at line 291 of file PhotonState.h.

◆ m_x

double m_x = 0

private

position in x

Definition at line 275 of file PhotonState.h.

◆ m_y

double m_y = 0

private

position in y

Definition at line 276 of file PhotonState.h.

◆ m_y0

double m_y0 = 0

private

origin in y for unfolding

Definition at line 288 of file PhotonState.h.

◆ m_yD

double m_yD = 0

private

unfolded prism detection position in y

Definition at line 289 of file PhotonState.h.

◆ m_z

double m_z = 0

private

position in z

Definition at line 277 of file PhotonState.h.

◆ m_zD

double m_zD = 0

private

unfolded prism detection position in z

Definition at line 290 of file PhotonState.h.

◆ s_maxLen

double s_maxLen = 10000

staticprivate

maximal allowed propagation length

Definition at line 294 of file PhotonState.h.

The documentation for this class was generated from the following files:

top/reconstruction_cpp/include/PhotonState.h
top/reconstruction_cpp/src/PhotonState.cc

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Enumeration Documentation

◆ EType

Constructor & Destructor Documentation

◆ PhotonState() [1/4]

◆ PhotonState() [2/4]

◆ PhotonState() [3/4]

◆ PhotonState() [4/4]

Member Function Documentation

◆ flipKy()

◆ getA()

◆ getB()

◆ getDirection()

◆ getKx()

◆ getKy()

◆ getKz()

◆ getNx()

◆ getNy()

◆ getPosition()

◆ getPropagationLen()

◆ getPropagationStatus()

◆ getSegmentType()

◆ getTotalReflStatus()

◆ getUnfoldedX()

◆ getUnfoldedY()

◆ getX()

◆ getXD()

◆ getY()

◆ getYD()

◆ getZ()

◆ getZD()

◆ isInside() [1/3]

◆ isInside() [2/3]

◆ isInside() [3/3]

◆ propagate() [1/2]

◆ propagate() [2/2]

◆ propagateExact()

◆ propagateSemiLinear()

◆ setMaxPropagationLen()

Member Data Documentation

◆ m_A

◆ m_B

◆ m_cosx

◆ m_cosy

◆ m_kx

◆ m_ky

◆ m_kz

◆ m_nx

◆ m_ny

◆ m_propLen

◆ m_status

◆ m_type

◆ m_x

◆ m_y

◆ m_y0

◆ m_yD

◆ m_z

◆ m_zD

◆ s_maxLen