Commit 0e5afe4b authored by Reza  ANSARI's avatar Reza ANSARI

Ajout commentaire avec valeur du coefficient de largeur pour beam gaussian...

Ajout commentaire avec valeur du coefficient de largeur pour beam gaussian determine par Jean-Eric en passant de SphericalBesselJ (bug) a BesselJ , Reza
parent 9d32df75
......@@ -141,14 +141,15 @@ public:
double Value(UnitVector const& uvo)
{
// circular beam response
double alp=acos(bdir_.Psc(uvo))*DoL_;
// remplacer acos(...) par sin(acos(...)) ds l'expression ci-dessus pour avoir
double alp=acos(bdir_.Psc(uvo))*DoL_;
// remplacer acos(...) par sin(acos(...)) ds l'expression ci-dessus pour avoir l'expression exacte qui est sin(theta)
if (fggauss_) {
// Reza / Sep 2019 : le facteur 2.61854 assure la meme largeur FWHM entre gaussian et Bessel
// 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...]
// Jean-Eric donne 0.28 pi^2 = 2.76 avec un fit de BesselJ
//---- 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)
}
......
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