Changeset 3799 in Sophya for trunk/Cosmo/SimLSS/genefluct3d.cc

Timestamp:

Jun 30, 2010, 6:23:54 PM (15 years ago)

Author:

cmv

Message:

Vlos * by (1+z), Dlos(com) do not * by (1+z), cmv 30/06/2010

File:

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

trunk/Cosmo/SimLSS/genefluct3d.cc

@@ -421 +421 @@
  for(int i=0;i<nptd;i++) {
    double z = zredmin_dpsd_ + i*dz;
-   double v = -cosmo_->H(z) * growth_->DsDz(z,good_dzinc_) / (*growth_)(z);
+   double om = cosmo_->Omatter(z);
+   double h = cosmo_->H(z);
+   double dsdz = growth_->DsDz(z,good_dzinc_);
+   double d = (*growth_)(z);
+   double v = -h * (1.+z) * dsdz / d;
    dpsdfrzred_.push_back(v);
-   if(lp_ && (i%nmod==0 || i==nptd-1)) cout<<"    z="<<z<<"  D'/D="<<v<<" km/s/Mpc"<<endl;
+   if(lp_ && (i%nmod==0 || i==nptd-1))
+     cout<<"    z="<<z<<" H="<<h<<" Beta="<<-(1.+z)*dsdz/d
+         <<" (Om^0.6="<<pow(om,0.6)
+         <<")  -H*(1+z)*D'/D="<<v<<" km/s/Mpc"<<endl;
  }
 
@@ -901 +908 @@
 {
  double zpk = compute_pk_redsh_ref_;
- double dpsd = -cosmo_->H(zpk) * growth_->DsDz(zpk,good_dzinc_) / (*growth_)(zpk);
- if(lp_>0) cout<<"--- ToVelLoS --- at z="<<zpk<<", D'/D="<<dpsd
+ double dpsd = -cosmo_->H(zpk) * (1.+zpk) * growth_->DsDz(zpk,good_dzinc_) / (*growth_)(zpk);
+ if(lp_>0) cout<<"--- ToVelLoS --- at z="<<zpk<<", -H*(1+z)*D'/D="<<dpsd
                <<" (km/s)/Mpc (comp. for dz="<<good_dzinc_<<")"<<endl;
  check_array_alloc();
@@ -980 +987 @@
 
  double zpk = compute_pk_redsh_ref_;
- double dpsd_orig = - cosmo_->H(zpk) * growth_->DsDz(zpk,good_dzinc_) / (*growth_)(zpk);
+ double dpsd_orig = - cosmo_->H(zpk) * (1.+zpk) * growth_->DsDz(zpk,good_dzinc_) / (*growth_)(zpk);
 
  InterpFunc interpinv(loscom2zred_min_,loscom2zred_max_,loscom2zred_);