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

#include <SurfaceIntersectionTest.h>

Inheritance diagram for Trk::SurfaceIntersectionTest:
Collaboration diagram for Trk::SurfaceIntersectionTest:

Public Member Functions

 SurfaceIntersectionTest (const std::string &name, ISvcLocator *pSvcLocator)
 Standard Athena-Algorithm Constructor. More...
 
StatusCode finalize ()
 
StatusCode runTest ()
 
StatusCode initializeTest ()
 initialize the test, i.e. More...
 
StatusCode initialize ()
 standard Athena-Algorithm method More...
 
StatusCode execute ()
 standard Athena-Algorithm method More...
 
virtual StatusCode bookTree ()
 
virtual StatusCode sysInitialize () override
 Override sysInitialize. More...
 
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies. 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 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
 

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

Protected Attributes

Rndm::Numbers * m_gaussDist
 Random Number setup. More...
 
Rndm::Numbers * m_flatDist
 
size_t m_numTests
 number of tests More...
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_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

ToolHandle< ISurfaceBuilderm_surfaceBuilder
 
size_t m_testsPerSurface
 
size_t m_eventCounter
 
std::vector< size_t > m_surfaces
 
std::vector< size_t > m_surfaceAttempts
 
std::vector< size_t > m_surfaceSuccessful
 
std::vector< std::string > m_surfaceNames
 
std::vector< bool > m_surfaceNamesSet
 
bool m_writeTTree
 
std::string m_treeFolder
 
DataObjIDColl m_extendedExtraObjects
 
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

Definition at line 29 of file SurfaceIntersectionTest.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ SurfaceIntersectionTest()

Trk::SurfaceIntersectionTest::SurfaceIntersectionTest ( const std::string &  name,
ISvcLocator *  pSvcLocator 
)

Standard Athena-Algorithm Constructor.

Definition at line 22 of file SurfaceIntersectionTest.cxx.

22  :
23  Trk::TrkDetDescrUnitTestBase(name, pSvcLocator),
25  m_testsPerSurface(100),
26  m_eventCounter(0),
27  m_surfaces(6,0),
28  m_surfaceAttempts(6,0),
30  m_surfaceNames(6,""),
31  m_surfaceNamesSet(6,false),
32  m_writeTTree(true),
33  m_treeFolder("/val/")
34 {
35  declareProperty("SurfaceBuilder", m_surfaceBuilder);
36  declareProperty("NumberOfTestsPerSurface", m_testsPerSurface);
37  declareProperty("WriteTTree", m_writeTTree);
38 
39 }

Member Function Documentation

◆ bookTree()

StatusCode Trk::TrkDetDescrUnitTestBase::bookTree ( )
virtualinherited

Reimplemented in Trk::MappingTest.

Definition at line 60 of file TrkDetDescrUnitTestBase.cxx.

60 { return StatusCode::SUCCESS; }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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; }

◆ execute()

StatusCode Trk::TrkDetDescrUnitTestBase::execute ( )
inherited

standard Athena-Algorithm method

Definition at line 48 of file TrkDetDescrUnitTestBase.cxx.

49 {
50  return runTest();
51 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 50 of file AthAlgorithm.cxx.

51 {
52  // If we didn't find any symlinks to add, just return the collection
53  // from the base class. Otherwise, return the extended collection.
54  if (!m_extendedExtraObjects.empty()) {
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ finalize()

StatusCode Trk::SurfaceIntersectionTest::finalize ( )

Definition at line 41 of file SurfaceIntersectionTest.cxx.

41  {
42 
43  ATH_MSG_INFO("Test statistics output ---------------------- (start) -------------");
44  for (size_t it = 0; it < m_surfaces.size(); ++it){
45  if (m_surfaceNamesSet[it]){
46  ATH_MSG_INFO("Tested " << m_surfaces[it] << " created "<< m_surfaceNames[it] << " objects.");
47  ATH_MSG_INFO(" -> Result [ attempts / successful (percentage)] : "
48  << m_surfaceAttempts[it] << " / " << m_surfaceSuccessful[it]
49  << " ( " << double(m_surfaceSuccessful[it])/m_surfaceAttempts[it] << " )" );
50  }
51  }
52  ATH_MSG_INFO("Test statistics output ----------------------- (end) --------------");
53 
55 }

◆ initialize()

StatusCode Trk::TrkDetDescrUnitTestBase::initialize ( )
inherited

standard Athena-Algorithm method

Definition at line 27 of file TrkDetDescrUnitTestBase.cxx.

28 {
29  ATH_MSG_INFO( "Creating random number services, call bookTree() and initializeTest()" );
30 
31  // intialize the random number generators
32  m_gaussDist = new Rndm::Numbers(randSvc(), Rndm::Gauss(0.,1.));
33  m_flatDist = new Rndm::Numbers(randSvc(), Rndm::Flat(0.,1.));
34 
35  if (bookTree().isFailure()){
36  ATH_MSG_FATAL( "Could not book the TTree object" );
37  return StatusCode::FAILURE;
38  }
39 
40  if (initializeTest().isFailure()){
41  ATH_MSG_FATAL( "Could not initialize the test" );
42  return StatusCode::FAILURE;
43  }
44 
45  return StatusCode::SUCCESS;
46 }

◆ initializeTest()

StatusCode Trk::SurfaceIntersectionTest::initializeTest ( )
virtual

initialize the test, i.e.

retrieve the TrackingGeometry Svc

Reimplemented from Trk::TrkDetDescrUnitTestBase.

Definition at line 57 of file SurfaceIntersectionTest.cxx.

58 {
59  if ( m_surfaceBuilder.retrieve().isFailure() ){
60  ATH_MSG_FATAL("Could not retrieve surface builder. Aborting.");
61  return StatusCode::FAILURE;
62  }
63  // success return
64  return StatusCode::SUCCESS;
65 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

MsgStream& AthCommonMsg< Algorithm >::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< Algorithm >::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< Algorithm > >::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< Algorithm > >::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< Algorithm > >::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  }

◆ runTest()

StatusCode Trk::SurfaceIntersectionTest::runTest ( )
virtual

Implements Trk::TrkDetDescrUnitTestBase.

Definition at line 70 of file SurfaceIntersectionTest.cxx.

71 {
72  ATH_MSG_VERBOSE("Running the SurfaceIntersectionTest Test");
73 
74  // ----------------- creation of the surfaces & track parameters -------------
75  for (size_t it = 0; it < TrkDetDescrUnitTestBase::m_numTests; ++it ){
76 
77  // get a surface from the surface builder
78  const Trk::Surface* sf = m_surfaceBuilder->surface();
79  if (!sf) continue;
80 
81  // get the type
82  Trk::SurfaceType sType = sf->type();
83  ++m_surfaces[static_cast<int>(sType)];
84 
85  ATH_MSG_VERBOSE("Surface of type " << sf->name() << " at postion " << Amg::toString(sf->center()));
86  if (!m_surfaceNamesSet[static_cast<int>(sType)]) {
87  m_surfaceNames[static_cast<int>(sType)] = sf->name();
88  m_surfaceNamesSet[static_cast<int>(sType)] = true;
89  }
90  // create the Tree for this TrackingVolume
91  TString treeName = "IntersectionTest_Event";
93  TTree* currentTree = new TTree(treeName.Data(),"Intersection");
94  treeName += "_Info";
95  TTree* currentInfoTree = new TTree(treeName.Data(),"Surface Information");
96 
97  SmartIF<ITHistSvc> tHistSvc;
98  if (m_writeTTree && (tHistSvc = service("THistSvc")).isValid()) {
99  ATH_MSG_FATAL( "initialize() Could not find Hist Service! Aborting." );
100  return StatusCode::FAILURE;
101  }
102  // create the name under which the tree is registered
103  TString fullTreeName = m_treeFolder+"Intersection_Event";
104  fullTreeName += m_eventCounter;
105  fullTreeName += "_Number";
106  fullTreeName += it;
107  fullTreeName += "_"+sf->name();
108 
109  if (m_writeTTree && tHistSvc->regTree(fullTreeName.Data(), currentTree).isFailure()) {
110  ATH_MSG_ERROR( "initialize() Could not register the validation Tree! Aborting." );
111  delete currentTree;
112  return StatusCode::FAILURE;
113  }
114 
115  TString fullInfoTreeName = m_treeFolder+"SurfaceInfo_Event";
116  fullInfoTreeName += m_eventCounter;
117  fullInfoTreeName += "_Number";
118  fullInfoTreeName += it;
119  fullInfoTreeName += "_"+sf->name();
120  if (m_writeTTree && tHistSvc->regTree(fullInfoTreeName.Data(), currentInfoTree).isFailure()) {
121  ATH_MSG_ERROR( "initialize() Could not register the validation Tree! Aborting." );
122  delete currentInfoTree;
123  return StatusCode::FAILURE;
124  }
125 
126  // the ntuple branches
127  Trk::SurfaceNtupleBranch surfaceBranch;
128  surfaceBranch.initForWrite(*currentInfoTree);
129  surfaceBranch.fill(*sf);
130 
131  Trk::VectorNtupleBranch intersectionBranch;
132  intersectionBranch.initForWrite(*currentTree,"Intersection_",3);
133 
134  {
135  for (size_t ist = 0; ist < m_testsPerSurface; ++ist){
136  // starting around the perigee
137  Amg::Vector3D startPosition;
138  startPosition.x() = TrkDetDescrUnitTestBase::m_gaussDist->shoot()* 50.;
139  startPosition.y() = TrkDetDescrUnitTestBase::m_gaussDist->shoot()* 50.;
140  startPosition.z() = TrkDetDescrUnitTestBase::m_gaussDist->shoot()* 150.;
141  // aiming towards the surface center
142  const Amg::Vector3D& sfCenter = sf->center();
143  Amg::Vector3D direction = sf->center();
144  direction.x() = sfCenter.x()*(1.+0.25*TrkDetDescrUnitTestBase::m_gaussDist->shoot());
145  direction.y() = sfCenter.y()*(1.+0.25*TrkDetDescrUnitTestBase::m_gaussDist->shoot());
146  direction.z() = sfCenter.z()*(1.+0.25*TrkDetDescrUnitTestBase::m_gaussDist->shoot());
147  direction = direction.normalized();
148  ATH_MSG_VERBOSE(" -> Line: " << Amg::toString(startPosition) << " + t * " << Amg::toString(direction));
149 
150  // now doing the distance estimation
151  Trk::Intersection intersection = sf->straightLineIntersection(startPosition,direction,false);
152  ++m_surfaceAttempts[static_cast<int>(sType)];
153  // now check if the intesection is valid or not
154  if (intersection.valid){
155  ++m_surfaceSuccessful[static_cast<int>(sType)];
156  Amg::Vector3D consistency = (intersection.position-intersection.pathLength*direction) - startPosition;
157  if (m_writeTTree){
158  intersectionBranch.fill(intersection.position);
159  currentTree->Fill();
160  }
161  ATH_MSG_VERBOSE( " -o Intersection at " << Amg::toString(intersection.position) );
162  ATH_MSG_VERBOSE( " -x Cross-checked at " << Amg::toString(consistency) );
163  }
164  }
165  // write the tree
166  if (m_writeTTree) currentInfoTree->Fill();
167  delete sf;
168  }
169  }
170  ++m_eventCounter;
171  return StatusCode::SUCCESS;
172 }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

Reimplemented in AthAnalysisAlgorithm, AthFilterAlgorithm, PyAthena::Alg, and AthHistogramAlgorithm.

Definition at line 66 of file AthAlgorithm.cxx.

66  {
68 
69  if (sc.isFailure()) {
70  return sc;
71  }
72  ServiceHandle<ICondSvc> cs("CondSvc",name());
73  for (auto h : outputHandles()) {
74  if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
75  // do this inside the loop so we don't create the CondSvc until needed
76  if ( cs.retrieve().isFailure() ) {
77  ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
78  return StatusCode::SUCCESS;
79  }
80  if (cs->regHandle(this,*h).isFailure()) {
81  sc = StatusCode::FAILURE;
82  ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
83  << " with CondSvc");
84  }
85  }
86  }
87  return sc;
88 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::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< Algorithm > >::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< Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eventCounter

size_t Trk::SurfaceIntersectionTest::m_eventCounter
private

Definition at line 49 of file SurfaceIntersectionTest.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthAlgorithm::m_extendedExtraObjects
privateinherited

Definition at line 79 of file AthAlgorithm.h.

◆ m_flatDist

Rndm::Numbers* Trk::TrkDetDescrUnitTestBase::m_flatDist
protectedinherited

Definition at line 58 of file TrkDetDescrUnitTestBase.h.

◆ m_gaussDist

Rndm::Numbers* Trk::TrkDetDescrUnitTestBase::m_gaussDist
protectedinherited

Random Number setup.

Definition at line 57 of file TrkDetDescrUnitTestBase.h.

◆ m_numTests

size_t Trk::TrkDetDescrUnitTestBase::m_numTests
protectedinherited

number of tests

Definition at line 61 of file TrkDetDescrUnitTestBase.h.

◆ m_surfaceAttempts

std::vector<size_t> Trk::SurfaceIntersectionTest::m_surfaceAttempts
private

Definition at line 52 of file SurfaceIntersectionTest.h.

◆ m_surfaceBuilder

ToolHandle<ISurfaceBuilder> Trk::SurfaceIntersectionTest::m_surfaceBuilder
private

Definition at line 46 of file SurfaceIntersectionTest.h.

◆ m_surfaceNames

std::vector<std::string> Trk::SurfaceIntersectionTest::m_surfaceNames
private

Definition at line 54 of file SurfaceIntersectionTest.h.

◆ m_surfaceNamesSet

std::vector<bool> Trk::SurfaceIntersectionTest::m_surfaceNamesSet
private

Definition at line 55 of file SurfaceIntersectionTest.h.

◆ m_surfaces

std::vector<size_t> Trk::SurfaceIntersectionTest::m_surfaces
private

Definition at line 51 of file SurfaceIntersectionTest.h.

◆ m_surfaceSuccessful

std::vector<size_t> Trk::SurfaceIntersectionTest::m_surfaceSuccessful
private

Definition at line 53 of file SurfaceIntersectionTest.h.

◆ m_testsPerSurface

size_t Trk::SurfaceIntersectionTest::m_testsPerSurface
private

Definition at line 48 of file SurfaceIntersectionTest.h.

◆ m_treeFolder

std::string Trk::SurfaceIntersectionTest::m_treeFolder
private

Definition at line 58 of file SurfaceIntersectionTest.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_writeTTree

bool Trk::SurfaceIntersectionTest::m_writeTTree
private

Definition at line 57 of file SurfaceIntersectionTest.h.


The documentation for this class was generated from the following files:
Trk::TrkDetDescrUnitTestBase::m_flatDist
Rndm::Numbers * m_flatDist
Definition: TrkDetDescrUnitTestBase.h:58
Trk::TrkDetDescrUnitTestBase::m_numTests
size_t m_numTests
number of tests
Definition: TrkDetDescrUnitTestBase.h:61
Trk::SurfaceNtupleBranch::fill
void fill(const Trk::Surface &surface)
fill a vector
Definition: SurfaceNtupleBranch.cxx:30
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
Trk::Intersection
Definition: Intersection.h:24
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::TrkDetDescrUnitTestBase::m_gaussDist
Rndm::Numbers * m_gaussDist
Random Number setup.
Definition: TrkDetDescrUnitTestBase.h:57
Trk::VectorNtupleBranch
This class provides a simple interface to write Vectors to a root tree The size of the vectors cannot...
Definition: VectorNtupleBranch.h:34
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
Trk::TrkDetDescrUnitTestBase
Definition: TrkDetDescrUnitTestBase.h:27
Trk::TrkDetDescrUnitTestBase::bookTree
virtual StatusCode bookTree()
Definition: TrkDetDescrUnitTestBase.cxx:60
skel.it
it
Definition: skel.GENtoEVGEN.py:396
Trk::SurfaceIntersectionTest::m_eventCounter
size_t m_eventCounter
Definition: SurfaceIntersectionTest.h:49
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
Trk::SurfaceType
SurfaceType
Definition: SurfaceTypes.h:17
Trk::SurfaceIntersectionTest::m_surfaceNamesSet
std::vector< bool > m_surfaceNamesSet
Definition: SurfaceIntersectionTest.h:55
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
Trk::TrkDetDescrUnitTestBase::runTest
virtual StatusCode runTest()=0
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
isValid
bool isValid(const T &p)
Definition: AtlasPID.h:225
intersection
std::vector< std::string > intersection(std::vector< std::string > &v1, std::vector< std::string > &v2)
Definition: compareFlatTrees.cxx:25
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
Trk::VectorNtupleBranch::initForWrite
bool initForWrite(TTree &tree, const std::string &varname, int nrow, const std::string &prefix="")
initialize class for writing
Definition: VectorNtupleBranch.cxx:10
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
Trk::VectorNtupleBranch::fill
bool fill(const Amg::VectorX &matrix)
fill a vector
Definition: VectorNtupleBranch.cxx:34
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
Amg::toString
std::string toString(const Translation3D &translation, int precision=4)
GeoPrimitvesToStringConverter.
Definition: GeoPrimitivesToStringConverter.h:40
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
Trk::TrkDetDescrUnitTestBase::finalize
StatusCode finalize()
standard Athena-Algorithm method
Definition: TrkDetDescrUnitTestBase.cxx:53
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
Trk::SurfaceIntersectionTest::m_surfaceSuccessful
std::vector< size_t > m_surfaceSuccessful
Definition: SurfaceIntersectionTest.h:53
Trk::SurfaceIntersectionTest::m_surfaceAttempts
std::vector< size_t > m_surfaceAttempts
Definition: SurfaceIntersectionTest.h:52
dumpFileToPlots.treeName
string treeName
Definition: dumpFileToPlots.py:20
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
Trk::SurfaceIntersectionTest::m_testsPerSurface
size_t m_testsPerSurface
Definition: SurfaceIntersectionTest.h:48
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
Trk::SurfaceNtupleBranch::initForWrite
void initForWrite(TTree &tree, const std::string &prefix="")
initialize class for writing
Definition: SurfaceNtupleBranch.cxx:11
Trk::SurfaceNtupleBranch
This class provides a simple interface to write Surfaces to a root tree.
Definition: SurfaceNtupleBranch.h:36
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
Trk::SurfaceIntersectionTest::m_writeTTree
bool m_writeTTree
Definition: SurfaceIntersectionTest.h:57
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
a
TList * a
Definition: liststreamerinfos.cxx:10
h
mapkey::sf
@ sf
Definition: TElectronEfficiencyCorrectionTool.cxx:38
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Trk::SurfaceIntersectionTest::m_surfaceBuilder
ToolHandle< ISurfaceBuilder > m_surfaceBuilder
Definition: SurfaceIntersectionTest.h:46
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
Trk::SurfaceIntersectionTest::m_surfaces
std::vector< size_t > m_surfaces
Definition: SurfaceIntersectionTest.h:51
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
Trk::TrkDetDescrUnitTestBase::initializeTest
virtual StatusCode initializeTest()
Definition: TrkDetDescrUnitTestBase.cxx:62
Trk::Surface
Definition: Tracking/TrkDetDescr/TrkSurfaces/TrkSurfaces/Surface.h:75
Trk::SurfaceIntersectionTest::m_treeFolder
std::string m_treeFolder
Definition: SurfaceIntersectionTest.h:58
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
Trk::SurfaceIntersectionTest::m_surfaceNames
std::vector< std::string > m_surfaceNames
Definition: SurfaceIntersectionTest.h:54
fitman.k
k
Definition: fitman.py:528
ServiceHandle< ICondSvc >