From 61461a8859f02afedebed496c9adba9dcfb257dd Mon Sep 17 00:00:00 2001 From: Wouter Deconinck <wdconinc@gmail.com> Date: Mon, 26 Sep 2022 18:16:52 -0500 Subject: [PATCH] fix: ddsim --gun.distribution eta should set eta_min, eta_max correctly --- DDG4/src/Geant4IsotropeGenerator.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/DDG4/src/Geant4IsotropeGenerator.cpp b/DDG4/src/Geant4IsotropeGenerator.cpp index f7d132c6c..b549442f7 100644 --- a/DDG4/src/Geant4IsotropeGenerator.cpp +++ b/DDG4/src/Geant4IsotropeGenerator.cpp @@ -76,8 +76,8 @@ void Geant4IsotropeGenerator::getParticleDirectionEta(int, ROOT::Math::XYZVector // See https://en.wikipedia.org/wiki/Pseudorapidity const double dmin = std::numeric_limits<double>::epsilon(); double phi = m_phiMin+(m_phiMax-m_phiMin)*rnd.rndm(); - double eta_min = Distribution::eta(m_thetaMin>dmin ? m_thetaMin : dmin); - double eta_max = Distribution::eta(m_thetaMax>(M_PI-dmin) ? m_thetaMax : M_PI-dmin); + double eta_max = Distribution::eta(m_thetaMin>dmin ? m_thetaMin : dmin); + double eta_min = Distribution::eta(m_thetaMax>(M_PI-dmin) ? M_PI-dmin : m_thetaMax); double eta = eta_min + (eta_max-eta_min)*rnd.rndm(); double x1 = std::cos(phi); double x2 = std::sin(phi); -- GitLab