From 4bcd7062d6d8f1e800e4b072c09dfa2aa571ee4e Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Fri, 24 Feb 2023 20:44:41 +0100
Subject: [PATCH] add hepevt to header
---
modules/fluka/FLUKA.hpp | 44 +++++++++++++++++++++++++++++++++++++
tests/modules/testFluka.cpp | 2 +-
2 files changed, 45 insertions(+), 1 deletion(-)
diff --git a/modules/fluka/FLUKA.hpp b/modules/fluka/FLUKA.hpp
index 76c542b03..b7f9c3e73 100644
--- a/modules/fluka/FLUKA.hpp
+++ b/modules/fluka/FLUKA.hpp
@@ -9,7 +9,51 @@
#pragma once
namespace fluka {
+ size_t constexpr nmxhep = 10000;
+ template <typename T>
+ using hepmc_array = std::array<T, nmxhep>;
+
extern "C" {
+ /*
+ c---------------------------------------------------------------------------
+ c hep standard event commonblock.
+ c---------------------------------------------------------------------------
+ c---------------------------------------------------------------------------
+ c
+ c nevhep - event number
+ c nhep - number of entries in the event record
+ c
+ c isthep(i) - status code
+ c idhep(i) - particle id (particle data group standard)
+ c
+ c jmohep(1,i) - position of mother particle in list
+ c jmohep(2,i) - position of second mother particle in list
+ c jdahep(1,i) - position of first daughter in list
+ c jdahep(2,i) - position of first daughter in list
+ c
+ c phep(1,i) - p_x momentum in gev/c
+ c phep(2,i) - p_y momentum in gev/c
+ c phep(3,i) - p_z momentum in gev/c
+ c phep(4,i) - energy in gev
+ c phep(5,i) - mass in gev/c**2
+ c
+ c vhep(1,i) - x position of production vertex in mm
+ c vhep(2,i) - y position of production vertex in mm
+ c vhep(3,i) - z position of production vertex in mm
+ c vhep(4,i) - time of production in mm/c
+ */
+
+ extern struct {
+ int nevhep; // event number
+ int nhep; // number of entries
+ hepmc_array<int> isthep; // status code
+ hepmc_array<int> idhep; // PDG particle id
+ hepmc_array<int[2]> jmohep; // position of first, second mother
+ hepmc_array<int[2]> jdahep; // position of first, last daughter
+ hepmc_array<double[5]> phep; // 4-momemtum, mass (GeV)
+ hepmc_array<double[4]> vhep; // vertex, production time in mm
+ } hepevt_;
+
/* Iflxyz must be consistent in all calls, to Stpxyz, Sgmxyz, Evtxyz
* Iflxyz = 1 -> only inelastic
* Iflxyz = 10 -> only elastic
diff --git a/tests/modules/testFluka.cpp b/tests/modules/testFluka.cpp
index c154cbbde..818ce0963 100644
--- a/tests/modules/testFluka.cpp
+++ b/tests/modules/testFluka.cpp
@@ -79,7 +79,7 @@ TEST_CASE("FLUKA") {
};
HEPEnergyType const p = 100_GeV;
- for (auto const& [projectileCode, targetCode] : combinations) {
+ for (auto const& [projectileCode, targetCode] : combinationsOK) {
auto const projectile4mom =
FourVector{calculate_total_energy(p, get_mass(projectileCode)),
MomentumVector{cs, 0_eV, 0_eV, p}};
--
GitLab