ATLAS Offline Software
Public Member Functions | Protected Member Functions | Private Types | Private Member Functions | Private Attributes | List of all members
Trk::LayerMaterialInspector Class Reference

#include <LayerMaterialInspector.h>

Inheritance diagram for Trk::LayerMaterialInspector:
Collaboration diagram for Trk::LayerMaterialInspector:

Public Member Functions

 LayerMaterialInspector (const std::string &, const std::string &, const IInterface *)
 Constructor. More...
 
virtual ~LayerMaterialInspector ()
 Destructor. More...
 
StatusCode initialize ()
 AlgTool initialize method. More...
 
StatusCode finalize ()
 AlgTool finalize method. More...
 
virtual StatusCode process (TrackingGeometry &tgeo) const
 Processor Action to work on TrackingGeometry& tgeo. More...
 
virtual StatusCode process (TrackingVolume &tvol, size_t level=0) const
 Processor Action to work on TrackingVolumes - the level is for the hierachy tree. More...
 
virtual StatusCode process (Layer &lay, size_t level=0) const
 Processor Action to work on Layers. More...
 
virtual StatusCode process (Surface &surf, size_t level=0) const
 Processor Action to work on Surfaces. 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 (IGeometryProcessor, 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< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

virtual StatusCode processNode (const TrackingVolume &tvol, size_t level=0) const
 Processor Action to work on TrackingVolumes - the level is for the hierachy tree. More...
 
virtual StatusCode processNode (const Layer &lay, size_t level=0) const
 Processor Action to work on Layers. More...
 
virtual StatusCode processNode (const Surface &surf, size_t level=0) const
 Processor Action to work on Surfaces. More...
 
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_treeFolder
 
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

This helper AlgTool can be called to record a dedicated tree that writes the LayerMaterialProperties into an ntuple.

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

Definition at line 35 of file LayerMaterialInspector.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

◆ LayerMaterialInspector()

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

Constructor.

Definition at line 27 of file LayerMaterialInspector.cxx.

27  :
29  m_treeFolder("/val/")
30 {}

◆ ~LayerMaterialInspector()

Trk::LayerMaterialInspector::~LayerMaterialInspector ( )
virtualdefault

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.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

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

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

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

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

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

200  {
201  return PBASE::declareProperty(t);
202  }

◆ DeclareInterfaceID()

Trk::IGeometryProcessor::DeclareInterfaceID ( IGeometryProcessor  ,
,
 
)
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
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ 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
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

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

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ 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
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

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

◆ 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
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

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

◆ declareProperty() [6/6]

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

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ 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::RecursiveGeometryProcessor::finalize ( )
inherited

AlgTool finalize method.

Definition at line 39 of file RecursiveGeometryProcessor.cxx.

40 {
41  ATH_MSG_INFO( "finalize() successful" );
42  return StatusCode::SUCCESS;
43 }

◆ initialize()

StatusCode Trk::RecursiveGeometryProcessor::initialize ( )
inherited

AlgTool initialize method.

Definition at line 32 of file RecursiveGeometryProcessor.cxx.

33 {
34  ATH_MSG_INFO( "initialize()" );
35  return StatusCode::SUCCESS;
36 }

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

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

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

◆ process() [1/4]

StatusCode Trk::RecursiveGeometryProcessor::process ( Trk::Layer lay,
size_t  level = 0 
) const
virtualinherited

Processor Action to work on Layers.

Implements Trk::IGeometryProcessor.

Definition at line 128 of file RecursiveGeometryProcessor.cxx.

128  {
129 
130  std::stringstream displayBuffer;
131  for (size_t il = 0; il < level; ++il) displayBuffer << " ";
132  ATH_MSG_VERBOSE(displayBuffer.str() << " processing Layer with Index: " << lay.layerIndex() );
133 
134  // process the node itself
135  if (processNode(lay, level).isFailure()){
136  ATH_MSG_FATAL("Failed to call processNode(const Layer&). Aborting.");
137  return StatusCode::FAILURE;
138  }
139  // get the subsurface array
140  Trk::SurfaceArray* surfArray = lay.surfaceArray();
141  if (surfArray) {
142  Trk::BinnedArraySpan<Trk::Surface * const > layerSurfaces = surfArray->arrayObjects();
143  ATH_MSG_VERBOSE(displayBuffer.str() << " ---> has " << layerSurfaces.size() << " surfaces on the layer.");
144 
145  auto laySurfIter = layerSurfaces.begin();
146  auto laySurfIterEnd = layerSurfaces.end();
147  // loop over the surfaces and draw them
148  for ( ; laySurfIter != laySurfIterEnd; ++laySurfIter) {
149  if (!(*laySurfIter))
150  ATH_MSG_WARNING("Zero-pointer found in SurfaceArray - indicates problem !");
151  if ((*laySurfIter) && process(**laySurfIter, level).isFailure()){
152  ATH_MSG_FATAL("Failed to call process(const Surface&) on confined layer surfaces. Aborting.");
153  return StatusCode::FAILURE;
154  }
155  }
156  }
157  // return SUCCESS
158  return StatusCode::SUCCESS;
159 }

◆ process() [2/4]

StatusCode Trk::RecursiveGeometryProcessor::process ( Trk::Surface surf,
size_t  level = 0 
) const
virtualinherited

Processor Action to work on Surfaces.

Implements Trk::IGeometryProcessor.

Definition at line 162 of file RecursiveGeometryProcessor.cxx.

162  {
163  return processNode(surf, level);
164 }

◆ process() [3/4]

StatusCode Trk::RecursiveGeometryProcessor::process ( Trk::TrackingGeometry tgeo) const
virtualinherited

Processor Action to work on TrackingGeometry& tgeo.

Implements Trk::IGeometryProcessor.

Definition at line 46 of file RecursiveGeometryProcessor.cxx.

46  {
47 
48  ATH_MSG_VERBOSE("Start processing the TrackingGeometry recursively");
49  // retrieve the highest tracking volume
50  Trk::TrackingVolume* worldVolume = tgeo.highestTrackingVolume();
51  if (worldVolume){
52  ATH_MSG_VERBOSE("TrackingVolume '" << worldVolume->volumeName() << "' retrieved as highest level node.");
53  return process(*worldVolume, 0);
54  }
55  // abort job
56  ATH_MSG_FATAL("No highest level TrackingVolume found. Stopping recursive parsing, abort job.");
57  return StatusCode::FAILURE;
58 }

◆ process() [4/4]

StatusCode Trk::RecursiveGeometryProcessor::process ( Trk::TrackingVolume tvol,
size_t  level = 0 
) const
virtualinherited

Processor Action to work on TrackingVolumes - the level is for the hierachy tree.

Implements Trk::IGeometryProcessor.

Definition at line 61 of file RecursiveGeometryProcessor.cxx.

61  {
62 
63  std::stringstream displayBuffer;
64  for (size_t il = 0; il < level; ++il) displayBuffer << " ";
65  // formatted screen output
66  ATH_MSG_VERBOSE(displayBuffer.str() << "TrackingVolume '" << tvol.volumeName() << "'");
67 
68  // create the action on the volume part of the TrackingVolume
69  if (processNode(tvol, level).isFailure() ){
70  ATH_MSG_FATAL("Failed to call processNode(const TrackingVolume&). Aborting.");
71  return StatusCode::FAILURE;
72  }
73 
74  // Process the contained layers if they exist
75  Trk::LayerArray* layerArray = tvol.confinedLayers();
76  if (layerArray) {
77  // display output
78  auto layers = layerArray->arrayObjects();
79  ATH_MSG_VERBOSE(displayBuffer.str() << "--> has " << layers.size() << " confined layers." );
80  for (const auto & layIter : layers){
81  if (!layIter)
82  ATH_MSG_WARNING("Zero-pointer found in LayerArray - indicates problem !");
83  if ((layIter) && process(*layIter, level).isFailure()){
84  ATH_MSG_FATAL("Failed to call process(const Layer&) on confined layers. Aborting.");
85  return StatusCode::FAILURE;
86  }
87  }
88  }
89 
90  // Process the boundary surface layers
91  auto bSurfaces = tvol.boundarySurfaces();
92  for (auto & bSurface : bSurfaces){
93  if (bSurface->surfaceRepresentation().associatedLayer()){
94  ATH_MSG_VERBOSE(displayBuffer.str() << "--> has a boundary layer." );
95  if (process(
96  const_cast<Trk::Layer&>(
97  *bSurface->surfaceRepresentation().associatedLayer()),level)
98  .isFailure()) {
99  ATH_MSG_FATAL("Failed to call process(const Layer&) on boundary "
100  "layer. Aborting.");
101  return StatusCode::FAILURE;
102  }
103  }
104  }
105 
106 
107  // Process the contained TrackingVolumes (recursively) if they exist
109  // register the next round
110  if (confinedVolumes) {
111  auto volumes = confinedVolumes->arrayObjects();
112  for (const auto & volumesIter : volumes){
113  if (!volumesIter)
114  ATH_MSG_WARNING("Zero-pointer found in VolumeArray - indicates problem !");
115  if (volumesIter && process(*volumesIter, ++level).isFailure() ){
116  ATH_MSG_FATAL("Failed to call process(const TrackingVolume&) on confined volumes. Aborting.");
117  return StatusCode::FAILURE;
118  }
119  }
120  }
121 
122  // return
123  return StatusCode::SUCCESS;
124 
125 }

◆ processNode() [1/3]

StatusCode Trk::LayerMaterialInspector::processNode ( const Layer lay,
size_t  level = 0 
) const
privatevirtual

Processor Action to work on Layers.

Reimplemented from Trk::RecursiveGeometryProcessor.

Definition at line 41 of file LayerMaterialInspector.cxx.

42 {
43 
44  const Trk::TrackingVolume* tvol = lay.enclosingTrackingVolume();
45  TString folderId = tvol ? tvol->volumeName() : std::string("Unknown");
46  TString folderName = m_treeFolder+folderId;
47  folderName.ReplaceAll("::","_");
48 
49  // prepare the information
50  int lIndex = lay.layerIndex().value();
51  const Amg::Vector3D& lCenter = lay.surfaceRepresentation().center();
52  const Trk::LayerMaterialProperties* lMaterial = lay.layerMaterialProperties();
53 
54  // skip navigation layers
55  if (!lMaterial) return StatusCode::SUCCESS;
56 
57  int lType = 0; // 1 - cylinder, 2 - disc
58  float dim0 = 0.;
59  float dim1 = 0.;
60 
61  int bins0 = 1;
62  int bins1 = 1;
63 
64  float centerX = lCenter.x();
65  float centerY = lCenter.y();
66  float centerZ = lCenter.z();
67 
68  // cylinder bounds
69  if ( lay.surfaceRepresentation().type() == Trk::SurfaceType::Cylinder ){
70  lType = 1;
71  // cylinder bounds
72  const Trk::CylinderBounds* cb = dynamic_cast<const Trk::CylinderBounds*>(&(lay.surfaceRepresentation().bounds()));
73  if (cb){
74  dim0 = cb->r();
75  dim1 = cb->halflengthZ();
76  }
77  } else if ( lay.surfaceRepresentation().type() == Trk::SurfaceType::Disc ) {
78  lType = 2;
79  // disc bounds
80  const Trk::DiscBounds* db = dynamic_cast<const Trk::DiscBounds*>(&(lay.surfaceRepresentation().bounds()));
81  if (db){
82  dim0 = db->rMin();
83  dim1 = db->rMax();
84  }
85  }
86 
87  if (lMaterial){
88  const Trk::BinUtility* lBinUtility = lMaterial->binUtility();
89  bins0 = lBinUtility ? lBinUtility->max(0)+1 : 1;
90  bins1 = lBinUtility ? lBinUtility->max(1)+1 : 1;
91  }
92 
93 
94  // -------------------------------------------------------------------------------------
95  TString hName = ( lType == 1 ) ? "CylinderLayer_" : "DiscLayer_";
96  hName += int(lIndex);
97 
98  TString pXo = "_pX0";
99  TString info = "_Information";
100 
101  TH2F* lMaterialHist = lType == 1 ?
102  new TH2F(hName+pXo, hName, bins0, -M_PI*dim0, M_PI*dim0, bins1, -dim1, dim1) :
103  new TH2F(hName+pXo, hName, bins0, dim0, dim1, bins1, -M_PI, M_PI);
104 
105  TTree* lTreeInformation = new TTree(hName+info,hName);
106  lTreeInformation->Branch("LayerCenterX", &centerX, "cX/F");
107  lTreeInformation->Branch("LayerCenterY", &centerY, "cY/F");
108  lTreeInformation->Branch("LayerCenterZ", &centerZ, "cZ/F");
109  lTreeInformation->Branch("LayerIndex", &lIndex, "idx/I");
110  if (lType == 1){
111  lTreeInformation->Branch("Radius", &dim0, "r/F");
112  lTreeInformation->Branch("HalflengthZ", &dim1, "hz/F");
113  lTreeInformation->Branch("BinsRPhi", &bins0, "bin0/I");
114  lTreeInformation->Branch("BinsHalfZ", &bins1, "bin1/I");
115  } else if (lType == 2){
116  lTreeInformation->Branch("InnerRadius", &dim0, "rmin/F");
117  lTreeInformation->Branch("OuterRadius", &dim1, "rmax/F");
118  lTreeInformation->Branch("BinsR", &bins0, "bin0/I");
119  lTreeInformation->Branch("BinsPhi", &bins1, "bin1/I");
120  }
121 
122  // --------------- registration for histogram and tree ---------------------------------
123  TString regHistName = folderName+"/"+hName+"/"+hName+pXo;
124  TString regTreeName = folderName+"/"+hName+"/"+hName+info;
125 
126  ITHistSvc* tHistSvc = nullptr;
127  if (service("THistSvc",tHistSvc).isFailure()) {
128  ATH_MSG_FATAL( "initialize() Could not find Hist Service! Aborting." );
129  return StatusCode::FAILURE;
130  }
131 
132  if ( (tHistSvc->regTree(std::string(regTreeName.Data()),lTreeInformation)).isFailure() ){
133  ATH_MSG_WARNING("Could not register TTree with name " << hName);
134  }
135  lTreeInformation->Fill();
136 
137  if ( (tHistSvc->regHist(std::string(regHistName.Data()),lMaterialHist)).isFailure() ){
138  ATH_MSG_WARNING("Could not register THist with name " << hName);
139  }
140 
141  if (lMaterial){
142  // fill the material into the bins
143  for ( int ib0 = 0; ib0 < bins0; ++ib0 )
144  for ( int ib1 = 0; ib1 < bins1; ++ib1 ){
145  // get the material
146  const Trk::MaterialProperties* mps = lMaterial->material(size_t(ib0),size_t(ib1));
147  if (mps)
148  lMaterialHist->SetBinContent(ib0+1,ib1+1,mps->thicknessInX0());
149  }
150  }
151 
152  return StatusCode::SUCCESS;
153 }

◆ processNode() [2/3]

StatusCode Trk::LayerMaterialInspector::processNode ( const Surface surf,
size_t  level = 0 
) const
privatevirtual

Processor Action to work on Surfaces.

Reimplemented from Trk::RecursiveGeometryProcessor.

Definition at line 155 of file LayerMaterialInspector.cxx.

156 {
157  return StatusCode::SUCCESS;
158 }

◆ processNode() [3/3]

StatusCode Trk::LayerMaterialInspector::processNode ( const TrackingVolume tvol,
size_t  level = 0 
) const
privatevirtual

Processor Action to work on TrackingVolumes - the level is for the hierachy tree.

Reimplemented from Trk::RecursiveGeometryProcessor.

Definition at line 36 of file LayerMaterialInspector.cxx.

37 {
38  return StatusCode::SUCCESS;
39 }

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

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

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

364  {
365  handlesArray.renounce();
366  }

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

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

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_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_treeFolder

std::string Trk::LayerMaterialInspector::m_treeFolder
private

Definition at line 56 of file LayerMaterialInspector.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:
grepfile.info
info
Definition: grepfile.py:38
PlotCalibFromCool.il
il
Definition: PlotCalibFromCool.py:381
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
Trk::RecursiveGeometryProcessor::processNode
virtual StatusCode processNode(const TrackingVolume &tvol, size_t level=0) const
Dedicated action for the different processors.
Definition: RecursiveGeometryProcessor.cxx:167
python.PerfMonSerializer.p
def p
Definition: PerfMonSerializer.py:743
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
CaloCellPos2Ntuple.int
int
Definition: CaloCellPos2Ntuple.py:24
TH2F
Definition: rootspy.cxx:420
Trk::BinUtility::max
size_t max(size_t ba=0) const
First bin maximal value.
Definition: BinUtility.h:215
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
Trk::TrackingVolume::boundarySurfaces
std::vector< SharedObject< BoundarySurface< TrackingVolume > > > & boundarySurfaces()
Method to return the BoundarySurfaces.
Definition: TrackingVolume.cxx:982
CaloCondBlobAlgs_fillNoiseFromASCII.db
db
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:43
Trk::MaterialProperties::thicknessInX0
float thicknessInX0() const
Return the radiationlength fraction.
CSV_InDetExporter.new
new
Definition: CSV_InDetExporter.py:145
module_driven_slicing.layers
layers
Definition: module_driven_slicing.py:114
M_PI
#define M_PI
Definition: ActiveFraction.h:11
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
Trk::RecursiveGeometryProcessor::process
virtual StatusCode process(TrackingGeometry &tgeo) const
Processor Action to work on TrackingGeometry& tgeo.
Definition: RecursiveGeometryProcessor.cxx:46
Trk::LayerMaterialInspector::m_treeFolder
std::string m_treeFolder
Definition: LayerMaterialInspector.h:56
Trk::LayerMaterialProperties::material
virtual const MaterialProperties * material(size_t ib0, size_t ib1) const =0
Direct access via bins to the MaterialProperties.
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
python.TrigEgammaMonitorHelper.TH2F
def TH2F(name, title, nxbins, bins_par2, bins_par3, bins_par4, bins_par5=None, bins_par6=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:45
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
Trk::TrackingVolume::confinedLayers
const LayerArray * confinedLayers() const
Return the subLayer array.
Trk::Layer::surfaceArray
const SurfaceArray * surfaceArray() const
Return the entire SurfaceArray, returns nullptr if no SurfaceArray.
Trk::TrackingGeometry::highestTrackingVolume
const TrackingVolume * highestTrackingVolume() const
return the world
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
beamspotman.n
n
Definition: beamspotman.py:731
Trk::CylinderBounds
Definition: CylinderBounds.h:46
Trk::LayerMaterialProperties
Definition: LayerMaterialProperties.h:62
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Trk::LayerMaterialProperties::binUtility
virtual const BinUtility * binUtility() const =0
Return the BinUtility.
TH2F::SetBinContent
void SetBinContent(int, double)
Definition: rootspy.cxx:426
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
Trk::BinUtility
Definition: BinUtility.h:39
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Trk::TrackingVolume::volumeName
const std::string & volumeName() const
Returns the VolumeName - for debug reason, might be depreciated later.
CaloCellTimeCorrFiller.folderName
string folderName
Definition: CaloCellTimeCorrFiller.py:20
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:195
Trk::BinnedArray::arrayObjects
virtual BinnedArraySpan< T *const > arrayObjects()=0
Return all objects of the Array non-const we can still modify the T.
Trk::RecursiveGeometryProcessor
Definition: RecursiveGeometryProcessor.h:39
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
Trk::MaterialProperties
Definition: MaterialProperties.h:40
a
TList * a
Definition: liststreamerinfos.cxx:10
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Trk::TrackingVolume::confinedVolumes
const TrackingVolumeArray * confinedVolumes() const
Return the subLayer array.
Trk::SurfaceType::Disc
@ Disc
Trk::SurfaceType::Cylinder
@ Cylinder
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:616
declareProperty
#define declareProperty(n, p, h)
Definition: BaseFakeBkgTool.cxx:15
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:790
Trk::BinnedArray
Definition: BinnedArray.h:38
Trk::TrackingVolume
Definition: TrackingVolume.h:121
Trk::BinnedArraySpan
std::span< T > BinnedArraySpan
Definition: BinnedArray.h:34
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
fitman.k
k
Definition: fitman.py:528
Trk::DiscBounds
Definition: DiscBounds.h:44
Trk::CylinderBounds::halflengthZ
double halflengthZ() const
This method returns the halflengthZ.
Trk::CylinderBounds::r
virtual double r() const override final
This method returns the radius.
Trk::Layer
Definition: Layer.h:73
Trk::Layer::layerIndex
const LayerIndex & layerIndex() const
get the layerIndex