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

#include <TRT_PrepDataToxAOD.h>

Inheritance diagram for TRT_PrepDataToxAOD:
Collaboration diagram for TRT_PrepDataToxAOD:

Public Types

enum  EGasType {
  kUnset =-128, kXenon =0, kArgon =1, kKrypton =2,
  kEmAr =6, kEmKr =7
}
 

Public Member Functions

 TRT_PrepDataToxAOD (const std::string &name, ISvcLocator *pSvcLocator)
 
virtual StatusCode initialize ()
 
virtual StatusCode execute ()
 
virtual StatusCode finalize ()
 
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< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

Gaudi::Property< bool > m_useTruthInfo {this, "UseTruthInfo",false,"Use data object"}
 
Gaudi::Property< bool > m_writeSDOs {this, "WriteSDOs", false,"Use data object"}
 
SG::ReadHandleKey< InDet::TRT_DriftCircleContainerm_driftcirclecontainer {this,"DriftCircleContainer","TRT_DriftCircles","RHK to retrieve the drift Circle Container"}
 
SG::ReadHandleKey< PRD_MultiTruthCollectionm_multiTruth {this,"PRD_MultiTruth","PRD_MultiTruthTRT","RHK to retrieve PRD multitruth"}
 
SG::ReadHandleKey< InDetSimDataCollectionm_SDOcontainer {this,"MC_TRTUncompressedHit","TRT_SDO_Map","RHK to retrieve TRT SDO map"}
 
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainerm_xAodContainer {this,"TRTxAodContainer","TRT_DriftCircles","WHK to write the TRT drift circle container"}
 
SG::WriteHandleKey< std::vector< unsigned int > > m_xAodOffset {this,"TRTxAodOffset","TRT_DriftCirclesOffsets","WHK to write the TRT drift circle offsets"}
 
ToolHandle< ITRT_DriftFunctionToolm_driftFunctionTool
 DriftFunctionTool. More...
 
ToolHandle< ITRT_CalDbToolm_trtcaldbTool
 
ServiceHandle< ITRT_StrawNeighbourSvcm_neighbourSvc
 
ToolHandle< ITRT_StrawStatusSummaryToolm_TRTStrawSummaryTool
 
const TRT_IDm_TRTHelper
 
const InDetDD::TRT_DetectorManagerm_trtman
 
bool m_firstEventWarnings
 
DataObjIDColl m_extendedExtraObjects
 
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 40 of file TRT_PrepDataToxAOD.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ EGasType

Enumerator
kUnset 
kXenon 
kArgon 
kKrypton 
kEmAr 
kEmKr 

Definition at line 43 of file TRT_PrepDataToxAOD.h.

43 {kUnset=-128,kXenon=0,kArgon=1,kKrypton=2, kEmAr=6,kEmKr=7};

Constructor & Destructor Documentation

◆ TRT_PrepDataToxAOD()

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

Definition at line 41 of file TRT_PrepDataToxAOD.cxx.

41  :
42  AthAlgorithm(name,pSvcLocator),
43  m_driftFunctionTool("TRT_DriftFunctionTool", this),
44  m_trtcaldbTool("TRT_CalDbTool", this),
45  m_neighbourSvc("TRT_StrawNeighbourSvc", name),
46  m_TRTStrawSummaryTool("TRT_StrawStatusSummaryTool",this),
47  m_TRTHelper(nullptr),
48  m_trtman(nullptr),
50 {
51 
52  // --- Services and Tools
53  declareProperty("TRTDriftFunctionTool", m_driftFunctionTool);
54  declareProperty("TRTCalDbTool", m_trtcaldbTool);
55  declareProperty("NeighbourSvc", m_neighbourSvc);
56  declareProperty("TRTStrawSummaryTool", m_TRTStrawSummaryTool);
57 
58 }

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;
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 TRT_PrepDataToxAOD::execute ( )
virtual

Definition at line 97 of file TRT_PrepDataToxAOD.cxx.

98 {
99  //This is needed for the algorithm. If not there, it fails
100 
102  if (not h_trtPrds.isValid()) {
103  ATH_MSG_ERROR("Cannot retrieve TRT PrepDataContainer " << m_driftcirclecontainer.key());
104  return StatusCode::FAILURE;
105  }
106 
107  //This is optional for the algorithm. If not there, just print a one-time warning
108  // On ESD
109 
110  const PRD_MultiTruthCollection* prdmtColl = nullptr; // to be used in the loop later
111  if (m_useTruthInfo && (!m_multiTruth.key().empty()) ) {
113  if (not h_prdmtColl.isValid()){
114  if (m_firstEventWarnings) {
115  ATH_MSG_WARNING("PRD MultiTruth collection not available (" << m_multiTruth.key() << "). Skipping this info although requested.");}
116  } else {
117  prdmtColl = h_prdmtColl.cptr();
118  }
119  }
120 
121  //This is optional for the algorithm. If not there, just print a one-time warning
122  // On RDO
123 
124  const InDetSimDataCollection* sdoCollection = nullptr; // to be used in the loop later
125  if (m_writeSDOs && m_useTruthInfo && (!m_SDOcontainer.key().empty()) ) {
127  if (not h_sdoCollection.isValid()) {
129  ATH_MSG_WARNING("SDO Collection not available (" << m_SDOcontainer.key() << "). Skipping this info although requested.");}
130  } else{
131  sdoCollection = h_sdoCollection.cptr();
132  }
133  }
134 
135 
136  // Create the xAOD container and its auxiliary store:
138  ATH_CHECK(xaod.record(std::make_unique<xAOD::TrackMeasurementValidationContainer>(),std::make_unique<xAOD::TrackMeasurementValidationAuxContainer>() ) );
139 
141  ATH_CHECK(offsets.record(std::make_unique<std::vector<unsigned int>>(m_TRTHelper->straw_layer_hash_max() , 0) ));
142 
143  InDet::TRT_DriftCircleContainer::const_iterator it = h_trtPrds->begin();
144  InDet::TRT_DriftCircleContainer::const_iterator it_end = h_trtPrds->end();
145  unsigned int counter(0);
146  for( ; it!=it_end; ++it ) {
147 
148  //Fill Offset container
149  if( m_TRTHelper->straw_layer_hash_max() <= (*it)->identifyHash() )
150  ATH_MSG_ERROR("My assumption about the maximum size of the hash was wrong");
151  (*offsets)[ (*it)->identifyHash() ] = counter;
152 
153  // skip empty collections
154  if( (*it)->empty() ) continue;
155 
156  // loop over collection and convert to xAOD
157  for( const auto prd : **it ){
158 
159  ++counter;
160 
161  // create and add xAOD object
163  xaod->push_back(xprd);
164 
165  Identifier surfaceID = prd->identify();
166 
167  xprd->setIdentifier(surfaceID.get_compact());
168 
169 
170  // Save ID info:
171  AUXDATA(xprd, int, bec) = m_TRTHelper->barrel_ec( surfaceID ) ;
172  AUXDATA(xprd, int, layer) = m_TRTHelper->layer_or_wheel(surfaceID ) ;
173  AUXDATA(xprd, int, phi_module) = m_TRTHelper->phi_module( surfaceID ) ;
174  AUXDATA(xprd, int, strawlayer) = m_TRTHelper->straw_layer( surfaceID ) ;
175  AUXDATA(xprd, int, strawnumber) = m_TRTHelper->straw( surfaceID ) ;
176  const InDetDD::TRT_BaseElement* element = m_trtman->getElement(surfaceID);
177  AUXDATA(xprd, float, strawphi) = element->center(surfaceID).phi();
178 
179  int chip=0;
180  int board=-1;
181  m_neighbourSvc->getChip(surfaceID,chip);
182  if(abs(m_TRTHelper->barrel_ec(surfaceID))<2)
183  board=m_neighbourSvc->chipToBoardBarrel(chip,m_TRTHelper->layer_or_wheel(surfaceID));
184  else if (chip<12)
185  board=0;
186  else {
187  chip=chip-20;
188  board=1;
189  }
190  AUXDATA(xprd, int, TRTboard) = board ;
191  AUXDATA(xprd, int, TRTchip) = chip ;
192 
193 
194  //Set Local Position
195  const Amg::Vector2D& locpos = prd->localPosition();
196  float locY(0.);
197  float locX = locpos.x();
198  if ( !(std::isinf(locpos.y()) || std::isnan(locpos.y())) ){
199  if (locpos.y()>=1e-07)
200  locY=locpos.y();
201  } else {
202  locY = -9999.;
203  }
204 
205  // Set local error matrix
206  xprd->setLocalPosition(locX,locY);
207 
208  const Amg::MatrixX& localCov = prd->localCovariance();
209  if(localCov.size() == 1){
210  xprd->setLocalPositionError( sqrt(localCov(0,0)), 0., 0. );
211  } else if(localCov.size() == 2){
212  xprd->setLocalPositionError( sqrt(localCov(0,0)), sqrt(localCov(1,1)), sqrt(localCov(0,1)) );
213  } else {
214  xprd->setLocalPositionError(0.,0.,0.);
215  }
216 
217 
218  //Set Global position
219  const Amg::Vector3D gpos = prd->detectorElement()->surface(surfaceID).localToGlobal(prd->localPosition());
220  xprd->setGlobalPosition(gpos.x(),gpos.y(),gpos.z());
221 
222 
223  //TRT hit bit word
224  unsigned int word = prd->getWord();
225 
226 
227  double tot = prd->timeOverThreshold();
228  bool isvalid=false;
229  AUXDATA(xprd, float, drifttime) = prd->driftTime(isvalid) ;
230  AUXDATA(xprd, int, status) = isvalid;
231  AUXDATA(xprd, float, tot) = tot ;
232  AUXDATA(xprd, char, isHT) = prd->highLevel() ;
233  AUXDATA(xprd, float, T0) = m_trtcaldbTool->getT0(surfaceID) ;
234 
235 
236  // Save time info:
237  AUXDATA(xprd, float, leadingEdge) = prd->rawDriftTime();
238 
239  //Drift Time corrections
240  AUXDATA(xprd, float, driftTimeToTCorrection) = m_driftFunctionTool->driftTimeToTCorrection(tot,surfaceID);
241  if(prd->highLevel()) {
242  AUXDATA(xprd, float, driftTimeHTCorrection) = m_driftFunctionTool->driftTimeHTCorrection(surfaceID);
243  }
244  else {
245  AUXDATA(xprd, float, driftTimeHTCorrection) = 0;
246  }
247 
248 
249  //HT Bit patterens
250  unsigned int theWord = word & 0x04020100; //HT bits mask
251  char highThreshold = 0;
252  //this is the agreed-upon binning for HT bit patterns...
253  if (theWord == 0x04000000) //pattern 1
254  highThreshold = 1;
255  else if (theWord == 0x00020000) //pattern 2
256  highThreshold = 2;
257  else if (theWord == 0x00000100) //pattern 3
258  highThreshold = 3;
259  else if (theWord == 0x04020000) //pattern 4
260  highThreshold = 4;
261  else if (theWord == 0x00020100) //pattern 5
262  highThreshold = 5;
263  else if (theWord == 0x04000100) //pattern 6
264  highThreshold = 6;
265  else if (theWord == 0x04020100) //pattern 7
266  highThreshold = 7;
267 
268  AUXDATA(xprd, char, highThreshold) = highThreshold;
269 
270  //Full bit patternword from the TRT
271  AUXDATA(xprd, unsigned int, bitPattern) = word;
272 
273  char gas_type = kUnset;
274  if (!m_TRTStrawSummaryTool.empty()) {
275  int stat = m_TRTStrawSummaryTool->getStatusHT(surfaceID, Gaudi::Hive::currentContext());
276 
277  if ( stat==1 || stat==4 ) { gas_type = kArgon; }
278  else if ( stat==5 ) { gas_type = kKrypton; }
279  else if ( stat==2 || stat==3 ) { gas_type = kXenon; }
280  else if ( stat==6 ) gas_type = kEmAr;
281  else if ( stat==7 ) gas_type = kEmKr;
282  }
283  AUXDATA(xprd, char, gasType) = gas_type;
284 
285  // Use the MultiTruth Collection to get a list of all true particle contributing to the DC
286  if (m_useTruthInfo){
287  if(prdmtColl){
288  std::vector<int> barcodes;
289  auto range = prdmtColl->equal_range(surfaceID);
290  for (auto i = range.first; i != range.second; ++i) {
291  barcodes.push_back( i->second.barcode() );
292  }
293  AUXDATA(xprd, std::vector<int> , truth_barcode) = barcodes;
294  }
295  }
296  if (m_writeSDOs) {
297  if(sdoCollection){
298  // find hit
299  auto pos = sdoCollection->find(surfaceID);
300  int sdo_word = -1000000;
301  if( pos != sdoCollection->end() ) {
302  sdo_word = pos->second.word();
303  }
304  AUXDATA(xprd, int, sdo_word) = sdo_word;
305  }
306  }
307 
308  }
309  }
310  ATH_MSG_DEBUG( " recorded TRT_PrepData obejcts: size " << xaod->size() );
311 
312 
313 
314  // Code to test that something was added to SG:
315  StatusCode sc = StatusCode::SUCCESS;
316 
317  // --- end of event. Disable one-time warnings
318  m_firstEventWarnings = false;
319 
320  return sc;
321 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

◆ extraOutputDeps()

const DataObjIDColl & AthAlgorithm::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 50 of file AthAlgorithm.cxx.

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

◆ finalize()

StatusCode TRT_PrepDataToxAOD::finalize ( )
virtual

Definition at line 328 of file TRT_PrepDataToxAOD.cxx.

329 {
330  return StatusCode::SUCCESS;
331 }

◆ initialize()

StatusCode TRT_PrepDataToxAOD::initialize ( )
virtual

Definition at line 65 of file TRT_PrepDataToxAOD.cxx.

66 {
67  ATH_MSG_DEBUG("Initialize");
68 
69  //Initialize VarHandleKeys
73  ATH_CHECK(m_xAodContainer.initialize());
75 
76  // --- Retrieve services and tools
77  CHECK ( detStore()->retrieve(m_TRTHelper, "TRT_ID") );
78 
79  CHECK ( detStore()->retrieve(m_trtman, "TRT") );
80 
81  CHECK ( m_neighbourSvc.retrieve() );
82 
83  CHECK ( m_trtcaldbTool.retrieve() );
84 
85  CHECK ( m_TRTStrawSummaryTool.retrieve() );
86 
87  CHECK( m_driftFunctionTool.retrieve() );
88 
89  return StatusCode::SUCCESS;
90 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::inputHandles ( ) const
overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ msg() [1/2]

MsgStream& AthCommonMsg< Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

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

◆ msg() [2/2]

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

Definition at line 27 of file AthCommonMsg.h.

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

◆ msgLvl()

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

Definition at line 30 of file AthCommonMsg.h.

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

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles ( ) const
overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< Algorithm > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::renounceArray ( SG::VarHandleKeyArray handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

StatusCode AthAlgorithm::sysInitialize ( )
overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

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

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

Reimplemented from AthCommonDataStore< AthCommonMsg< Algorithm > >.

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

Definition at line 66 of file AthAlgorithm.cxx.

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

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

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

Member Data Documentation

◆ m_detStore

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

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_driftcirclecontainer

SG::ReadHandleKey<InDet::TRT_DriftCircleContainer> TRT_PrepDataToxAOD::m_driftcirclecontainer {this,"DriftCircleContainer","TRT_DriftCircles","RHK to retrieve the drift Circle Container"}
private

Definition at line 59 of file TRT_PrepDataToxAOD.h.

◆ m_driftFunctionTool

ToolHandle< ITRT_DriftFunctionTool > TRT_PrepDataToxAOD::m_driftFunctionTool
private

DriftFunctionTool.

Definition at line 67 of file TRT_PrepDataToxAOD.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_firstEventWarnings

bool TRT_PrepDataToxAOD::m_firstEventWarnings
private

Definition at line 75 of file TRT_PrepDataToxAOD.h.

◆ m_multiTruth

SG::ReadHandleKey<PRD_MultiTruthCollection> TRT_PrepDataToxAOD::m_multiTruth {this,"PRD_MultiTruth","PRD_MultiTruthTRT","RHK to retrieve PRD multitruth"}
private

Definition at line 60 of file TRT_PrepDataToxAOD.h.

◆ m_neighbourSvc

ServiceHandle<ITRT_StrawNeighbourSvc> TRT_PrepDataToxAOD::m_neighbourSvc
private

Definition at line 69 of file TRT_PrepDataToxAOD.h.

◆ m_SDOcontainer

SG::ReadHandleKey<InDetSimDataCollection> TRT_PrepDataToxAOD::m_SDOcontainer {this,"MC_TRTUncompressedHit","TRT_SDO_Map","RHK to retrieve TRT SDO map"}
private

Definition at line 61 of file TRT_PrepDataToxAOD.h.

◆ m_trtcaldbTool

ToolHandle<ITRT_CalDbTool> TRT_PrepDataToxAOD::m_trtcaldbTool
private

Definition at line 68 of file TRT_PrepDataToxAOD.h.

◆ m_TRTHelper

const TRT_ID* TRT_PrepDataToxAOD::m_TRTHelper
private

Definition at line 71 of file TRT_PrepDataToxAOD.h.

◆ m_trtman

const InDetDD::TRT_DetectorManager* TRT_PrepDataToxAOD::m_trtman
private

Definition at line 72 of file TRT_PrepDataToxAOD.h.

◆ m_TRTStrawSummaryTool

ToolHandle<ITRT_StrawStatusSummaryTool> TRT_PrepDataToxAOD::m_TRTStrawSummaryTool
private

Definition at line 70 of file TRT_PrepDataToxAOD.h.

◆ m_useTruthInfo

Gaudi::Property<bool> TRT_PrepDataToxAOD::m_useTruthInfo {this, "UseTruthInfo",false,"Use data object"}
private

Definition at line 55 of file TRT_PrepDataToxAOD.h.

◆ m_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

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

Definition at line 398 of file AthCommonDataStore.h.

◆ m_writeSDOs

Gaudi::Property<bool> TRT_PrepDataToxAOD::m_writeSDOs {this, "WriteSDOs", false,"Use data object"}
private

Definition at line 56 of file TRT_PrepDataToxAOD.h.

◆ m_xAodContainer

SG::WriteHandleKey<xAOD::TrackMeasurementValidationContainer> TRT_PrepDataToxAOD::m_xAodContainer {this,"TRTxAodContainer","TRT_DriftCircles","WHK to write the TRT drift circle container"}
private

Definition at line 62 of file TRT_PrepDataToxAOD.h.

◆ m_xAodOffset

SG::WriteHandleKey<std::vector<unsigned int> > TRT_PrepDataToxAOD::m_xAodOffset {this,"TRTxAodOffset","TRT_DriftCirclesOffsets","WHK to write the TRT drift circle offsets"}
private

Definition at line 63 of file TRT_PrepDataToxAOD.h.


The documentation for this class was generated from the following files:
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
AllowedVariables::e
e
Definition: AsgElectronSelectorTool.cxx:37
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
Trk::locX
@ locX
Definition: ParamDefs.h:37
Trk::locY
@ locY
local cartesian
Definition: ParamDefs.h:38
TRT_PrepDataToxAOD::m_driftcirclecontainer
SG::ReadHandleKey< InDet::TRT_DriftCircleContainer > m_driftcirclecontainer
Definition: TRT_PrepDataToxAOD.h:59
ParticleGun_SamplingFraction.bec
int bec
Definition: ParticleGun_SamplingFraction.py:89
make_unique
std::unique_ptr< T > make_unique(Args &&... args)
Definition: SkimmingToolEXOT5.cxx:23
Amg::Vector2D
Eigen::Matrix< double, 2, 1 > Vector2D
Definition: GeoPrimitives.h:48
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
AthCommonDataStore< AthCommonMsg< Algorithm > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
xAOD::TrackMeasurementValidation
TrackMeasurementValidation_v1 TrackMeasurementValidation
Reference the current persistent version:
Definition: TrackMeasurementValidation.h:13
TRT_PrepDataToxAOD::m_firstEventWarnings
bool m_firstEventWarnings
Definition: TRT_PrepDataToxAOD.h:75
PRD_MultiTruthCollection
A PRD is mapped onto all contributing particles.
Definition: PRD_MultiTruthCollection.h:24
TRT_ID::straw_layer_hash_max
size_type straw_layer_hash_max(void) const
Definition: TRT_ID.h:920
TRT_PrepDataToxAOD::m_trtman
const InDetDD::TRT_DetectorManager * m_trtman
Definition: TRT_PrepDataToxAOD.h:72
AUXDATA
#define AUXDATA(OBJ, TYP, NAME)
Definition: TRT_PrepDataToxAOD.cxx:33
skel.it
it
Definition: skel.GENtoEVGEN.py:396
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
Identifier::get_compact
value_type get_compact() const
Get the compact id.
InDetSimDataCollection
Definition: InDetSimDataCollection.h:25
TRT_PrepDataToxAOD::m_TRTHelper
const TRT_ID * m_TRTHelper
Definition: TRT_PrepDataToxAOD.h:71
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
TRT_PrepDataToxAOD::m_multiTruth
SG::ReadHandleKey< PRD_MultiTruthCollection > m_multiTruth
Definition: TRT_PrepDataToxAOD.h:60
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
AthCommonDataStore< AthCommonMsg< Algorithm > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
TRT_PrepDataToxAOD::m_SDOcontainer
SG::ReadHandleKey< InDetSimDataCollection > m_SDOcontainer
Definition: TRT_PrepDataToxAOD.h:61
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:29
TRT_PrepDataToxAOD::m_neighbourSvc
ServiceHandle< ITRT_StrawNeighbourSvc > m_neighbourSvc
Definition: TRT_PrepDataToxAOD.h:69
xAOD::TrackMeasurementValidation_v1
Class describing a TrackMeasurementValidation.
Definition: TrackMeasurementValidation_v1.h:27
TRT_PrepDataToxAOD::kXenon
@ kXenon
Definition: TRT_PrepDataToxAOD.h:43
TRT_ID::straw
int straw(const Identifier &id) const
Definition: TRT_ID.h:902
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
AthAlgorithm::sysInitialize
virtual StatusCode sysInitialize() override
Override sysInitialize.
Definition: AthAlgorithm.cxx:66
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
AthCommonDataStore< AthCommonMsg< Algorithm > >::outputHandles
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
xAOD::TrackMeasurementValidation_v1::setLocalPosition
void setLocalPosition(float localX, float localY)
Sets the local position.
python.Dumpers.barcodes
def barcodes(beg, end, sz)
Definition: Dumpers.py:2800
lumiFormat.i
int i
Definition: lumiFormat.py:85
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
TRT_PrepDataToxAOD::m_TRTStrawSummaryTool
ToolHandle< ITRT_StrawStatusSummaryTool > m_TRTStrawSummaryTool
Definition: TRT_PrepDataToxAOD.h:70
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
TRT_PrepDataToxAOD::kEmKr
@ kEmKr
Definition: TRT_PrepDataToxAOD.h:43
plotBeamSpotVxVal.range
range
Definition: plotBeamSpotVxVal.py:195
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
AthCommonDataStore< AthCommonMsg< Algorithm > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
TRT_ID::barrel_ec
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: TRT_ID.h:866
TRT_ID::straw_layer
int straw_layer(const Identifier &id) const
Definition: TRT_ID.h:893
xAOD::TrackMeasurementValidation_v1::setIdentifier
void setIdentifier(uint64_t identifier)
Sets the identifier.
TRT_ID::layer_or_wheel
int layer_or_wheel(const Identifier &id) const
Definition: TRT_ID.h:884
beamspotman.stat
stat
Definition: beamspotman.py:266
TRT_PrepDataToxAOD::m_useTruthInfo
Gaudi::Property< bool > m_useTruthInfo
Definition: TRT_PrepDataToxAOD.h:55
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
TRT_PrepDataToxAOD::m_xAodOffset
SG::WriteHandleKey< std::vector< unsigned int > > m_xAodOffset
Definition: TRT_PrepDataToxAOD.h:63
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
xAOD::TrackMeasurementValidation_v1::setGlobalPosition
void setGlobalPosition(float globalX, float globalY, float globalZ)
Sets the global position.
Definition: TrackMeasurementValidation_v1.cxx:50
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
InDetDD::TRT_DetectorManager::getElement
const TRT_BaseElement * getElement(Identifier id) const
Access Elements Generically---------------------------------------------—.
Definition: TRT_DetectorManager.cxx:148
xAOD::TrackMeasurementValidation_v1::setLocalPositionError
void setLocalPositionError(float localXError, float localYError, float localXYCorrelation)
Sets the local position error.
Definition: TrackMeasurementValidation_v1.cxx:37
Amg::Vector3D
Eigen::Matrix< double, 3, 1 > Vector3D
Definition: GeoPrimitives.h:47
TRT_ID::phi_module
int phi_module(const Identifier &id) const
Definition: TRT_ID.h:875
TRT_PrepDataToxAOD::kEmAr
@ kEmAr
Definition: TRT_PrepDataToxAOD.h:43
python.LumiBlobConversion.pos
pos
Definition: LumiBlobConversion.py:18
TRT_PrepDataToxAOD::m_trtcaldbTool
ToolHandle< ITRT_CalDbTool > m_trtcaldbTool
Definition: TRT_PrepDataToxAOD.h:68
SG::WriteHandle
Definition: StoreGate/StoreGate/WriteHandle.h:76
AthAlgorithm::m_extendedExtraObjects
DataObjIDColl m_extendedExtraObjects
Definition: AthAlgorithm.h:79
a
TList * a
Definition: liststreamerinfos.cxx:10
h
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
InDetDD::TRT_BaseElement::center
virtual const Amg::Vector3D & center() const override final
Element Surface: center of a straw layer.
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
TRT_PrepDataToxAOD::m_xAodContainer
SG::WriteHandleKey< xAOD::TrackMeasurementValidationContainer > m_xAodContainer
Definition: TRT_PrepDataToxAOD.h:62
TRT_PrepDataToxAOD::kArgon
@ kArgon
Definition: TRT_PrepDataToxAOD.h:43
AthAlgorithm::AthAlgorithm
AthAlgorithm()
Default constructor:
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
merge.status
status
Definition: merge.py:17
LHEF::Writer
Pythia8::Writer Writer
Definition: Prophecy4fMerger.cxx:12
TRT_PrepDataToxAOD::m_driftFunctionTool
ToolHandle< ITRT_DriftFunctionTool > m_driftFunctionTool
DriftFunctionTool.
Definition: TRT_PrepDataToxAOD.h:67
test_pyathena.counter
counter
Definition: test_pyathena.py:15
TRT_PrepDataToxAOD::m_writeSDOs
Gaudi::Property< bool > m_writeSDOs
Definition: TRT_PrepDataToxAOD.h:56
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
TRT_PrepDataToxAOD::kUnset
@ kUnset
Definition: TRT_PrepDataToxAOD.h:43
fitman.k
k
Definition: fitman.py:528
InDetDD::TRT_BaseElement
Definition: TRT_BaseElement.h:57
ServiceHandle< ICondSvc >
Identifier
Definition: IdentifierFieldParser.cxx:14
TRT_PrepDataToxAOD::kKrypton
@ kKrypton
Definition: TRT_PrepDataToxAOD.h:43