Process.hpp

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

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

#include <cmath>

namespace Fatras {

template <typename physics_t, typename selector_in_t, typename selector_out_t,
          typename selector_child_t>

struct Process {

  physics_t process;

  selector_in_t selectorIn;
  selector_out_t selectorOut;
  selector_child_t selectorChild;

  template <typename generator_t, typename detector_t, typename particle_t>
  bool operator()(generator_t &gen, const detector_t &det, particle_t &in,
                  std::vector<particle_t> &out) const {
    // check if the process applies
    if (selectorIn(det, in)) {
      // apply energy loss and get eventual children
      auto children = process(gen, det, in);
      if (children.size()) {
        // copy the children that comply with the child selector
        std::copy_if(
            children.begin(), children.end(), std::back_inserter(out),
            [this, det](const particle_t &p) { return selectorChild(det, p); });
      }
    }
    // check if this killed the particle,
    // or pushed below threshold
    return (!selectorOut(det, in));
  }
};

} // namespace Fatras