Commit b1f4993a authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

correction du facteur determinant la largeur du beam gaussian,...

correction du facteur determinant la largeur du beam gaussian, exp(-2.61854*a^2) au lieu de exp(-2.1767*a^2) , Reza
parent 9065c96c
...@@ -141,10 +141,16 @@ public: ...@@ -141,10 +141,16 @@ public:
double Value(UnitVector const& uvo) double Value(UnitVector const& uvo)
{ {
// circular beam response // circular beam response
double alp=acos(bdir_.Psc(uvo))*DoL_; double alp=acos(bdir_.Psc(uvo))*DoL_;
// remplacer acos(...) par sin(acos(...)) ds l'expression ci-dessus pour avoir
if (fggauss_) { if (fggauss_) {
// JEC , wiki BAORadio 2.1767 = 0.22054*Pi^2 // Reza / Sep 2019 : le facteur 2.61854 assure la meme largeur FWHM entre gaussian et Bessel
return ( exp(-2.1767*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp) // On pose : a = D/lambda * gamma ( ou D/lambda * sin(gamma)
// Beam Bessel : beamJ(a) = ( 2 J1(pi a) / (pi a) )^2 Beam Gaussian : beamG(a) = Exp[ - K a^2 ]
// BeamJ(a=1.238098) = 0.5 et beamG(a=1.238098) = 0.5 pour K = 2.61854
// K=2.61854 au lieu de K=2.1767 (valeur de JEC , mais avec SphericalBesselJ[1...]
//---- return ( exp(-2.1767*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp)
return ( exp(-2.61854*alp*alp)*NormFac_ ); // Jiao : Check the factor 2 in the exp( ) -> exp(-2.*alp*alp)
} }
else { else {
if (alp<1.e-19) return NormFac_; if (alp<1.e-19) return NormFac_;
......
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