HGTD::TimeCompatibilityCheckAlg Class Reference

#include <TimeCompatibilityCheckAlg.h>

Inheritance diagram for HGTD::TimeCompatibilityCheckAlg:

Collaboration diagram for HGTD::TimeCompatibilityCheckAlg:

Classes
struct	Hit

Public Member Functions
	TimeCompatibilityCheckAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~TimeCompatibilityCheckAlg ()
virtual StatusCode	initialize () override final
virtual StatusCode	execute (const EventContext &ctx) const override final
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
	Get filter decision:
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
	Set filter decision:
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
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.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
float	calculateChi2 (const std::vector< Hit > &hits) const
std::vector< TimeCompatibilityCheckAlg::Hit >	getValidHits (const xAOD::TrackParticle *track_particle) const
bool	passesDeltaT (const std::vector< Hit > &hits) const
std::vector< TimeCompatibilityCheckAlg::Hit >	getTimeCompatibleHits (const xAOD::TrackParticle *track_particle) const
bool	lastHitIsOnLastSurface (const xAOD::TrackParticle &track_particle) const
const Trk::TrackParameters *	getLastHitOnTrack (const Trk::Track &track) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
ToolHandle< IHGTD_HolesITkTool >	m_extensionTool {this, "ITkHoles", "HGTD_IterativeExtensionTool/HGTD_IterativeExtensionTool", "Tool for calculating the holes in"}
SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_trackParticleContainerKey {this, "TrackParticleContainerName", "InDetTrackParticles", "Name of the TrackParticle container"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerHasExtensionKey {this, "HGTD_has_extension", m_trackParticleContainerKey, "HGTD_has_extension", "deco with a handle for an extension"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerExtensionChi2Key {this, "HGTD_extension_chi2", m_trackParticleContainerKey, "HGTD_extension_chi2", "deco with a handle for a ch2 of extension"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerClusterRawTimeKey {this, "HGTD_cluster_raw_time", m_trackParticleContainerKey, "HGTD_cluster_raw_time", "deco with a handle for layer cluster raw time"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerClusterTimeKey {this, "HGTD_cluster_time", m_trackParticleContainerKey, "HGTD_cluster_time", "deco with a handle for cluster time"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerClusterTruthClassKey {this, "HGTD_cluster_truth_class", m_trackParticleContainerKey, "HGTD_cluster_truth_class", "deco with a handle for a truth time"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerClusterShadowedKey {this, "HGTD_cluster_shadowed", m_trackParticleContainerKey, "HGTD_cluster_shadowed", "deco with a handle for a shadowed cluster"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerClusterMergedKey {this, "HGTD_cluster_merged", m_trackParticleContainerKey, "HGTD_cluster_merged", "deco with a handle for a merged cluster"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_layerPrimaryExpectedKey {this, "HGTD_primary_expected", m_trackParticleContainerKey, "HGTD_primary_expected", "deco with a handle for an expected primary"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_extrapXKey {this, "HGTD_extrap_x",m_trackParticleContainerKey, "HGTD_extrap_x", "deco with a handle for an x of extrap"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_extrapYKey {this, "HGTD_extrap_y", m_trackParticleContainerKey, "HGTD_extrap_y", "deco with a handle for an y of extrap"}
SG::ReadDecorHandleKey< xAOD::TrackParticleContainer >	m_holesHGTDKey {this, "HGTD_holes", m_trackParticleContainerKey, "HGTD_holes", "deco with the holes on track in HGTD"}
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer >	m_compatibleHitsTimesKey {this, "HGTD_times_of_compatible_hits",m_trackParticleContainerKey, "HGTD_times_of_compatible_hits", "deco with a handle for the time compatible hits' times"}
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer >	m_lastHitInITkCutKey {this, "HGTD_last_hit_in_ITk_cut", m_trackParticleContainerKey, "HGTD_last_hit_in_ITk_cut", "deco with a handle for the last hit to be close to HGTDrequirement"}
SG::WriteDecorHandleKey< xAOD::TrackParticleContainer >	m_holesITkKey {this, "HGTD_holes_in_ITk", m_trackParticleContainerKey, "HGTD_holes_in_ITk", "deco with a handle for the number of holes on track in ITk between the last hit on track and HGTD"}
FloatProperty	m_chi2_threshold
FloatProperty	m_delta_cut
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default).
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default).
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 44 of file TimeCompatibilityCheckAlg.h.

Member Typedef Documentation

StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

TimeCompatibilityCheckAlg()

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

Definition at line 20 of file TimeCompatibilityCheckAlg.cxx.

        : AthReentrantAlgorithm(name, pSvcLocator){}

~TimeCompatibilityCheckAlg()

virtual HGTD::TimeCompatibilityCheckAlg::~TimeCompatibilityCheckAlg ( )

inlinevirtual

Definition at line 48 of file TimeCompatibilityCheckAlg.h.

{}

Member Function Documentation

calculateChi2()

float HGTD::TimeCompatibilityCheckAlg::calculateChi2 ( const std::vector< Hit > & hits ) const

private

Definition at line 117 of file TimeCompatibilityCheckAlg.cxx.

                                                                              {
 
    //Calculation of the mean
    float sum = 0.;
    for (const auto& hit : hits){
        sum += hit.m_time;
    }
    float mean = sum / (float)hits.size();
 
    //Calculation of the Chi2
    float chi2 = 0.;
    for (size_t i = 0; i < hits.size(); i++) {
        chi2 += (hits.at(i).m_time - mean) * (hits.at(i).m_time - mean) /
                (hits.at(i).m_resolution * hits.at(i).m_resolution);
    }
    return chi2; 
 }

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode HGTD::TimeCompatibilityCheckAlg::execute ( const EventContext & ctx ) const

finaloverridevirtual

This is a first attempt to design a bitfield with an HGTD track summary

uint32_t HGTD_track_summary = 0;

Field will contain the following values: HGTD_acceptance : 1; HGTD_has_extension : 3; (max 4 extensions) HGTD_extension_chi2 : 12; (each layer (4) has 1 option (3 bits): <1 / <2 / <3 / <4) HGTD_number_of_holes : 3; (max 4 holes) HGTD_hits_compatible_in_time : 3 (max 4)

Definition at line 49 of file TimeCompatibilityCheckAlg.cxx.

                                                                           {
 
    SG::ReadHandle<xAOD::TrackParticleContainer> trk_ptkl_container_handle(
        m_trackParticleContainerKey, ctx);
    const xAOD::TrackParticleContainer* track_particles = trk_ptkl_container_handle.cptr();
if (not track_particles) {
    ATH_MSG_ERROR("[TimeCompatibilityCheckAlg] TrackParticleContainer not found, "
                   "aborting execute!");
     return StatusCode::FAILURE;
}
 
    SG::WriteDecorHandle<xAOD::TrackParticleContainer, std::vector<float>> compatibleHitsTimesHandle(
        m_compatibleHitsTimesKey, ctx);  
    SG::WriteDecorHandle<xAOD::TrackParticleContainer, char> lastHitInITkCutHandle(
        m_lastHitInITkCutKey, ctx);
    SG::WriteDecorHandle<xAOD::TrackParticleContainer, int> holesITkHandle(
        m_holesITkKey, ctx);
        
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<float>> layerExtensionChi2Handle(
        m_layerExtensionChi2Key, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<float>> layerClusterRawTimeHandle(
        m_layerClusterRawTimeKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<int>> layerClusterTruthClassHandle(
        m_layerClusterTruthClassKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerClusterShadowedHandle(
        m_layerClusterShadowedKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerClusterMergedHandle(
        m_layerClusterMergedKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerPrimaryExpectedHandle(
        m_layerPrimaryExpectedKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, float> extrapXHandle(
        m_extrapXKey, ctx);
    SG::ReadDecorHandle<xAOD::TrackParticleContainer, float> extrapYHandle(
        m_extrapYKey, ctx);
 
  for (const auto* track_ptkl : *track_particles) {
 
// Filling here the decoration with times of the compatible hits per track
 
    for(const auto hit : getTimeCompatibleHits(track_ptkl)){
        compatibleHitsTimesHandle(*track_ptkl).push_back(hit.m_time);
    }
    lastHitInITkCutHandle(*track_ptkl) = lastHitIsOnLastSurface(*track_ptkl);
 
    holesITkHandle(*track_ptkl) = (m_extensionTool->getHolesITk(ctx, *track_ptkl)).size();
 

 
  }
  return StatusCode::SUCCESS;
}

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::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 & AthCommonAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 89 of file AthCommonAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

filterPassed()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

getLastHitOnTrack()

const Trk::TrackParameters * HGTD::TimeCompatibilityCheckAlg::getLastHitOnTrack ( const Trk::Track & track ) const

private

Definition at line 257 of file TimeCompatibilityCheckAlg.cxx.

                                                                        {
 
  const DataVector<const Trk::TrackStateOnSurface>* tsos =
      track.trackStateOnSurfaces();
  if (not tsos) {
    ATH_MSG_ERROR("Failed to retrieve track state on surfaces");
    return nullptr;
  }
  // loop over the associated hits in ITk in reverse order, since we want to
  // select the one closest to HGTD to start the extrapolation
  for (auto i = tsos->rbegin(); i != tsos->rend(); ++i) {
    const auto* curr_last_tsos = *i;
    if (not curr_last_tsos) {
      continue;
    }
    if (curr_last_tsos->type(Trk::TrackStateOnSurface::Measurement) and
        curr_last_tsos->trackParameters() and
        curr_last_tsos->measurementOnTrack()) {
      return curr_last_tsos->trackParameters();
    }
  }
  return nullptr;
}

getTimeCompatibleHits()

std::vector< TimeCompatibilityCheckAlg::Hit > HGTD::TimeCompatibilityCheckAlg::getTimeCompatibleHits ( const xAOD::TrackParticle * track_particle ) const

private

Definition at line 174 of file TimeCompatibilityCheckAlg.cxx.

                                                   {
// get all available hits (see the struct Hit) in a first step
auto valid_hits = getValidHits(track_particle);
 
size_t vts = valid_hits.size();
 
// if there is only one hit, not time consistency check can be done
if (vts <= 1) {
    return valid_hits;
}
 
// in case of two hits, check for time consistencyf
if (vts == 2) {
    if (passesDeltaT(valid_hits)) {
        return valid_hits;
    } else {
        // if times are too far away form each other, do not accept time
        return {};
    }
}
 
// if there are 3 or 4 hits, perfome chi2 outlier removal
// calculate the chi2 value of the available hits in a first step
float chi2 = calculateChi2(valid_hits);
 
// if the chi2 value does not surpass the set threshold, the hits are accepted as compatible in time
if (chi2 < m_chi2_threshold) {
    return valid_hits;
}
 
std::vector<Hit> time_candidates_copy = valid_hits; // TODO do I need this copy?
  bool searching = true;
  while (searching) {
    // calculate chi2 contribution of each value
    std::vector<float> chi2_contributions(time_candidates_copy.size(), 0.0);
    for (size_t i = 0; i < time_candidates_copy.size(); i++) {
      std::vector<Hit> buff = time_candidates_copy;
      buff.erase(buff.begin() + i);
 
      // calculate the chi2 value we would get when removing the i-th hit
      double local_chi2 = calculateChi2(buff);
 
      chi2_contributions.at(i) = local_chi2;
    }
    // if removing one of the hits gives a much smaller chi2, it should be
    // removed so find the position where the "local chi2" is the smallest, and
    // this is the hit that should be removed (since it gave a big
    // contribution)]
 
    // find minimum local chi2
    int position = std::distance(
        chi2_contributions.begin(),
        std::min_element(chi2_contributions.begin(), chi2_contributions.end()));
 
    // and remove it from the hits
    time_candidates_copy.erase(time_candidates_copy.begin() + position);
 
    // recompute chi2 value
    chi2 = calculateChi2(time_candidates_copy);
 
    // check for accepted chi2
    if (chi2 < m_chi2_threshold) {
      // if the threshold is now fulfilled, break out of the while loop
      searching = false;
    } 
    // if everything except 2 values has been removed, stop algo
    if (time_candidates_copy.size() == 2) {
      if (passesDeltaT(time_candidates_copy)) {
        return time_candidates_copy;
      } else {
       ATH_MSG_DEBUG("[TimeCompatibilityCheckAlg::getValidHits] times of hits are too far away");
        // if times are too far away, don't accept any TODO maybe accept one,
        // can the spatial chi2 be used?
       return {};
      }
    }
  }
  return time_candidates_copy;
}

getValidHits()

std::vector< TimeCompatibilityCheckAlg::Hit > HGTD::TimeCompatibilityCheckAlg::getValidHits ( const xAOD::TrackParticle * track_particle ) const

private

Definition at line 136 of file TimeCompatibilityCheckAlg.cxx.

                                                                                     {
 
  SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<float>> layerClusterTimeHandle(m_layerClusterTimeKey);
  std::vector<float> times = layerClusterTimeHandle(*track_particle);
 
  SG::ReadDecorHandle<xAOD::TrackParticleContainer, std::vector<char>> layerHasExtensionHandle(m_layerHasExtensionKey);
  std::vector<char> has_clusters = layerHasExtensionHandle(*track_particle);
    
  std::vector<Hit> valid_hits;
  valid_hits.reserve(4);
 
for (size_t i = 0; i < has_clusters.size(); i++){
    if (not has_clusters.at(i)) {
       ATH_MSG_DEBUG("[TimeCompatibilityCheckAlg::getValidHits] NO clusters found");
       continue;  
    }
    Hit newhit;
 
    newhit.m_time = times.at(i);
 
    valid_hits.push_back(newhit);
}
return valid_hits;
 
}

initialize()

StatusCode HGTD::TimeCompatibilityCheckAlg::initialize ( )

finaloverridevirtual

Definition at line 24 of file TimeCompatibilityCheckAlg.cxx.

                                                 {
 
  ATH_CHECK(m_extensionTool.retrieve());
  ATH_CHECK(m_trackParticleContainerKey.initialize());
  ATH_CHECK(m_layerHasExtensionKey.initialize());
  ATH_CHECK(m_layerExtensionChi2Key.initialize());
  ATH_CHECK(m_layerClusterRawTimeKey.initialize());
  ATH_CHECK(m_layerClusterTimeKey.initialize());
  ATH_CHECK(m_layerClusterTruthClassKey.initialize());
  ATH_CHECK(m_layerClusterShadowedKey.initialize());
  ATH_CHECK(m_layerClusterMergedKey.initialize());
  ATH_CHECK(m_layerPrimaryExpectedKey.initialize());
  ATH_CHECK(m_extrapXKey.initialize());
  ATH_CHECK(m_extrapYKey.initialize());
  ATH_CHECK(m_compatibleHitsTimesKey.initialize ());
  ATH_CHECK(m_lastHitInITkCutKey.initialize());
  ATH_CHECK(m_holesHGTDKey.initialize());
  ATH_CHECK(m_holesITkKey.initialize());
 
  return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::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.

isClonable()

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

                                           {
    return true;
  }

lastHitIsOnLastSurface()

bool HGTD::TimeCompatibilityCheckAlg::lastHitIsOnLastSurface ( const xAOD::TrackParticle & track_particle ) const

private

Definition at line 281 of file TimeCompatibilityCheckAlg.cxx.

                                                                                                    {
  const Trk::Track* track = track_particle.track();
  const Trk::TrackParameters* last_hit_param = getLastHitOnTrack(*track);
  double radius = hypot(last_hit_param->position().x(), last_hit_param->position().y());
  double abs_z = fabs(last_hit_param->position().z());
 
  if (abs_z > 2700) {
    return true;
  }
  if (radius < 350 and abs_z > 2400) {
    return true;
  }
  // region 2
  if (radius > 205 and radius < 350 and abs_z > 2100) {
    return true;
  }
  // region 3
  if (radius < 220 and abs_z > 2200) {
    return true;
  }
 
  if (radius < 140 and abs_z > 1890) {
    return true;
  }
 
  return false;
}

msg()

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

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

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::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.

passesDeltaT()

bool HGTD::TimeCompatibilityCheckAlg::passesDeltaT ( const std::vector< Hit > & hits ) const

private

Definition at line 162 of file TimeCompatibilityCheckAlg.cxx.

                                                                             {
  if (not(hits.size()==2)){
    return false;
  }
// pass if the distance in units of the resolution passes the cut
if (fabs(hits.at(0).m_time - hits.at(1).m_time) <
    m_delta_cut * hypot(hits.at(0).m_resolution, hits.at(1).m_resolution)) {
        return true;
    }
    return false;
}

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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::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.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonAlgorithm< Gaudi::Algorithm >::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 AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

                                                      {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::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< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_chi2_threshold

FloatProperty HGTD::TimeCompatibilityCheckAlg::m_chi2_threshold

private

Initial value:

{this, "Chi2Threshold", 1.5,
                        "Quality cut for decision to keep hits compatible in time"}

Definition at line 82 of file TimeCompatibilityCheckAlg.h.

                              {this, "Chi2Threshold", 1.5,
                        "Quality cut for decision to keep hits compatible in time"};

m_compatibleHitsTimesKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_compatibleHitsTimesKey {this, "HGTD_times_of_compatible_hits",m_trackParticleContainerKey, "HGTD_times_of_compatible_hits", "deco with a handle for the time compatible hits' times"}

private

Definition at line 71 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_times_of_compatible_hits",m_trackParticleContainerKey, "HGTD_times_of_compatible_hits", "deco with a handle for the time compatible hits' times"};

m_delta_cut

FloatProperty HGTD::TimeCompatibilityCheckAlg::m_delta_cut

private

Initial value:

{this, "DeltaTCut", 2.0,
                        "Upper limit for a cluster delta t cut"}

Definition at line 85 of file TimeCompatibilityCheckAlg.h.

                         {this, "DeltaTCut", 2.0,
                        "Upper limit for a cluster delta t cut"};

m_detStore

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

privateinherited

Pointer to StoreGate (detector store by default).

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

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

privateinherited

Pointer to StoreGate (event store by default).

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

m_extensionTool

ToolHandle<IHGTD_HolesITkTool> HGTD::TimeCompatibilityCheckAlg::m_extensionTool {this, "ITkHoles", "HGTD_IterativeExtensionTool/HGTD_IterativeExtensionTool", "Tool for calculating the holes in"}

private

Definition at line 53 of file TimeCompatibilityCheckAlg.h.

{this, "ITkHoles", "HGTD_IterativeExtensionTool/HGTD_IterativeExtensionTool", "Tool for calculating the holes in"};

m_extrapXKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_extrapXKey {this, "HGTD_extrap_x",m_trackParticleContainerKey, "HGTD_extrap_x", "deco with a handle for an x of extrap"}

private

Definition at line 66 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_extrap_x",m_trackParticleContainerKey, "HGTD_extrap_x", "deco with a handle for an x of extrap"};

m_extrapYKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_extrapYKey {this, "HGTD_extrap_y", m_trackParticleContainerKey, "HGTD_extrap_y", "deco with a handle for an y of extrap"}

private

Definition at line 67 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_extrap_y", m_trackParticleContainerKey, "HGTD_extrap_y", "deco with a handle for an y of extrap"};

m_holesHGTDKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_holesHGTDKey {this, "HGTD_holes", m_trackParticleContainerKey, "HGTD_holes", "deco with the holes on track in HGTD"}

private

Definition at line 68 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_holes", m_trackParticleContainerKey,  "HGTD_holes", "deco with the holes on track in HGTD"};

m_holesITkKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_holesITkKey {this, "HGTD_holes_in_ITk", m_trackParticleContainerKey, "HGTD_holes_in_ITk", "deco with a handle for the number of holes on track in ITk between the last hit on track and HGTD"}

private

Definition at line 73 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_holes_in_ITk", m_trackParticleContainerKey, "HGTD_holes_in_ITk", "deco with a handle for the number of holes on track in ITk between the last hit on track and HGTD"};

m_lastHitInITkCutKey

SG::WriteDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_lastHitInITkCutKey {this, "HGTD_last_hit_in_ITk_cut", m_trackParticleContainerKey, "HGTD_last_hit_in_ITk_cut", "deco with a handle for the last hit to be close to HGTDrequirement"}

private

Definition at line 72 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_last_hit_in_ITk_cut", m_trackParticleContainerKey, "HGTD_last_hit_in_ITk_cut", "deco with a handle for the last hit to be close to HGTDrequirement"};

m_layerClusterMergedKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerClusterMergedKey {this, "HGTD_cluster_merged", m_trackParticleContainerKey, "HGTD_cluster_merged", "deco with a handle for a merged cluster"}

private

Definition at line 64 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_cluster_merged", m_trackParticleContainerKey, "HGTD_cluster_merged", "deco with a handle for a merged cluster"};

m_layerClusterRawTimeKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerClusterRawTimeKey {this, "HGTD_cluster_raw_time", m_trackParticleContainerKey, "HGTD_cluster_raw_time", "deco with a handle for layer cluster raw time"}

private

Definition at line 60 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_cluster_raw_time", m_trackParticleContainerKey, "HGTD_cluster_raw_time", "deco with a handle for layer cluster raw time"};

m_layerClusterShadowedKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerClusterShadowedKey {this, "HGTD_cluster_shadowed", m_trackParticleContainerKey, "HGTD_cluster_shadowed", "deco with a handle for a shadowed cluster"}

private

Definition at line 63 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_cluster_shadowed", m_trackParticleContainerKey, "HGTD_cluster_shadowed", "deco with a handle for a shadowed cluster"};

m_layerClusterTimeKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerClusterTimeKey {this, "HGTD_cluster_time", m_trackParticleContainerKey, "HGTD_cluster_time", "deco with a handle for cluster time"}

private

Definition at line 61 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_cluster_time", m_trackParticleContainerKey, "HGTD_cluster_time", "deco with a handle for cluster time"};

m_layerClusterTruthClassKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerClusterTruthClassKey {this, "HGTD_cluster_truth_class", m_trackParticleContainerKey, "HGTD_cluster_truth_class", "deco with a handle for a truth time"}

private

Definition at line 62 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_cluster_truth_class", m_trackParticleContainerKey,  "HGTD_cluster_truth_class", "deco with a handle for a truth time"};

m_layerExtensionChi2Key

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerExtensionChi2Key {this, "HGTD_extension_chi2", m_trackParticleContainerKey, "HGTD_extension_chi2", "deco with a handle for a ch2 of extension"}

private

Definition at line 59 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_extension_chi2", m_trackParticleContainerKey, "HGTD_extension_chi2", "deco with a handle for a ch2 of extension"};

m_layerHasExtensionKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerHasExtensionKey {this, "HGTD_has_extension", m_trackParticleContainerKey, "HGTD_has_extension", "deco with a handle for an extension"}

private

Definition at line 58 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_has_extension", m_trackParticleContainerKey, "HGTD_has_extension", "deco with a handle for an extension"};

m_layerPrimaryExpectedKey

SG::ReadDecorHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_layerPrimaryExpectedKey {this, "HGTD_primary_expected", m_trackParticleContainerKey, "HGTD_primary_expected", "deco with a handle for an expected primary"}

private

Definition at line 65 of file TimeCompatibilityCheckAlg.h.

{this, "HGTD_primary_expected", m_trackParticleContainerKey, "HGTD_primary_expected", "deco with a handle for an expected primary"};

m_trackParticleContainerKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> HGTD::TimeCompatibilityCheckAlg::m_trackParticleContainerKey {this, "TrackParticleContainerName", "InDetTrackParticles", "Name of the TrackParticle container"}

private

Definition at line 57 of file TimeCompatibilityCheckAlg.h.

{this, "TrackParticleContainerName", "InDetTrackParticles", "Name of the TrackParticle container"};

m_varHandleArraysDeclared

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

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

The documentation for this class was generated from the following files:

• TimeCompatibilityCheckAlg.h
• TimeCompatibilityCheckAlg.cxx