ParticleSmearing.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2019 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/Units.hpp"
#include "Acts/EventData/ParticleHypothesis.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/EventData/SimParticle.hpp"
#include "ActsExamples/EventData/Track.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/IAlgorithm.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/Framework/RandomNumbers.hpp"

#include <array>
#include <limits>
#include <memory>
#include <optional>
#include <string>

namespace ActsExamples {
class RandomNumbers;
struct AlgorithmContext;

class ParticleSmearing final : public IAlgorithm {
 public:
  struct Config {
    std::string inputParticles;
    std::string outputTrackParameters;
    double sigmaD0 = 20 * Acts::UnitConstants::um;
    double sigmaD0PtA = 30 * Acts::UnitConstants::um;
    double sigmaD0PtB = 0.3 / Acts::UnitConstants::GeV;
    double sigmaZ0 = 20 * Acts::UnitConstants::um;
    double sigmaZ0PtA = 30 * Acts::UnitConstants::um;
    double sigmaZ0PtB = 0.3 / Acts::UnitConstants::GeV;
    double sigmaT0 = 1 * Acts::UnitConstants::ns;
    double sigmaPhi = 1 * Acts::UnitConstants::degree;
    double sigmaTheta = 1 * Acts::UnitConstants::degree;
    double sigmaPRel = 0.05;
    std::optional<std::array<double, 6>> initialSigmas = std::nullopt;
    std::array<double, 6> initialVarInflation = {1., 1., 1., 1., 1., 1.};
    std::shared_ptr<const RandomNumbers> randomNumbers = nullptr;
    std::optional<Acts::ParticleHypothesis> particleHypothesis = std::nullopt;
  };

  ParticleSmearing(const Config& config, Acts::Logging::Level level);

  ProcessCode execute(const AlgorithmContext& ctx) const override;

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

 private:
  Config m_cfg;

  ReadDataHandle<SimParticleContainer> m_inputParticles{this, "InputParticles"};

  WriteDataHandle<TrackParametersContainer> m_outputTrackParameters{
      this, "OutputTrackParameters"};
};

}  // namespace ActsExamples