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

Timestamp:

Dec 21, 2006, 4:45:09 PM (19 years ago)

Author:

cmv

Message:

Suite de la mise dans la base cvs cmv 21/12/2006

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -12 +12 @@
 #include "srandgen.h"
 #include "perandom.h"
+#include "fabtwriter.h"
 
 #include "arrctcast.h"
@@ -31 +32 @@
      <<" -z nz,dz : taille en z (npix,Mpc)"<<endl
      <<" -Z zref : redshift median"<<endl
-     <<" -s scale : sigma du bruit en Msol"<<endl
+     <<" -s snoise : sigma du bruit en Msol"<<endl
      <<endl;
 }
@@ -163 +164 @@
  pkz.SetZ(zref);
 
- Histo hpkz (lkmin,lkmax,npt); hpkz.ReCenterBin();
+ Histo hpkz(lkmin,lkmax,npt); hpkz.ReCenterBin();
  FuncToHisto(pkz,hpkz,true);
  {
@@ -330 +331 @@
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
 
+ if(1) {
+   cout<<"\n---Writing Cube to Fits file"<<endl;
+   FitsImg3DWriter fwrt("!cmvobserv3dr.fits",FLOAT_IMG,5);
+   fwrt.WriteKey("ZREF",zref," redshift");
+   vector<long> n = fluct3d.GetNpix();
+   fwrt.WriteKey("NX",n[0]," axe transverse 1");
+   fwrt.WriteKey("NY",n[1]," axe transverse 2");
+   fwrt.WriteKey("NZ",n[2]," axe longitudinal (redshift)");
+   vector<r_8> d;
+   d = fluct3d.GetDinc();
+   fwrt.WriteKey("DX",d[0]," Mpc");
+   fwrt.WriteKey("DY",d[1]," Mpc");
+   fwrt.WriteKey("DZ",d[2]," Mpc");
+   d = fluct3d.GetKinc();
+   fwrt.WriteKey("DKX",d[0]," Mpc^-1");
+   fwrt.WriteKey("DKY",d[1]," Mpc^-1");
+   fwrt.WriteKey("DKZ",d[2]," Mpc^-1");
+   fwrt.WriteKey("SNOISE",snoise," Msol");
+   fwrt.Write(rgen);
+ }
+
  cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<endl;
  fluct3d.AddNoise2Real(snoise);
@@ -364 +386 @@
 objaoper pkgen sliceyz 0
 
-n/plot hpkz.val%x x>0 ! "connectpoints"
-n/plot hpkgen.val%log10(x) x>0 ! "same red connectpoints"
-n/plot hpkgenf.val%log10(x) x>0 ! "same pink connectpoints"
-n/plot hpkrec.val%log10(x) x>0 ! "same blue connectpoints"
+n/plot hpkz.val%x ! ! "nsta connectpoints"
+n/plot hpkgen.val%log10(x) x>0 ! "nsta same red connectpoints"
+n/plot hpkgenf.val%log10(x) x>0 ! "nsta same orange connectpoints"
+n/plot hpkrec.val%log10(x) x>0 ! "nsta same blue connectpoints"
 
 set k pow(10.,x)
@@ -379 +401 @@
 rgensl 0
 
+disp hmdndm
+disp tirhmdndm
+addline 0 1 20 1 "red"
+
  */