TrackTimeExtensionAlg.cxx

Go to the documentation of this file.

 
#include "TrackTimeExtensionAlg.h"
 
#include "HGTD_RIO_OnTrack/HGTD_ClusterOnTrack.h"
#include "StoreGate/WriteDecorHandle.h"
#include "xAODTruth/TruthParticleContainer.h"
 
namespace {
static constexpr unsigned short n_hgtd_layers = 4;
}
 
namespace HGTD {
 
TrackTimeExtensionAlg::TrackTimeExtensionAlg(const std::string& name,
                                             ISvcLocator* pSvcLocator)
    : AthReentrantAlgorithm(name, pSvcLocator) {}
 
StatusCode TrackTimeExtensionAlg::initialize() {
 
  ATH_CHECK(m_extension_tool.retrieve());
 
  ATH_CHECK(m_truth_tool.retrieve());
 
  ATH_CHECK(m_clustercont_rh_key.initialize());
  ATH_CHECK(m_sdo_coll_rh_key.initialize());
  ATH_CHECK(m_mc_coll_rh_key.initialize());
  ATH_CHECK(m_trk_ptkl_rh_key.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());
 
 
  return StatusCode::SUCCESS;
}

StatusCode TrackTimeExtensionAlg::execute(const EventContext& ctx) const {
 
  ATH_MSG_DEBUG("Start event");
 
  SG::ReadHandle<HGTD_ClusterContainer> cluster_container_handle(
      m_clustercont_rh_key, ctx);
  const HGTD_ClusterContainer* cluster_container =
      cluster_container_handle.cptr();
  if (not cluster_container) {
    ATH_MSG_ERROR("[TrackTimeExtensionAlg] HGTD_ClusterContainer not found, "
                  "aborting execute!");
    return StatusCode::FAILURE;
  }
 
  SG::ReadHandle<InDetSimDataCollection> sdo_collection_handle(
      m_sdo_coll_rh_key, ctx);
  const InDetSimDataCollection* sdo_collection = sdo_collection_handle.cptr();
  if (not sdo_collection) {
    ATH_MSG_WARNING("[TrackTimeExtensionAlg] SDO Collection not found, no "
                    "truth info available!");
  }
 
  SG::ReadHandle<McEventCollection> mc_collection_handle(m_mc_coll_rh_key, ctx);
  const McEventCollection* mc_coll = mc_collection_handle.cptr();
  if (not mc_coll) {
    ATH_MSG_ERROR("[TrackTimeExtensionAlg] McEventCollection not found, "
                  "aborting execute!");
    return StatusCode::FAILURE;
  }
 
  const HepMC::GenEvent* hs_event = nullptr;
  if (mc_coll->size() > 0) {
    hs_event = mc_coll->at(0);
  } else {
    hs_event = nullptr;
  }
 
  // for each track, run the extension if the track is in HGTD acceptance
  SG::ReadHandle<xAOD::TrackParticleContainer> trk_ptkl_container_handle(
      m_trk_ptkl_rh_key, ctx);
  const xAOD::TrackParticleContainer* track_particles =
      trk_ptkl_container_handle.cptr();
  if (not track_particles) {
    ATH_MSG_ERROR("[TrackTimeExtensionAlg] TrackParticleContainer not found, "
                  "aborting execute!");
    return StatusCode::FAILURE;
  }
 
  DecorHandles dh (*this, ctx);
 
  for (const auto* track_ptkl : *track_particles) {
 
    ATH_MSG_DEBUG("Track eta: " << track_ptkl->eta()
                                << " pt: " << track_ptkl->pt());
 
    HGTD::ExtensionObject extension;
 
    if (std::abs(track_ptkl->eta()) < m_eta_cut) {
      ATH_MSG_DEBUG("Track out of acceptance");
      // decorate all track particle objects to avoid issues with the
      // decorations
      ATH_CHECK(decorateTrackParticle(dh, track_ptkl, extension, sdo_collection,
                                      hs_event, true));
      continue;
    }
 
    // this should not happen?
    if (track_ptkl->track() == nullptr) {
      ATH_MSG_DEBUG("There is no Trk::Track");
      ATH_CHECK(decorateTrackParticle(dh, track_ptkl, extension, sdo_collection,
                                      hs_event, true));
      continue;
    }
 
    // return 4 track states on surface objects as a result of the extension
    extension = m_extension_tool->extendTrackToHGTD(ctx,
                                                    *track_ptkl,
                                                    cluster_container,
                                                    hs_event,
                                                    sdo_collection);
 
    // TODO here:
    // retrieve truth info for associated and non-associated HGTD hits, merging
    // or shadowing info
 
    // decorate the track
    ATH_CHECK(decorateTrackParticle(dh, track_ptkl, extension, sdo_collection,
                                    hs_event, false));
 
  } // END LOOP over tracks
 
  return StatusCode::SUCCESS;
}

 
TrackTimeExtensionAlg::DecorHandles::DecorHandles
  (const TrackTimeExtensionAlg& tool, const EventContext& ctx)
    : layerHasExtensionHandle (tool.m_layerHasExtensionKey, ctx),
      layerExtensionChi2Handle (tool.m_layerExtensionChi2Key, ctx),
      layerClusterRawTimeHandle (tool.m_layerClusterRawTimeKey, ctx),
      layerClusterTimeHandle (tool.m_layerClusterTimeKey, ctx),
      layerClusterTruthClassHandle (tool.m_layerClusterTruthClassKey, ctx),
      layerClusterShadowedHandle (tool.m_layerClusterShadowedKey, ctx),
      layerClusterMergedHandle (tool.m_layerClusterMergedKey, ctx),
      layerPrimaryExpectedHandle (tool.m_layerPrimaryExpectedKey, ctx),
      extrapXHandle (tool.m_extrapXKey, ctx),
      extrapYHandle (tool.m_extrapYKey, ctx)
{
}
 
StatusCode TrackTimeExtensionAlg::decorateTrackParticle(
    DecorHandles& dh,
    const xAOD::TrackParticle* track_ptkl, const HGTD::ExtensionObject& extension,
    const InDetSimDataCollection* sdo_collection,
    const HepMC::GenEvent* hs_event, bool skip_deco) const {
 
  std::vector<bool> has_cluster_vec;
  has_cluster_vec.reserve(n_hgtd_layers);
  std::vector<float> chi2_vec;
  chi2_vec.reserve(n_hgtd_layers);
  std::vector<float> raw_time_vec;
  raw_time_vec.reserve(n_hgtd_layers);
  std::vector<float> time_vec;
  time_vec.reserve(n_hgtd_layers);
  std::vector<int> truth_vec;
  truth_vec.reserve(n_hgtd_layers);
  std::vector<bool> is_shadowed_vec;
  is_shadowed_vec.reserve(n_hgtd_layers);
  std::vector<bool> is_merged_vec;
  is_merged_vec.reserve(n_hgtd_layers);
  std::vector<bool> primary_exists_vec;
  primary_exists_vec.reserve(n_hgtd_layers);
 
  for (unsigned short i = 0; i < n_hgtd_layers; i++) {
 
    const std::unique_ptr<const Trk::TrackStateOnSurface>& trk_state =
        extension.m_hits.at(i);
    const HGTD_Cluster* primary_cluster = extension.m_truth_primary_hits.at(i);
 
    primary_exists_vec.push_back(primary_cluster != nullptr);
 
    if (trk_state) {
      ATH_MSG_DEBUG("[decorateTrackParticle] extension found");
      has_cluster_vec.emplace_back(true);
 
      chi2_vec.emplace_back(
          trk_state->fitQualityOnSurface().chiSquared() /
          trk_state->fitQualityOnSurface().doubleNumberDoF());
 
      const HGTD_ClusterOnTrack* cot =
          dynamic_cast<const HGTD_ClusterOnTrack*>(trk_state->measurementOnTrack());
 
      time_vec.emplace_back(cot->time());
 
      // get the cluster
      const HGTD_Cluster* cluster = cot->prepRawData();
 
      raw_time_vec.emplace_back(cluster->time());
 
      // get the truth particle
      static const SG::AuxElement::Accessor<ElementLink<xAOD::TruthParticleContainer>>
          acc_tpl("truthParticleLink");
      const xAOD::TruthParticle* truth_particle = nullptr;
      if (acc_tpl.isAvailable(*track_ptkl)) {
        const auto& truth_match_link = acc_tpl(*track_ptkl);
        if (truth_match_link.isValid()) {
          truth_particle = *truth_match_link;
        }
      }
 
      ClusterTruthOrigin truth_origin = ClusterTruthOrigin::UNIDENTIFIED;
      bool is_shadowed = false;
      bool is_merged = false;
 
      if (truth_particle) {
        auto truth_info = m_truth_tool->classifyCluster(
            cluster, truth_particle, sdo_collection, hs_event);
        truth_origin = truth_info.origin;
        is_shadowed = truth_info.is_shadowed;
        is_merged = truth_info.is_merged;
      }
 
      ATH_MSG_DEBUG("Truth origin: " << (int)truth_origin);
 
      truth_vec.emplace_back((int)truth_origin);
      is_shadowed_vec.emplace_back(is_shadowed);
      is_merged_vec.emplace_back(is_merged);
    } else {
      ATH_MSG_DEBUG("[decorateTrackParticle] NO extension found");
      has_cluster_vec.emplace_back(false);
      if (not skip_deco) {
        chi2_vec.emplace_back(-1.);
        raw_time_vec.emplace_back(-1.);
        time_vec.emplace_back(-1.);
        truth_vec.emplace_back(-1);
        is_shadowed_vec.emplace_back(false);
        is_merged_vec.emplace_back(false);
      }
    }
 
  } // END LOOP over TrackStateOnSurface
 
  dh.layerHasExtensionHandle(*track_ptkl) = has_cluster_vec;
  dh.layerExtensionChi2Handle(*track_ptkl) = chi2_vec;
  dh.layerClusterRawTimeHandle(*track_ptkl) = raw_time_vec;
  dh.layerClusterTimeHandle(*track_ptkl) = time_vec;
  dh.layerClusterTruthClassHandle(*track_ptkl) = truth_vec;
  dh.layerClusterShadowedHandle(*track_ptkl) = is_shadowed_vec;
  dh.layerClusterMergedHandle(*track_ptkl) = is_merged_vec;
  dh.layerPrimaryExpectedHandle(*track_ptkl) = primary_exists_vec;
  dh.extrapXHandle(*track_ptkl) = extension.m_extrap_x;
  dh.extrapYHandle(*track_ptkl) = extension.m_extrap_y;
 
  return StatusCode::SUCCESS;
}
 
} // namespace HGTD