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

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -11 +11 @@
 #include "cosmocalc.h"
 #include "geneutils.h"
-#include "integfunc.h"
 #include "schechter.h"
 #include "planckspectra.h"
 
 /* --- Check Peterson at al. astro-ph/0606104 v1
-cmvdefsurv -z 0.0025 -x 1 -U 0.75,0.3,0.7,-1,1 -V 300 -O 400000,6000 -A 75 -M 6.156e9 -F 3 -2 1.5
+cmvdefsurv -z 0.0025 -x 1 -U 0.75,0.3,0.7,-1,1 -V 300 -O 400000,6000 -N 75 -M 6.156e9 -F 3 -2 1.5
  --- */
 
@@ -35 +34 @@
       <<" -L lobewidth : taille du lobe d\'observation (FWHM) en arcmin (def= 1\')"<<endl
       <<" -O surf,tobs : surface effective (m^2) et temps d\'observation (s)"<<endl
-      <<" -A Tsys : temperature du system (K)"<<endl
+      <<" -N Tsys : temperature du system (K)"<<endl
       <<" -S Tsynch,indnu : temperature (K) synch a 408 Mhz, index d\'evolution"<<endl
       <<"                   (indnu==0 no evolution with freq.)"<<endl
@@ -43 +42 @@
       <<" -U h100,om0,ol0,w0,or0,flat : cosmology"<<endl
       <<" -2 : two polarisations measured"<<endl
+      <<" -A <log10(S_agn)> : moyenne du flux AGN en Jy dans le pixel"<<endl
       <<" redshift : redshift moyen du survey"<<endl
       <<endl;
@@ -59 +59 @@
  double alpha = -1.37;
  cout<<"nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;
-
 
  // --- Arguments
@@ -76 +75 @@
  double vrot = 300.; // largeur en vitesse (km/s) pour elargissement doppler
  double facpolar = 0.5; // si on ne mesure les 2 polars -> 1.0
+ double lflux_agn = -3.;
 
  // --- Decodage arguments
  char c;
-  while((c = getopt(narg,arg,"hP2x:y:z:A:S:O:M:F:V:U:L:")) != -1) {
+  while((c = getopt(narg,arg,"hP2x:y:z:N:S:O:M:F:V:U:L:A:")) != -1) {
   switch (c) {
   case 'P' :
@@ -99 +99 @@
     sscanf(optarg,"%lf",&lobewidth);
     break;
-  case 'A' :
+  case 'N' :
     sscanf(optarg,"%lf",&Tsys);
     break;
@@ -119 +119 @@
   case '2' :
     facpolar = 1.0;
+    break;
+  case 'A' :
+    sscanf(optarg,"%lf",&lflux_agn);
     break;
   case 'h' :
@@ -281 +284 @@
  double lobecyl = sqrt(8.)*slobe; // diametre du lobe cylindrique equiv en arcmin
  double lobearea = M_PI*lobecyl*lobecyl/4.;  // en arcmin^2 (hypothese lobe gaussien)
+ double nlobes = rad2min(adtx)*rad2min(adty)/lobearea;
+ if(lobewidth<=0.) nlobes = 1.;
  cout<<"\nBeam FWHM = "<<lobewidth<<"\' -> sigma = "<<slobe<<"\' -> "
      <<" Dcyl = "<<lobecyl<<"\' -> area = "<<lobearea<<" arcmin^2"<<endl;
- double nlobes = rad2min(adtx)*rad2min(adty)/lobearea;
  cout<<"Number of beams in one transversal pixel = "<<nlobes<<endl;
 
@@ -290 +294 @@
  cout<<"flux factor = "<<hifactor<<" at redshift = "<<redshift<<endl;
 
- double fhi = hifactor*FluxHI(mhiref,dlum);
+ double fhi = hifactor*Msol2FluxHI(mhiref,dlum);
  cout<<"FluxHI("<<dlum<<" Mpc) all polar:"<<endl
      <<"  Flux= "<<fhi<<" W/m^2 = "<<fhi/Jansky2Watt_cst<<" Jy.Hz"<<endl
-     <<"      in ["<<hifactor*FluxHI(mhiref,dlumlim[0])
-     <<","<<hifactor*FluxHI(mhiref,dlumlim[1])<<"] W/m^2"<<endl;
+     <<"      in ["<<hifactor*Msol2FluxHI(mhiref,dlumlim[0])
+     <<","<<hifactor*Msol2FluxHI(mhiref,dlumlim[1])<<"] W/m^2"<<endl;
  double sfhi = fhi / (dnuhiz*1e9) / Jansky2Watt_cst;
  cout<<"If spread over "<<dnuhiz<<" GHz, flux density = "<<sfhi<<" Jy"<<endl;
@@ -342 +346 @@
  cout<<"                       flux density = "<<scmb<<" Jy"<<endl;
 
+ // AGN
+ double flux_agn = pow(10.,lflux_agn);
+ cout<<"AGN: log10(S_agn)="<<lflux_agn<<" -> S_agn="<<flux_agn<<" Jy"<<endl;
+ double flux_agn_pix = flux_agn*(dnuhiz*1e9);
+ double mass_agn_pix = FluxHI2Msol(flux_agn_pix*Jansky2Watt_cst,dlum);
+ double lmass_agn_pix = log10(mass_agn_pix);
+ cout<<"...pixel: f="<<flux_agn_pix<<" 10^-26 W/m^2"
+     <<" -> "<<mass_agn_pix<<" Msol -> log10 = "<<lmass_agn_pix<<endl;
+
  // --- Noise analysis
  cout<<"\n--- Noise analysis"<<endl;