CsvTrackParameterReader.cpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CsvTrackParameterReader.cpp

// This file is part of the Acts project.
//
// Copyright (C) 2022 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#include "ActsExamples/Io/Csv/CsvTrackParameterReader.hpp"

#include "Acts/Definitions/Algebra.hpp"
#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/Surfaces/PerigeeSurface.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "ActsExamples/EventData/Track.hpp"
#include "ActsExamples/Framework/AlgorithmContext.hpp"
#include "ActsExamples/Utilities/Paths.hpp"

#include <algorithm>
#include <stdexcept>
#include <string>

#include <dfe/dfe_io_dsv.hpp>

#include "CsvOutputData.hpp"

ActsExamples::CsvTrackParameterReader::CsvTrackParameterReader(
    const ActsExamples::CsvTrackParameterReader::Config& config,
    Acts::Logging::Level level)
    : m_cfg(config),
      m_eventsRange(
          determineEventFilesRange(m_cfg.inputDir, m_cfg.inputStem + ".csv")),
      m_logger(Acts::getDefaultLogger("CsvTrackParameterReader", level)) {
  if (m_cfg.inputStem.empty()) {
    throw std::invalid_argument("Missing input filename stem");
  }
  if (m_cfg.outputTrackParameters.empty()) {
    throw std::invalid_argument("Missing output collection");
  }

  m_outputTrackParameters.initialize(m_cfg.outputTrackParameters);
}

std::string
ActsExamples::CsvTrackParameterReader::CsvTrackParameterReader::name() const {
  return "CsvTrackParameterReader";
}

std::pair<size_t, size_t>
ActsExamples::CsvTrackParameterReader::availableEvents() const {
  return m_eventsRange;
}

ActsExamples::ProcessCode ActsExamples::CsvTrackParameterReader::read(
    const ActsExamples::AlgorithmContext& ctx) {
  TrackParametersContainer trackParameters;

  auto surface = Acts::Surface::makeShared<Acts::PerigeeSurface>(
      Acts::Vector3(m_cfg.beamspot[0], m_cfg.beamspot[1], m_cfg.beamspot[2]));

  auto path = perEventFilepath(m_cfg.inputDir, m_cfg.inputStem + ".csv",
                               ctx.eventNumber);
  dfe::NamedTupleCsvReader<TrackParameterData> reader(path);
  TrackParameterData d{};

  while (reader.read(d)) {
    Acts::BoundVector params = Acts::BoundVector::Zero();
    params[Acts::eBoundLoc0] = d.d0;
    params[Acts::eBoundLoc1] = d.z0;
    params[Acts::eBoundPhi] = d.phi;
    params[Acts::eBoundTheta] = d.theta;
    params[Acts::eBoundQOverP] = d.qop;

    Acts::BoundSquareMatrix cov = Acts::BoundSquareMatrix::Zero();
    cov(Acts::eBoundLoc0, Acts::eBoundLoc0) = d.var_d0;
    cov(Acts::eBoundLoc1, Acts::eBoundLoc1) = d.var_z0;
    cov(Acts::eBoundPhi, Acts::eBoundPhi) = d.var_phi;
    cov(Acts::eBoundTheta, Acts::eBoundTheta) = d.var_theta;
    cov(Acts::eBoundQOverP, Acts::eBoundQOverP) = d.var_qop;
    cov(Acts::eBoundTime, Acts::eBoundTime) = 1;

    cov(Acts::eBoundLoc0, Acts::eBoundLoc1) = d.cov_d0z0;
    cov(Acts::eBoundLoc0, Acts::eBoundPhi) = d.cov_d0phi;
    cov(Acts::eBoundLoc0, Acts::eBoundTheta) = d.cov_d0theta;
    cov(Acts::eBoundLoc0, Acts::eBoundQOverP) = d.cov_d0qop;

    cov(Acts::eBoundLoc1, Acts::eBoundLoc0) = d.cov_z0d0;
    cov(Acts::eBoundLoc1, Acts::eBoundPhi) = d.cov_z0phi;
    cov(Acts::eBoundLoc1, Acts::eBoundTheta) = d.cov_z0theta;
    cov(Acts::eBoundLoc1, Acts::eBoundQOverP) = d.cov_z0qop;

    cov(Acts::eBoundPhi, Acts::eBoundLoc0) = d.cov_phid0;
    cov(Acts::eBoundPhi, Acts::eBoundLoc1) = d.cov_phiz0;
    cov(Acts::eBoundPhi, Acts::eBoundTheta) = d.cov_phitheta;
    cov(Acts::eBoundPhi, Acts::eBoundQOverP) = d.cov_phiqop;

    cov(Acts::eBoundTheta, Acts::eBoundLoc0) = d.cov_thetad0;
    cov(Acts::eBoundTheta, Acts::eBoundLoc1) = d.cov_thetaz0;
    cov(Acts::eBoundTheta, Acts::eBoundPhi) = d.cov_thetaphi;
    cov(Acts::eBoundTheta, Acts::eBoundQOverP) = d.cov_thetaqop;

    cov(Acts::eBoundQOverP, Acts::eBoundLoc0) = d.cov_qopd0;
    cov(Acts::eBoundQOverP, Acts::eBoundLoc1) = d.cov_qopz0;
    cov(Acts::eBoundQOverP, Acts::eBoundPhi) = d.cov_qopphi;
    cov(Acts::eBoundQOverP, Acts::eBoundTheta) = d.cov_qoptheta;

    // TODO we do not have a hypothesis at hand here. defaulting to pion
    trackParameters.emplace_back(surface, params, cov,
                                 Acts::ParticleHypothesis::pion());
  }

  m_outputTrackParameters(ctx, std::move(trackParameters));

  return ProcessCode::SUCCESS;
}