Implementing second PDA parametrization for calculation of D-terms

235a22d4 · Jose Manuel Morgado Chávez · d2146832 · 235a22d4
Commit 235a22d4 authored Jun 25, 2021 by Jose Manuel Morgado Chávez
--- a/src/NumA/RadonTransform/RadonTransform.cpp
+++ b/src/NumA/RadonTransform/RadonTransform.cpp
@@ -562,6 +562,7 @@ double RadonTransform::computeDterm( const Eigen::VectorXd &DD, const double x,

    // PDA \phi(x)=6x(1-x)
    double PDA = 0.75*(1-pow(alpha,2))/xi;
+    // double PDA = 1.445*(1-pow(alpha,2))*sqrt(2.74589-pow(alpha,2)-2.56139*sqrt(1-pow(alpha,2)))/xi;

    // Contribution
    Dterm += PDA;