Skip to content
Snippets Groups Projects
Commit 61461a88 authored by Wouter Deconinck's avatar Wouter Deconinck Committed by Andre Sailer
Browse files

fix: ddsim --gun.distribution eta should set eta_min, eta_max correctly

parent a9b6a1bf
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
DDG4/src/Geant4IsotropeGenerator.cpp
+ 2
− 2
View file @ 61461a88
...@@ -76,8 +76,8 @@ void Geant4IsotropeGenerator::getParticleDirectionEta(int, ROOT::Math::XYZVector ...@@ -76,8 +76,8 @@ void Geant4IsotropeGenerator::getParticleDirectionEta(int, ROOT::Math::XYZVector
// See https://en.wikipedia.org/wiki/Pseudorapidity // See https://en.wikipedia.org/wiki/Pseudorapidity
const double dmin = std::numeric_limits<double>::epsilon(); const double dmin = std::numeric_limits<double>::epsilon();
double phi = m_phiMin+(m_phiMax-m_phiMin)*rnd.rndm(); 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_thetaMin>dmin ? m_thetaMin : dmin);
double eta_max = Distribution::eta(m_thetaMax>(M_PI-dmin) ? m_thetaMax : M_PI-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 eta = eta_min + (eta_max-eta_min)*rnd.rndm();
double x1 = std::cos(phi); double x1 = std::cos(phi);
double x2 = std::sin(phi); double x2 = std::sin(phi);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment