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

Timestamp:

Oct 23, 2007, 4:59:17 PM (18 years ago)

Author:

cmv

Message:

add possibility to just multiply dRho.Rho , cmv 23/10/2007

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -28 +28 @@
      <<"      (default: no, spectrum Pk(z=z_median) for all cube)"<<endl
      <<" -F : filter spectrum by pixel shape (0=no 1=yes(default)"<<endl
-     <<" -U : do not poisson fluctuate with Ngal, convert directly to HI mass"<<endl
+     <<" -U typ : do not poisson fluctuate with Ngal, convert directly to HI mass"<<endl
+     <<"          (typ<0 just multiply dRho/Rho by mass_by_pixel)"<<endl
      <<" -x nx,dx : size along x axis (npix,Mpc)"<<endl
      <<" -y ny,dy : size along y axis (npix,Mpc)"<<endl
@@ -100 +101 @@
  string schmassdistfile = "";
  bool no_poisson = false;
+ int no_poisson_type = 1;
 
  double scalemass = -1.;
@@ -133 +135 @@
 
  char c;
- while((c = getopt(narg,arg,"ha0PWSV2UG:F:x:y:z:s:Z:M:A:T:N:Q:R:8:")) != -1) {
+ while((c = getopt(narg,arg,"ha0PWSV2U:G:F:x:y:z:s:Z:M:A:T:N:Q:R:8:")) != -1) {
   int nth = 0;
   switch (c) {
@@ -147 +149 @@
   case 'U' :
     no_poisson = true;
+    sscanf(optarg,"%d",&no_poisson_type);
     break;
   case 'F' :
@@ -230 +233 @@
      <<"), schmax="<<schmax<<" ("<<lschmax<<") Msol"
      <<", schnpt="<<schnpt<<endl;
- if(no_poisson) cout<<"No poisson fluctuation, direct conversion to HI mass"<<endl;
+ if(no_poisson) cout<<"No poisson fluctuation, direct conversion to HI mass, typ="
+                    <<no_poisson_type<<endl;
  cout<<"snoise="<<snoise<<" equivalent Msol, evolution="<<noise_evol<<endl;
  cout<<"scalemass="<<scalemass<<endl;
@@ -291 +295 @@
  cout<<"Spectrum normalisation = "<<pkz.GetScale()<<endl;
 
+ {
+ Histo hpkz0(lkmin,lkmax,npt); hpkz0.ReCenterBin();
+ FuncToHisto(pkz,hpkz0,true);
+ tagobs = "hpkz0"; posobs.PutObject(hpkz0,tagobs);
+ }
+
  pkz.SetZ(zref);
 
+ {
  Histo hpkz(lkmin,lkmax,npt); hpkz.ReCenterBin();
  FuncToHisto(pkz,hpkz,true);
- {
  tagobs = "hpkz"; posobs.PutObject(hpkz,tagobs);
  }
@@ -542 +552 @@
    cout<<"\n--- Converting !!!DIRECTLY!!! mass into HI mass: mass per pixel ="
        <<mass_by_pixel<<endl;
-   rm = fluct3d.TurnFluct2MeanNumber(mass_by_mpc3); // ici on doit donner Msol/Mpc^3
+   if(no_poisson_type<0) {
+     cout<<"!!!WARNING!!! dRho/Rho is just multiplied by "<<mass_by_pixel<<endl;
+     fluct3d.ScaleOffset(mass_by_pixel);
+   } else {
+     rm = fluct3d.TurnFluct2MeanNumber(mass_by_mpc3); // ici on doit donner Msol/Mpc^3
+   }
  } else {
    cout<<"\n--- Converting mass into galaxy number: gal per pixel ="