TrackParamsEstimationAlgorithm.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2021 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/.

#pragma once

#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/Definitions/Units.hpp"
#include "Acts/EventData/ParticleHypothesis.hpp"
#include "Acts/Geometry/TrackingGeometry.hpp"
#include "Acts/MagneticField/ConstantBField.hpp"
#include "Acts/MagneticField/InterpolatedBFieldMap.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/EventData/ProtoTrack.hpp"
#include "ActsExamples/EventData/SimSeed.hpp"
#include "ActsExamples/EventData/Track.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/IAlgorithm.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/MagneticField/MagneticField.hpp"

#include <algorithm>
#include <array>
#include <functional>
#include <memory>
#include <string>
#include <vector>

namespace Acts {
class TrackingGeometry;
class MagneticFieldProvider;
}  // namespace Acts

namespace ActsExamples {
struct AlgorithmContext;

class TrackParamsEstimationAlgorithm final : public IAlgorithm {
 public:
  struct Config {
    std::string inputSeeds;
    std::string inputProtoTracks;
    std::string outputTrackParameters;
    std::string outputSeeds;
    std::string outputProtoTracks;
    std::shared_ptr<const Acts::TrackingGeometry> trackingGeometry;
    std::shared_ptr<const Acts::MagneticFieldProvider> magneticField;
    double bFieldMin = 0.1 * Acts::UnitConstants::T;
    std::array<double, 6> initialSigmas = {
        25 * Acts::UnitConstants::um,       100 * Acts::UnitConstants::um,
        0.02 * Acts::UnitConstants::degree, 0.02 * Acts::UnitConstants::degree,
        0.1 / Acts::UnitConstants::GeV,     10 * Acts::UnitConstants::ns};
    std::array<double, 6> initialVarInflation = {1., 1., 1., 1., 1., 1.};
    Acts::ParticleHypothesis particleHypothesis =
        Acts::ParticleHypothesis::pion();
  };

  TrackParamsEstimationAlgorithm(Config cfg, Acts::Logging::Level lvl);

  ProcessCode execute(const AlgorithmContext& ctx) const override;

  const Config& config() const { return m_cfg; }

 private:
  Config m_cfg;

  Acts::BoundSquareMatrix m_covariance = Acts::BoundSquareMatrix::Zero();

  ReadDataHandle<SimSeedContainer> m_inputSeeds{this, "InputSeeds"};
  ReadDataHandle<ProtoTrackContainer> m_inputTracks{this, "InputTracks"};

  WriteDataHandle<TrackParametersContainer> m_outputTrackParameters{
      this, "OutputTrackParameters"};
  WriteDataHandle<SimSeedContainer> m_outputSeeds{this, "OutputSeeds"};
  WriteDataHandle<ProtoTrackContainer> m_outputTracks{this, "OutputTracks"};
};

}  // namespace ActsExamples