Changeset 3283 in Sophya for trunk/Cosmo/SimLSS/cmvobserv3d.cc

Timestamp:

Jul 25, 2007, 2:41:39 PM (18 years ago)

Author:

cmv

Message:

scale du spectre final cmv 24/07/2007

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -208 +208 @@
  GrowthFactor growth(om0,ol0);
  // GrowthFactor growth(1.,0.); // D(z) = 1/(1+z)
+ double growth_at_z = growth(zref);
+ cout<<"...Growth factor at z="<<zref<<" = "<<growth_at_z<<endl;
 
  PkSpectrum0 pk0(pkini,tf);
@@ -439 +441 @@
    double varr;
    fluct3d.VarianceFrReal(R,varr);
+   cout<<"...Computed variance = "<<varr
+       <<" , Theorical variance at (z=0) = "<<pow(sigmaR,1.)
+       <<" , at (z="<<zref<<") = "<<pow(sigmaR*growth_at_z,1.)<<endl;
    PrtTim(">>>> End Check variance sigmaR in real space");
  }
@@ -511 +516 @@
  }
 
+ double scalecube = -1., offsetcube=0.;
  if(snoise>0.) {
    cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<", evolution="<<isnoise_evol<<endl;
    fluct3d.AddNoise2Real(snoise,(isnoise_evol>0? true:false));
      nm = fluct3d.MeanSigma2(rm,rs2);
+     if(rs2>0. && sr2int>0.) scalecube = sqrt(sr2int)/sqrt(rs2);
+     offsetcube = -rm;
    PrtTim(">>>> End Add noise");
  }
 
+ if(scalecube>0.) {
+   cout<<"\n--- Scale and offset scale="<<scalecube
+       <<" off="<<offsetcube<<endl;
+   fluct3d.ScaleOffset(scalecube,offsetcube);
+ }
+ 
  if(wfits) {
    fluct3d.WriteFits("!cmvobserv3d_rf.fits");