Trk::LayerMaterialAnalyser Class Reference

#include <LayerMaterialAnalyser.h>

Inheritance diagram for Trk::LayerMaterialAnalyser:

Collaboration diagram for Trk::LayerMaterialAnalyser:

Public Member Functions
	LayerMaterialAnalyser (const std::string &, const std::string &, const IInterface *)
	Constructor. More...
	~LayerMaterialAnalyser ()
	Destructor. More...
StatusCode	initialize ()
	AlgTool initialize method. More...
StatusCode	finalize ()
	AlgTool finalize method. More...
StatusCode	analyseLayerMaterial (const Layer &lay) const
	process the layer - after material creation and loading More...
StatusCode	analyseLayerMaterial (const Layer &lay, const LayerMaterialProperties &lmp) const
	process the layer material proerties - after material creation and before loading More...
StatusCode	analyseLayerMaterial (const Layer &lay, const MaterialPropertiesMatrix &lmr) const
	process the layer material record - before material creation More...
StatusCode	analyseLayerMaterial (const Layer &lay, const LayerMaterialRecord &lmr) const
	process the layer material record - before material creation 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
	DeclareInterfaceID (ILayerMaterialAnalyser, 1, 0)
	Creates the InterfaceID and interfaceID() method. More...

Protected Member Functions
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...

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	analyse (const Layer &lay, const MaterialPropertiesMatrix &lmr, const std::vector< std::vector< unsigned int > > *bCounter=0) const
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
std::string	m_layerMaterialName
TTree *	m_validationTree
	The validation tree. More...
std::string	m_validationTreeName
	validation tree name - to be accessed by this from root More...
std::string	m_validationTreeDescription
	validation tree description - second argument in TTree More...
std::string	m_validationTreeFolder
	stream/folder to for the TTree to be written out More...
int	m_layerIndex
	the layer index given by the TrackingGeometry More...
int	m_layerType
	the type of the layer 1 - cylinder, 2 - disk More...
std::vector< float > *	m_layerTranslation
	center of the transform More...
std::vector< float > *	m_layerRotation
	orientation of the layer More...
float	m_layerDimension0
	dimension 0 : cylinder r, disk r_min More...
float	m_layerDimension1
	dimension 1 : cylinder z, disk r_max More...
int	m_layerBins
	total number of bins - loc0 * loc 1 More...
int	m_layerBins0
	total number of bins - loc 0 More...
int	m_layerBins1
	total number of bins - loc 0 More...
std::vector< int > *	m_bin0
	bin 0 More...
std::vector< int > *	m_bin1
	bin 1 More...
std::vector< float > *	m_thickness
	gathered thickness from material mapping/material properties More...
std::vector< float > *	m_X0
	gathered X0 from material mapping/material properties More...
std::vector< float > *	m_L0
	gathered L0 from material mapping/material properties More...
std::vector< float > *	m_A
	gathered A from material mapping/material properties More...
std::vector< float > *	m_Z
	gathered Z from material mapping/material properties More...
std::vector< float > *	m_Rho
	gathered rho from material mapping/material properties More...
std::vector< int > *	m_elements
	gathered number of elements from material mapping/material properties More...
std::vector< int > *	m_binCounter
	how often was this bin hit / used More...
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

LayerMaterialProperties creator for CompoundLayerMaterial

Author

Andre.nosp@m.as.S.nosp@m.alzbu.nosp@m.rger.nosp@m.@cern.nosp@m..ch

Definition at line 38 of file LayerMaterialAnalyser.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

LayerMaterialAnalyser()

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

Constructor.

Definition at line 33 of file LayerMaterialAnalyser.cxx.

: AthAlgTool(t,n,p),
  m_layerMaterialName("UnspecifiedLayerMaterialMap"),
  m_validationTree(nullptr),
  m_validationTreeName("LayerMaterialAnalyser"),
  m_validationTreeDescription("LayerMaterialAnalyser information"),
  m_validationTreeFolder("/val/LayerMaterialAnalyser"),
  m_layerIndex(0),        
  m_layerType(0),
  m_layerTranslation(nullptr),
  m_layerRotation(nullptr),
  m_layerDimension0(0.),   
  m_layerDimension1(0.),   
  m_layerBins(0),         
  m_layerBins0(0),         
  m_layerBins1(0),
  m_bin0(nullptr),
  m_bin1(nullptr),
  m_thickness(nullptr),
  m_X0(nullptr),
  m_L0(nullptr),
  m_A(nullptr),
  m_Z(nullptr),
  m_Rho(nullptr),
  m_elements(nullptr),
  m_binCounter(nullptr)
{
    declareInterface<Trk::ILayerMaterialAnalyser>(this);
    // give the map a name
    declareProperty("LayerMaterialName",          m_layerMaterialName);
    declareProperty("ValidationTreeName",         m_validationTreeName);
    declareProperty("ValidationTreeDescription",  m_validationTreeDescription);
    declareProperty("ValidationTreeFolder",       m_validationTreeFolder);
}

~LayerMaterialAnalyser()

Trk::LayerMaterialAnalyser::~LayerMaterialAnalyser ( )

default

Destructor.

Member Function Documentation

analyse()

StatusCode Trk::LayerMaterialAnalyser::analyse ( const Layer &  lay, const MaterialPropertiesMatrix &  lmr, const std::vector< std::vector< unsigned int > > *  bCounter = 0  ) const

private

Definition at line 201 of file LayerMaterialAnalyser.cxx.

{
  
    // general layer information
    m_layerIndex = layer.layerIndex().value();
    const Trk::Surface& lSurface = layer.surfaceRepresentation();
    m_layerTranslation->at(0) = lSurface.center().x();
    m_layerTranslation->at(1) = lSurface.center().y();
    m_layerTranslation->at(2) = lSurface.center().z();
    
    AmgMatrix(3,3) rMatrix = lSurface.transform().rotation();
    m_layerRotation->at(0) = rMatrix(0,0);
    m_layerRotation->at(1) = rMatrix(1,0);
    m_layerRotation->at(2) = rMatrix(2,0);
    m_layerRotation->at(3) = rMatrix(0,1);
    m_layerRotation->at(4) = rMatrix(1,1);
    m_layerRotation->at(5) = rMatrix(2,1);
    m_layerRotation->at(6) = rMatrix(0,2);
    m_layerRotation->at(7) = rMatrix(1,2);
    m_layerRotation->at(8) = rMatrix(2,2);
 
    // cylinder bounds 
    if ( lSurface.type() == Trk::SurfaceType::Cylinder ){
        m_layerType = 1;
        // cylinder bounds
        const Trk::CylinderBounds* cb = dynamic_cast<const Trk::CylinderBounds*>(&(lSurface.bounds()));
        if (cb){
            m_layerDimension0 = cb->r();
            m_layerDimension1 = cb->halflengthZ();
        }
    } else if ( lSurface.type() == Trk::SurfaceType::Disc ) {
        m_layerType = 2;
        // disc bounds
        const Trk::DiscBounds* db = dynamic_cast<const Trk::DiscBounds*>(&(lSurface.bounds()));
        if (db){
            m_layerDimension0 = db->rMin();
            m_layerDimension1 = db->rMax();
        }
    }
 
    // now get the material matrix and record all single bins;
    m_layerBins0  = 0;
    m_layerBins1  = 0;
    m_layerBins   = 0; 
    int bin1      = 0;
    for (const auto & outerIter : mpMatrix){
        int bin0 = 0;
        for (const auto & innerIter : outerIter ){
            m_bin0->at(m_layerBins) = bin0;
            m_bin1->at(m_layerBins) = bin1;
            // get the material
            const Trk::MaterialProperties* mProperties = innerIter;
            if (mProperties){
                m_thickness->at(m_layerBins) = mProperties->thickness();
                m_X0->at(m_layerBins)        = mProperties->x0();
                m_L0->at(m_layerBins)        = mProperties->l0();
                m_A->at(m_layerBins)         = mProperties->averageA();
                m_Z->at(m_layerBins)         = mProperties->averageZ();
                m_Rho->at(m_layerBins)       = mProperties->averageRho();
                m_elements->at(m_layerBins)  = mProperties->material().composition ? mProperties->material().composition->size() : 0;    
            } else  {
                m_thickness->at(m_layerBins) = 0.;
                m_X0->at(m_layerBins)        = 0.;
                m_L0->at(m_layerBins)        = 0.;
                m_A->at(m_layerBins)         = 0.;
                m_Z->at(m_layerBins)         = 0.;
                m_Rho->at(m_layerBins)       = 0.;
                m_elements->at(m_layerBins)  = 0.;                
            }
            // set the bin Counter 
            m_binCounter->at(m_layerBins)    = bCounter ? (*bCounter)[bin1][bin0] : 1;
            //     
            ++bin0;
            if (!bin1) ++m_layerBins0;
            ++m_layerBins;    
        }
        ++bin1;
        ++m_layerBins1;
    }
    m_validationTree->Fill();
    
    // return 0 - since this is only an analyser
    return StatusCode::SUCCESS;
}    

analyseLayerMaterial() [1/4]

StatusCode Trk::LayerMaterialAnalyser::analyseLayerMaterial ( const Layer &  lay ) const

virtual

process the layer - after material creation and loading

Implements Trk::ILayerMaterialAnalyser.

Definition at line 152 of file LayerMaterialAnalyser.cxx.

{
    
    // get the LayerMaterial 
    const Trk::LayerMaterialProperties* lMaterial = layer.layerMaterialProperties();
    if (!lMaterial) return StatusCode::SUCCESS;
 
    return analyseLayerMaterial(layer, *lMaterial);
}

analyseLayerMaterial() [2/4]

StatusCode Trk::LayerMaterialAnalyser::analyseLayerMaterial ( const Layer &  lay, const LayerMaterialProperties &  lmp  ) const

virtual

process the layer material proerties - after material creation and before loading

Implements Trk::ILayerMaterialAnalyser.

Definition at line 163 of file LayerMaterialAnalyser.cxx.

{
    
    // binned material can directly call the filling method
    const Trk::BinnedLayerMaterial* blMaterial = dynamic_cast<const Trk::BinnedLayerMaterial*>(&lMaterial);
    if (blMaterial) {
        ATH_MSG_DEBUG( "Recieved BinnedLayerMaterial - analyzing it." );
        return analyseLayerMaterial(layer, blMaterial->fullMaterial());
    }
    
    // we need to create a MaterialProperties matrix for this
    const Trk::BinUtility* bUtility = lMaterial.binUtility();
    size_t mBins0 = bUtility ? bUtility->max(0)+1 : 1;
    size_t mBins1 = bUtility ? bUtility->max(1)+1 : 1;
    
    Trk::MaterialPropertiesMatrix mpMatrix(mBins1, std::vector< const Trk::MaterialProperties*>(mBins0, nullptr));    
    for (size_t ibin1 = 0; ibin1 < mBins1; ++ ibin1){
        for (size_t ibin0 = 0; ibin0 < mBins0; ++ibin0)
            mpMatrix[ibin1][ibin0] = lMaterial.material(ibin0, ibin1);
    }    
    // now send it to the analyser
    return analyse(layer, mpMatrix);
}

analyseLayerMaterial() [3/4]

StatusCode Trk::LayerMaterialAnalyser::analyseLayerMaterial ( const Layer &  lay, const LayerMaterialRecord &  lmr  ) const

virtual

process the layer material record - before material creation

Implements Trk::ILayerMaterialAnalyser.

Definition at line 195 of file LayerMaterialAnalyser.cxx.

{
    ATH_MSG_DEBUG( "Recieved LayerMaterialRecord  - analyzing it." );
    return analyse(layer, lmRecord.associatedLayerMaterial(), (const std::vector< std::vector< unsigned int > >*)&lmRecord.binCounts());
}

analyseLayerMaterial() [4/4]

StatusCode Trk::LayerMaterialAnalyser::analyseLayerMaterial ( const Layer &  lay, const MaterialPropertiesMatrix &  lmr  ) const

virtual

process the layer material record - before material creation

Implements Trk::ILayerMaterialAnalyser.

Definition at line 188 of file LayerMaterialAnalyser.cxx.

{
    ATH_MSG_DEBUG( "Recieved MaterialPropertyMatrix  - analyzing it." );
    return analyse(layer, mpMatrix);
}

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);
  }

DeclareInterfaceID()

Trk::ILayerMaterialAnalyser::DeclareInterfaceID ( ILayerMaterialAnalyser  , 1  , 0    )

inherited

Creates the InterfaceID and interfaceID() method.

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.

{ 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.

{ 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.

{ 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::LayerMaterialAnalyser::finalize

(

)

AlgTool finalize method.

Definition at line 133 of file LayerMaterialAnalyser.cxx.

{
    delete m_layerTranslation ;
    delete m_layerRotation    ;
    delete m_bin0             ;
    delete m_bin1             ;
    delete m_thickness        ;
    delete m_X0               ;
    delete m_L0               ;
    delete m_A                ;
    delete m_Z                ;
    delete m_Rho              ;
    delete m_elements         ;
    delete m_binCounter       ;
    return StatusCode::SUCCESS;
}

initialize()

StatusCode Trk::LayerMaterialAnalyser::initialize

(

)

AlgTool initialize method.

Definition at line 74 of file LayerMaterialAnalyser.cxx.

{
  
    m_layerTranslation = new std::vector<float>(3, 0.);
    m_layerRotation    = new std::vector<float>(9, 0.);
    m_bin0             = new std::vector<int>(LAYERMAXBINS, 0);
    m_bin1             = new std::vector<int>(LAYERMAXBINS, 0);
    m_thickness        = new std::vector<float>(LAYERMAXBINS, 0.);
    m_X0               = new std::vector<float>(LAYERMAXBINS, 0.);
    m_L0               = new std::vector<float>(LAYERMAXBINS, 0.);
    m_A                = new std::vector<float>(LAYERMAXBINS, 0.);
    m_Z                = new std::vector<float>(LAYERMAXBINS, 0.);
    m_Rho              = new std::vector<float>(LAYERMAXBINS, 0.);
    m_elements         = new std::vector<int>(LAYERMAXBINS, 0);
    m_binCounter       = new std::vector<int>(LAYERMAXBINS, 0);
    
    // now register the Tree
 
    // ------------- validation section ------------------------------------------
    m_validationTree = new TTree(m_validationTreeName.c_str(), m_validationTreeDescription.c_str());
 
    // position coordinates of the update
    m_validationTree->Branch("LayerIndex",            &m_layerIndex         );
    m_validationTree->Branch("LayerType",             &m_layerType          );
    m_validationTree->Branch("LayerTranslation",      &m_layerTranslation   );
    m_validationTree->Branch("LayerRotation",         &m_layerRotation      );
    m_validationTree->Branch("LayerDimension0",       &m_layerDimension0    );
    m_validationTree->Branch("LayerDimension1",       &m_layerDimension1    );
    m_validationTree->Branch("LayerBins",             &m_layerBins          );
    m_validationTree->Branch("LayerBins0",            &m_layerBins0         );
    m_validationTree->Branch("LayerBins1",            &m_layerBins1         );
    m_validationTree->Branch("LayerBin0",             &m_bin0               );
    m_validationTree->Branch("LayerBin1",             &m_bin1               );
    m_validationTree->Branch("LayerBinCounter",       &m_binCounter         );
    m_validationTree->Branch("LayerThickness",        &m_thickness          );
    m_validationTree->Branch("LayerX0",               &m_X0                 );
    m_validationTree->Branch("LayerL0",               &m_L0                 );
    m_validationTree->Branch("LayerA",                &m_A                  );
    m_validationTree->Branch("LayerZ",                &m_Z                  );
    m_validationTree->Branch("LayerRo",               &m_Rho                );
    m_validationTree->Branch("LayerElements",         &m_elements           );
    
    // now register the Tree
    SmartIF<ITHistSvc> tHistSvc{service("THistSvc")};
    if (!tHistSvc) {
        ATH_MSG_ERROR("initialize() Could not find Hist Service -> Switching ValidationMode Off !" );
        delete m_validationTree; m_validationTree = nullptr;
        return StatusCode::SUCCESS;
    }
    if ((tHistSvc->regTree(m_validationTreeFolder.c_str(), m_validationTree)).isFailure()) {
        ATH_MSG_ERROR("initialize() Could not register the validation Tree -> Switching ValidationMode Off !" );
        delete m_validationTree; m_validationTree = nullptr;
        return StatusCode::SUCCESS;
    }
 
    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.

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.

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_A

std::vector<float>* Trk::LayerMaterialAnalyser::m_A

private

gathered A from material mapping/material properties

Definition at line 93 of file LayerMaterialAnalyser.h.

m_bin0

std::vector<int>* Trk::LayerMaterialAnalyser::m_bin0

private

bin 0

Definition at line 88 of file LayerMaterialAnalyser.h.

m_bin1

std::vector<int>* Trk::LayerMaterialAnalyser::m_bin1

private

bin 1

Definition at line 89 of file LayerMaterialAnalyser.h.

m_binCounter

std::vector<int>* Trk::LayerMaterialAnalyser::m_binCounter

private

how often was this bin hit / used

Definition at line 97 of file LayerMaterialAnalyser.h.

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_elements

std::vector<int>* Trk::LayerMaterialAnalyser::m_elements

private

gathered number of elements from material mapping/material properties

Definition at line 96 of file LayerMaterialAnalyser.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_L0

std::vector<float>* Trk::LayerMaterialAnalyser::m_L0

private

gathered L0 from material mapping/material properties

Definition at line 92 of file LayerMaterialAnalyser.h.

m_layerBins

int Trk::LayerMaterialAnalyser::m_layerBins

mutableprivate

total number of bins - loc0 * loc 1

Definition at line 85 of file LayerMaterialAnalyser.h.

m_layerBins0

int Trk::LayerMaterialAnalyser::m_layerBins0

mutableprivate

total number of bins - loc 0

Definition at line 86 of file LayerMaterialAnalyser.h.

m_layerBins1

int Trk::LayerMaterialAnalyser::m_layerBins1

mutableprivate

total number of bins - loc 0

Definition at line 87 of file LayerMaterialAnalyser.h.

m_layerDimension0

float Trk::LayerMaterialAnalyser::m_layerDimension0

mutableprivate

dimension 0 : cylinder r, disk r_min

Definition at line 83 of file LayerMaterialAnalyser.h.

m_layerDimension1

float Trk::LayerMaterialAnalyser::m_layerDimension1

mutableprivate

dimension 1 : cylinder z, disk r_max

Definition at line 84 of file LayerMaterialAnalyser.h.

m_layerIndex

int Trk::LayerMaterialAnalyser::m_layerIndex

mutableprivate

the layer index given by the TrackingGeometry

Definition at line 79 of file LayerMaterialAnalyser.h.

m_layerMaterialName

std::string Trk::LayerMaterialAnalyser::m_layerMaterialName

private

Definition at line 72 of file LayerMaterialAnalyser.h.

m_layerRotation

std::vector<float>* Trk::LayerMaterialAnalyser::m_layerRotation

private

orientation of the layer

Definition at line 82 of file LayerMaterialAnalyser.h.

m_layerTranslation

std::vector<float>* Trk::LayerMaterialAnalyser::m_layerTranslation

private

center of the transform

Definition at line 81 of file LayerMaterialAnalyser.h.

m_layerType

int Trk::LayerMaterialAnalyser::m_layerType

mutableprivate

the type of the layer 1 - cylinder, 2 - disk

Definition at line 80 of file LayerMaterialAnalyser.h.

m_Rho

std::vector<float>* Trk::LayerMaterialAnalyser::m_Rho

private

gathered rho from material mapping/material properties

Definition at line 95 of file LayerMaterialAnalyser.h.

m_thickness

std::vector<float>* Trk::LayerMaterialAnalyser::m_thickness

private

gathered thickness from material mapping/material properties

Definition at line 90 of file LayerMaterialAnalyser.h.

m_validationTree

TTree* Trk::LayerMaterialAnalyser::m_validationTree

private

The validation tree.

Definition at line 74 of file LayerMaterialAnalyser.h.

m_validationTreeDescription

std::string Trk::LayerMaterialAnalyser::m_validationTreeDescription

private

validation tree description - second argument in TTree

Definition at line 76 of file LayerMaterialAnalyser.h.

m_validationTreeFolder

std::string Trk::LayerMaterialAnalyser::m_validationTreeFolder

private

stream/folder to for the TTree to be written out

Definition at line 77 of file LayerMaterialAnalyser.h.

m_validationTreeName

std::string Trk::LayerMaterialAnalyser::m_validationTreeName

private

validation tree name - to be accessed by this from root

Definition at line 75 of file LayerMaterialAnalyser.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.

m_X0

std::vector<float>* Trk::LayerMaterialAnalyser::m_X0

private

gathered X0 from material mapping/material properties

Definition at line 91 of file LayerMaterialAnalyser.h.

m_Z

std::vector<float>* Trk::LayerMaterialAnalyser::m_Z

private

gathered Z from material mapping/material properties

Definition at line 94 of file LayerMaterialAnalyser.h.

---

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

• LayerMaterialAnalyser.h
• LayerMaterialAnalyser.cxx