sts/RecoSetupUnit.cxx

Go to the documentation of this file.

/* Copyright (C) 2025 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergei Zharko [committer] */

 
#include "algo/detectors/sts/RecoSetupUnit.h"
 
#include <limits>
#include <unordered_set>
 
using cbm::algo::sts::RecoSetupUnit;
 
// ---------------------------------------------------------------------------------------------------------------------
//
RecoSetupUnit::RecoSetupUnit(std::pair<std::vector<GeoVolume>, std::vector<GeoVolume>>&& stationVolumes,
                             StationIdMap_t&& stationIdMap)
  : TrackingSetupIfs<RecoSetupUnit>(std::move(stationVolumes))
  , fvStationIdMap(std::move(stationIdMap))
{
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
bool RecoSetupUnit::IsInitialized() const
{
  // STS tracking detector is initialized, if:
  // - the fvStationIdMap has elements
  // - the fvStationIdMap contains valid station indices in a range [0, GetNofTrackingStations() - 1]
 
  bool res = true;
  if (!fvStationIdMap.size()) {
    LOG(error) << "\t - station ID map is empty";
    res = false;
  }
 
 
  // Collect the indices of tracking stations, which are stored in fvStationIdMap, in a set; check, if the set
  // has a size == GetNofTrackingStations(), and if the all the station IDs are in the range
  // [0, GetNofTrackingStations() - 1].
  std::unordered_set<int> knownStationIds;
  for (auto item : fvStationIdMap) {
    knownStationIds.insert(item.second);
    if (item.first < std::numeric_limits<decltype(item.first)>::max()) {
      knownStationIds.insert(item.second + 1);
    }
  }
 
  if (knownStationIds.size() != static_cast<size_t>(GetNofTrackingStations())) {
    LOG(error) << "\t - number of known tracking station IDs (" << knownStationIds.size() << ") differs from the "
               << "number of stations, provided by geometry (" << GetNofTrackingStations() << ')';
    res = false;
  }
 
  for (int stationId : knownStationIds) {
    if (stationId < 0 || stationId >= GetNofTrackingStations()) {
      LOG(error) << "\t - found a tracking station < 0 or >= " << GetNofTrackingStations()
                 << " (geometry): " << stationId;
      res = false;
    }
  }
  return res;
}