RecParticle.hpp

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <cstdint>
#include <vector>
 
#include "SHiP/RecHit.hpp"
#include "SHiP/SimParticle.hpp"
 
namespace SHiP {
 
struct RecParticle {
  std::int32_t track_id{0};
  std::int32_t parent_id{0};
  std::int32_t pdg_code{0};
  std::array<double, 3> vertex{0, 0, 0};    
  std::array<double, 3> endpoint{0, 0, 0};  
  std::array<double, 3> momentum{0, 0, 0};  
  std::vector<RecHit> hits{};
  double energy{0};  
  double time{0};    
  std::int32_t creator_process{0};
  double ip_pv{0};  
 
  bool operator==(RecParticle const&) const = default;
};
 
inline RecParticle from_sim_particle(SimParticle const& sp) {
  return {
      .track_id = sp.track_id,
      .parent_id = sp.parent_id,
      .pdg_code = sp.pdg_code,
      .vertex = sp.vertex,
      .endpoint = sp.endpoint,
      .momentum = sp.momentum,
      .energy = sp.energy,
      .time = sp.time,
      .creator_process = sp.creator_process,
      .ip_pv = 0.0,
  };
}
 
}  // namespace SHiP