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 17 of file GsfMaterialEffectsUpdator.h.

Member Typedef Documentation

typedef materialEffect::Covariance GsfMaterialEffectsUpdator::Covariance

Definition at line 19 of file GsfMaterialEffectsUpdator.h.

typedef materialEffect::CovIndex GsfMaterialEffectsUpdator::CovIndex

Definition at line 21 of file GsfMaterialEffectsUpdator.h.

typedef materialEffect::Effect GsfMaterialEffectsUpdator::Effect

Definition at line 20 of file GsfMaterialEffectsUpdator.h.

Constructor & Destructor Documentation

GsfMaterialEffectsUpdator::GsfMaterialEffectsUpdator ( float  mass, uint32_t  is  )

inline

Constructor with explicit mass hypothesis

Definition at line 26 of file GsfMaterialEffectsUpdator.h.

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

virtual GsfMaterialEffectsUpdator::~GsfMaterialEffectsUpdator ( )

inlinevirtual

Definition at line 29 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, GsfMaterialEffectsAdapter, GsfMultipleScatteringUpdator, and GsfCombinedMaterialEffectsUpdator.

Referenced by GsfCombinedMaterialEffectsUpdator::compute().

float GsfMaterialEffectsUpdator::mass ( ) const

inline

Particle mass assigned at construction.

Definition at line 39 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 50 of file GsfMaterialEffectsUpdator.h.

References m_size.

Referenced by Vispa.Gui.TextDialog.TextDialog::__init__(), Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::__init__(), Vispa.Main.MainWindow.MainWindow::_loadIni(), GsfBetheHeitlerUpdator::GsfBetheHeitlerUpdator(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{ m_size=is;}

size_t GsfMaterialEffectsUpdator::size ( void  ) const

inline

Definition at line 48 of file GsfMaterialEffectsUpdator.h.

References m_size.

Referenced by 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 15 of file GsfMaterialEffectsUpdator.cc.

References MultiTrajectoryStateAssembler::addState(), SurfaceSideDefinition::afterSurface, alongMomentum, SurfaceSideDefinition::beforeSurface, MultiTrajectoryStateAssembler::combinedState(), TrajectoryStateOnSurface::components(), bookConverter::compute(), ProxyBase11< T >::data(), TrajectoryStateOnSurface::globalParameters(), TrajectoryStateOnSurface::hasError(), MediumProperties::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), GlobalTrajectoryParameters::magneticField(), LocalTrajectoryError::matrix(), Surface::mediumProperties(), query::result, findQualityFiles::size, 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.components().size()>1 )
    throw cms::Exception("LogicError") << "GsfMaterialEffectsUpdator::updateState used with MultiTSOS";
  double 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
  //
  //   edm::LogDebug("GsfMaterialEffectsUpdator") << "found " << Ws.size() << " components\n"
  //                         << "  input state has weight " << TSoS.weight();
  for ( auto const & effect : effects ) {
    //      edm::LogDebug("GsfMaterialEffectsUpdator") << "w, dp, sigp = "
    //       << Ws[ic] << ", "
    //       << dPs[ic] << ", "
    //       << sqrt((deltaErrors[ic])[0][0]);
    //
    // 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));

      //       edm::LogDebug("GsfMaterialEffectsUpdator") 
      //    << "adding state with weight " << weight*Ws[ic];
    }
    else {
      result.addState(TrajectoryStateOnSurface(lp,surface,
                           &(TSoS.globalParameters().magneticField()),
                           side));
    }
  }
  //   edm::LogDebug("GsfMaterialEffectsUpdator") 
  //     << "  output state has weight " << result.combinedState().weight();
  return result.combinedState();
}

Member Data Documentation

uint32_t GsfMaterialEffectsUpdator::m_size

private

Definition at line 54 of file GsfMaterialEffectsUpdator.h.

Referenced by resize(), and size().

float GsfMaterialEffectsUpdator::theMass

private

Definition at line 53 of file GsfMaterialEffectsUpdator.h.

Referenced by mass().