Changeset 3352 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

Oct 16, 2007, 7:25:51 PM (18 years ago)

Author:

cmv

Message:

ajour calcul OmegaHI equiv, cmv 16/10/2007

File:

• trunk/Cosmo/SimLSS/cmvdefsurv.cc (modified) (13 diffs)

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -43 +43 @@
       <<" -K k,dk,pk : k(Mpc^-1) dk(Mpc^-1) pk(a z=0 en Mpc^-3) pour estimer la variance cosmique"<<endl
       <<"----------------"<<endl
-      <<" -O surf,tobs : surface effective (m^2) et temps d\'observation (s)"<<endl
+      <<" -O surf,tobs,eta : surface effective (m^2), temps d\'observation (s), efficacite d\'ouverture"<<endl
       <<" -N Tsys : temperature du system (K)"<<endl
       <<" -L lobewidth,freqlob : taille du lobe d\'observation (FWHM) en arcmin (def= 1\')"<<endl
@@ -86 +86 @@
    int nz=0; char unit_z = 'Z';
  double redshift = 0.;
- double tobs = 6000., surfeff = 400000.;
+ double tobs = 6000., surftot = 400000., eta_eff  = 1.;
  // variance cosmique (default = standard SDSSII)
  double kcosm = 0.05, dkcosm = -1., pkcosm = 40000.;
@@ -115 +115 @@
     break;
   case 'O' :
-    sscanf(optarg,"%lf,%lf",&surfeff,&tobs);
+    sscanf(optarg,"%lf,%lf,%lf",&surftot,&tobs,&eta_eff);
     break;
   case 'L' :
@@ -169 +169 @@
  univ.Print(redshift);
  GrowthFactor growth(om0,ol0);
- double grothfac = growth(redshift);
- cout<<"Facteur de croissance lineaire = "<<grothfac
+ double growthfac = growth(redshift);
+ cout<<"Facteur de croissance lineaire = "<<growthfac
      <<" ^2 , ( 1/(1+z)="<<1./(1.+redshift)<<" )"<<endl;
+ double rhoc0 = univ.Rhoc(0.)*GCm3toMsolMpc3_Cst;
+ double rhocz = univ.Rhoc(redshift)*GCm3toMsolMpc3_Cst;
+ cout<<"Rho_c a z=0: "<<rhoc0<<", a z="<<redshift<<": "<<rhocz<<" Msol/Mpc^3"<<endl;
 
  double dang    = univ.Dang(redshift);
@@ -313 +316 @@
  // nombre de mode = 1/2 * Vsurvey/(2Pi)^3 * 4Pi*k^2*dk
  if(kcosm>0. && pkcosm>0.) {
-   double pk = pkcosm*pow(grothfac,2.);
+   double pk = pkcosm*pow(growthfac,2.);
    cout<<"\n>>>>\n>>>> variance cosmique pour k="<<kcosm<<" Mpc^-1, pk(z=0)="
        <<pkcosm<<", pk(z="<<redshift<<")="<<pk<<"\n>>>>"<<endl;
@@ -346 +349 @@
  cout<<"Pixel mass = "<<masshipix<<" Msol"<<endl
      <<"Total mass in survey = "<<masshitot<<" Msol"<<endl;
+ cout<<"OmegaHI a z=0: "<<masshimpc3/rhoc0
+     <<", a z="<<redshift<<": "<<masshimpc3/rhocz<<endl;
  if(mhiref<=0.)  mhiref = masshipix;
 
@@ -361 +366 @@
  // --- Survey values
  cout<<"\n>>>>\n>>>> Observations\n>>>>"<<endl;
+ double surfeff = surftot * eta_eff;
  double unplusz = 1.+redshift;
  double nuhiz = Fr_HyperFin_Par / unplusz;  // GHz
@@ -368 +374 @@
  double dnuhiz = Fr_HyperFin_Par *dred/(unplusz*unplusz)
                / (1.- pow(dred/.2/unplusz,2.));
- cout<<"      surf_eff="<<surfeff<<" m^2, tobs="<<tobs<<" s"<<endl
+ cout<<"      surf="<<surftot<<" m^2, eta="<<eta_eff<<"  surf_eff="<<surfeff<<" m^2, tobs="<<tobs<<" s"<<endl
      <<"      nu="<<nuhiz<<" GHz, dnu="<<dnuhiz*1.e3<<" Mhz"<<endl;
  cout<<"dang lumi = "<<dlum<<" in ["<<dlumlim[0]<<","<<dlumlim[1]<<"] Mpc"<<endl;
@@ -399 +405 @@
      <<","<<hifactor*Msol2FluxHI(mhiref,dlumlim[1])<<"] W/m^2"<<endl;
  double sfhi = fhi / (dnuhiz*1e9) / Jansky2Watt_cst;
- cout<<"If spread over pixel depth ("<<dnuhiz<<" GHz), flux density = "<<sfhi<<" Jy"<<endl;
+ cout<<"If spread over pixel depth ("<<dnuhiz<<" GHz), flux density = "<<sfhi*1.e6<<" mu_Jy"<<endl;
 
  // --- Signal analysis
@@ -421 +427 @@
      <<"     P="<<psig<<" W"<<endl
      <<"     flux density = "<<ssig*1.e6<<" mu_Jy  (for Dnu="<<dnuhiz<<" GHz)"<<endl
+     <<"                   ("<<ssig_2polar<<" mu_Jy for 2 polars)"<<endl
      <<"     Antenna temperature: tsig="<<tsig<<" K"<<endl;
 
@@ -474 +481 @@
  //     Pb = k * Tsys * dNu  (W)
  // Pour une source (non-polarisee) de densite de flux (totale 2 polars)
- //     St (exprimee en Jy=W/m^2/Hz)
+ //     St (exprimee en Jy = 10^-26 W/m^2/Hz)
  //     Pt = St * Ae * dNu   (puissance totale emise en W pour 2 polars)
  //     P1 = 1/2 * St * Ae * dNu   (puissance emise en W pour une polar)
@@ -565 +572 @@
    smass = mhiref / ssig_2polar * slim;
    cout<<"for 1 lobe:"<<endl
-       <<"   Slim    = "<<slim<<" Jy"<<endl
+       <<"   Slim    = "<<slim*1.e6<<" mu_Jy"<<endl
        <<"   S/N     = "<<SsN<<endl
        <<"   Mass HI = "<<smass<<" Msol"<<endl;
@@ -573 +580 @@
    smass_nl = mhiref / ssig_2polar * slim_nl;
    cout<<"for "<<nlobes<<" lobes:"<<endl
-       <<"   Flux    = "<<slim_nl<<" Jy"<<endl
+       <<"   Flux    = "<<slim_nl*1.e6<<" mu_Jy"<<endl
        <<"   S/N     = "<<SsN_nl<<endl
        <<"   Mass HI = "<<smass_nl<<" Msol"<<endl;