SteppingLogger.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2016-2018 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/Algebra.hpp"
#include "Acts/Definitions/Common.hpp"
#include "Acts/Definitions/Direction.hpp"
#include "Acts/Geometry/GeometryIdentifier.hpp"
#include "Acts/Propagator/ConstrainedStep.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Logger.hpp"

#include <memory>
#include <vector>

namespace Acts {

class Surface;

namespace detail {

struct Step {
  ConstrainedStep stepSize;
  Direction navDir;
  Vector3 position = Vector3(0., 0., 0.);
  Vector3 momentum = Vector3(0., 0., 0.);
  std::shared_ptr<const Surface> surface = nullptr;
  GeometryIdentifier geoID = 0;
};

struct SteppingLogger {
  struct this_result {
    std::vector<Step> steps;
  };

  using result_type = this_result;

  bool sterile = false;

  template <typename propagator_state_t, typename stepper_t,
            typename navigator_t>
  void operator()(propagator_state_t& state, const stepper_t& stepper,
                  const navigator_t& navigator, result_type& result,
                  const Logger& /*logger*/) const {
    // Don't log if you have reached the target or are sterile
    if (sterile or navigator.targetReached(state.navigation)) {
      return;
    }
    // Record the propagation state
    Step step;
    step.stepSize = state.stepping.stepSize;
    step.navDir = state.options.direction;
    step.position = stepper.position(state.stepping);
    step.momentum = stepper.momentum(state.stepping);

    // Record the information about the surface
    if (navigator.currentSurface(state.navigation) != nullptr) {
      step.surface = navigator.currentSurface(state.navigation)->getSharedPtr();
      step.geoID = step.surface->geometryId();
    } else if (navigator.currentVolume(state.navigation) != nullptr) {
      // If there's no surface but a volume, this sets the geoID
      step.geoID = navigator.currentVolume(state.navigation)->geometryId();
    }
    result.steps.push_back(std::move(step));
  }
};

}  // namespace detail
}  // namespace Acts