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

Timestamp:

Jun 4, 2007, 7:02:03 PM (18 years ago)

Author:

cmv

Message:

add sigma computation with bad pixel cmv 04/06/2007

File:

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

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -435 +435 @@
  cout<<"1+rgen: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200,true,0.);
+ cout<<"1+rgen with_neg_a_zero: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+     <<rs2<<" -> "<<sqrt(rs2)<<endl;
  PrtTim(">>>> End Converting fluctuations into mass");
 
@@ -441 +444 @@
  cout<<rm<<" galaxies put into survey"<<endl;
  nm = fluct3d.MeanSigma2(rm,rs2,0.);
- cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+ cout<<"galaxy mass: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+     <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200,true,0.);
+ cout<<"galaxy mass with_neg_a_zero: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
  PrtTim(">>>> End Converting mass into galaxy number");
 
- cout<<"\n--- Set negative pixels to BAD"<<endl;
+ cout<<"\n--- Set negative and null pixels to BAD"<<endl;
  nm = fluct3d.SetToVal(0.,1e+200,-999.);
  cout<<nm<<" negative in survey set to BAD"<<endl;
  nm = fluct3d.MeanSigma2(rm,rs2,-998.);
- cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+ cout<<"galaxy mass: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
  PrtTim(">>>> End Set negative pixels to BAD etc...");
@@ -457 +463 @@
  cout<<nm<<" galaxies into survey after poisson"<<endl;
  nm = fluct3d.MeanSigma2(rm,rs2,-998.);
- cout<<"galaxy: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+ cout<<"galaxy number: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+     <<rs2<<" -> "<<sqrt(rs2)<<endl;
+ nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200,true,0.);
+ cout<<"galaxy number with_neg_a_zero: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
  PrtTim(">>>> End Apply poisson on galaxy number");
@@ -474 +483 @@
  double mhi = fluct3d.TurnNGal2Mass(tirhmdndm,true);
  cout<<mhi<<" MSol in survey / "<<mass_by_mpc3*fluct3d.GetVol()<<endl;
- nm = fluct3d.MeanSigma2(rm,rs2,0.);
+ nm = fluct3d.MeanSigma2(rm,rs2,-998.);
  cout<<"HI mass: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
      <<rs2<<" -> "<<sqrt(rs2)<<endl;
  cout<<"Equivalent: "<<rm*nm/fluct3d.NPix()<<" Msol / pixels"<<endl;
+ nm = fluct3d.MeanSigma2(rm,rs2,0.,1e200,true,0.);
+ cout<<"HI mass with_neg_a_zero: ("<<nm<<") Mean = "<<rm<<", Sigma^2 = "
+     <<rs2<<" -> "<<sqrt(rs2)<<endl;
  PrtTim(">>>> End Convert Galaxy number to HI mass");