MaterialEffectsUpdator Class Reference

#include <MaterialEffectsUpdator.h>

Inheritance diagram for MaterialEffectsUpdator:

Public Types
typedef materialEffect::Covariance	Covariance
typedef materialEffect::CovIndex	CovIndex
typedef materialEffect::Effect	Effect

Public Member Functions
virtual MaterialEffectsUpdator *	clone () const =0
virtual void	compute (const TrajectoryStateOnSurface &, const PropagationDirection, Effect &effect) const =0
float	mass () const
	MaterialEffectsUpdator (float mass)
virtual TrajectoryStateOnSurface	updateState (const TrajectoryStateOnSurface &TSoS, const PropagationDirection propDir) const
virtual bool	updateStateInPlace (TrajectoryStateOnSurface &TSoS, const PropagationDirection propDir) const
virtual	~MaterialEffectsUpdator ()

Private Attributes
float	theMass

Detailed Description

Interface for adding material effects during propagation. Updates to TrajectoryStateOnSurface are implemented in this class. Ported from ORCA.

Moved "state" into an independent struct "Effect"

Definition at line 56 of file MaterialEffectsUpdator.h.

Member Typedef Documentation

typedef materialEffect::Covariance MaterialEffectsUpdator::Covariance

Definition at line 58 of file MaterialEffectsUpdator.h.

typedef materialEffect::CovIndex MaterialEffectsUpdator::CovIndex

Definition at line 60 of file MaterialEffectsUpdator.h.

typedef materialEffect::Effect MaterialEffectsUpdator::Effect

Definition at line 59 of file MaterialEffectsUpdator.h.

Constructor & Destructor Documentation

MaterialEffectsUpdator::MaterialEffectsUpdator

(

float

mass

)

Constructor with explicit mass hypothesis

Definition at line 8 of file MaterialEffectsUpdator.cc.

: theMass(mass) {}

MaterialEffectsUpdator::~MaterialEffectsUpdator ( )

virtual

Definition at line 10 of file MaterialEffectsUpdator.cc.

{}

Member Function Documentation

virtual MaterialEffectsUpdator* MaterialEffectsUpdator::clone ( ) const

pure virtual

Implemented in EnergyLossUpdator, CombinedMaterialEffectsUpdator, and MultipleScatteringUpdator.

virtual void MaterialEffectsUpdator::compute ( const TrajectoryStateOnSurface &  , const PropagationDirection  , Effect &  effect  ) const

pure virtual

Implemented in CombinedMaterialEffectsUpdator, EnergyLossUpdator, and MultipleScatteringUpdator.

Referenced by GsfMaterialEffectsAdapter::compute(), and updateStateInPlace().

float MaterialEffectsUpdator::mass ( ) const

inline

Particle mass assigned at construction.

Definition at line 83 of file MaterialEffectsUpdator.h.

References clone(), and bookConverter::compute().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), MultipleScatteringUpdator::compute(), EnergyLossUpdator::compute(), EnergyLossUpdator::computeBetheBloch(), oldComputeBetheBloch(), and oldMUcompute().

{ return theMass; }

TrajectoryStateOnSurface MaterialEffectsUpdator::updateState ( const TrajectoryStateOnSurface &  TSoS, const PropagationDirection  propDir  ) const

virtual

Updates TrajectoryStateOnSurface with material effects (momentum and covariance matrix are potentially affected.

Definition at line 15 of file MaterialEffectsUpdator.cc.

References updateStateInPlace().

                                                                                                       {
  TrajectoryStateOnSurface shallowCopy = TSoS;
  // A TSOS is a proxy. Its contents will be really copied only if/when the updateStateInPlace attempts to change them
  return updateStateInPlace(shallowCopy, propDir) ? shallowCopy : TrajectoryStateOnSurface();
}

bool MaterialEffectsUpdator::updateStateInPlace ( TrajectoryStateOnSurface &  TSoS, const PropagationDirection  propDir  ) const

virtual

Updates in place TrajectoryStateOnSurface with material effects (momentum and covariance matrix are potentially affected) Will return 'false' if the 'updateState' would have returned an invalid TSOS Note that the TSoS might be very well unchanged from this method (just like 'updateState' can return the same TSOS)

Definition at line 26 of file MaterialEffectsUpdator.cc.

References materialEffect::Covariance::add(), SurfaceSideDefinition::afterSurface, alongMomentum, anyDirection, SurfaceSideDefinition::atCenterOfSurface, SurfaceSideDefinition::beforeSurface, compute(), materialEffect::Effect::deltaCov, materialEffect::Effect::deltaP, TrajectoryStateOnSurface::hasError(), MediumProperties::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryError::matrix(), Surface::mediumProperties(), oppositeToMomentum, TrajectoryStateOnSurface::surface(), TrajectoryStateOnSurface::surfaceSide(), TrajectoryStateOnSurface::update(), and LocalTrajectoryParameters::updateP().

Referenced by PropagatorWithMaterial::propagateWithPath(), and updateState().

                                                                                          {
  //
  // Check if
  // - material is associated to surface
  // - propagation direction is not anyDirection
  // - side of surface is not atCenterOfSurface (could be handled with 50% material?)
  //
  const Surface& surface = TSoS.surface();
  if (!surface.mediumProperties().isValid() || propDir == anyDirection || TSoS.surfaceSide() == atCenterOfSurface)
    return true;
  //
  // Check, if already on right side of surface
  //
  if ((propDir == alongMomentum && TSoS.surfaceSide() == afterSurface) ||
      (propDir == oppositeToMomentum && TSoS.surfaceSide() == beforeSurface))
    return true;
  //
  // Update momentum. In case of failure: return invalid state
  //
  LocalTrajectoryParameters lp = TSoS.localParameters();
  Effect effect;
  compute(TSoS, propDir, effect);
  if (!lp.updateP(effect.deltaP))
    return false;
  //
  // Update covariance matrix?
  //
  SurfaceSide side = propDir == alongMomentum ? afterSurface : beforeSurface;
  if (TSoS.hasError()) {
    AlgebraicSymMatrix55 eloc = TSoS.localError().matrix();
    effect.deltaCov.add(eloc);
    //TSoS = TrajectoryStateOnSurface(lp,LocalTrajectoryError(eloc),surface, &(TSoS.globalParameters().magneticField()),side);
    //TSoS.update(lp,LocalTrajectoryError(eloc),side);
    TSoS.update(lp, eloc, side);
  } else {
    TSoS.update(lp, side);
    //TSoS = TrajectoryStateOnSurface(lp,surface,&(TSoS.globalParameters().magneticField()),side);
  }
  return true;
}

Member Data Documentation

float MaterialEffectsUpdator::theMass

private

Definition at line 91 of file MaterialEffectsUpdator.h.