MaterialAccountingStep Class Reference

#include <MaterialAccountingStep.h>

Public Member Functions
void	clear (void)
double	energyLoss (void) const
const GlobalPoint &	in (void) const
double	length (void) const
	MaterialAccountingStep (void)
	MaterialAccountingStep (double position, double radlen, double loss, const GlobalPoint &in, const GlobalPoint &out)
MaterialAccountingStep &	operator*= (const MaterialAccountingStep &step)
	multiply two steps, usefull to implement (co)variance More...
MaterialAccountingStep &	operator*= (double x)
	multiply by a scalar More...
MaterialAccountingStep &	operator+= (const MaterialAccountingStep &step)
	add a step More...
MaterialAccountingStep &	operator-= (const MaterialAccountingStep &step)
	subtract a step More...
MaterialAccountingStep &	operator/= (double x)
	divide by a scalar More...
MaterialAccountingStep &	operator= (const MaterialAccountingStep &step)
	assignement operator More...
const GlobalPoint &	out (void) const
double	radiationLengths (void) const
std::pair< MaterialAccountingStep, MaterialAccountingStep >	split (double fraction) const
	split the step (0..1) in (0..f) + (f..1) using linear interpolation More...

Private Attributes
double	m_energyLoss
GlobalPoint	m_in
double	m_length
GlobalPoint	m_out
double	m_radiationLengths

Detailed Description

Definition at line 9 of file MaterialAccountingStep.h.

Constructor & Destructor Documentation

MaterialAccountingStep() [1/2]

MaterialAccountingStep::MaterialAccountingStep ( void  )

inline

Definition at line 11 of file MaterialAccountingStep.h.

: m_length(0.), m_radiationLengths(0.), m_energyLoss(0.), m_in(), m_out() {}

MaterialAccountingStep() [2/2]

MaterialAccountingStep::MaterialAccountingStep ( double  position, double  radlen, double  loss, const GlobalPoint &  in, const GlobalPoint &  out  )

inline

Definition at line 13 of file MaterialAccountingStep.h.

      : m_length(position), m_radiationLengths(radlen), m_energyLoss(loss), m_in(in), m_out(out) {}

Member Function Documentation

clear()

void MaterialAccountingStep::clear ( void  )

inline

Definition at line 16 of file MaterialAccountingStep.h.

References m_energyLoss, m_in, m_length, m_out, and m_radiationLengths.

Referenced by MaterialAccountingDetector::clear(), and MaterialAccountingTrack::reset().

                   {
    m_length = 0.;
    m_radiationLengths = 0.;
    m_energyLoss = 0.;
    m_in = GlobalPoint();
    m_out = GlobalPoint();
  }

energyLoss()

double MaterialAccountingStep::energyLoss ( void  ) const

inline

Definition at line 36 of file MaterialAccountingStep.h.

References m_energyLoss.

Referenced by DD4hep_MaterialAccountingGroup::averageEnergyLoss(), MaterialAccountingGroup::averageEnergyLoss(), DD4hep_MaterialAccountingGroup::endOfTrack(), MaterialAccountingGroup::endOfTrack(), DD4hep_MaterialAccountingGroup::sigmaEnergyLoss(), and MaterialAccountingGroup::sigmaEnergyLoss().

{ return m_energyLoss; }

in()

const GlobalPoint& MaterialAccountingStep::in ( void  ) const

inline

Definition at line 38 of file MaterialAccountingStep.h.

References m_in.

Referenced by DD4hep_MaterialAccountingGroup::endOfTrack(), and MaterialAccountingGroup::endOfTrack().

{ return m_in; }

length()

double MaterialAccountingStep::length ( void  ) const

inline

Definition at line 32 of file MaterialAccountingStep.h.

References m_length.

Referenced by DD4hep_MaterialAccountingGroup::averageLength(), MaterialAccountingGroup::averageLength(), MaterialAccountingTrack::enterDetector(), MaterialAccountingTrack::leaveDetector(), DD4hep_MaterialAccountingGroup::sigmaLength(), and MaterialAccountingGroup::sigmaLength().

{ return m_length; }

operator*=() [1/2]

MaterialAccountingStep& MaterialAccountingStep::operator*= ( const MaterialAccountingStep &  step )

inline

multiply two steps, usefull to implement (co)variance

Definition at line 100 of file MaterialAccountingStep.h.

References m_energyLoss, m_length, and m_radiationLengths.

                                                                         {
    m_length *= step.m_length;
    m_radiationLengths *= step.m_radiationLengths;
    m_energyLoss *= step.m_energyLoss;
    return *this;
  }

operator*=() [2/2]

MaterialAccountingStep& MaterialAccountingStep::operator*= ( double  x )

inline

multiply by a scalar

Definition at line 108 of file MaterialAccountingStep.h.

References m_energyLoss, m_length, m_radiationLengths, and x.

                                               {
    m_length *= x;
    m_radiationLengths *= x;
    m_energyLoss *= x;
    return *this;
  }

operator+=()

MaterialAccountingStep& MaterialAccountingStep::operator+= ( const MaterialAccountingStep &  step )

inline

add a step

Definition at line 68 of file MaterialAccountingStep.h.

References m_energyLoss, m_in, m_length, m_out, m_radiationLengths, or, and PV3DBase< T, PVType, FrameType >::perp2().

                                                                         {
    m_length += step.m_length;
    m_radiationLengths += step.m_radiationLengths;
    m_energyLoss += step.m_energyLoss;

    // assume that perp2 is 0 only for uninitialized steps
    if ((m_in.perp2() == 0.0) or (step.m_in.perp2() < m_in.perp2()))
      m_in = step.m_in;

    if ((m_out.perp2() == 0.0) or (step.m_out.perp2() > m_out.perp2()))
      m_out = step.m_out;

    return *this;
  }

operator-=()

MaterialAccountingStep& MaterialAccountingStep::operator-= ( const MaterialAccountingStep &  step )

inline

subtract a step

Definition at line 84 of file MaterialAccountingStep.h.

References m_energyLoss, m_in, m_length, m_out, m_radiationLengths, and PV3DBase< T, PVType, FrameType >::perp2().

                                                                         {
    m_length -= step.m_length;
    m_radiationLengths -= step.m_radiationLengths;
    m_energyLoss -= step.m_energyLoss;

    // can anything more sensible be done for m_in and/or m_out ?
    if ((step.m_in.perp2() <= m_in.perp2()) and (step.m_out.perp2() >= m_in.perp2()))
      m_in = step.m_out;

    if ((step.m_out.perp2() >= m_out.perp2()) and (step.m_in.perp2() <= m_out.perp2()))
      m_out = step.m_in;

    return *this;
  }

operator/=()

MaterialAccountingStep& MaterialAccountingStep::operator/= ( double  x )

inline

divide by a scalar

Definition at line 116 of file MaterialAccountingStep.h.

References m_energyLoss, m_length, m_radiationLengths, and x.

                                               {
    m_length /= x;
    m_radiationLengths /= x;
    m_energyLoss /= x;
    return *this;
  }

operator=()

MaterialAccountingStep& MaterialAccountingStep::operator= ( const MaterialAccountingStep &  step )

inline

assignement operator

Definition at line 58 of file MaterialAccountingStep.h.

References m_energyLoss, m_in, m_length, m_out, and m_radiationLengths.

                                                                        {
    m_length = step.m_length;
    m_radiationLengths = step.m_radiationLengths;
    m_energyLoss = step.m_energyLoss;
    m_in = step.m_in;
    m_out = step.m_out;
    return *this;
  }

out()

const GlobalPoint& MaterialAccountingStep::out ( void  ) const

inline

Definition at line 40 of file MaterialAccountingStep.h.

References m_out.

Referenced by DD4hep_MaterialAccountingGroup::endOfTrack(), and MaterialAccountingGroup::endOfTrack().

{ return m_out; }

radiationLengths()

double MaterialAccountingStep::radiationLengths ( void  ) const

inline

Definition at line 34 of file MaterialAccountingStep.h.

References m_radiationLengths.

Referenced by DD4hep_MaterialAccountingGroup::averageRadiationLengths(), MaterialAccountingGroup::averageRadiationLengths(), DD4hep_MaterialAccountingGroup::endOfTrack(), MaterialAccountingGroup::endOfTrack(), DD4hep_MaterialAccountingGroup::sigmaRadiationLengths(), and MaterialAccountingGroup::sigmaRadiationLengths().

{ return m_radiationLengths; }

split()

std::pair<MaterialAccountingStep, MaterialAccountingStep> MaterialAccountingStep::split ( double  fraction ) const

inline

split the step (0..1) in (0..f) + (f..1) using linear interpolation

Definition at line 43 of file MaterialAccountingStep.h.

References HLT_2022v15_cff::fraction, remoteMonitoring_LASER_era2018_cfg::limit, m_energyLoss, m_in, m_length, m_out, m_radiationLengths, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                       {
    // no check is done to ensure that 0 <= f <= 1 !
    GlobalPoint limit(m_in.x() * fraction + m_out.x() * (1. - fraction),
                      m_in.y() * fraction + m_out.y() * (1. - fraction),
                      m_in.z() * fraction + m_out.z() * (1. - fraction));

    MaterialAccountingStep part1(
        fraction * m_length, fraction * m_radiationLengths, fraction * m_energyLoss, m_in, limit);

    MaterialAccountingStep part2(
        (1 - fraction) * m_length, (1 - fraction) * m_radiationLengths, (1 - fraction) * m_energyLoss, limit, m_out);
    return std::make_pair(part1, part2);
  }

Member Data Documentation

m_energyLoss

double MaterialAccountingStep::m_energyLoss

private

Definition at line 27 of file MaterialAccountingStep.h.

Referenced by clear(), energyLoss(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), and split().

m_in

GlobalPoint MaterialAccountingStep::m_in

private

Definition at line 28 of file MaterialAccountingStep.h.

Referenced by clear(), in(), operator+=(), operator-=(), operator=(), and split().

m_length

double MaterialAccountingStep::m_length

private

Definition at line 25 of file MaterialAccountingStep.h.

Referenced by clear(), length(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), and split().

m_out

GlobalPoint MaterialAccountingStep::m_out

private

Definition at line 29 of file MaterialAccountingStep.h.

Referenced by clear(), operator+=(), operator-=(), operator=(), out(), and split().

m_radiationLengths

double MaterialAccountingStep::m_radiationLengths

private

Definition at line 26 of file MaterialAccountingStep.h.

Referenced by clear(), operator*=(), operator+=(), operator-=(), operator/=(), operator=(), radiationLengths(), and split().