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

#include <CETmaterial.h>

Inheritance diagram for Trk::CETmaterial:
Collaboration diagram for Trk::CETmaterial:

Public Member Functions

 CETmaterial (const std::string &name, ISvcLocator *pSvcLocator)
 Standard Athena-Algorithm Constructor. More...
 
 ~CETmaterial ()
 Default Destructor. More...
 
StatusCode initialize ()
 standard Athena-Algorithm method More...
 
StatusCode execute ()
 standard Athena-Algorithm method More...
 
StatusCode finalize ()
 standard Athena-Algorithm method More...
 
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...
 

Private Types

typedef ServiceHandle< IChronoStatSvc > IChronoStatSvc_t
 
typedef ServiceHandle< StoreGateSvc > StoreGateSvc_t
 

Private Member Functions

void printMat (double th, double ph, double mat, double dtheta=0., double dphi=0.) const
 
void printMatScan (double theta, double phi, double r, double z, double mat, const std::string &name) const
 
void printMatPrec (double theta, double phi, const Trk::TrackParameters *, const Trk::TrackParameters *, double mat, int id, const std::string &name)
 
void printMatComp (double theta, double phi, const Trk::TrackParameters *currPar, const std::string &name, double mat, double matApp, double dx, double dy) 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

ToolHandle< IExtrapolator > m_extrapolator
 The Extrapolator(s) to be retrieved. More...
 
ToolHandle< IExtrapolator > m_extraprec
 
Rndm::Numbers * m_gaussDist
 Random Number setup. More...
 
Rndm::Numbers * m_flatDist
 
double m_sigmaD0
 Sigma of distribution for D0. More...
 
double m_minZ0
 min Z0 More...
 
double m_maxZ0
 max Z0 More...
 
double m_minEta
 Minimal eta value. More...
 
double m_maxEta
 Maximal eta value. More...
 
double m_minTheta
 Minimal theta value. More...
 
double m_maxTheta
 Maximal theta value. More...
 
double m_minPhi
 Minimal phi value. More...
 
double m_maxPhi
 Maximal phi value. More...
 
double m_minP
 Minimal p value. More...
 
double m_maxP
 Maximal p value. More...
 
double m_charge
 charge More...
 
unsigned int m_numScan
 number of scan tracks More...
 
bool m_checkStepWise
 
bool m_printMaterial
 
bool m_printActive
 
const char * m_matTotFile
 
const char * m_matScanFile
 
const char * m_matActiveFile
 
const char * m_matCompFile
 
bool m_backward
 
bool m_domsentry
 
bool m_doprecision
 
double m_th
 
double m_ph
 
int m_id
 
double m_matSaved
 
Trk::TrackParameters * m_next
 
Amg::MatrixX * m_err
 
const Trk::Surface * m_outerBoundary
 
const Trk::TrackingGeometry * m_trackingGeometry
 
const Trk::TrackingVolume * m_msentry
 
int m_particleType
 the particle typre for the extrap. More...
 
IChronoStatSvc_t m_chronoStatSvc
 
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

The ExtrapolatorTest Algorithm runs a number of n test extrapolations from randomly distributed Track Parameters to geometry outer boundary and back to perigee

Author
Sarka Todorova sarka.nosp@m..tod.nosp@m.orova.nosp@m.@cer.nosp@m.n.ch

Definition at line 38 of file CETmaterial.h.

Member Typedef Documentation

◆ IChronoStatSvc_t

typedef ServiceHandle<IChronoStatSvc> Trk::CETmaterial::IChronoStatSvc_t
private

Definition at line 105 of file CETmaterial.h.

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ CETmaterial()

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

Standard Athena-Algorithm Constructor.

Definition at line 26 of file CETmaterial.cxx.

27  :
28  AthAlgorithm(name,pSvcLocator),
29  m_extrapolator("Trk::Extrapolator/AtlasExtrapolator"),
30  m_extraprec("Trk::Extrapolator/MuonExtrapolator"),
31  m_gaussDist(nullptr),
32  m_flatDist(nullptr),
33  m_sigmaD0(17.*Gaudi::Units::micrometer),
34  m_minZ0(0.),
35  m_maxZ0(0.),
36  m_minEta(-3.),
37  m_maxEta(3.),
38  m_minTheta(0.),
39  m_maxTheta(M_PI),
40  m_minPhi(-M_PI),
41  m_maxPhi( M_PI),
42  m_minP(50000*Gaudi::Units::GeV),
43  m_maxP(50000*Gaudi::Units::GeV),
44  m_charge(1.),
45  m_numScan(10),
46  m_checkStepWise(false),
47  m_printMaterial(false),
48  m_printActive(false),
49  m_matTotFile("material.txt"),
50  m_matScanFile("material_scan.txt"),
51  m_matActiveFile("mat_active.txt"),
52  m_matCompFile("material_comp.txt"),
53  m_backward(false),
54  m_domsentry(false),
55  m_doprecision(false),
56  m_th(0.),
57  m_ph(0.),
58  m_id(0),
59  m_matSaved(0.),
60  m_next(nullptr),
61  m_err(nullptr),
62  m_outerBoundary(nullptr),
63  m_trackingGeometry(nullptr),
64  m_msentry(nullptr),
65  m_particleType(Trk::muon),
66  m_chronoStatSvc( "ChronoStatSvc", name )
67 {
68  // used algorithms and alg tools
69  declareProperty("Extrapolator" , m_extrapolator);
70  declareProperty("ExtraExtrapolator" , m_extraprec);
71 
72  // algorithm steering
73  declareProperty("StartPerigeeSigmaD0" , m_sigmaD0);
74  declareProperty("StartPerigeeMinZ0" , m_minZ0);
75  declareProperty("StartPerigeeMaxZ0" , m_maxZ0);
76  declareProperty("StartPerigeeMinEta" , m_minEta);
77  declareProperty("StartPerigeeMaxEta" , m_maxEta);
78  declareProperty("StartPerigeeMinTheta" , m_minTheta);
79  declareProperty("StartPerigeeMaxTheta" , m_maxTheta);
80  declareProperty("StartPerigeeMinPhi" , m_minPhi);
81  declareProperty("StartPerigeeMaxPhi" , m_maxPhi);
82  declareProperty("StartPerigeeMinP" , m_minP);
83  declareProperty("StartPerigeeMaxP" , m_maxP);
84  declareProperty("StartPerigeeCharge" , m_charge);
85  declareProperty("ParticleType" , m_particleType);
86  declareProperty("NumberOfScanTracks" , m_numScan);
87  declareProperty("CheckActiveLayers" , m_checkStepWise);
88  declareProperty("PrintMaterial" , m_printMaterial);
89  declareProperty("PrintActivePos" , m_printActive);
90  //declareProperty("TotalMaterialOutput" , m_matTotFile);
91  //declareProperty("ScanMaterialOutput" , m_matScanFile);
92  declareProperty("CheckBackward" , m_backward);
93  declareProperty("CheckMSentry" , m_domsentry);
94  declareProperty("CheckPrecision" , m_doprecision);
95 
96 }

◆ ~CETmaterial()

Trk::CETmaterial::~CETmaterial ( )

Default Destructor.

Definition at line 100 of file CETmaterial.cxx.

101 {
102  delete m_gaussDist;
103  delete m_flatDist;
104  delete m_err;
105 }

Member Function Documentation

◆ 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;
147  return AthCommonDataStore<PBASE>::declareGaudiProperty(t, 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::CETmaterial::execute ( )

standard Athena-Algorithm method

Definition at line 150 of file CETmaterial.cxx.

151 {
152  ATH_MSG_INFO( "execute()" );
153  const EventContext& ctx = Gaudi::Hive::currentContext();
154  // retrieve outer boundary cylinder surface
155  if (!m_outerBoundary) {
156  m_trackingGeometry = m_extrapolator->trackingGeometry();
157  m_outerBoundary = &(m_trackingGeometry->highestTrackingVolume()->boundarySurfaces()[2]->surfaceRepresentation());
158  if (!m_outerBoundary) {
159  ATH_MSG_FATAL( "Could not retrieve cylinder boundary from " << m_extrapolator << ". Exiting." );
160  return StatusCode::FAILURE;
161  }
162  ATH_MSG_INFO( "boundary retrieved " );
163  }
164 
165  // generate with random number generator
166  //double d0 = m_gaussDist->shoot() * m_sigmaD0;
167  //double z0 = m_minZ0 + m_flatDist->shoot() * (m_maxZ0-m_minZ0);
168  //double phi = 2.*M_PI * m_flatDist->shoot() - M_PI;
169  //double eta = m_minEta + m_flatDist->shoot()*(m_maxEta-m_minEta);
170  //double ctheta = -1. + 2* m_flatDist->shoot();
171  //double theta = acos(ctheta);
172  //double p = m_minP + m_flatDist->shoot()*(m_maxP-m_minP);
173  //double charge = (m_flatDist->shoot() > 0.5 ) ? -1. : 1.;
174  //double qOverP = charge/(p);
175 
176  if (m_chronoStatSvc) m_chronoStatSvc->chronoStart("MS::scan");
177 
178  // scan
179  std::vector<const TrackStateOnSurface*> material;
180  std::vector<const TrackStateOnSurface*> matPrec;
181  double phi = -M_PI;
182  for ( unsigned int it = 0; it < m_numScan+1; it++) {
183  // the initial perigee
184  double z0 = m_minZ0 + (m_maxZ0-m_minZ0)/m_numScan *it;
185  phi += 1*Gaudi::Units::deg; if (phi>M_PI) phi -=2*M_PI;
186 
187  double theta = m_minTheta + (m_maxTheta-m_minTheta)/m_numScan*it;
188  double p = m_minP + (m_maxP-m_minP)/m_numScan *it;
189  Trk::PerigeeSurface surface( Amg::Vector3D(0.,0.,0.));
190  Trk::Perigee initialPerigee(0., z0, phi, theta, m_charge/p, surface);
191 
192  const Trk::TrackParameters* seed = initialPerigee.clone();
193  const Trk::PerigeeSurface& pSf = initialPerigee.associatedSurface();
194  material.clear(); matPrec.clear();
195  const Trk::TrackParameters* currPar = seed;
196  const Trk::TrackParameters* precPar = seed;
197  if (m_domsentry) {
198  if (!m_msentry) {
199  m_msentry = m_trackingGeometry->trackingVolume("Calo::Containers::Calorimeter");
200  }
201  if (m_msentry) {
202  const Trk::TrackParameters* msEntry =
203  m_extrapolator->extrapolateToVolume(
204  ctx,
205  *currPar,
206  *m_msentry,
207  Trk::alongMomentum,
208  (Trk::ParticleHypothesis)m_particleType).release();
209  if (msEntry) {
210  printMat(
211  theta,
212  phi,
213  currPar->momentum().mag() - msEntry->momentum().mag(),
214  Amg::error(msEntry->covariance()->inverse().eval(), Trk::theta),
215  Amg::error(msEntry->covariance()->inverse().eval(), Trk::phi));
216 
217  const std::vector<const Trk::TrackStateOnSurface*>* mmsentry = m_extrapolator->extrapolateM(ctx,
218  *currPar,
219  msEntry->associatedSurface(),
220  Trk::alongMomentum,
221  false,
222  (Trk::ParticleHypothesis)m_particleType);
223  if (mmsentry ) {
224  for (unsigned int i=0; i< mmsentry->size(); i++)
225  if ((*mmsentry)[i])
226  ATH_MSG_DEBUG( "position:eloss:" << i << "," << (*mmsentry)[i]->trackParameters()->position() << ":"
227  << (*mmsentry)[i]->trackParameters()->momentum().mag()-currPar->momentum().mag());
228  currPar = (mmsentry->back()) ? mmsentry->back()->trackParameters() : msEntry;
229 
230  const std::vector<const Trk::TrackStateOnSurface*>* peri = m_extrapolator->extrapolateM(ctx,
231  *currPar,
232  pSf,
233  Trk::oppositeMomentum,
234  false,
235  (Trk::ParticleHypothesis)m_particleType);
236  ATH_MSG_INFO ( "material scan:backward:" );
237  if (peri){
238  ATH_MSG_DEBUG ("trPar vector size:" << peri->size() );
239  } else {
240  ATH_MSG_ERROR ("Perigee pointer is null in CETmaterial.cxx");
241  delete msEntry;
242  return StatusCode::FAILURE;
243  }
244  for (unsigned int i=0; i< peri->size(); i++)
245  if ((*peri)[i] && (*peri)[i]->trackParameters())
246  ATH_MSG_DEBUG( "position:eloss:" << i << "," << (*peri)[i]->trackParameters()->position() << ":"
247  << (*peri)[i]->trackParameters()->momentum().mag()-msEntry->momentum().mag() );
248 
249  if (peri->back() && peri->back()->trackParameters()) {
250  ATH_MSG_INFO( "extrapolation to perigee:input: "
251  << initialPerigee.parameters()[0] << ","
252  << initialPerigee.parameters()[1] << ","
253  << initialPerigee.parameters()[2] << ","
254  << initialPerigee.parameters()[3] << ","
255  << initialPerigee.momentum().mag() );
256  ATH_MSG_INFO( "extrapolation to perigee:output: "
257  << peri->back()->trackParameters()->parameters()[0] << ","
258  << peri->back()->trackParameters()->parameters()[1] << ","
259  << peri->back()->trackParameters()->parameters()[2] << ","
260  << peri->back()->trackParameters()->parameters()[3] << ","
261  << peri->back()->trackParameters()->momentum().mag() );
262  } else {
263  ATH_MSG_ERROR( "extrapolation to perigee failed for input parameters: " << msEntry->parameters() );
264  }
265  delete peri;
266  delete msEntry;
267  } else {
268  ATH_MSG_ERROR( "extrapolation to MSentry failed for input parameters: " << currPar->parameters() );
269  printMat(theta,phi,0.);
270  }
271  }
272  }
273  delete currPar;
274  continue;
275  }
276  if (m_checkStepWise) {
277  double matApp = 0.;
278  while (currPar) {
279  std::pair<std::unique_ptr<Trk::TrackParameters>,const Trk::Layer*> next = m_extrapolator->extrapolateToNextActiveLayerM(
280  ctx,
281  *currPar,
282  Trk::alongMomentum,
283  true,
284  material,
285  (Trk::ParticleHypothesis)m_particleType);
286 
287 
288  const Trk::TrackParameters* nextPar = next.first.release();
289  const Trk::Layer* lay = next.second;
290  currPar = nextPar;
291 
292  if (m_doprecision && precPar && currPar ) {
293  // try to extrapolate to the same surface
294  const std::vector<const Trk::TrackStateOnSurface*>* nextPrec = m_extraprec->extrapolateM(
295  ctx,
296  *precPar,currPar->associatedSurface(),
297  Trk::alongMomentum,
298  false,
299  (Trk::ParticleHypothesis)m_particleType);
300  delete precPar;
301  // collect material
302  if (nextPrec) {
303  for (const auto *i : *nextPrec) {
304  const Trk::MaterialEffectsBase* mEff = i->materialEffectsOnTrack();
305  const Trk::TrackParameters* trPar = i->trackParameters();
306  if (mEff && trPar) {
307  matApp += mEff->thicknessInX0();
308  }
309  }
310  }
311  // stop of extrapolation failed
312  if (!lay || !nextPrec || nextPrec->empty() || !nextPrec->back() ) break;
313  precPar = nextPrec->back()->trackParameters();
314  double mat=0.;
315  if (!material.empty()) for (auto & i : material) {
316  if (i->materialEffectsOnTrack()) mat += i->materialEffectsOnTrack()->thicknessInX0();
317  }
318  if ( precPar ) printMatComp(theta,phi,currPar,lay->enclosingDetachedTrackingVolume()->name(),mat,matApp,currPar->parameters()[0]-precPar->parameters()[0],
319  currPar->parameters()[1]-precPar->parameters()[1]);
320  else if (currPar) {
321  //precPar is nullptr here
322  ATH_MSG_INFO( "expected layer not reached:" << currPar->position() );
323  }
324  }
325  if (nextPar && m_printActive) {
326  int id = 0;
327  if (lay) id = lay->layerType();
328  double matc=0.;
329  if (!material.empty()) for (auto & i : material) {
330  if (i->materialEffectsOnTrack()) matc += i->materialEffectsOnTrack()->thicknessInX0();
331  }
332  else ATH_MSG_INFO( "mat & error:" << theta << "," << phi << "," << matc << ","
333  << Amg::error(nextPar->covariance()->inverse().eval(),Trk::theta) << ","
334  << Amg::error(nextPar->covariance()->inverse().eval(),Trk::phi) );
335 
336  printMatPrec(theta,phi,nextPar,nextPar,matc,id,"unknown");
337  }
338  if (!lay) break;
339  }
340  if (m_printMaterial) {
341  double mat=0.;
342  if (!material.empty()) for (auto & i : material) {
343  if (i->materialEffectsOnTrack()) {
344  mat += i->materialEffectsOnTrack()->thicknessInX0();
345  }
346  }
347  printMat(theta,phi,mat);
348  }
349  } else {
350  const std::vector<const Trk::TrackStateOnSurface*>* destParameters = m_extrapolator->extrapolateM(
351  ctx,
352  *currPar,
353  *m_outerBoundary,
354  Trk::alongMomentum,
355  false,
356  (Trk::ParticleHypothesis)m_particleType);
357 
358  if (m_printMaterial) {
359  double mat=0.;
360  if (destParameters) for (const auto *destParameter : *destParameters) {
361  const Trk::MaterialEffectsBase* mEff = destParameter->materialEffectsOnTrack();
362  const Trk::TrackParameters* trPar = destParameter->trackParameters();
363  if (trPar) {
364  //const Trk::MeasuredTrackParameters* mdest = dynamic_cast<const Trk::MeasuredTrackParameters*> (trPar);
365  //if (mdest) ATH_MSG_INFO( "radiation thickness and errors(theta,phi):" << theta << "," << phi << "," << mat << "," <<
366  // mdest->localErrorMatrix().error(Trk::theta) << "," << mdest->localErrorMatrix().error(Trk::phi) );
367  }
368  if (mEff && trPar) {
369  mat += mEff->thicknessInX0();
370  // find volume
371  std::vector<const Trk::DetachedTrackingVolume*> detVols = m_extrapolator->trackingGeometry()->lowestDetachedTrackingVolumes(trPar->position());
372  if (!detVols.empty()) printMatScan(theta,phi,trPar->position().perp(),trPar->position().z(),mEff->thicknessInX0(),(detVols)[0]->name());
373  else printMatScan(theta,phi,trPar->position().perp(),trPar->position().z(),mEff->thicknessInX0(),m_extrapolator->trackingGeometry()->lowestStaticTrackingVolume(trPar->position())->volumeName());
374  }
375  }
376  if (destParameters) {
377  //const Trk::MeasuredTrackParameters* mdest = dynamic_cast<const Trk::MeasuredTrackParameters*> ((*destParameters).back()->trackParameters());
378  //if (mdest) ATH_MSG_INFO( "radiation thickness and errors(theta,phi):" << theta << "," << phi << "," << mat << "," <<
379  //mdest->localErrorMatrix().error(Trk::theta) << "," << mdest->localErrorMatrix().error(Trk::phi) );
380  }
381  printMat(theta,phi,mat);
382  }
383 
384  if (!destParameters || destParameters->empty() ) {
385  ATH_MSG_ERROR( "extrapolation to outer boundary failed for input parameters: " << initialPerigee.parameters() );
386  } else if (destParameters->back()->trackParameters()) {
387  // forward extrapolation ok
388  //ATH_MSG_INFO( "last material reached at: " << destParameters->back()->trackParameters()->position().perp() <<","<<destParameters->back()->position().z() );
389 
390  if (m_backward) {
391  material.clear();
392  const std::vector<const Trk::TrackStateOnSurface*>* peri = m_extrapolator->extrapolateM(
393  ctx, *(destParameters->back()->trackParameters()),
394  pSf,
395  Trk::oppositeMomentum,
396  false,
397  (Trk::ParticleHypothesis)m_particleType);
398 
399  if (peri) {
400  ATH_MSG_INFO( "trPar vector size:" << peri->size() );
401  for (unsigned int i=0; i< peri->size(); i++)
402  ATH_MSG_INFO( "position:" << i << "," << (*peri)[i]->trackParameters()->position() );
403  ATH_MSG_INFO( "extrapolation to perigee:input: " << initialPerigee.parameters() );
404  ATH_MSG_INFO( "extrapolation to perigee:output: " << peri->back()->trackParameters()->parameters() );
405  } else {
406  ATH_MSG_ERROR( "extrapolation to perigee failed for input parameters: " << destParameters->back()->trackParameters()->parameters() );
407  }
408  delete peri;
409  }
410  }
411 
412  delete destParameters;
413  }
414  }
415 
416  if (m_chronoStatSvc) m_chronoStatSvc->chronoStop("MS::scan");
417 
418  return StatusCode::SUCCESS;
419 }

◆ 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()) {
55  return m_extendedExtraObjects;
56  }
57  return Algorithm::extraOutputDeps();
58 }

◆ finalize()

StatusCode Trk::CETmaterial::finalize ( )

standard Athena-Algorithm method

Definition at line 139 of file CETmaterial.cxx.

140 {
141  // Code entered here will be executed once at the end of the program run.
142 
143  if (m_chronoStatSvc) m_chronoStatSvc->chronoPrint("MS::scan");
144 
145  return StatusCode::SUCCESS;
146 }

◆ initialize()

StatusCode Trk::CETmaterial::initialize ( )

standard Athena-Algorithm method

Definition at line 110 of file CETmaterial.cxx.

111 {
112  // Code entered here will be executed once at program start.
113 
114  ATH_MSG_INFO( "initialize()" );
115 
116  // Get Extrapolator from ToolService
117  if (m_extrapolator.retrieve().isFailure()) {
118  ATH_MSG_FATAL( "Could not retrieve Tool " << m_extrapolator << ". Exiting." );
119  return StatusCode::FAILURE;
120  }
121  if (m_extraprec.retrieve().isFailure()) {
122  ATH_MSG_FATAL( "Could not retrieve Tool " << m_extraprec << ". Exiting." );
123  return StatusCode::FAILURE;
124  }
125 
126  //if (m_chronoStatSvc.retrieve().isFailure()) {
127  // ATH_MSG_WARNING << "Could not retrieve Tool " << m_chronoStatSvc << ". Exiting." );
128  //}
129 
130  m_gaussDist = new Rndm::Numbers(randSvc(), Rndm::Gauss(0.,1.));
131  m_flatDist = new Rndm::Numbers(randSvc(), Rndm::Flat(0.,1.));
132 
133  ATH_MSG_INFO( "initialize() successful" );
134  return StatusCode::SUCCESS;
135 }

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

◆ printMat()

void Trk::CETmaterial::printMat ( double  th,
double  ph,
double  mat,
double  dtheta = 0.,
double  dphi = 0. 
) const
private

Definition at line 423 of file CETmaterial.cxx.

423  {
424 
425  std::ofstream myfilemat;
426  myfilemat.open (m_matTotFile,std::ios::app);
427  myfilemat<<theta<<" "<<phi<<" "<<mat<<" "<<dtheta<<" "<<dphi<<std::endl;
428 }

◆ printMatComp()

void Trk::CETmaterial::printMatComp ( double  theta,
double  phi,
const Trk::TrackParameters *  currPar,
const std::string &  name,
double  mat,
double  matApp,
double  dx,
double  dy 
) const
private

Definition at line 504 of file CETmaterial.cxx.

505 {
506  std::ofstream myfilemat;
507  myfilemat.open(m_matCompFile,std::ios::app);
508  myfilemat << theta << " " << phi << " " << currPar->position().perp() << " " << currPar->position().z() << " " << name.substr(0,2)
509  << " " << mat << " " << matApp << " " << dx << " " << dy << std::endl;
510 }

◆ printMatPrec()

void Trk::CETmaterial::printMatPrec ( double  theta,
double  phi,
const Trk::TrackParameters *  nextPar,
const Trk::TrackParameters *  mdest,
double  mat,
int  id,
const std::string &  name 
)
private

Definition at line 438 of file CETmaterial.cxx.

438  {
439 
440  if (name.empty()) {}; // dummy to get rid of warning message (unused variable name)
441  std::ofstream myfilemat;
442  myfilemat.open(m_matActiveFile,std::ios::app);
443 
444  if (!m_th && !m_ph) {
445  m_th = theta;
446  m_ph = phi;
447  m_id = id;
448  m_matSaved = mat;
449  delete m_next; m_next=nextPar->clone();
450  delete m_err;
451  m_err=nullptr;
452  if (mdest) {
453  m_err=new Amg::MatrixX;
454  *m_err = mdest->covariance()->inverse().eval();
455  }
456  return;
457  }
458 
459  if ( theta!=m_th || phi!=m_ph ) {
460 
461  if (m_err && m_id>0) {
462  myfilemat << m_th << " " << m_ph << " " << 1 << " " << m_id << " " << m_matSaved << std::endl;
463  myfilemat << m_next->parameters()[Trk::locX] << " " << m_next->parameters()[Trk::locY] << " " << m_next->parameters()[Trk::phi]
464  << " " << m_next->parameters()[Trk::theta] << " " << m_next->parameters()[Trk::qOverP] << std::endl;
465  myfilemat << Amg::error(*m_err,Trk::locX) << " " << Amg::error(*m_err,Trk::locY)
466  << " " << Amg::error(*m_err,Trk::phi) << " " << Amg::error(*m_err,Trk::theta)
467  << " " << Amg::error(*m_err,Trk::qOverP) << std::endl;
468  } else {
469  myfilemat << m_th << " " << m_ph << " " << 0 << " " << m_id << std::endl;
470  myfilemat << m_next->parameters()[Trk::locX] << " " << m_next->parameters()[Trk::locY] << " " << m_next->parameters()[Trk::phi]
471  << " " << m_next->parameters()[Trk::theta] << " " << m_next->parameters()[Trk::qOverP] << std::endl;
472  }
473  m_th = theta;
474  m_ph = phi;
475  m_id = id;
476  m_matSaved = mat;
477  delete m_next; m_next=nextPar->clone();
478  delete m_err;
479  m_err=nullptr;
480 
481  if (mdest) {
482  m_err=new Amg::MatrixX;
483  *m_err = mdest->covariance()->inverse().eval();
484  }
485  return;
486  }
487 
488  // update data
489  if (id>1) {
490  m_th = theta;
491  m_ph = phi;
492  m_id = id;
493  m_matSaved = mat;
494  delete m_next; m_next=nextPar->clone();
495  delete m_err;
496  m_err=nullptr;
497  if (mdest) {
498  m_err=new Amg::MatrixX;
499  *m_err = mdest->covariance()->inverse().eval();
500  }
501  }
502  }

◆ printMatScan()

void Trk::CETmaterial::printMatScan ( double  theta,
double  phi,
double  r,
double  z,
double  mat,
const std::string &  name 
) const
private

Definition at line 431 of file CETmaterial.cxx.

431  {
432 
433  std::ofstream myfilemat;
434  myfilemat.open(m_matScanFile,std::ios::app);
435  myfilemat << theta << " " << phi << " " << r << " " << z << " " << mat << " " << name << std::endl;
436 }

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

◆ 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  {
67  StatusCode sc=AthCommonDataStore<AthCommonMsg<Algorithm>>::sysInitialize();
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_backward

bool Trk::CETmaterial::m_backward
private

Definition at line 89 of file CETmaterial.h.

◆ m_charge

double Trk::CETmaterial::m_charge
private

charge

Definition at line 80 of file CETmaterial.h.

◆ m_checkStepWise

bool Trk::CETmaterial::m_checkStepWise
private

Definition at line 82 of file CETmaterial.h.

◆ m_chronoStatSvc

IChronoStatSvc_t Trk::CETmaterial::m_chronoStatSvc
private

Definition at line 106 of file CETmaterial.h.

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

bool Trk::CETmaterial::m_domsentry
private

Definition at line 90 of file CETmaterial.h.

◆ m_doprecision

bool Trk::CETmaterial::m_doprecision
private

Definition at line 91 of file CETmaterial.h.

◆ m_err

Amg::MatrixX* Trk::CETmaterial::m_err
private

Definition at line 97 of file CETmaterial.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_extrapolator

ToolHandle<IExtrapolator> Trk::CETmaterial::m_extrapolator
private

The Extrapolator(s) to be retrieved.

Definition at line 62 of file CETmaterial.h.

◆ m_extraprec

ToolHandle<IExtrapolator> Trk::CETmaterial::m_extraprec
private

Definition at line 63 of file CETmaterial.h.

◆ m_flatDist

Rndm::Numbers* Trk::CETmaterial::m_flatDist
private

Definition at line 67 of file CETmaterial.h.

◆ m_gaussDist

Rndm::Numbers* Trk::CETmaterial::m_gaussDist
private

Random Number setup.

Definition at line 66 of file CETmaterial.h.

◆ m_id

int Trk::CETmaterial::m_id
private

Definition at line 94 of file CETmaterial.h.

◆ m_matActiveFile

const char* Trk::CETmaterial::m_matActiveFile
private

Definition at line 87 of file CETmaterial.h.

◆ m_matCompFile

const char* Trk::CETmaterial::m_matCompFile
private

Definition at line 88 of file CETmaterial.h.

◆ m_matSaved

double Trk::CETmaterial::m_matSaved
private

Definition at line 95 of file CETmaterial.h.

◆ m_matScanFile

const char* Trk::CETmaterial::m_matScanFile
private

Definition at line 86 of file CETmaterial.h.

◆ m_matTotFile

const char* Trk::CETmaterial::m_matTotFile
private

Definition at line 85 of file CETmaterial.h.

◆ m_maxEta

double Trk::CETmaterial::m_maxEta
private

Maximal eta value.

Definition at line 73 of file CETmaterial.h.

◆ m_maxP

double Trk::CETmaterial::m_maxP
private

Maximal p value.

Definition at line 79 of file CETmaterial.h.

◆ m_maxPhi

double Trk::CETmaterial::m_maxPhi
private

Maximal phi value.

Definition at line 77 of file CETmaterial.h.

◆ m_maxTheta

double Trk::CETmaterial::m_maxTheta
private

Maximal theta value.

Definition at line 75 of file CETmaterial.h.

◆ m_maxZ0

double Trk::CETmaterial::m_maxZ0
private

max Z0

Definition at line 71 of file CETmaterial.h.

◆ m_minEta

double Trk::CETmaterial::m_minEta
private

Minimal eta value.

Definition at line 72 of file CETmaterial.h.

◆ m_minP

double Trk::CETmaterial::m_minP
private

Minimal p value.

Definition at line 78 of file CETmaterial.h.

◆ m_minPhi

double Trk::CETmaterial::m_minPhi
private

Minimal phi value.

Definition at line 76 of file CETmaterial.h.

◆ m_minTheta

double Trk::CETmaterial::m_minTheta
private

Minimal theta value.

Definition at line 74 of file CETmaterial.h.

◆ m_minZ0

double Trk::CETmaterial::m_minZ0
private

min Z0

Definition at line 70 of file CETmaterial.h.

◆ m_msentry

const Trk::TrackingVolume* Trk::CETmaterial::m_msentry
private

Definition at line 101 of file CETmaterial.h.

◆ m_next

Trk::TrackParameters* Trk::CETmaterial::m_next
private

Definition at line 96 of file CETmaterial.h.

◆ m_numScan

unsigned int Trk::CETmaterial::m_numScan
private

number of scan tracks

Definition at line 81 of file CETmaterial.h.

◆ m_outerBoundary

const Trk::Surface* Trk::CETmaterial::m_outerBoundary
private

Definition at line 99 of file CETmaterial.h.

◆ m_particleType

int Trk::CETmaterial::m_particleType
private

the particle typre for the extrap.

Definition at line 103 of file CETmaterial.h.

◆ m_ph

double Trk::CETmaterial::m_ph
private

Definition at line 93 of file CETmaterial.h.

◆ m_printActive

bool Trk::CETmaterial::m_printActive
private

Definition at line 84 of file CETmaterial.h.

◆ m_printMaterial

bool Trk::CETmaterial::m_printMaterial
private

Definition at line 83 of file CETmaterial.h.

◆ m_sigmaD0

double Trk::CETmaterial::m_sigmaD0
private

Sigma of distribution for D0.

Definition at line 69 of file CETmaterial.h.

◆ m_th

double Trk::CETmaterial::m_th
private

Definition at line 92 of file CETmaterial.h.

◆ m_trackingGeometry

const Trk::TrackingGeometry* Trk::CETmaterial::m_trackingGeometry
private

Definition at line 100 of file CETmaterial.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.

The documentation for this class was generated from the following files:
beamspotman.r
def r
Definition: beamspotman.py:676
Trk::CETmaterial::m_id
int m_id
Definition: CETmaterial.h:94
GeV
#define GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:17
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
Trk::CETmaterial::m_ph
double m_ph
Definition: CETmaterial.h:93
Trk::z
@ z
global position (cartesian)
Definition: ParamDefs.h:57
Trk::CETmaterial::m_flatDist
Rndm::Numbers * m_flatDist
Definition: CETmaterial.h:67
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::locX
@ locX
Definition: ParamDefs.h:37
Trk::locY
@ locY
local cartesian
Definition: ParamDefs.h:38
Trk::CETmaterial::m_charge
double m_charge
charge
Definition: CETmaterial.h:80
Trk::CETmaterial::printMatComp
void printMatComp(double theta, double phi, const Trk::TrackParameters *currPar, const std::string &name, double mat, double matApp, double dx, double dy) const
Definition: CETmaterial.cxx:504
Trk::PerigeeSurface
Definition: PerigeeSurface.h:43
Trk::next
@ next
Definition: BinningData.h:33
Trk::ParametersBase::position
const Amg::Vector3D & position() const
Access method for the position.
Trk::oppositeMomentum
@ oppositeMomentum
Definition: PropDirection.h:21
Trk::ParametersBase::associatedSurface
virtual const Surface & associatedSurface() const override=0
Access to the Surface associated to the Parameters.
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
Trk::ParametersT
Dummy class used to allow special convertors to be called for surfaces owned by a detector element.
Definition: EMErrorDetail.h:25
Trk::TrackingVolume::boundarySurfaces
std::vector< SharedObject< BoundarySurface< TrackingVolume > > > & boundarySurfaces()
Method to return the BoundarySurfaces.
Definition: TrackingVolume.cxx:982
mat
GeoMaterial * mat
Definition: LArDetectorConstructionTBEC.cxx:55
Trk::CETmaterial::m_gaussDist
Rndm::Numbers * m_gaussDist
Random Number setup.
Definition: CETmaterial.h:66
skel.it
it
Definition: skel.GENtoEVGEN.py:396
M_PI
#define M_PI
Definition: ActiveFraction.h:11
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::CETmaterial::m_next
Trk::TrackParameters * m_next
Definition: CETmaterial.h:96
deg
#define deg
Definition: SbPolyhedron.cxx:17
Trk::z0
@ z0
Definition: ParamDefs.h:64
Trk::CETmaterial::m_matScanFile
const char * m_matScanFile
Definition: CETmaterial.h:86
Trk::CETmaterial::m_maxTheta
double m_maxTheta
Maximal theta value.
Definition: CETmaterial.h:75
Trk::Layer::enclosingDetachedTrackingVolume
const DetachedTrackingVolume * enclosingDetachedTrackingVolume() const
get the confining DetachedTrackingVolume
Trk::MaterialEffectsBase::thicknessInX0
double thicknessInX0() const
returns the actually traversed material .
Trk::alongMomentum
@ alongMomentum
Definition: PropDirection.h:20
Trk::CETmaterial::m_domsentry
bool m_domsentry
Definition: CETmaterial.h:90
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
Trk::MaterialEffectsBase
base class to integrate material effects on Trk::Track in a flexible way.
Definition: MaterialEffectsBase.h:35
Trk::CETmaterial::m_numScan
unsigned int m_numScan
number of scan tracks
Definition: CETmaterial.h:81
Trk::CETmaterial::printMatPrec
void printMatPrec(double theta, double phi, const Trk::TrackParameters *, const Trk::TrackParameters *, double mat, int id, const std::string &name)
Definition: CETmaterial.cxx:438
Trk::Layer::layerType
int layerType() const
get the Layer coding
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
Trk::CETmaterial::m_minZ0
double m_minZ0
min Z0
Definition: CETmaterial.h:70
Trk::ParticleHypothesis
ParticleHypothesis
Definition: ParticleHypothesis.h:25
Trk::CETmaterial::printMat
void printMat(double th, double ph, double mat, double dtheta=0., double dphi=0.) const
Definition: CETmaterial.cxx:423
Trk::CETmaterial::m_trackingGeometry
const Trk::TrackingGeometry * m_trackingGeometry
Definition: CETmaterial.h:100
Trk::CETmaterial::m_outerBoundary
const Trk::Surface * m_outerBoundary
Definition: CETmaterial.h:99
Trk::TrackingGeometry::highestTrackingVolume
const TrackingVolume * highestTrackingVolume() const
return the world
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
Generate_dsid_ranseed.seed
seed
Definition: Generate_dsid_ranseed.py:10
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::CETmaterial::m_maxP
double m_maxP
Maximal p value.
Definition: CETmaterial.h:79
ParticleGun_EoverP_Config.momentum
momentum
Definition: ParticleGun_EoverP_Config.py:63
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
Trk::CETmaterial::m_matCompFile
const char * m_matCompFile
Definition: CETmaterial.h:88
lumiFormat.i
int i
Definition: lumiFormat.py:85
Trk::CETmaterial::m_doprecision
bool m_doprecision
Definition: CETmaterial.h:91
Trk::theta
@ theta
Definition: ParamDefs.h:66
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
Trk::CETmaterial::printMatScan
void printMatScan(double theta, double phi, double r, double z, double mat, const std::string &name) const
Definition: CETmaterial.cxx:431
Trk::CETmaterial::m_err
Amg::MatrixX * m_err
Definition: CETmaterial.h:97
Trk::CETmaterial::m_minTheta
double m_minTheta
Minimal theta value.
Definition: CETmaterial.h:74
Trk::CETmaterial::m_maxPhi
double m_maxPhi
Maximal phi value.
Definition: CETmaterial.h:77
Trk::CETmaterial::m_chronoStatSvc
IChronoStatSvc_t m_chronoStatSvc
Definition: CETmaterial.h:106
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
Trk::ParametersBase
Definition: ParametersBase.h:55
python.SystemOfUnits.micrometer
int micrometer
Definition: SystemOfUnits.py:71
Trk::CETmaterial::m_msentry
const Trk::TrackingVolume * m_msentry
Definition: CETmaterial.h:101
Trk::muon
@ muon
Definition: ParticleHypothesis.h:28
Trk::CETmaterial::m_backward
bool m_backward
Definition: CETmaterial.h:89
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::CETmaterial::m_th
double m_th
Definition: CETmaterial.h:92
Trk::CETmaterial::m_maxEta
double m_maxEta
Maximal eta value.
Definition: CETmaterial.h:73
Trk::CETmaterial::m_sigmaD0
double m_sigmaD0
Sigma of distribution for D0.
Definition: CETmaterial.h:69
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
id
SG::auxid_t id
Definition: Control/AthContainers/Root/debug.cxx:220
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
Amg::error
double error(const Amg::MatrixX &mat, int index)
return diagonal error of the matrix caller should ensure the matrix is symmetric and the index is in ...
Definition: EventPrimitivesHelpers.h:40
Trk::CETmaterial::m_minEta
double m_minEta
Minimal eta value.
Definition: CETmaterial.h:72
Trk::CETmaterial::m_printMaterial
bool m_printMaterial
Definition: CETmaterial.h:83
Trk::DetachedTrackingVolume::name
const std::string & name() const
returns the Name
Definition: DetachedTrackingVolume.h:134
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
makeTRTBarrelCans.dy
tuple dy
Definition: makeTRTBarrelCans.py:21
Trk::ParametersBase::momentum
const Amg::Vector3D & momentum() const
Access method for the momentum.
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
Trk::CETmaterial::m_matSaved
double m_matSaved
Definition: CETmaterial.h:95
a
TList * a
Definition: liststreamerinfos.cxx:10
h
Trk::CETmaterial::m_maxZ0
double m_maxZ0
max Z0
Definition: CETmaterial.h:71
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Trk::CETmaterial::m_checkStepWise
bool m_checkStepWise
Definition: CETmaterial.h:82
Trk::CETmaterial::m_printActive
bool m_printActive
Definition: CETmaterial.h:84
Trk::CETmaterial::m_minPhi
double m_minPhi
Minimal phi value.
Definition: CETmaterial.h:76
makeTRTBarrelCans.dx
tuple dx
Definition: makeTRTBarrelCans.py:20
Trk::qOverP
@ qOverP
perigee
Definition: ParamDefs.h:67
Trk::CETmaterial::m_extrapolator
ToolHandle< IExtrapolator > m_extrapolator
The Extrapolator(s) to be retrieved.
Definition: CETmaterial.h:62
Trk::TrackingGeometry::trackingVolume
const TrackingVolume * trackingVolume(const std::string &name) const
return the tracking Volume by name, 0 if it doesn't exist
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
Trk::CETmaterial::m_matActiveFile
const char * m_matActiveFile
Definition: CETmaterial.h:87
AthAlgorithm::AthAlgorithm
AthAlgorithm()
Default constructor:
Trk::phi
@ phi
Definition: ParamDefs.h:75
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
Trk::CETmaterial::m_minP
double m_minP
Minimal p value.
Definition: CETmaterial.h:78
Trk::CETmaterial::m_matTotFile
const char * m_matTotFile
Definition: CETmaterial.h:85
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
mag
Scalar mag() const
mag method
Definition: AmgMatrixBasePlugin.h:26
Trk::CETmaterial::m_particleType
int m_particleType
the particle typre for the extrap.
Definition: CETmaterial.h:103
fitman.k
k
Definition: fitman.py:528
Trk::ParametersBase::clone
virtual ParametersBase< DIM, T > * clone() const override=0
clone method for polymorphic deep copy
Trk::Layer
Definition: Layer.h:73
ServiceHandle< ICondSvc >
Trk::CETmaterial::m_extraprec
ToolHandle< IExtrapolator > m_extraprec
Definition: CETmaterial.h:63
Generated on Thu Nov 7 2024 22:28:06 for ATLAS Offline Software by   doxygen 1.8.18