GsfMaterialEffectsUpdator Class Reference

#include <GsfMaterialEffectsUpdator.h>

Inheritance diagram for GsfMaterialEffectsUpdator:

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

Public Member Functions
virtual GsfMaterialEffectsUpdator *	clone () const =0
virtual void	compute (const TrajectoryStateOnSurface &, const PropagationDirection, Effect effects[]) const =0
	GsfMaterialEffectsUpdator (float mass, uint32_t is)
float	mass () const
size_t	size () const
virtual TrajectoryStateOnSurface	updateState (const TrajectoryStateOnSurface &TSoS, const PropagationDirection propDir) const
virtual	~GsfMaterialEffectsUpdator ()

Protected Member Functions
void	resize (size_t is)

Private Attributes
uint32_t	m_size
float	theMass

Detailed Description

Interface for adding material effects during propagation as a Gaussian mixture. Similar to MaterialEffectsUpdator, but returns MultiTrajectoryState.

Definition at line 15 of file GsfMaterialEffectsUpdator.h.

Member Typedef Documentation

typedef materialEffect::Covariance GsfMaterialEffectsUpdator::Covariance

Definition at line 17 of file GsfMaterialEffectsUpdator.h.

typedef materialEffect::CovIndex GsfMaterialEffectsUpdator::CovIndex

Definition at line 19 of file GsfMaterialEffectsUpdator.h.

typedef materialEffect::Effect GsfMaterialEffectsUpdator::Effect

Definition at line 18 of file GsfMaterialEffectsUpdator.h.

Constructor & Destructor Documentation

GsfMaterialEffectsUpdator::GsfMaterialEffectsUpdator ( float  mass, uint32_t  is  )

inline

Constructor with explicit mass hypothesis

Definition at line 23 of file GsfMaterialEffectsUpdator.h.

: theMass(mass), m_size(is) {}

virtual GsfMaterialEffectsUpdator::~GsfMaterialEffectsUpdator ( )

inlinevirtual

Definition at line 25 of file GsfMaterialEffectsUpdator.h.

{}

Member Function Documentation

virtual GsfMaterialEffectsUpdator* GsfMaterialEffectsUpdator::clone ( ) const

pure virtual

Implemented in GsfBetheHeitlerUpdator, GsfMultipleScatteringUpdator, GsfCombinedMaterialEffectsUpdator, and GsfMaterialEffectsAdapter.

virtual void GsfMaterialEffectsUpdator::compute ( const TrajectoryStateOnSurface &  , const PropagationDirection  , Effect  effects[]  ) const

pure virtual

Implemented in GsfBetheHeitlerUpdator, GsfMultipleScatteringUpdator, GsfMaterialEffectsAdapter, and GsfCombinedMaterialEffectsUpdator.

Referenced by GsfCombinedMaterialEffectsUpdator::compute().

float GsfMaterialEffectsUpdator::mass ( ) const

inline

Particle mass assigned at construction.

Definition at line 35 of file GsfMaterialEffectsUpdator.h.

References theMass.

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), and GsfMultipleScatteringUpdator::compute().

{ return theMass; }

void GsfMaterialEffectsUpdator::resize ( size_t  is )

inlineprotected

Definition at line 45 of file GsfMaterialEffectsUpdator.h.

References m_size.

Referenced by GsfBetheHeitlerUpdator::GsfBetheHeitlerUpdator().

{ m_size = is; }

size_t GsfMaterialEffectsUpdator::size ( void  ) const

inline

Definition at line 42 of file GsfMaterialEffectsUpdator.h.

References m_size.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and GsfCombinedMaterialEffectsUpdator::compute().

{ return m_size; }

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

virtual

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

Definition at line 14 of file GsfMaterialEffectsUpdator.cc.

References MultiTrajectoryStateAssembler::addState(), SurfaceSideDefinition::afterSurface, alongMomentum, SurfaceSideDefinition::beforeSurface, MultiTrajectoryStateAssembler::combinedState(), bookConverter::compute(), ProxyBase11< T >::data(), materialEffect::elos, TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::hasError(), MediumProperties::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), LogDebug, GlobalTrajectoryParameters::magneticField(), LocalTrajectoryError::matrix(), Surface::mediumProperties(), mps_fire::result, TrajectoryStateOnSurface::singleState(), findQualityFiles::size, mathSSE::sqrt(), TrajectoryStateOnSurface::surface(), LocalTrajectoryParameters::updateP(), histoStyle::weight, and TrajectoryStateOnSurface::weight().

Referenced by FullConvolutionWithMaterial::operator()().

                                                                                                          {
  //
  // get components of input state and check if material is associated to surface
  //
  const Surface& surface = TSoS.surface();
  if (!surface.mediumProperties().isValid())
    return TSoS;
  SurfaceSide side = propDir == alongMomentum ? afterSurface : beforeSurface;
  // single input state?
  if (!TSoS.singleState())
    throw cms::Exception("LogicError") << "GsfMaterialEffectsUpdator::updateState used with MultiTSOS";
  auto weight = TSoS.weight();
//
// Get components (will force recalculation, if necessary)
//
#if __clang__
  std::vector<Effect> effects(size());
  compute(TSoS, propDir, effects.data());
#else
  Effect effects[size()];
  compute(TSoS, propDir, effects);
#endif

  //
  // prepare output vector
  //
  MultiTrajectoryStateAssembler result;
  //
  // loop over components
  //
  LogDebug("GsfMaterialEffectsUpdator") << "found " << size() << " components "
                                        << "  input state has weight " << TSoS.weight();
  for (auto const& effect : effects) {
    LogDebug("GsfMaterialEffectsUpdatorDETAIL") << "w, dp, sigp = " << effect.weight << ", " << effect.deltaP << ", "
                                                << std::sqrt(effect.deltaCov[materialEffect::elos]);
    //
    // Update momentum. In case of failure: return invalid state.
    // Use deltaP method to ensure update of cache, if necessary!
    //
    LocalTrajectoryParameters lp = TSoS.localParameters();
    if (!lp.updateP(effect.deltaP))
      return TrajectoryStateOnSurface();
    //
    // Update covariance matrix?
    //
    if (TSoS.hasError()) {
      AlgebraicSymMatrix55 eloc = TSoS.localError().matrix();
      effect.deltaCov.add(eloc);
      result.addState(TrajectoryStateOnSurface(weight * effect.weight,
                                               lp,
                                               LocalTrajectoryError(eloc),
                                               surface,
                                               &(TSoS.globalParameters().magneticField()),
                                               side));
      LogDebug("GsfMaterialEffectsUpdatorDETAIL") << "adding state with weight " << weight * effect.weight;
    } else {
      result.addState(TrajectoryStateOnSurface(lp, surface, &(TSoS.globalParameters().magneticField()), side));
    }
  }
  LogDebug("GsfMaterialEffectsUpdator") << "  output state has weight " << result.combinedState().weight();
  return result.combinedState();
}

Member Data Documentation

uint32_t GsfMaterialEffectsUpdator::m_size

private

Definition at line 49 of file GsfMaterialEffectsUpdator.h.

Referenced by resize(), and size().

float GsfMaterialEffectsUpdator::theMass

private

Definition at line 48 of file GsfMaterialEffectsUpdator.h.

Referenced by mass().