Loading [MathJax]/jax/output/SVG/config.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
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 {this, "Extrapolator", "Trk::Extrapolator/AtlasExtrapolator"}
 The Extrapolator(s) to be retrieved. More...
 
ToolHandle< IExtrapolator > m_extraprec {this, "ExtraExtrapolator", "Trk::Extrapolator/MuonExtrapolator"}
 
DoubleProperty m_minZ0 {this, "StartPerigeeMinZ0", 0.}
 
DoubleProperty m_maxZ0 {this, "StartPerigeeMaxZ0", 0.}
 
DoubleProperty m_minTheta {this, "StartPerigeeMinTheta", 0.}
 
DoubleProperty m_maxTheta {this, "StartPerigeeMaxTheta", M_PI}
 
DoubleProperty m_minP {this, "StartPerigeeMinP", 50000*Gaudi::Units::GeV}
 
DoubleProperty m_maxP {this, "StartPerigeeMaxP", 50000*Gaudi::Units::GeV}
 
DoubleProperty m_charge {this, "StartPerigeeCharge", 1.}
 
UnsignedIntegerProperty m_numScan {this, "NumberOfScanTracks", 10}
 
BooleanProperty m_checkStepWise {this, "CheckActiveLayers", false}
 
BooleanProperty m_printMaterial {this, "PrintMaterial", false}
 
BooleanProperty m_printActive {this, "PrintActivePos", false}
 
const char * m_matTotFile = "material.txt"
 
const char * m_matScanFile = "material_scan.txt"
 
const char * m_matActiveFile = "mat_active.txt"
 
const char * m_matCompFile = "material_comp.txt"
 
BooleanProperty m_backward {this, "CheckBackward", false}
 
BooleanProperty m_domsentry {this, "CheckMSentry", false}
 
BooleanProperty m_doprecision {this, "CheckPrecision", false}
 
double m_th = 0.
 
double m_ph = 0.
 
int m_id = 0
 
double m_matSaved = 0.
 
Trk::TrackParameters * m_next = nullptr
 
Amg::MatrixX * m_err = nullptr
 
const Trk::Surface * m_outerBoundary = nullptr
 
const Trk::TrackingGeometry * m_trackingGeometry = nullptr
 
const Trk::TrackingVolume * m_msentry = nullptr
 
IntegerProperty m_particleType
 
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 39 of file CETmaterial.h.

Member Typedef Documentation

◆ IChronoStatSvc_t

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

Definition at line 103 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 24 of file CETmaterial.cxx.

25  :
26  AthAlgorithm(name,pSvcLocator),
27  m_chronoStatSvc( "ChronoStatSvc", name )
28 { }

◆ ~CETmaterial()

Trk::CETmaterial::~CETmaterial ( )

Default Destructor.

Definition at line 32 of file CETmaterial.cxx.

33 {
34  delete m_err;
35 }

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 73 of file CETmaterial.cxx.

74 {
75  ATH_MSG_INFO( "execute()" );
76  const EventContext& ctx = Gaudi::Hive::currentContext();
77  // retrieve outer boundary cylinder surface
78  if (!m_outerBoundary) {
79  m_trackingGeometry = m_extrapolator->trackingGeometry();
80  m_outerBoundary = &(m_trackingGeometry->highestTrackingVolume()->boundarySurfaces()[2]->surfaceRepresentation());
81  if (!m_outerBoundary) {
82  ATH_MSG_FATAL( "Could not retrieve cylinder boundary from " << m_extrapolator << ". Exiting." );
83  return StatusCode::FAILURE;
84  }
85  ATH_MSG_INFO( "boundary retrieved " );
86  }
87 
88  if (m_chronoStatSvc) m_chronoStatSvc->chronoStart("MS::scan");
89 
90  // scan
91  std::vector<const TrackStateOnSurface*> material;
92  std::vector<const TrackStateOnSurface*> matPrec;
93  double phi = -M_PI;
94  for ( unsigned int it = 0; it < m_numScan+1; it++) {
95  // the initial perigee
96  double z0 = m_minZ0 + (m_maxZ0-m_minZ0)/m_numScan *it;
97  phi += 1*Gaudi::Units::deg; if (phi>M_PI) phi -=2*M_PI;
98 
99  double theta = m_minTheta + (m_maxTheta-m_minTheta)/m_numScan*it;
100  double p = m_minP + (m_maxP-m_minP)/m_numScan *it;
101  Trk::PerigeeSurface surface( Amg::Vector3D(0.,0.,0.));
102  Trk::Perigee initialPerigee(0., z0, phi, theta, m_charge/p, surface);
103 
104  const Trk::TrackParameters* seed = initialPerigee.clone();
105  const Trk::PerigeeSurface& pSf = initialPerigee.associatedSurface();
106  material.clear(); matPrec.clear();
107  const Trk::TrackParameters* currPar = seed;
108  const Trk::TrackParameters* precPar = seed;
109  if (m_domsentry) {
110  if (!m_msentry) {
111  m_msentry = m_trackingGeometry->trackingVolume("Calo::Containers::Calorimeter");
112  }
113  if (m_msentry) {
114  const Trk::TrackParameters* msEntry =
115  m_extrapolator->extrapolateToVolume(
116  ctx,
117  *currPar,
118  *m_msentry,
119  Trk::alongMomentum,
120  static_cast<Trk::ParticleHypothesis>(m_particleType.value())).release();
121  if (msEntry) {
122  printMat(
123  theta,
124  phi,
125  currPar->momentum().mag() - msEntry->momentum().mag(),
126  Amg::error(msEntry->covariance()->inverse().eval(), Trk::theta),
127  Amg::error(msEntry->covariance()->inverse().eval(), Trk::phi));
128 
129  const std::vector<const Trk::TrackStateOnSurface*>* mmsentry = m_extrapolator->extrapolateM(ctx,
130  *currPar,
131  msEntry->associatedSurface(),
132  Trk::alongMomentum,
133  false,
134  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
135  if (mmsentry ) {
136  for (const auto& entry : *mmsentry) {
137  if (entry) {
138  ATH_MSG_DEBUG("position:eloss:"
139  << entry->trackParameters()->position() << ":"
140  << entry->trackParameters()->momentum().mag() - currPar->momentum().mag());
141  }
142  }
143 
144  currPar = (mmsentry->back()) ? mmsentry->back()->trackParameters() : msEntry;
145 
146  const std::vector<const Trk::TrackStateOnSurface*>* peri = m_extrapolator->extrapolateM(ctx,
147  *currPar,
148  pSf,
149  Trk::oppositeMomentum,
150  false,
151  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
152  ATH_MSG_INFO ( "material scan:backward:" );
153  if (peri){
154  ATH_MSG_DEBUG ("trPar vector size:" << peri->size() );
155  } else {
156  ATH_MSG_ERROR ("Perigee pointer is null in CETmaterial.cxx");
157  delete msEntry;
158  return StatusCode::FAILURE;
159  }
160  for (const auto& entry : *peri) {
161  if (entry && entry->trackParameters()) {
162  ATH_MSG_DEBUG("position:eloss:"
163  << entry->trackParameters()->position() << ":"
164  << entry->trackParameters()->momentum().mag() - msEntry->momentum().mag());
165  }
166  }
167 
168  if (peri->back() && peri->back()->trackParameters()) {
169  ATH_MSG_INFO( "extrapolation to perigee:input: "
170  << initialPerigee.parameters()[0] << ","
171  << initialPerigee.parameters()[1] << ","
172  << initialPerigee.parameters()[2] << ","
173  << initialPerigee.parameters()[3] << ","
174  << initialPerigee.momentum().mag() );
175  ATH_MSG_INFO( "extrapolation to perigee:output: "
176  << peri->back()->trackParameters()->parameters()[0] << ","
177  << peri->back()->trackParameters()->parameters()[1] << ","
178  << peri->back()->trackParameters()->parameters()[2] << ","
179  << peri->back()->trackParameters()->parameters()[3] << ","
180  << peri->back()->trackParameters()->momentum().mag() );
181  } else {
182  ATH_MSG_ERROR( "extrapolation to perigee failed for input parameters: " << msEntry->parameters() );
183  }
184  delete peri;
185  delete msEntry;
186  } else {
187  ATH_MSG_ERROR( "extrapolation to MSentry failed for input parameters: " << currPar->parameters() );
188  printMat(theta,phi,0.);
189  }
190  }
191  }
192  delete currPar;
193  continue;
194  }
195  if (m_checkStepWise) {
196  double matApp = 0.;
197  while (currPar) {
198  std::pair<std::unique_ptr<Trk::TrackParameters>,const Trk::Layer*> next = m_extrapolator->extrapolateToNextActiveLayerM(
199  ctx,
200  *currPar,
201  Trk::alongMomentum,
202  true,
203  material,
204  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
205 
206 
207  const Trk::TrackParameters* nextPar = next.first.release();
208  const Trk::Layer* lay = next.second;
209  currPar = nextPar;
210 
211  if (m_doprecision && precPar && currPar ) {
212  // try to extrapolate to the same surface
213  const std::vector<const Trk::TrackStateOnSurface*>* nextPrec = m_extraprec->extrapolateM(
214  ctx,
215  *precPar,currPar->associatedSurface(),
216  Trk::alongMomentum,
217  false,
218  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
219  delete precPar;
220  // collect material
221  if (nextPrec) {
222  for (const auto *i : *nextPrec) {
223  const Trk::MaterialEffectsBase* mEff = i->materialEffectsOnTrack();
224  const Trk::TrackParameters* trPar = i->trackParameters();
225  if (mEff && trPar) {
226  matApp += mEff->thicknessInX0();
227  }
228  }
229  }
230  // stop of extrapolation failed
231  if (!lay || !nextPrec || nextPrec->empty() || !nextPrec->back() ) break;
232  precPar = nextPrec->back()->trackParameters();
233  double mat=0.;
234  if (!material.empty()) for (auto & i : material) {
235  if (i->materialEffectsOnTrack()) mat += i->materialEffectsOnTrack()->thicknessInX0();
236  }
237  if ( precPar ) printMatComp(theta,phi,currPar,lay->enclosingDetachedTrackingVolume()->name(),mat,matApp,currPar->parameters()[0]-precPar->parameters()[0],
238  currPar->parameters()[1]-precPar->parameters()[1]);
239  else if (currPar) {
240  //precPar is nullptr here
241  ATH_MSG_INFO( "expected layer not reached:" << currPar->position() );
242  }
243  }
244  if (nextPar && m_printActive) {
245  int id = 0;
246  if (lay) id = lay->layerType();
247  double matc=0.;
248  if (!material.empty()) for (auto & i : material) {
249  if (i->materialEffectsOnTrack()) matc += i->materialEffectsOnTrack()->thicknessInX0();
250  }
251  else ATH_MSG_INFO( "mat & error:" << theta << "," << phi << "," << matc << ","
252  << Amg::error(nextPar->covariance()->inverse().eval(),Trk::theta) << ","
253  << Amg::error(nextPar->covariance()->inverse().eval(),Trk::phi) );
254 
255  printMatPrec(theta,phi,nextPar,nextPar,matc,id,"unknown");
256  }
257  if (!lay) break;
258  }
259  if (m_printMaterial) {
260  double mat=0.;
261  if (!material.empty()) for (auto & i : material) {
262  if (i->materialEffectsOnTrack()) {
263  mat += i->materialEffectsOnTrack()->thicknessInX0();
264  }
265  }
266  printMat(theta,phi,mat);
267  }
268  } else {
269  const std::vector<const Trk::TrackStateOnSurface*>* destParameters = m_extrapolator->extrapolateM(
270  ctx,
271  *currPar,
272  *m_outerBoundary,
273  Trk::alongMomentum,
274  false,
275  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
276 
277  if (m_printMaterial) {
278  double mat=0.;
279  if (destParameters) for (const auto *destParameter : *destParameters) {
280  const Trk::MaterialEffectsBase* mEff = destParameter->materialEffectsOnTrack();
281  const Trk::TrackParameters* trPar = destParameter->trackParameters();
282  if (trPar) {
283  //const Trk::MeasuredTrackParameters* mdest = dynamic_cast<const Trk::MeasuredTrackParameters*> (trPar);
284  //if (mdest) ATH_MSG_INFO( "radiation thickness and errors(theta,phi):" << theta << "," << phi << "," << mat << "," <<
285  // mdest->localErrorMatrix().error(Trk::theta) << "," << mdest->localErrorMatrix().error(Trk::phi) );
286  }
287  if (mEff && trPar) {
288  mat += mEff->thicknessInX0();
289  // find volume
290  std::vector<const Trk::DetachedTrackingVolume*> detVols = m_extrapolator->trackingGeometry()->lowestDetachedTrackingVolumes(trPar->position());
291  if (!detVols.empty()) printMatScan(theta,phi,trPar->position().perp(),trPar->position().z(),mEff->thicknessInX0(),(detVols)[0]->name());
292  else printMatScan(theta,phi,trPar->position().perp(),trPar->position().z(),mEff->thicknessInX0(),m_extrapolator->trackingGeometry()->lowestStaticTrackingVolume(trPar->position())->volumeName());
293  }
294  }
295  if (destParameters) {
296  //const Trk::MeasuredTrackParameters* mdest = dynamic_cast<const Trk::MeasuredTrackParameters*> ((*destParameters).back()->trackParameters());
297  //if (mdest) ATH_MSG_INFO( "radiation thickness and errors(theta,phi):" << theta << "," << phi << "," << mat << "," <<
298  //mdest->localErrorMatrix().error(Trk::theta) << "," << mdest->localErrorMatrix().error(Trk::phi) );
299  }
300  printMat(theta,phi,mat);
301  }
302 
303  if (!destParameters || destParameters->empty() ) {
304  ATH_MSG_ERROR( "extrapolation to outer boundary failed for input parameters: " << initialPerigee.parameters() );
305  } else if (destParameters->back()->trackParameters()) {
306  // forward extrapolation ok
307  if (m_backward) {
308  material.clear();
309  const std::vector<const Trk::TrackStateOnSurface*>* peri = m_extrapolator->extrapolateM(
310  ctx, *(destParameters->back()->trackParameters()),
311  pSf,
312  Trk::oppositeMomentum,
313  false,
314  static_cast<Trk::ParticleHypothesis>(m_particleType.value()));
315 
316  if (peri) {
317  ATH_MSG_INFO( "trPar vector size:" << peri->size() );
318  for (unsigned int i=0; i< peri->size(); i++)
319  ATH_MSG_INFO( "position:" << i << "," << (*peri)[i]->trackParameters()->position() );
320  ATH_MSG_INFO( "extrapolation to perigee:input: " << initialPerigee.parameters() );
321  ATH_MSG_INFO( "extrapolation to perigee:output: " << peri->back()->trackParameters()->parameters() );
322  } else {
323  ATH_MSG_ERROR( "extrapolation to perigee failed for input parameters: " << destParameters->back()->trackParameters()->parameters() );
324  }
325  delete peri;
326  }
327  }
328 
329  delete destParameters;
330  }
331  }
332 
333  if (m_chronoStatSvc) m_chronoStatSvc->chronoStop("MS::scan");
334 
335  return StatusCode::SUCCESS;
336 }

◆ 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 62 of file CETmaterial.cxx.

63 {
64  // Code entered here will be executed once at the end of the program run.
65 
66  if (m_chronoStatSvc) m_chronoStatSvc->chronoPrint("MS::scan");
67 
68  return StatusCode::SUCCESS;
69 }

◆ initialize()

StatusCode Trk::CETmaterial::initialize ( )

standard Athena-Algorithm method

Definition at line 40 of file CETmaterial.cxx.

41 {
42  // Code entered here will be executed once at program start.
43 
44  ATH_MSG_INFO( "initialize()" );
45 
46  // Get Extrapolator from ToolService
47  if (m_extrapolator.retrieve().isFailure()) {
48  ATH_MSG_FATAL( "Could not retrieve Tool " << m_extrapolator << ". Exiting." );
49  return StatusCode::FAILURE;
50  }
51  if (m_extraprec.retrieve().isFailure()) {
52  ATH_MSG_FATAL( "Could not retrieve Tool " << m_extraprec << ". Exiting." );
53  return StatusCode::FAILURE;
54  }
55 
56  ATH_MSG_INFO( "initialize() successful" );
57  return StatusCode::SUCCESS;
58 }

◆ 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 340 of file CETmaterial.cxx.

340  {
341 
342  std::ofstream myfilemat;
343  myfilemat.open (m_matTotFile,std::ios::app);
344  myfilemat<<theta<<" "<<phi<<" "<<mat<<" "<<dtheta<<" "<<dphi<<std::endl;
345 }

◆ 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 421 of file CETmaterial.cxx.

422 {
423  std::ofstream myfilemat;
424  myfilemat.open(m_matCompFile,std::ios::app);
425  myfilemat << theta << " " << phi << " " << currPar->position().perp() << " " << currPar->position().z() << " " << name.substr(0,2)
426  << " " << mat << " " << matApp << " " << dx << " " << dy << std::endl;
427 }

◆ 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 355 of file CETmaterial.cxx.

355  {
356 
357  if (name.empty()) {}; // dummy to get rid of warning message (unused variable name)
358  std::ofstream myfilemat;
359  myfilemat.open(m_matActiveFile,std::ios::app);
360 
361  if (!m_th && !m_ph) {
362  m_th = theta;
363  m_ph = phi;
364  m_id = id;
365  m_matSaved = mat;
366  delete m_next; m_next=nextPar->clone();
367  delete m_err;
368  m_err=nullptr;
369  if (mdest) {
370  m_err=new Amg::MatrixX;
371  *m_err = mdest->covariance()->inverse().eval();
372  }
373  return;
374  }
375 
376  if ( theta!=m_th || phi!=m_ph ) {
377 
378  if (m_err && m_id>0) {
379  myfilemat << m_th << " " << m_ph << " " << 1 << " " << m_id << " " << m_matSaved << std::endl;
380  myfilemat << m_next->parameters()[Trk::locX] << " " << m_next->parameters()[Trk::locY] << " " << m_next->parameters()[Trk::phi]
381  << " " << m_next->parameters()[Trk::theta] << " " << m_next->parameters()[Trk::qOverP] << std::endl;
382  myfilemat << Amg::error(*m_err,Trk::locX) << " " << Amg::error(*m_err,Trk::locY)
383  << " " << Amg::error(*m_err,Trk::phi) << " " << Amg::error(*m_err,Trk::theta)
384  << " " << Amg::error(*m_err,Trk::qOverP) << std::endl;
385  } else {
386  myfilemat << m_th << " " << m_ph << " " << 0 << " " << m_id << std::endl;
387  myfilemat << m_next->parameters()[Trk::locX] << " " << m_next->parameters()[Trk::locY] << " " << m_next->parameters()[Trk::phi]
388  << " " << m_next->parameters()[Trk::theta] << " " << m_next->parameters()[Trk::qOverP] << std::endl;
389  }
390  m_th = theta;
391  m_ph = phi;
392  m_id = id;
393  m_matSaved = mat;
394  delete m_next; m_next=nextPar->clone();
395  delete m_err;
396  m_err=nullptr;
397 
398  if (mdest) {
399  m_err=new Amg::MatrixX;
400  *m_err = mdest->covariance()->inverse().eval();
401  }
402  return;
403  }
404 
405  // update data
406  if (id>1) {
407  m_th = theta;
408  m_ph = phi;
409  m_id = id;
410  m_matSaved = mat;
411  delete m_next; m_next=nextPar->clone();
412  delete m_err;
413  m_err=nullptr;
414  if (mdest) {
415  m_err=new Amg::MatrixX;
416  *m_err = mdest->covariance()->inverse().eval();
417  }
418  }
419  }

◆ printMatScan()

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

Definition at line 348 of file CETmaterial.cxx.

348  {
349 
350  std::ofstream myfilemat;
351  myfilemat.open(m_matScanFile,std::ios::app);
352  myfilemat << theta << " " << phi << " " << r << " " << z << " " << mat << " " << name << std::endl;
353 }

◆ 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

BooleanProperty Trk::CETmaterial::m_backward {this, "CheckBackward", false}
private

Definition at line 85 of file CETmaterial.h.

◆ m_charge

DoubleProperty Trk::CETmaterial::m_charge {this, "StartPerigeeCharge", 1.}
private

Definition at line 74 of file CETmaterial.h.

◆ m_checkStepWise

BooleanProperty Trk::CETmaterial::m_checkStepWise {this, "CheckActiveLayers", false}
private

Definition at line 76 of file CETmaterial.h.

◆ m_chronoStatSvc

IChronoStatSvc_t Trk::CETmaterial::m_chronoStatSvc
private

Definition at line 104 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

BooleanProperty Trk::CETmaterial::m_domsentry {this, "CheckMSentry", false}
private

Definition at line 86 of file CETmaterial.h.

◆ m_doprecision

BooleanProperty Trk::CETmaterial::m_doprecision {this, "CheckPrecision", false}
private

Definition at line 87 of file CETmaterial.h.

◆ m_err

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

Definition at line 94 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 {this, "Extrapolator", "Trk::Extrapolator/AtlasExtrapolator"}
private

The Extrapolator(s) to be retrieved.

Definition at line 63 of file CETmaterial.h.

◆ m_extraprec

ToolHandle<IExtrapolator> Trk::CETmaterial::m_extraprec {this, "ExtraExtrapolator", "Trk::Extrapolator/MuonExtrapolator"}
private

Definition at line 65 of file CETmaterial.h.

◆ m_id

int Trk::CETmaterial::m_id = 0
private

Definition at line 91 of file CETmaterial.h.

◆ m_matActiveFile

const char* Trk::CETmaterial::m_matActiveFile = "mat_active.txt"
private

Definition at line 82 of file CETmaterial.h.

◆ m_matCompFile

const char* Trk::CETmaterial::m_matCompFile = "material_comp.txt"
private

Definition at line 83 of file CETmaterial.h.

◆ m_matSaved

double Trk::CETmaterial::m_matSaved = 0.
private

Definition at line 92 of file CETmaterial.h.

◆ m_matScanFile

const char* Trk::CETmaterial::m_matScanFile = "material_scan.txt"
private

Definition at line 81 of file CETmaterial.h.

◆ m_matTotFile

const char* Trk::CETmaterial::m_matTotFile = "material.txt"
private

Definition at line 80 of file CETmaterial.h.

◆ m_maxP

DoubleProperty Trk::CETmaterial::m_maxP {this, "StartPerigeeMaxP", 50000*Gaudi::Units::GeV}
private

Definition at line 73 of file CETmaterial.h.

◆ m_maxTheta

DoubleProperty Trk::CETmaterial::m_maxTheta {this, "StartPerigeeMaxTheta", M_PI}
private

Definition at line 71 of file CETmaterial.h.

◆ m_maxZ0

DoubleProperty Trk::CETmaterial::m_maxZ0 {this, "StartPerigeeMaxZ0", 0.}
private

Definition at line 69 of file CETmaterial.h.

◆ m_minP

DoubleProperty Trk::CETmaterial::m_minP {this, "StartPerigeeMinP", 50000*Gaudi::Units::GeV}
private

Definition at line 72 of file CETmaterial.h.

◆ m_minTheta

DoubleProperty Trk::CETmaterial::m_minTheta {this, "StartPerigeeMinTheta", 0.}
private

Definition at line 70 of file CETmaterial.h.

◆ m_minZ0

DoubleProperty Trk::CETmaterial::m_minZ0 {this, "StartPerigeeMinZ0", 0.}
private

Definition at line 68 of file CETmaterial.h.

◆ m_msentry

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

Definition at line 98 of file CETmaterial.h.

◆ m_next

Trk::TrackParameters* Trk::CETmaterial::m_next = nullptr
private

Definition at line 93 of file CETmaterial.h.

◆ m_numScan

UnsignedIntegerProperty Trk::CETmaterial::m_numScan {this, "NumberOfScanTracks", 10}
private

Definition at line 75 of file CETmaterial.h.

◆ m_outerBoundary

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

Definition at line 96 of file CETmaterial.h.

◆ m_particleType

IntegerProperty Trk::CETmaterial::m_particleType
private
Initial value:
{this, "ParticleType", Trk::muon,
"the particle type for the extrap."}

Definition at line 100 of file CETmaterial.h.

◆ m_ph

double Trk::CETmaterial::m_ph = 0.
private

Definition at line 90 of file CETmaterial.h.

◆ m_printActive

BooleanProperty Trk::CETmaterial::m_printActive {this, "PrintActivePos", false}
private

Definition at line 78 of file CETmaterial.h.

◆ m_printMaterial

BooleanProperty Trk::CETmaterial::m_printMaterial {this, "PrintMaterial", false}
private

Definition at line 77 of file CETmaterial.h.

◆ m_th

double Trk::CETmaterial::m_th = 0.
private

Definition at line 89 of file CETmaterial.h.

◆ m_trackingGeometry

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

Definition at line 97 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:91
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:90
Trk::CETmaterial::m_printMaterial
BooleanProperty m_printMaterial
Definition: CETmaterial.h:77
Trk::z
@ z
global position (cartesian)
Definition: ParamDefs.h:57
Trk::CETmaterial::m_charge
DoubleProperty m_charge
Definition: CETmaterial.h:74
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::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:421
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::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
skel.it
it
Definition: skel.GENtoEVGEN.py:407
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:93
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:81
Trk::Layer::enclosingDetachedTrackingVolume
const DetachedTrackingVolume * enclosingDetachedTrackingVolume() const
get the confining DetachedTrackingVolume
Trk::CETmaterial::m_domsentry
BooleanProperty m_domsentry
Definition: CETmaterial.h:86
Trk::MaterialEffectsBase::thicknessInX0
double thicknessInX0() const
returns the actually traversed material .
Trk::alongMomentum
@ alongMomentum
Definition: PropDirection.h:20
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::printMatPrec
void printMatPrec(double theta, double phi, const Trk::TrackParameters *, const Trk::TrackParameters *, double mat, int id, const std::string &name)
Definition: CETmaterial.cxx:355
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:29
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:340
Trk::CETmaterial::m_trackingGeometry
const Trk::TrackingGeometry * m_trackingGeometry
Definition: CETmaterial.h:97
Trk::CETmaterial::m_outerBoundary
const Trk::Surface * m_outerBoundary
Definition: CETmaterial.h:96
Trk::TrackingGeometry::highestTrackingVolume
const TrackingVolume * highestTrackingVolume() const
return the world
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
Trk::CETmaterial::m_minTheta
DoubleProperty m_minTheta
Definition: CETmaterial.h:70
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_maxTheta
DoubleProperty m_maxTheta
Definition: CETmaterial.h:71
Trk::CETmaterial::m_minZ0
DoubleProperty m_minZ0
Definition: CETmaterial.h:68
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:83
lumiFormat.i
int i
Definition: lumiFormat.py:85
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::m_numScan
UnsignedIntegerProperty m_numScan
Definition: CETmaterial.h:75
Trk::CETmaterial::printMatScan
void printMatScan(double theta, double phi, double r, double z, double mat, const std::string &name) const
Definition: CETmaterial.cxx:348
Trk::CETmaterial::m_err
Amg::MatrixX * m_err
Definition: CETmaterial.h:94
Trk::CETmaterial::m_chronoStatSvc
IChronoStatSvc_t m_chronoStatSvc
Definition: CETmaterial.h:104
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
Trk::CETmaterial::m_msentry
const Trk::TrackingVolume * m_msentry
Definition: CETmaterial.h:98
Trk::muon
@ muon
Definition: ParticleHypothesis.h:28
Trk::CETmaterial::m_particleType
IntegerProperty m_particleType
Definition: CETmaterial.h:100
Trk::CETmaterial::m_checkStepWise
BooleanProperty m_checkStepWise
Definition: CETmaterial.h:76
Trk::CETmaterial::m_minP
DoubleProperty m_minP
Definition: CETmaterial.h:72
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
Trk::CETmaterial::m_th
double m_th
Definition: CETmaterial.h:89
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
GetAllXsec.entry
list entry
Definition: GetAllXsec.py:132
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
id
SG::auxid_t id
Definition: Control/AthContainers/Root/debug.cxx:239
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
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::DetachedTrackingVolume::name
const std::string & name() const
returns the Name
Definition: DetachedTrackingVolume.h:134
Trk::CETmaterial::m_printActive
BooleanProperty m_printActive
Definition: CETmaterial.h:78
Trk::CETmaterial::m_doprecision
BooleanProperty m_doprecision
Definition: CETmaterial.h:87
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
Trk::CETmaterial::m_maxP
DoubleProperty m_maxP
Definition: CETmaterial.h:73
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:92
a
TList * a
Definition: liststreamerinfos.cxx:10
Trk::CETmaterial::m_backward
BooleanProperty m_backward
Definition: CETmaterial.h:85
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
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:64
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:629
Trk::CETmaterial::m_matActiveFile
const char * m_matActiveFile
Definition: CETmaterial.h:82
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_matTotFile
const char * m_matTotFile
Definition: CETmaterial.h:80
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
Trk::CETmaterial::m_maxZ0
DoubleProperty m_maxZ0
Definition: CETmaterial.h:69
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::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:66
Generated on Fri Apr 11 2025 22:14:58 for ATLAS Offline Software by   doxygen 1.8.18