CbmMvdRecoSetupUnitFactory.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 "CbmMvdRecoSetupUnitFactory.h"
 
#include "CbmMvdAddress.h"
#include "CbmMvdDetector.h"
#include "CbmMvdSensor.h"
#include "CbmMvdStationPar.h"
 
#include <Logger.h>
 
#include <TString.h>
 
#include <set>
#include <sstream>
#include <utility>
 
using cbm::algo::GeoVolume;
using cbm::algo::mvd::RecoSetupUnit;
using cbm::mvd::RecoSetupUnitFactory;
 
 
// ---------------------------------------------------------------------------------------------------------------------
//
RecoSetupUnit RecoSetupUnitFactory::Create() const
{
  RecoSetupUnit res{CreateStationVolumes(), CreateStationIdMap()};
  if (!res.IsValid()) {
    LOG(fatal) << "mvd::RecoSetupUnitFactory: creation of a valid MVD reco-setup unit failed";
  }
  return res;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
RecoSetupUnitFactory::StationIdMap_t RecoSetupUnitFactory::CreateStationIdMap() const
{
  StationIdMap_t res;
  const auto& mvdSensorMap = CbmMvdDetector::Instance()->GetSensorMap();
  res.resize(mvdSensorMap.size(), -1);
  for (auto [dummy, pSensor] : mvdSensorMap) {
    res[pSensor->GetSensorNr()] = pSensor->GetStationNr();
  }
  return res;
}
 
// ---------------------------------------------------------------------------------------------------------------------
//
std::pair<std::vector<GeoVolume>, std::vector<GeoVolume>> RecoSetupUnitFactory::CreateStationVolumes() const
{
  // FIXME: Properly get full and active volumes from geometry (see TRD and TOF setup unit factories for example)
  auto* mvdStationPar = CbmMvdDetector::Instance()->GetParameterFile();
  assert(mvdStationPar);
 
  int nStations = mvdStationPar->GetStationCount();
  std::vector<GeoVolume> fullVolumes;
  fullVolumes.reserve(nStations);
  for (int iSt = 0; iSt < nStations; ++iSt) {
    fullVolumes.emplace_back(-mvdStationPar->GetWidth(iSt), +mvdStationPar->GetWidth(iSt),
                             -mvdStationPar->GetHeight(iSt), +mvdStationPar->GetHeight(iSt),
                             mvdStationPar->GetZPosition(iSt) - 0.5 * mvdStationPar->GetZThickness(iSt),
                             mvdStationPar->GetZPosition(iSt) + 0.5 * mvdStationPar->GetZThickness(iSt));
  }
  std::vector<GeoVolume> activeVolumes = fullVolumes;
 
  return {std::move(fullVolumes), std::move(activeVolumes)};
}