|
ATLAS Offline Software
|
#include <TrackTimeDefAndQualityAlg.h>
|
| TrackTimeDefAndQualityAlg (const std::string &name, ISvcLocator *pSvcLocator) |
|
virtual | ~TrackTimeDefAndQualityAlg () |
|
virtual StatusCode | initialize () override final |
|
virtual StatusCode | execute (const EventContext &ctx) const override final |
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual bool | isClonable () const override |
| Specify if the algorithm is clonable. More...
|
|
virtual unsigned int | cardinality () const override |
| Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
|
|
virtual StatusCode | sysExecute (const EventContext &ctx) override |
| Execute an algorithm. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
virtual bool | filterPassed (const EventContext &ctx) const |
|
virtual void | setFilterPassed (bool state, const EventContext &ctx) const |
|
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 |
|
|
CleaningResult | runTimeConsistencyCuts (const xAOD::TrackParticle *track_particle) const |
|
std::array< Hit, n_hgtd_layers > | getValidHits (const xAOD::TrackParticle *track_particle) const |
|
float | calculateChi2 (const std::array< Hit, n_hgtd_layers > &hits) const |
| Calculates the chi2 of the hit times given their resolution. More...
|
|
bool | passesDeltaT (const std::array< Hit, n_hgtd_layers > &hits) const |
| Checks two hits for time compatibility. More...
|
|
float | meanTime (const std::array< Hit, n_hgtd_layers > &hits) const |
| Calculates the arithmetic mean of the valid hit times;. More...
|
|
float | trackTimeResolution (const std::array< Hit, n_hgtd_layers > &hits) const |
| Calculates the combined resolution. More...
|
|
short | findLayerWithBadChi2 (std::array< Hit, n_hgtd_layers > hits) const |
| Identifies time outliers by finding the layer within which a hit contributes negatively to the overall chi2 value and returns the layer number. More...
|
|
void | setLayerAsInvalid (std::array< Hit, n_hgtd_layers > &hits, short layer) const |
| Given a layer number, the hit sitting on this layer is flagged as invalid. More...
|
|
short | getValidPattern (const std::array< Hit, n_hgtd_layers > &hits) const |
| Returns the pattern of valid hits in HGTD as a 4-bit bitfield, where a 1 encodes that a valid hit was found, while a 0 means no valid hit found. More...
|
|
bool | lastHitIsOnLastSurface (const xAOD::TrackParticle &track_particle) const |
| Checks if the last hit on track was found on a pre-specified set of Pixel and Strip layers close to the HGTD surfaces. More...
|
|
const Trk::TrackParameters * | getLastHitOnTrack (const Trk::Track &track) 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...
|
|
Definition at line 57 of file TrackTimeDefAndQualityAlg.h.
◆ StoreGateSvc_t
◆ TrackTimeDefAndQualityAlg()
HGTD::TrackTimeDefAndQualityAlg::TrackTimeDefAndQualityAlg |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ ~TrackTimeDefAndQualityAlg()
virtual HGTD::TrackTimeDefAndQualityAlg::~TrackTimeDefAndQualityAlg |
( |
| ) |
|
|
inlinevirtual |
◆ calculateChi2()
float HGTD::TrackTimeDefAndQualityAlg::calculateChi2 |
( |
const std::array< Hit, n_hgtd_layers > & |
hits | ) |
const |
|
private |
Calculates the chi2 of the hit times given their resolution.
- Parameters
-
- Returns
- Chi2 value.
Definition at line 227 of file TrackTimeDefAndQualityAlg.cxx.
233 for (
const auto& hit :
hits) {
236 (hit.m_resolution * hit.m_resolution);
◆ cardinality()
unsigned int AthReentrantAlgorithm::cardinality |
( |
| ) |
const |
|
overridevirtualinherited |
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 55 of file AthReentrantAlgorithm.cxx.
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation 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.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode HGTD::TrackTimeDefAndQualityAlg::execute |
( |
const EventContext & |
ctx | ) |
const |
|
finaloverridevirtual |
Definition at line 40 of file TrackTimeDefAndQualityAlg.cxx.
45 trk_ptkl_container_handle.cptr();
46 if (not track_particles) {
48 "[TrackTimeDefAndQualityAlg] TrackParticleContainer not found, "
50 return StatusCode::FAILURE;
62 for (
const auto* track_ptkl : *track_particles) {
76 short prime_pattern = 0x0;
77 for (
short i = 0;
i < n_hgtd_layers;
i++) {
78 if (
res.m_hits.at(
i).m_isprime) {
79 prime_pattern |= (1 <<
i);
85 time_handle(*track_ptkl) =
res.m_time;
86 timeres_handle(*track_ptkl) =
res.m_resolution;
87 hasValidTime_handle(*track_ptkl) =
res.m_hasValidTime;
88 summary_handle(*track_ptkl) =
res.m_field;
90 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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 & AthReentrantAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 79 of file AthReentrantAlgorithm.cxx.
86 return Algorithm::extraOutputDeps();
◆ filterPassed()
virtual bool AthReentrantAlgorithm::filterPassed |
( |
const EventContext & |
ctx | ) |
const |
|
inlinevirtualinherited |
◆ findLayerWithBadChi2()
short HGTD::TrackTimeDefAndQualityAlg::findLayerWithBadChi2 |
( |
std::array< Hit, n_hgtd_layers > |
hits | ) |
const |
|
private |
Identifies time outliers by finding the layer within which a hit contributes negatively to the overall chi2 value and returns the layer number.
Copies the hits, since they have to be invalidated one by one.
- Parameters
-
- Returns
- Layer with outlier time, value between 0 and 3.
Definition at line 293 of file TrackTimeDefAndQualityAlg.cxx.
295 short remove_layer = -1;
296 float local_min_chi2 = 999999;
297 for (
auto& hit :
hits) {
299 bool validbuff = hit.m_isvalid;
300 hit.m_isvalid =
false;
302 hit.m_isvalid = validbuff;
303 if (local_chi2 < local_min_chi2) {
304 local_min_chi2 = local_chi2;
305 remove_layer = hit.m_layer;
◆ getLastHitOnTrack()
Definition at line 322 of file TrackTimeDefAndQualityAlg.cxx.
325 track.trackStateOnSurfaces();
332 for (
auto i = tsos->rbegin();
i != tsos->rend(); ++
i) {
333 const auto* curr_last_tsos = *
i;
334 if (not curr_last_tsos) {
338 curr_last_tsos->trackParameters() and
339 curr_last_tsos->measurementOnTrack()) {
340 return curr_last_tsos->trackParameters();
◆ getValidHits()
Definition at line 183 of file TrackTimeDefAndQualityAlg.cxx.
188 std::vector<float>
times = layerClusterTimeHandle(*track_particle);
192 std::vector<bool> has_clusters = layerHasExtensionHandle(*track_particle);
196 std::vector<int> hit_classification =
197 layerClusterTruthClassHandle(*track_particle);
199 std::array<Hit, n_hgtd_layers> valid_hits;
201 for (
size_t i = 0;
i < n_hgtd_layers;
i++) {
203 if (has_clusters.at(
i)) {
204 newhit.m_time =
times.at(
i);
205 newhit.m_isprime = hit_classification.at(
i) == 1;
206 newhit.m_isvalid =
true;
209 valid_hits.at(
i) = newhit;
◆ getValidPattern()
short HGTD::TrackTimeDefAndQualityAlg::getValidPattern |
( |
const std::array< Hit, n_hgtd_layers > & |
hits | ) |
const |
|
private |
Returns the pattern of valid hits in HGTD as a 4-bit bitfield, where a 1 encodes that a valid hit was found, while a 0 means no valid hit found.
The LSB carries the hit information for the layer closest to the IP.
- Parameters
-
- Returns
- Bit pattern encoding valid hits in layers.
Definition at line 215 of file TrackTimeDefAndQualityAlg.cxx.
219 for (
short i = 0;
i < n_hgtd_layers;
i++) {
220 if (
hits.at(
i).m_isvalid) {
◆ initialize()
StatusCode HGTD::TrackTimeDefAndQualityAlg::initialize |
( |
| ) |
|
|
finaloverridevirtual |
◆ inputHandles()
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.
◆ isClonable()
bool AthReentrantAlgorithm::isClonable |
( |
| ) |
const |
|
overridevirtualinherited |
Specify if the algorithm is clonable.
Reentrant algorithms are clonable.
Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, InDet::SiSPGNNTrackMaker, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, ITkStripCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.
Definition at line 44 of file AthReentrantAlgorithm.cxx.
◆ lastHitIsOnLastSurface()
Checks if the last hit on track was found on a pre-specified set of Pixel and Strip layers close to the HGTD surfaces.
FIXME: should not be hardcoded. To be exchanged by a hole search at a later point.
- Parameters
-
[in] | track_particle | Track reconstructed in ITk. |
- Returns
- True if the last hit on track was within the specified volume.
Definition at line 346 of file TrackTimeDefAndQualityAlg.cxx.
352 float abs_z = std::abs(last_hit_param->
position().z());
357 if (radius < 350 and abs_z > 2400) {
361 if (
radius > 205 and radius < 350 and abs_z > 2100) {
365 if (radius < 220 and abs_z > 2200) {
369 if (radius < 140 and abs_z > 1890) {
◆ meanTime()
float HGTD::TrackTimeDefAndQualityAlg::meanTime |
( |
const std::array< Hit, n_hgtd_layers > & |
hits | ) |
const |
|
private |
Calculates the arithmetic mean of the valid hit times;.
- Parameters
-
- Returns
- Returns true if the times were in agreement.
Definition at line 265 of file TrackTimeDefAndQualityAlg.cxx.
270 for (
const auto& hit :
hits) {
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ outputHandles()
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.
◆ passesDeltaT()
bool HGTD::TrackTimeDefAndQualityAlg::passesDeltaT |
( |
const std::array< Hit, n_hgtd_layers > & |
hits | ) |
const |
|
private |
Checks two hits for time compatibility.
- Parameters
-
- Returns
- Returns true if the times were in agreement.
Definition at line 242 of file TrackTimeDefAndQualityAlg.cxx.
246 short n_valid = std::count_if(
hits.begin(),
hits.end(),
247 [](
const Hit& hit) { return hit.m_isvalid; });
252 std::vector<float>
times;
253 std::vector<float>
res;
254 for (
const auto& hit :
hits) {
256 times.push_back(hit.m_time);
257 res.push_back(hit.m_resolution);
◆ renounce()
◆ renounceArray()
◆ runTimeConsistencyCuts()
Definition at line 97 of file TrackTimeDefAndQualityAlg.cxx.
101 std::array<Hit, n_hgtd_layers> valid_hits =
getValidHits(track_particle);
104 result.m_hits = valid_hits;
108 result.m_hasValidTime = 0;
114 short nhits = std::count_if(valid_hits.begin(), valid_hits.end(),
115 [](
const Hit& hit) { return hit.m_isvalid; });
121 result.m_hasValidTime = recoed_pattern ? 1 : 0;
123 }
else if (nhits == 2) {
131 result.m_hasValidTime = 1;
149 nhits = std::count_if(valid_hits.begin(), valid_hits.end(),
150 [](
const Hit& hit) { return hit.m_isvalid; });
163 result.m_hasValidTime = 1;
176 result.m_hasValidTime = 1;
◆ setFilterPassed()
virtual void AthReentrantAlgorithm::setFilterPassed |
( |
bool |
state, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
inlinevirtualinherited |
◆ setLayerAsInvalid()
void HGTD::TrackTimeDefAndQualityAlg::setLayerAsInvalid |
( |
std::array< Hit, n_hgtd_layers > & |
hits, |
|
|
short |
layer |
|
) |
| const |
|
private |
Given a layer number, the hit sitting on this layer is flagged as invalid.
Used for outlier removal.
- Parameters
-
[in] | hits | Array of hits. |
[in] | layer | The layer that should be masked. |
Definition at line 311 of file TrackTimeDefAndQualityAlg.cxx.
314 for (
auto& hit :
hits) {
315 if (hit.m_layer ==
layer) {
316 hit.m_isvalid =
false;
◆ sysExecute()
StatusCode AthReentrantAlgorithm::sysExecute |
( |
const EventContext & |
ctx | ) |
|
|
overridevirtualinherited |
Execute an algorithm.
We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.
Definition at line 67 of file AthReentrantAlgorithm.cxx.
69 return Gaudi::Algorithm::sysExecute (ctx);
◆ sysInitialize()
StatusCode AthReentrantAlgorithm::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< Gaudi::Algorithm > >.
Reimplemented in InputMakerBase, and HypoBase.
Definition at line 96 of file AthReentrantAlgorithm.cxx.
107 if ( cs.retrieve().isFailure() ) {
109 return StatusCode::SUCCESS;
111 if (cs->regHandle(
this,*
h).isFailure()) {
112 sc = StatusCode::FAILURE;
113 ATH_MSG_ERROR(
"unable to register WriteCondHandle " <<
h->fullKey()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ trackTimeResolution()
float HGTD::TrackTimeDefAndQualityAlg::trackTimeResolution |
( |
const std::array< Hit, n_hgtd_layers > & |
hits | ) |
const |
|
private |
Calculates the combined resolution.
Definition at line 279 of file TrackTimeDefAndQualityAlg.cxx.
284 for (
const auto& hit :
hits) {
286 sum += 1. / (hit.m_resolution * hit.m_resolution);
290 :
static_cast<float>(std::sqrt(1. /
sum));
◆ updateVHKA()
◆ m_chi2_threshold
FloatProperty HGTD::TrackTimeDefAndQualityAlg::m_chi2_threshold |
|
private |
Initial value:{
this, "Chi2Threshold", 1.5,
"Quality cut for decision to keep hits compatible in time"}
Definition at line 120 of file TrackTimeDefAndQualityAlg.h.
◆ m_comp_ptrn_sft
const short HGTD::TrackTimeDefAndQualityAlg::m_comp_ptrn_sft = 8 |
|
private |
◆ m_default_time
FloatProperty HGTD::TrackTimeDefAndQualityAlg::m_default_time |
|
private |
Initial value:{
this, "DefaultTime", 0.0,
"Default time used for tracks without HGTD timing info"}
Definition at line 126 of file TrackTimeDefAndQualityAlg.h.
◆ m_default_time_res
FloatProperty HGTD::TrackTimeDefAndQualityAlg::m_default_time_res |
|
private |
Initial value:{
this, "DefaultTimeRes", 50. / std::sqrt(12.),
"Default time resolution used for tracks without HGTD timing info"}
Definition at line 131 of file TrackTimeDefAndQualityAlg.h.
◆ m_deltat_cut
FloatProperty HGTD::TrackTimeDefAndQualityAlg::m_deltat_cut |
|
private |
◆ m_detStore
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects |
|
privateinherited |
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Empty if no symlinks were found.
Definition at line 153 of file AthReentrantAlgorithm.h.
◆ m_hasValidTime_dec_key
Initial value:{
this, "hasValidTime", "InDetTrackParticles.hasValidTime", "Time assigned to this track is valid"}
Definition at line 95 of file TrackTimeDefAndQualityAlg.h.
◆ m_holes_ptrn_sft
const short HGTD::TrackTimeDefAndQualityAlg::m_holes_ptrn_sft = 12 |
|
private |
◆ m_layerClusterTimeKey
Initial value:{
this, "HGTD_cluster_time", "InDetTrackParticles.HGTD_cluster_time",
"deco with a handle for cluster time"}
Definition at line 81 of file TrackTimeDefAndQualityAlg.h.
◆ m_layerClusterTruthClassKey
Initial value:{
this, "HGTD_cluster_truth_class",
"InDetTrackParticles.HGTD_cluster_truth_class",
"deco with a handle for a truth time"}
Definition at line 85 of file TrackTimeDefAndQualityAlg.h.
◆ m_layerHasExtensionKey
Initial value:{
this, "HGTD_has_extension", "InDetTrackParticles.HGTD_has_extension",
"deco with a handle for an extension"}
Definition at line 78 of file TrackTimeDefAndQualityAlg.h.
◆ m_primes_ptrn_sft
const short HGTD::TrackTimeDefAndQualityAlg::m_primes_ptrn_sft = 16 |
|
private |
◆ m_recoed_ptrn_sft
const short HGTD::TrackTimeDefAndQualityAlg::m_recoed_ptrn_sft = 0 |
|
private |
◆ m_summarypattern_dec_key
Initial value:{this, "HGTD_summaryinfo",
"InDetTrackParticles.HGTD_summaryinfo",
"Bitfield for working point definition"}
Definition at line 98 of file TrackTimeDefAndQualityAlg.h.
◆ m_time_dec_key
Initial value:{
this, "time", "InDetTrackParticles.time", "Time assigned to this track"}
Definition at line 90 of file TrackTimeDefAndQualityAlg.h.
◆ m_time_res_dec_key
Initial value:{
this, "timeResolution", "InDetTrackParticles.timeResolution",
"Time resolution assigned to this track"}
Definition at line 92 of file TrackTimeDefAndQualityAlg.h.
◆ m_trackParticleContainerKey
Initial value:{
this, "TrackParticleContainerName", "InDetTrackParticles",
"Name of the TrackParticle container"}
Definition at line 74 of file TrackTimeDefAndQualityAlg.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
void mean(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
bool passesDeltaT(const std::array< Hit, n_hgtd_layers > &hits) const
Checks two hits for time compatibility.
void setLayerAsInvalid(std::array< Hit, n_hgtd_layers > &hits, short layer) const
Given a layer number, the hit sitting on this layer is flagged as invalid.
const Amg::Vector3D & position() const
Access method for the position.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
bool lastHitIsOnLastSurface(const xAOD::TrackParticle &track_particle) const
Checks if the last hit on track was found on a pre-specified set of Pixel and Strip layers close to t...
short findLayerWithBadChi2(std::array< Hit, n_hgtd_layers > hits) const
Identifies time outliers by finding the layer within which a hit contributes negatively to the overal...
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)
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer > m_layerClusterTruthClassKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_hasValidTime_dec_key
const Trk::TrackParameters * getLastHitOnTrack(const Trk::Track &track) const
virtual void setOwner(IDataHandleHolder *o)=0
AthReentrantAlgorithm()
Default constructor:
Handle class for reading a decoration on an object.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
float trackTimeResolution(const std::array< Hit, n_hgtd_layers > &hits) const
Calculates the combined resolution.
::StatusCode StatusCode
StatusCode definition for legacy code.
Handle class for adding a decoration to an object.
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_summarypattern_dec_key
short getValidPattern(const std::array< Hit, n_hgtd_layers > &hits) const
Returns the pattern of valid hits in HGTD as a 4-bit bitfield, where a 1 encodes that a valid hit was...
std::pair< std::vector< unsigned int >, bool > res
double chi2(TH1 *h0, TH1 *h1)
std::array< Hit, n_hgtd_layers > getValidHits(const xAOD::TrackParticle *track_particle) const
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
float calculateChi2(const std::array< Hit, n_hgtd_layers > &hits) const
Calculates the chi2 of the hit times given their resolution.
virtual void renounce()=0
CleaningResult runTimeConsistencyCuts(const xAOD::TrackParticle *track_particle) const
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
float meanTime(const std::array< Hit, n_hgtd_layers > &hits) const
Calculates the arithmetic mean of the valid hit times;.
FloatProperty m_default_time_res
FloatProperty m_deltat_cut
virtual StatusCode sysInitialize() override
Override sysInitialize.
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer > m_layerClusterTimeKey
const short m_comp_ptrn_sft
FloatProperty m_default_time
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticleContainerKey
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_time_dec_key
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
#define ATH_MSG_WARNING(x)
const short m_holes_ptrn_sft
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer > m_time_res_dec_key
const short m_primes_ptrn_sft
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer > m_layerHasExtensionKey
const short m_recoed_ptrn_sft
const Trk::Track * track() const
Returns a pointer (which can be NULL) to the Trk::Track which was used to make this TrackParticle.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
@ Measurement
This is a measurement, and will at least contain a Trk::MeasurementBase.
FloatProperty m_chi2_threshold