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

Timestamp:

Sep 22, 2008, 4:07:02 PM (17 years ago)

Author:

cmv

Message:

intro du remplissage du NTuple de debug cmv 22/09/2008

File:

• trunk/Cosmo/SimLSS/cmvobserv3d.cc (modified) (18 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 typ : do not poisson fluctuate with Ngal, convert directly to HI mass"<<endl
-     <<"          (typ<0 just multiply dRho/Rho by mass_by_pixel)"<<endl
+     <<" -U typ =0 compute <NGal>, poisson fluctuate then convert to HI mass (def)"<<endl
+     <<"        >0 compute directly <HI mass>, do NOT poisson fluctuate with <Ngal>"<<endl
+     <<"        <0 just multiply dRho/Rho by mass_by_pixel (possible negative pixel values)"<<endl
      <<" -x nx,dx : size along x axis (npix,Mpc)"<<endl
      <<" -y ny,dy : size along y axis (npix,Mpc)"<<endl
@@ -100 +101 @@
  bool recompute_schmassdist = true;
  string schmassdistfile = "";
- bool no_poisson = false;
- int no_poisson_type = 1;
+ int no_poisson_type = 0;
 
  double scalemass = -1.;
@@ -126 +126 @@
  bool wslice = false;
  bool compvarreal = false;
+ unsigned long ntnent = 10000;  // 0 = do not fill NTuple
 
  // --- Decodage arguments
@@ -148 +149 @@
     break;
   case 'U' :
-    no_poisson = true;
     sscanf(optarg,"%d",&no_poisson_type);
     break;
@@ -233 +233 @@
      <<"), schmax="<<schmax<<" ("<<lschmax<<") Msol"
      <<", schnpt="<<schnpt<<endl;
- if(no_poisson) cout<<"No poisson fluctuation, direct conversion to HI mass, typ="
-                    <<no_poisson_type<<endl;
+ if(no_poisson_type!=0) 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;
@@ -419 +419 @@
  if(computefourier0) fluct3d.ComputeFourier0(pkz);
    else              fluct3d.ComputeFourier(pkz);
+ fluct3d.NTupleCheck(posobs,string("ntpkgen"),ntnent);
  PrtTim(">>>> End Computing a realization in Fourier space");
 
@@ -446 +447 @@
    cout<<"\n--- Computing convolution by pixel shape"<<endl;
    fluct3d.FilterByPixel();
+   fluct3d.NTupleCheck(posobs,string("ntpkgenf"),ntnent);
    PrtTim(">>>> End Computing convolution by pixel shape");
 
@@ -501 +503 @@
  double rmin,rmax; fluct3d.MinMax(rmin,rmax);
  cout<<"rgen.Min = "<<rmin<<" , Max="<<rmax<<endl;
+ fluct3d.NTupleCheck(posobs,string("ntreal"),ntnent);
  PrtTim(">>>> End Computing a realization in real space");
 
@@ -509 +512 @@
    fluct3d.MinMax(rmin,rmax);
    cout<<"rgen.Min = "<<rmin<<" , Max="<<rmax<<endl;
+   fluct3d.NTupleCheck(posobs,string("ntgrow"),ntnent);
    PrtTim(">>>> End Applying growth factor");
  }
@@ -549 +553 @@
 
  //-----------------------------------------------------------------
- if(no_poisson) {
+ if(no_poisson_type!=0) {
    cout<<"\n--- Converting !!!DIRECTLY!!! mass into HI mass: mass per pixel ="
        <<mass_by_pixel<<endl;
@@ -558 +562 @@
      rm = fluct3d.TurnFluct2MeanNumber(mass_by_mpc3); // ici on doit donner Msol/Mpc^3
    }
+   fluct3d.NTupleCheck(posobs,string("ntmhi"),ntnent);
  } else {
    cout<<"\n--- Converting mass into galaxy number: gal per pixel ="
        <<ngal_by_mpc3*fluct3d.GetDVol()<<endl;
    rm = fluct3d.TurnFluct2MeanNumber(ngal_by_mpc3); // ici on doit donner Ngal/Mpc^3
+   fluct3d.NTupleCheck(posobs,string("ntmeang"),ntnent);
  }
  nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200);
@@ -568 +574 @@
  PrtTim(">>>> End Converting mass into galaxy number or mass");
 
- if( !no_poisson ) {
+ if( no_poisson_type==0 ) {
 
    cout<<"\n--- Apply poisson on galaxy number"<<endl;
@@ -575 +581 @@
      nm = fluct3d.MeanSigma2(rm,rs2);
    double xgalmin,xgalmax; fluct3d.MinMax(xgalmin,xgalmax,0.1,1.e50);
+   fluct3d.NTupleCheck(posobs,string("ntpois"),ntnent);
    PrtTim(">>>> End Apply poisson on galaxy number");
    if(wslice) {
@@ -598 +605 @@
      cout<<"Equivalent: "<<rm*nm/fluct3d.NPix()<<" Msol / pixels"<<endl;
      nm = fluct3d.MeanSigma2(rm,rs2);
+   fluct3d.NTupleCheck(posobs,string("ntmhi"),ntnent);
    PrtTim(">>>> End Convert Galaxy number to HI mass");
 
@@ -622 +630 @@
  ////  fluct3d.AddAGN(lfjy_agn,lsigma_agn,powlaw_agn);
  ////    nm = fluct3d.MeanSigma2(rm,rs2);
- ////   PrtTim(">>>> End Add AGN");
+ ////  fluct3d.NTupleCheck(posobs,string("ntagn"),ntnent);
+ ////  PrtTim(">>>> End Add AGN");
  ////}
 
@@ -632 +641 @@
    snoisesave = snoise;
      nm = fluct3d.MeanSigma2(rm,rs2);
+   fluct3d.NTupleCheck(posobs,string("ntnois"),ntnent);
    PrtTim(">>>> End Add noise");
  }
@@ -644 +654 @@
    PrtTim(">>>> End Scale cube");
  }
+ fluct3d.NTupleCheck(posobs,string("ntfin"),ntnent);
 
  //----------------------------------------------------------------- 
@@ -664 +675 @@
  cout<<endl<<"\n--- ReComputing spectrum from real space"<<endl;
  fluct3d.ReComputeFourier();
+ fluct3d.NTupleCheck(posobs,string("ntpkrec"),ntnent);
  PrtTim(">>>> End ReComputing spectrum");