Changeset 3196 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.cc

Timestamp:

Apr 3, 2007, 4:06:50 PM (18 years ago)

Author:

cmv

Message:

les AGN selon C.Jackson, une premiere approche simplifiee, recodage from Jim Rich. cmv 03/04/2007

File:

• trunk/Cosmo/SimLSS/genefluct3d.cc (modified) (2 diffs)

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -20 +20 @@
 
 #include "constcosmo.h"
-#include "integfunc.h"
 #include "geneutils.h"
 
@@ -970 +969 @@
 }
 
+void GeneFluct3D::AddAGN(double lmsol,double lsigma)
+// Add AGN flux into simulation:
+// --- Procedure:
+// 1. lancer "cmvdefsurv" avec les parametres du survey
+//    et recuperer l'angle solide "angsol sr" du pixel elementaire
+//    au centre du cube.
+// 2. lancer "cmvtstagn" pour cet angle solide -> cmvtstagn.ppf
+// 3. regarder l'histo "hlfang" et en deduire un equivalent gaussienne
+//    cad une moyenne <log10(S)> et un sigma "sig"
+//    Attention: la distribution n'etant pas gaussienne
+// 4. re-lancer "cmvdefsurv" en ajoutant l'info sur les AGN
+//    "cmvdefsurv ... -G <log10(S)> ..."
+//    et recuperer le log10(masse solaire equivalente)
+// --- Limitations actuelle du code:
+// a. l'angle solide du pixel est pris au centre du cube
+//    et ne varie pas avec la distance a l'interieur du cube
+// b. la taille en dNu des pixels (selon z) est supposee constante
+//    et egale a celle calculee pour le centre du cube
+// c. les AGN sont supposes plats en flux
+// d. le flux des AGN est mis dans une colonne Oz (indice k)
+//    et pas sur la ligne de visee
+// e. la distribution est approximee a une gaussienne
+// .. C'est une approximation pour un observateur loin du centre du cube
+//    et pour un cube peu epais / distance observateur 
+// --- Parametres de la routine:
+// lmsol : c'est le <log10(Msol)> qui est la conversion en masse solaire
+//         du flux depose dans un pixel par les AGN
+// lsigma : c'est le sigma de la distribution
+// - Comme on est en echelle log10():
+//   on tire log10(Msol) + X
+//   ou X est une realisation sur une gaussienne de variance "sig^2"
+//   La masse realisee est donc: Msol*10^X
+// - Pas de probleme de pixel negatif car on a une multiplication!
+{
+  if(lp_>0) cout<<"--- AddAGN: lmsol = "<<lmsol<<" , sigma = "<<lsigma<<endl;
+  check_array_alloc();
+
+  double m = pow(10.,lmsol);
+
+  double sum=0., sum2=0., nsum=0.;
+  for(long l=0;l<Nz_;l++) {
+    double a = lsigma*NorRand();
+    a = m*pow(10.,a);
+    // On met le meme tirage le long de Oz (indice k)
+    for(long i=0;i<Nx_;i++) for(long j=0;j<Ny_;j++) {
+      int_8 ip = IndexR(i,j,l);
+      data_[ip] += a;
+    }
+    sum += a; sum2 += a*a; nsum += 1.;
+  }
+
+ if(nsum>1.) {
+   sum /= nsum;
+   sum2 = sum2/nsum - sum*sum;
+   cout<<"...Mean mass="<<sum<<" Msol , s^2="<<sum2<<" s="<<sqrt(fabs(sum2))<<endl;
+ }
+ 
+}
+
 void GeneFluct3D::AddNoise2Real(double snoise)
 // add noise to every pixels (meme les <=0 !)