Changeset 3790 in Sophya for trunk/Cosmo/SimLSS/cmvtuniv.cc

Timestamp:

Jun 28, 2010, 2:12:24 PM (15 years ago)

Author:

cmv

Message:

etude des cubes generes par le GSM, cmv 28/06/2010

File:

• trunk/Cosmo/SimLSS/cmvtuniv.cc (modified) (3 diffs)

trunk/Cosmo/SimLSS/cmvtuniv.cc

@@ -10 +10 @@
 #include "ntuple.h"
 
+#include "constcosmo.h"
 #include "cosmocalc.h"
 #include "geneutils.h"
@@ -124 +125 @@
 {
  if(dz>z2-z1) dz = z2-z1;
+ double fz1 = 1000.*Fr_HyperFin_Par/(1.+z1);
+ double fz1pdz = 1000.*Fr_HyperFin_Par/(1.+z1+dz);
+ double fz2 = 1000.*Fr_HyperFin_Par/(1.+z2);
+ double fz2pdz = 1000.*Fr_HyperFin_Par/(1.+z2+dz);
  cout<<"\nTSTPRINT(): z1="<<z1<<" z2="<<z2<<" dz="<<dz<<endl;
- for(double z=z1;z<z2+dz/2.;z+=dz) {
-   cout<<"-- z = "<<z<<endl;
+ cout<<"  "<<fz1<<" -> "<<fz2<<" MHz,    df: "<<fz1-fz1pdz<<" , "<<fz2-fz2pdz<<" MHz"<<endl;
+ for(double z=z1;z<z2+dz/2.;z+=dz) {
+   cout<<"\n-- z = "<<z<<"      ("<<1000.*Fr_HyperFin_Par/(1.+z)<<" MHz)"<<endl;
    univ.Print(z);
    cout<<"Volume comoving in [z,z+"<<dz<<"] for 4Pi sr: "
@@ -138 +144 @@
    cout<<"1d -> "<<dang*a<<" Mpc = "<<dang*a*(1.+z)<<" Mpc com"<<endl;
    double dloscom = univ.Dhubble() / univ.E(z);
-   cout<<"dz=1 -> "<<dloscom<<" Mpc com"<<endl;
+   cout<<"dLOScom/dz -> "<<dloscom<<" Mpc com"<<endl;
  }