Changeset 3271 in Sophya for trunk/Cosmo/SimLSS/cmvdefsurv.cc

Timestamp:

Jun 20, 2007, 6:50:48 PM (18 years ago)

Author:

cmv

Message:

intro evolution bruit avec redshift cmv 20/06/2007

File:

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

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -185 +185 @@
  if(redlim[0]<=0.)
    {cout<<"Lower redshift limit "<<redlim[0]<<" should be >0"<<endl; return -3;}
- double dtrlim[2] = {univ.Dtrcom(redlim[0]),univ.Dtrcom(redlim[1])};
- double loslim[2] = {univ.Dloscom(redlim[0]), univ.Dloscom(redlim[1])};
+ double dtrlim[2]  = {univ.Dtrcom(redlim[0]),univ.Dtrcom(redlim[1])};
+ double loslim[2]  = {univ.Dloscom(redlim[0]), univ.Dloscom(redlim[1])};
  double dlumlim[2] = {univ.Dlum(redlim[0]),univ.Dlum(redlim[1])};
 
@@ -192 +192 @@
  cout<<"---- Line of Sight: Redshift = "<<redshift<<endl
      <<"dred = "<<dred<<"  redlarg = "<<redlarg<<endl
-     <<"  dz = "<<dz<<" Mpc  redlarg = "<<tzlarg<<" Mpc, nz = "<<nz<<" pix"<<endl;
+     <<"  dz = "<<dz<<" Mpc  redlarg = "<<tzlarg<<" Mpc com, nz = "<<nz<<" pix"<<endl;
  cout<<"---- Transverse X:"<<endl
      <<"adtx = "<<rad2min(adtx)<<"\',  atxlarg = "<<rad2deg(atxlarg)<<" d"<<endl
-     <<"  dx = "<<dx<<" Mpc,  txlarg = "<<txlarg<<" Mpc, nx = "<<nx<<" pix"<<endl;
+     <<"  dx = "<<dx<<" Mpc,  txlarg = "<<txlarg<<" Mpc com, nx = "<<nx<<" pix"<<endl;
  cout<<"---- Transverse Y:"<<endl
      <<"adty = "<<rad2min(adty)<<"\',  atylarg = "<<rad2deg(atylarg)<<" d"<<endl
-     <<"  dy = "<<dy<<" Mpc,  tylarg = "<<tylarg<<" Mpc, ny = "<<ny<<" pix"<<endl;
+     <<"  dy = "<<dy<<" Mpc,  tylarg = "<<tylarg<<" Mpc com, ny = "<<ny<<" pix"<<endl;
  cout<<"---- Npix total = "<<Npix<<" -> "<<Npix*sizeof(double)/1.e6<<" Mo"<<endl;
 
@@ -275 +275 @@
  // dnu = NuHi/(1.+z0-dz/2) - NuHi/(1.+z0+dz/2)
  //     = NuHi*dz/(1.+z0)^2 * 1/[1-(dz/(2*(1+z0)))^2]
+ //    ~= NuHi*dz/(1.+z0)^2
  double dnuhiz = Fr_HyperFin_Par *dred/(unplusz*unplusz)
                / (1.- pow(dred/.2/unplusz,2.));