#include <MaterialEffectsEngine.h>

Inheritance diagram for Trk::MaterialEffectsEngine:

Collaboration diagram for Trk::MaterialEffectsEngine:

Public Member Functions
	MaterialEffectsEngine (const std::string &, const std::string &, const IInterface *)
	Constructor. More...

	~MaterialEffectsEngine ()
	Destructor. More...

StatusCode	initialize ()
	AlgTool initialize method. More...

StatusCode	finalize ()
	AlgTool finalize method. More...

virtual ExtrapolationCode	handleMaterial (ExCellCharged &ecCharged, PropDirection dir=alongMomentum, MaterialUpdateStage matupstage=fullUpdate) const
	charged extrapolation More...

virtual ExtrapolationCode	handleMaterial (ExCellNeutral &ecNeutral, PropDirection dir=alongMomentum, MaterialUpdateStage matupstage=fullUpdate) const
	neutral extrapolation - only for Fatras, dummy implementation here More...

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) const

bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()
	AlgTool interface methods. More...

Protected Member Functions
TrackParameters *	updateTrackParameters (Trk::TrackParameters &parameters, Trk::ExCellCharged &eCell, Trk::PropDirection dir, Trk::MaterialUpdateStage matupstage) const
	charged extrapolation More...

void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)

void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Protected Attributes
MaterialInteraction	m_interactionFormulae
	the formulas concentrated More...

bool	m_eLossCorrection
	apply the energy loss correction More...

bool	m_eLossMpv
	apply the energy loss correction as most probable value More...

bool	m_mscCorrection
	apply the multiple (coulomb) scattering correction More...

std::string	m_sopPrefix
	prefix for screen output More...

std::string	m_sopPostfix
	prefix for screen output More...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

Material effects engine interface for charged and neutral (fast track simulation) , the update is alwyas on the:

eCell.leadParmaeters && eCell.leadLayer
if eCell.leadPameters == eCell.startParamters clone to new parameters else : update the new parameters

Author: Andreas Salzburger -at - cern.ch

Definition at line 38 of file MaterialEffectsEngine.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ MaterialEffectsEngine()

Trk::MaterialEffectsEngine::MaterialEffectsEngine	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Constructor.

Definition at line 16 of file MaterialEffectsEngine.cxx.

 : AthAlgTool(t,n,p),
   m_eLossCorrection(true),
   m_eLossMpv(true),
   m_mscCorrection(true)
 {
     declareInterface<Trk::IMaterialEffectsEngine>(this);
     // steering of the screen outoput (SOP)
     declareProperty("OutputPrefix"                          , m_sopPrefix);
     declareProperty("OutputPostfix"                         , m_sopPostfix);
     // steering of the material effects engine behaviour
     declareProperty("EnergyLossCorrection"                  , m_eLossCorrection);
     declareProperty("MostProbableEnergyLoss"                , m_eLossMpv);
     declareProperty("MultipleScatteringCorrection"          , m_mscCorrection);
 }

◆ ~MaterialEffectsEngine()

Trk::MaterialEffectsEngine::~MaterialEffectsEngine ( )

default

Destructor.

Member Function Documentation

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyArrayType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleKeyType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
  
   }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	hndl,
		const SG::VarHandleType &
	)

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

   {
     return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                        hndl.value(), 
                                                        hndl.documentation());
   }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                  {
     typedef typename SG::HandleClassifier<T>::type htype;
     return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
   }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ finalize()

StatusCode Trk::MaterialEffectsEngine::finalize ( )

AlgTool finalize method.

Definition at line 45 of file MaterialEffectsEngine.cxx.

 {
     EX_MSG_DEBUG( "", "finalize","",  "successful" );
     return StatusCode::SUCCESS;
 }

◆ handleMaterial() [1/2]

Trk::ExtrapolationCode Trk::MaterialEffectsEngine::handleMaterial	(	Trk::ExCellCharged &	eCell,
		Trk::PropDirection	dir = `alongMomentum`,
		Trk::MaterialUpdateStage	matupstage = `fullUpdate`
	)		const

virtual

charged extrapolation

Implements Trk::IMaterialEffectsEngine.

Definition at line 98 of file MaterialEffectsEngine.cxx.

 {
  
     // the Extrapolator made sure that the layer is the lead layer && the parameters are the lead parameters
     if (eCell.leadLayer && eCell.leadLayer->layerMaterialProperties()){
         EX_MSG_DEBUG( ++eCell.navigationStep, "layer",  eCell.leadLayer->layerIndex().value(), "handleMaterial for charged parameters called.");
         // update the track parameters
         eCell.leadParameters = updateTrackParameters(*eCell.leadParameters,eCell,dir,matupstage);
     }
     // only in case of post update it should not return InProgress
     return Trk::ExtrapolationCode::InProgress;
 }

◆ handleMaterial() [2/2]

Trk::ExtrapolationCode Trk::MaterialEffectsEngine::handleMaterial	(	Trk::ExCellNeutral &	eCell,
		Trk::PropDirection	dir = `alongMomentum`,
		Trk::MaterialUpdateStage	matupstage = `fullUpdate`
	)		const

virtual

neutral extrapolation - only for Fatras, dummy implementation here

neutral extrapolation - just collect material

Implements Trk::IMaterialEffectsEngine.

Definition at line 52 of file MaterialEffectsEngine.cxx.

 {
  
     // the Extrapolator made sure that the layer is the lead layer && the parameters are the lead parameters
     if (eCell.leadLayer && eCell.leadLayer->layerMaterialProperties()){
         EX_MSG_DEBUG( ++eCell.navigationStep, "layer",  eCell.leadLayer->layerIndex().value(), "handleMaterial for neutral parameters called - collect material.");
         // now calculate the pathCorrection from the layer surface - it is signed, gives you the relative direction to the layer
         const Trk::Layer* layer = eCell.leadLayer;
         // path correction
         double pathCorrection = layer->surfaceRepresentation().pathCorrection(eCell.leadParameters->position(),dir*(eCell.leadParameters->momentum()));
         // the relative direction wrt with the layer
         Trk::PropDirection rlDir = (pathCorrection >= 0. ? Trk::alongMomentum : Trk::oppositeMomentum);
         // multiply by the pre-and post-update factor
         double mFactor = layer->layerMaterialProperties()->factor(rlDir, matupstage);
         if (mFactor == 0.){
             EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material collection with "  << (matupstage > 0. ? "pre " : "post ")  << "factor 0.");
             // return the parameters untouched -
             return Trk::ExtrapolationCode::InProgress;
         }
         pathCorrection = mFactor*pathCorrection;
         // screen output
         EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material update with corr factor = " << pathCorrection);
         // get the actual material bin
         const Trk::MaterialProperties* materialProperties = layer->layerMaterialProperties()->fullMaterial(eCell.leadParameters->position());
         // and let's check if there's acutally something to do
         if (materialProperties && std::abs(pathCorrection)>0.){
             // thickness in X0
             double thicknessInX0          = materialProperties->thicknessInX0();
             // check if material filling was requested
             if (eCell.checkConfigurationMode(Trk::ExtrapolationMode::CollectMaterial)){
                 EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "collecting material of [t/X0] = " << thicknessInX0);
                 eCell.stepMaterial(eCell.leadParameters->associatedSurface(), layer, eCell.leadParameters->position(), pathCorrection, materialProperties);
             } else {
                 EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "adding material of [t/X0] = " << thicknessInX0);
                 eCell.addMaterial(pathCorrection, materialProperties);
             }
         }
     }
     // only in case of post update it should not return InProgress
     return Trk::ExtrapolationCode::InProgress;
 }

◆ initialize()

StatusCode Trk::MaterialEffectsEngine::initialize ( )

AlgTool initialize method.

Definition at line 38 of file MaterialEffectsEngine.cxx.

 {
     EX_MSG_DEBUG( "", "initialize","",  "successful" );
     return StatusCode::SUCCESS;
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ interfaceID()

static const InterfaceID& Trk::IMaterialEffectsEngine::interfaceID ( )

inlinestaticinherited

AlgTool interface methods.

Definition at line 47 of file IMaterialEffectsEngine.h.

47 { return IID_IMaterialEffectsEngine; }

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                                {
     return this->msgLevel(lvl);
   }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

   {
     h.renounce();
     PBASE::renounce (h);
   }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateTrackParameters()

Trk::TrackParameters * Trk::MaterialEffectsEngine::updateTrackParameters	(	Trk::TrackParameters &	parameters,
		Trk::ExCellCharged &	eCell,
		Trk::PropDirection	dir,
		Trk::MaterialUpdateStage	matupstage
	)		const

protected

charged extrapolation

dE/dl ionization energy loss per path unit

multiple scattering as function of dInX0

Definition at line 114 of file MaterialEffectsEngine.cxx.

 {
     // now calculate the pathCorrection from the layer surface - it is signed, gives you the relative direction to the layer
     const Trk::Layer* layer = eCell.leadLayer;
     // path correction
     double pathCorrection = layer->surfaceRepresentation().pathCorrection(parameters.position(),dir*(parameters.momentum()));
     // the relative direction wrt with the layer
     Trk::PropDirection rlDir = (pathCorrection >= 0. ? Trk::alongMomentum : Trk::oppositeMomentum);
     // multiply by the pre-and post-update factor
     double mFactor = layer->layerMaterialProperties()->factor(rlDir, matupstage);
     if (mFactor == 0.){
         EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material update with "  << (matupstage > 0. ? "pre " : "post ")  << "factor 0. No update done.");
         // return the parameters untouched -
         return (&parameters);
     }
     pathCorrection = mFactor*pathCorrection;
     // screen output
     EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material update with corr factor = " << pathCorrection);
     // get the actual material bin
     const Trk::MaterialProperties* materialProperties = layer->layerMaterialProperties()->fullMaterial(parameters.position());
     // and let's check if there's acutally something to do
     if (materialProperties && std::abs(pathCorrection)>0. &&
     ( m_eLossCorrection || m_mscCorrection || eCell.checkConfigurationMode(Trk::ExtrapolationMode::CollectMaterial)) ){
         // and add them
         int sign = int(eCell.materialUpdateMode);
         // a simple cross-check if the parameters are the initial ones
         AmgVector(5)      uParameters = parameters.parameters();
         std::unique_ptr<AmgSymMatrix(5)> uCovariance =
           parameters.covariance()
             ? std::make_unique<AmgSymMatrix(5)>(*parameters.covariance())
             : nullptr;
         // get the material itself & its parameters
         const Trk::Material& material = materialProperties->material();
         double thicknessInX0          = materialProperties->thicknessInX0();
         double thickness              = materialProperties->thickness();
         // calculate energy loss and multiple scattering
         double p      = parameters.momentum().mag();
         double m      = Trk::ParticleMasses::mass[eCell.pHypothesis];
         double E      = sqrt(p*p+m*m);
         double beta   = p/E;
         // (A) - energy loss correction
         if (m_eLossCorrection){
             double sigmaP = 0.;
             double kazl   = 0.;
             double dEdl = sign*dir*Trk::MaterialInteraction::dEdl_ionization(p, material, eCell.pHypothesis, sigmaP, kazl);
             double dE   = thickness*pathCorrection*dEdl;
             sigmaP *= thickness*pathCorrection;
             // calcuate the new momentum
             double newP = sqrt((E+dE)*(E+dE)-m*m);
             uParameters[Trk::qOverP] = parameters.charge()/newP;
             double sigmaDeltaE = thickness*pathCorrection*sigmaP;
             double sigmaQoverP = sigmaDeltaE/std::pow(beta*p,2);
             // update the covariance if needed
             if (uCovariance)
                (*uCovariance)(Trk::qOverP, Trk::qOverP) += sign*sigmaQoverP*sigmaQoverP;
     }
         // (B) - update the covariance if needed
         if (uCovariance && m_mscCorrection){
             double sigmaMS = Trk::MaterialInteraction::sigmaMS(thicknessInX0*pathCorrection, p, beta);
             double sinTheta = sin(parameters.parameters()[Trk::theta]);
             double sigmaDeltaPhiSq = sigmaMS*sigmaMS/(sinTheta*sinTheta);
             double sigmaDeltaThetaSq = sigmaMS*sigmaMS;
             // add or remove @TODO implement check for covariance matrix -> 0
             (*uCovariance)(Trk::phi,Trk::phi)      += sign*sigmaDeltaPhiSq;
             (*uCovariance)(Trk::theta, Trk::theta) += sign*sigmaDeltaThetaSq;
         }
         // check if material filling was requested
         if (eCell.checkConfigurationMode(Trk::ExtrapolationMode::CollectMaterial)){
             EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "collecting material of [t/X0] = " << thicknessInX0);
             eCell.stepMaterial(parameters.associatedSurface(), layer, parameters.position(), pathCorrection, materialProperties);
         } else {
             EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "adding material of [t/X0] = " << thicknessInX0);
             eCell.addMaterial(pathCorrection, materialProperties);
         }
         // now either create new ones or update - only start parameters can not be updated
         if (eCell.leadParameters != eCell.startParameters ){
             EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material update on non-initial parameters.");
             if (uCovariance)
               parameters.updateParameters(uParameters,*uCovariance);
             else
               parameters.updateParameters(uParameters);
         } else {
             EX_MSG_VERBOSE(eCell.navigationStep, "layer",  layer->layerIndex().value(), "material update on initial parameters, creating new ones.");
             // create new parameters
             const Trk::Surface& tSurface = parameters.associatedSurface();
             Trk::TrackParameters* tParameters = tSurface.createUniqueTrackParameters(uParameters[Trk::loc1],
                                                                                      uParameters[Trk::loc2],
                                                                                      uParameters[Trk::phi],
                                                                                      uParameters[Trk::theta],
                                                                                      uParameters[Trk::qOverP],
                                                                                      *uCovariance).release();
             // these are newly created
             return tParameters;
         }
     }
     //note aliasing input ...
     return (&parameters);
 }

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                                {
     // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
     //         << "  size: " << m_vhka.size() << endmsg;
     for (auto &a : m_vhka) {
       std::vector<SG::VarHandleKey*> keys = a->keys();
       for (auto k : keys) {
         k->setOwner(this);
       }
     }
   }

Member Data Documentation

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eLossCorrection

bool Trk::MaterialEffectsEngine::m_eLossCorrection

protected

apply the energy loss correction

Definition at line 71 of file MaterialEffectsEngine.h.

◆ m_eLossMpv

bool Trk::MaterialEffectsEngine::m_eLossMpv

protected

apply the energy loss correction as most probable value

Definition at line 72 of file MaterialEffectsEngine.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_interactionFormulae

MaterialInteraction Trk::MaterialEffectsEngine::m_interactionFormulae

protected

the formulas concentrated

Definition at line 70 of file MaterialEffectsEngine.h.

◆ m_mscCorrection

bool Trk::MaterialEffectsEngine::m_mscCorrection

protected

apply the multiple (coulomb) scattering correction

Definition at line 73 of file MaterialEffectsEngine.h.

◆ m_sopPostfix

std::string Trk::IMaterialEffectsEngine::m_sopPostfix

protectedinherited

prefix for screen output

Definition at line 62 of file IMaterialEffectsEngine.h.

◆ m_sopPrefix

std::string Trk::IMaterialEffectsEngine::m_sopPrefix

protectedinherited

prefix for screen output

Definition at line 61 of file IMaterialEffectsEngine.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ MaterialEffectsEngine()

◆ ~MaterialEffectsEngine()

Member Function Documentation

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ finalize()

◆ handleMaterial() [1/2]

◆ handleMaterial() [2/2]

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateTrackParameters()

◆ updateVHKA()

Member Data Documentation

◆ m_detStore

◆ m_eLossCorrection

◆ m_eLossMpv

◆ m_evtStore

◆ m_interactionFormulae

◆ m_mscCorrection

◆ m_sopPostfix

◆ m_sopPrefix

◆ m_varHandleArraysDeclared

◆ m_vhka