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

Timestamp:

Jan 19, 2007, 5:19:49 PM (19 years ago)

Author:

cmv

Message:

modifs d'options de cmvobserv3d.cc
creation de init_fftw() et deplacement a l'initialisation
creation de SetObservator pour positionner l'observateur + ecriture FITS et PPF

(mais le codage des distance/redshift reste a coder)

cmv 19/01/2007

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -28 +28 @@
      <<" -Z zref : redshift median"<<endl
      <<" -s snoise : sigma du bruit en Msol"<<endl
+     <<" -2 : compute 2D spectrum"<<endl
+     <<" -F : write cube in FITS format"<<endl
+     <<" -P : write cube in PPF format"<<endl
+     <<" -V : compute variance from real space"<<endl
      <<endl;
 }
@@ -73 +77 @@
  unsigned short nthread=4;
 
+ // *** What to do
+ bool comp2dspec = false;
+ bool wfits = false;
+ bool wppf = false;
+ bool compvarreal = false;
+
  // --- Decodage arguments
 
  char c;
- while((c = getopt(narg,arg,"ha0x:y:z:s:Z:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWV2x:y:z:s:Z:")) != -1) {
   switch (c) {
   case 'a' :
@@ -98 +108 @@
   case 'Z' :
     sscanf(optarg,"%lf",&zref);
+    break;
+  case '2' :
+    comp2dspec = true;
+    break;
+  case 'V' :
+    compvarreal = true;
+    break;
+  case 'W' :
+    wfits = true;
+    break;
+  case 'P' :
+    wppf = true;
     break;
   case 'h' :
@@ -185 +207 @@
 
  //-----------------------------------------------------------------
- cout<<endl<<"\n--- Compute galaxies density number"<<endl;
+ cout<<endl<<"\n--- Compute galaxy density number"<<endl;
 
  sch.SetOutValue(0);
@@ -249 +271 @@
  }
 
- if(1) {
+ if(comp2dspec) {
    cout<<"\n--- Checking realization 2D spectra"<<endl;
    Histo2DErr hpkgen2(0.,ktnyqmax,nherrt,0.,kznyqmax,nherrz);
@@ -265 +287 @@
  }
 
- if(1) fluct3d.WriteFits("!cmvobserv3d_k.fits");
- if(1) fluct3d.WritePPF("cmvobserv3d_k.ppf",false);
+ if(wfits) fluct3d.WriteFits("!cmvobserv3d_k0.fits");
+ if(wppf) fluct3d.WritePPF("cmvobserv3d_k0.ppf",false);
 
  if(1) {
@@ -279 +301 @@
  }
 
- if(1) {
+ if(comp2dspec) {
    cout<<"\n--- Checking realization 2D spectra after pixel shape convol."<<endl;
    Histo2DErr hpkgenf2(0.,ktnyqmax,nherrt,0.,kznyqmax,nherrz);
@@ -295 +317 @@
  cout<<"rgen.Min = "<<rmin<<" , Max="<<rmax<<endl;
 
- if(1) fluct3d.WriteFits("!cmvobserv3d_r0.fits");
- if(1) fluct3d.WritePPF("cmvobserv3d_r0.ppf",true);
+ if(wfits) fluct3d.WriteFits("!cmvobserv3d_r0.fits");
+ if(wppf) fluct3d.WritePPF("cmvobserv3d_r0.ppf",true);
 
  int_8 nm;
@@ -308 +330 @@
  }
 
- if(1) {
+ if(compvarreal) {
    cout<<"\n--- Check variance sigmaR in real space"<<endl;
    double varr;
@@ -336 +358 @@
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
 
- cout<<"\n--- Apply poisson on galaxie number"<<endl;
+ cout<<"\n--- Apply poisson on galaxy number"<<endl;
  nm = fluct3d.ApplyPoisson();
  cout<<nm<<" galaxies into survey after poisson"<<endl;
@@ -343 +365 @@
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
 
- cout<<"\n--- Convert Galaxies number to HI mass"<<endl;
+ cout<<"\n--- Convert Galaxy number to HI mass"<<endl;
  long nhmdndm = long( (lschmax-lschmin+1.)*100. + 0.5);
  Histo hmdndm(lschmin,lschmax,nhmdndm);
@@ -368 +390 @@
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
 
- if(1) fluct3d.WriteFits("!cmvobserv3d_r.fits");
- if(1) fluct3d.WritePPF("cmvobserv3d_r.ppf",true);
+ if(wfits) fluct3d.WriteFits("!cmvobserv3d_r.fits");
+ if(wppf) fluct3d.WritePPF("cmvobserv3d_r.ppf",true);
 
  cout<<"\n--- Add noise to HI Flux snoise="<<snoise<<endl;
@@ -384 +406 @@
    fluct3d.ReComputeFourier();
  }
+
+ if(wfits) fluct3d.WriteFits("!cmvobserv3d_k.fits");
+ if(wppf) fluct3d.WritePPF("cmvobserv3d_k.ppf",false);
 
  if(1) {
@@ -394 +419 @@
  }
 
- if(1) {
+ if(comp2dspec) {
    cout<<"\n--- Computing final 2D spectrum"<<endl;
    Histo2DErr hpkrec2(0.,ktnyqmax,nherrt,0.,kznyqmax,nherrz);
@@ -410 +435 @@
 /*
 ######################################################
+readfits cmvobserv3d_k0.fits
 readfits cmvobserv3d_k.fits
 readfits cmvobserv3d_r0.fits
 readfits cmvobserv3d_r.fits
 
+openppf cmvobserv3d_k0.ppf
 openppf cmvobserv3d_k.ppf
 openppf cmvobserv3d_r0.ppf