Changeset 3365 in Sophya for trunk/Cosmo/SimLSS/geneutils.cc

Timestamp:

Oct 29, 2007, 7:21:41 PM (18 years ago)

Author:

cmv

Message:

1./ grosses modifs de structure pour cmvdefsurv.cc
2./ PoissRandLimit enleve de geneutils.{h,cc} et mis dans srandgen.{h,cc}

cmv 29/10/2007

File:

• trunk/Cosmo/SimLSS/geneutils.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/geneutils.cc

@@ -278 +278 @@
 //               dtheta in [0,Pi]
 // Return: l'angle solide en steradian
-// Approx pour theta petit: PI theta^2
+// Approx pour theta petit: PI * theta^2
 {
   return   2.*M_PI * (1.-cos(dtheta));
 }
 
-//-------------------------------------------------------------------
-unsigned long PoissRandLimit(double mu,double mumax)
-{
-  double pp,ppi,x;
-  unsigned long n;
-
-  if(mu>=mumax) {
-    pp = sqrt(mu);
-    while( (x=pp*NorRand()) < -mu );
-    return (unsigned long)(mu+x+0.5);
-  }
-
-  ppi = pp = exp(-mu);
-  x = drand01();
-  n = 0;
-  while (x > ppi) {
-    n++;
-    pp = mu*pp/(double)n;
-    ppi += pp;
-  }
-  return n;
+double FrAngSol(double angsol)
+// Retourne la demi-ouverture "dtheta" d'une calotte spherique d'angle solide "angsol"
+// Input: angle solide de la calotte spherique en steradians
+// Return: demi-ouverture de la calotte spherique en radians
+{
+  angsol = 1. - angsol/(2.*M_PI);
+  if(angsol<-1. || angsol>1.) return -1.;
+  return acos(angsol);
 }