Commit cfd53652 authored by frey_m's avatar frey_m

modified: src/Distribution/ClosedOrbitFinder.h

Matthias: Just divide by "acon" to get \beta\gamma units.
parent 59002f18
......@@ -402,10 +402,9 @@ inline typename ClosedOrbitFinder<Value_type, Size_type, Stepper>::container_typ
*
* [a] = m --> [p] = m
*
* The momentum in \beta * \gamma is obtained by dividing by (\beta * \gamma * a)
* The momentum in \beta * \gamma is obtained by dividing by "a"
*/
value_type beta = std::sqrt(1.0 - 1.0 / ( gamma_m * gamma_m ) );
value_type factor = 1.0 / ( beta * gamma_m * acon_m(wo_m) );
value_type factor = 1.0 / acon_m(wo_m);
std::for_each(pr.begin(), pr.end(), [factor](value_type p) { return p * factor; });
return pr;
......
Markdown is supported
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