Changeset 3343 in Sophya

Timestamp:

Oct 8, 2007, 6:17:46 PM (18 years ago)

Author:

cmv

Message:

IMPORTANT: corrcetion bug Mstar article Zwaan h75^{2 , cmv 08/10/2007}

Location:

trunk/Cosmo/SimLSS

Files:

• cmvdefsurv.cc (modified) (6 diffs)
• cmvobserv3d.cc (modified) (1 diff)
• cmvschdist.cc (modified) (1 diff)
• cmvtstsch.cc (modified) (1 diff)

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -69 +69 @@
  double h75 = h100 / 0.75;
  double nstar = 0.006*pow(h75,3.);  //
- double mstar = pow(10.,9.8/(h75*h75));  // MSol
+ double mstar = pow(10.,9.8);  // MSol
  double alpha = -1.37;
  cout<<"nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;
@@ -299 +299 @@
  cout<<"mstar*sch(mstar) = "<<sch(mstar)<<" Msol/Mpc^3/Msol"<<endl;
  int npt = 10000;
- double lnx1=log10(1e-6), lnx2=log10(1e+14), dlnx=(lnx2-lnx1)/npt;
+ double lnx1=log10(1e+6), lnx2=log10(1e+14), dlnx=(lnx2-lnx1)/npt;
  double masshimpc3 = IntegrateFuncLog(sch,lnx1,lnx2,0.001,dlnx,10.*dlnx,6);
  cout<<"Mass density: "<<masshimpc3<<" Msol/Mpc^3"<<endl;
@@ -308 +308 @@
      <<"Total mass in survey = "<<masshitot<<" Msol"<<endl;
  if(mhiref<=0.)  mhiref = masshipix;
+
+ sch.SetOutValue(0);
+ cout<<"\nsch(mstar) = "<<sch(mstar)<<" /Mpc^3/Msol"<<endl;
+ cout<<"Galaxy number density:"<<endl;
+ for(double x=lnx1; x<10.5; x+=1.) {
+   double n = IntegrateFuncLog(sch,x,lnx2,0.001,dlnx,10.*dlnx,6);
+   cout<<"  m>"<<x<<" Msol: "<<n<<" /Mpc^3, "<<n*dvol<<" /pixel, "
+       <<n*vol<<" in survey"<<endl;
+ }
+ sch.SetOutValue(1);
+
 
  // --- Survey values
@@ -368 +379 @@
  double tsig = facpolar * tsig_2polar;
  double ssig = facpolar * ssig_2polar;
- cout<<"\nSignal("<<mhiref<<" Msol): P="<<psig<<" W"<<endl;
- cout<<"     flux density = "<<ssig<<" Jy  (for Dnu="<<dnuhiz<<" GHz)"<<endl;
- cout<<"     Antenna temperature: tsig="<<tsig<<" K"<<endl;
+ cout<<"\nSignal("<<mhiref<<" Msol):"<<endl
+     <<"     P="<<psig<<" W"<<endl
+     <<"     flux density = "<<ssig*1.e6<<" mu_Jy  (for Dnu="<<dnuhiz<<" GHz)"<<endl
+     <<"     Antenna temperature: tsig="<<tsig<<" K"<<endl;
 
  // Elargissement doppler de la raie a 21cm: dNu = vrot/C * Nu(21cm) / (1+z)
@@ -390 +402 @@
  cout<<"\nSynchrotron: T="<<Tsynch408<<" K ("<<nuhaslam<<" GHz), "
      <<tsynch<<" K ("<<nuhiz<<" GHz)"<<endl
-     <<"             P="<<psynch<<" W for pixel"<<endl;
- cout<<"             flux density = "<<ssynch<<" Jy for pixel solid angle"<<endl;
+     <<"             P="<<psynch<<" W for pixel"<<endl
+     <<"             flux density = "<<ssynch<<" Jy for pixel solid angle"<<endl;
 
  // CMB
@@ -400 +412 @@
  double pcmb = facpolar * pcmb_2polar;
  double scmb = facpolar * scmb_2polar;
- cout<<"\nCMB: T="<<tcmb<<" K -> P="<<pcmb<<" W for pixel"<<endl;
- cout<<"                       flux density = "<<scmb<<" Jy for pixel solid angle"<<endl;
+ cout<<"\nCMB: T="<<tcmb<<" K"<<endl
+     <<"     P="<<pcmb<<" W for pixel"<<endl
+     <<"     flux density = "<<scmb<<" Jy for pixel solid angle"<<endl;
 
  // AGN

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -86 +86 @@
  double h75 = h100 / 0.75;
  double nstar = 0.006*pow(h75,3.);
- double mstar = pow(10.,9.8/(h75*h75));  // MSol
+ double mstar = pow(10.,9.8);  // MSol
  double alpha = -1.37;
 

trunk/Cosmo/SimLSS/cmvschdist.cc

@@ -48 +48 @@
  double h75 = 0.71 / 0.75;
  double nstar = 0.006*pow(h75,3.);  //  
- double mstar = pow(10.,9.8/(h75*h75));  // MSol
+ double mstar = pow(10.,9.8);  // MSol
  double alpha = -1.37;
  cout<<"h75= "<<h75<<" nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;

trunk/Cosmo/SimLSS/cmvtstsch.cc

@@ -37 +37 @@
  double h75 = 0.71 / 0.75;
  double nstar = 0.006*pow(h75,3.);  //  
- double mstar = pow(10.,9.8/(h75*h75));  // MSol
+ double mstar = pow(10.,9.8);  // MSol
  double alpha = -1.37;
  cout<<"h75= "<<h75<<" nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;