Commit f9936cf1 authored by perdereau's avatar perdereau
Browse files
parents 6e5960df 0e5afe4b
...@@ -141,14 +141,15 @@ public: ...@@ -141,14 +141,15 @@ 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 // remplacer acos(...) par sin(acos(...)) ds l'expression ci-dessus pour avoir l'expression exacte qui est sin(theta)
if (fggauss_) { if (fggauss_) {
// Reza / Sep 2019 : le facteur 2.61854 assure la meme largeur FWHM entre gaussian et Bessel // 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) // 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 ] // 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 // 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...] // 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.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) 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