Changeset 3196 in Sophya for trunk/Cosmo/SimLSS/cmvobserv3d.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/cmvobserv3d.cc (modified) (7 diffs)

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -15 +15 @@
 #include "schechter.h"
 #include "geneutils.h"
-#include "integfunc.h"
 #include "genefluct3d.h"
 
@@ -32 +31 @@
      <<" -2 : compute 2D spectrum"<<endl
      <<" -M schmin,schmax,nsch : min,max mass and nb points for schechter HI"<<endl
+     <<" -A <log10(S_agn en Msol)>,sigma : distribution des AGN par pixel"<<endl
      <<" -W : write cube in FITS format"<<endl
      <<" -P : write cube in PPF format"<<endl
@@ -77 +77 @@
  double snoise= 0.;   // en equivalent MSol
 
+ // *** AGN
+ bool do_agn = false;
+ double lmsol_agn=-99., lsigma_agn=0.;   // en equivalent MSol
+
  // *** type de generation
  bool computefourier0=false;
@@ -91 +95 @@
 
  char c;
- while((c = getopt(narg,arg,"ha0PWV2Gx:y:z:s:Z:M:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWV2Gx:y:z:s:Z:M:A:")) != -1) {
   switch (c) {
   case 'a' :
@@ -122 +126 @@
   case 'M' :
     sscanf(optarg,"%lf,%lf,%d",&schmin,&schmax,&schnpt);
+    break;
+  case 'A' :
+    do_agn = true;
+    sscanf(optarg,"%lf,%lf",&lmsol_agn,&lsigma_agn);
     break;
   case 'V' :
@@ -153 +161 @@
      <<", schnpt="<<schnpt<<endl;
  cout<<"snoise="<<snoise<<" equivalent Msol"<<endl;
+ if(do_agn) cout<<"AGN: <log10(Msol)>="<<lmsol_agn<<" , sigma="<<lsigma_agn<<endl;
 
  //-----------------------------------------------------------------
@@ -476 +485 @@
  }
 
+ if(do_agn) {
+   cout<<"\n--- Add AGN: <mass>="<<lmsol_agn<<" , sigma="<<lsigma_agn<<endl;
+   fluct3d.AddAGN(lmsol_agn,lsigma_agn);
+   nm = fluct3d.MeanSigma2(rm,rs2);
+   cout<<"HI mass with AGN: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+       <<rs2<<" -> "<<sqrt(rs2)<<endl;
+   PrtTim(">>>> End Add AGN");
+ }
+
  if(snoise>0.) {
    cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<endl;