Changeset 3344 in Sophya for trunk/Cosmo/SimLSS

Timestamp:

Oct 9, 2007, 10:43:47 AM (18 years ago)

Author:

cmv

Message:

set schechter mmin=1e8 -> 1e7 cmv 09/10/2007

Location:

trunk/Cosmo/SimLSS

Files:

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

trunk/Cosmo/SimLSS/cmvdefsurv.cc

@@ -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(1.e+6), lnx2=log10(1.e+13), dlnx=(lnx2-lnx1)/npt;
  double masshimpc3 = IntegrateFuncLog(sch,lnx1,lnx2,0.001,dlnx,10.*dlnx,6);
  cout<<"Mass density: "<<masshimpc3<<" Msol/Mpc^3"<<endl;
@@ -312 +312 @@
  cout<<"\nsch(mstar) = "<<sch(mstar)<<" /Mpc^3/Msol"<<endl;
  cout<<"Galaxy number density:"<<endl;
- for(double x=lnx1; x<10.5; x+=1.) {
+ for(double x=lnx1; x<lnx2-0.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, "

trunk/Cosmo/SimLSS/cmvobserv3d.cc

@@ -89 +89 @@
  double alpha = -1.37;
 
- double schmin=1e8, schmax=1e13;
+ double schmin=1.e7, schmax=1.e13;
  int schnpt = -100;
  bool use_schmassdist = false;

trunk/Cosmo/SimLSS/cmvschdist.cc

@@ -6 +6 @@
 // ex: on a un pixel qui a N galaxies, quelle est la distribution
 //     statistique de masse dans un tel pixel ?
-// > cmvschdist -a -v -m 1e+8,1e+13 -n -100 -r 200 -N 100000
+// > cmvschdist -a -v -m 1e+7,1e+13 -n -100 -r 200 -N 100000
 // compare with:
-// > cmvschdist -a -v -m 1e+8,1e+13 -n -100 -r 200 -N 100000 -0
+// > cmvschdist -a -v -m 1e+7,1e+13 -n -100 -r 200 -N 100000 -0
 // Fabriquer un fichier pour la prod
-// > cmvschdist -a -m 1e+8,1e+13 -n -100 -r 2000 -N 5000000 -W schdist_2000.ppf
+// > cmvschdist -a -m 1e+7,1e+13 -n -100 -r 2000 -N 5000000 -W schdist_2000.ppf
 #include <iostream>
 #include <stdlib.h>
@@ -52 +52 @@
  cout<<"h75= "<<h75<<" nstar= "<<nstar<<"  mstar="<<mstar<<"  alpha="<<alpha<<endl;
 
- double xmin=1e8, xmax=1e13;
+ double xmin=1.e7, xmax=1.e13;
  int npt = -100;
  int ngmax = 200, ngmin = 1;