Changeset 3193 in Sophya for trunk/Cosmo/SimLSS/schechter.cc

Timestamp:

Mar 21, 2007, 7:18:25 PM (19 years ago)

Author:

cmv

Message:

petis changements cmv 21/03/2007

File:

• trunk/Cosmo/SimLSS/schechter.cc (modified) (2 diffs)

trunk/Cosmo/SimLSS/schechter.cc

@@ -54 +54 @@
 }
 
+void Schechter::Print(void)
+{
+  cout<<"Schechter::Print: nstar="<<nstar_<<" Mpc^-3"
+      <<"  mstar="<<mstar_<<" MSol"
+      <<"  alpha="<<alpha_
+      <<"  (outvalue="<<outvalue_<<" -> return ";
+  if(outvalue_) cout<<"m*dn/dm)"; else cout<<"dn/dm)";
+  cout<<endl;
+}
+
 ///////////////////////////////////////////////////////////
 //******************* Le Flux a 21cm ********************//
@@ -64 +74 @@
 // Return:
 //    le flux total emis a 21 cm en W/m^2
-// Ref: Binney & Merrifield, Galactic Astronomy p474 (ed:1998)
-//      J.Peterson & K.Bandura, astroph-0606104  (eq 7)
-//      F.Abdalla & Rawlings, astroph-0411342 (eq 7 may pb de d'unites?)
+// Ref:
+// -- Binney & Merrifield, Galactic Astronomy p474 (ed:1998)
+// S(W/m^2) = 1e-26 * Nu_21cm(Hz) * m(en masse solaire) /(2.35e5 * C(km/s) * Dlum^2)
+//          = 2.0e-28 * m / Dlum^2
+// -- J.Peterson & K.Bandura, astroph-0606104  (eq 7)
+//    F.Abdalla & Rawlings, astroph-0411342 (eq 7 mais pb de d'unites?)
+//          (A_21cm = 2.86888e-15 s^-1)
+// S(W/m^2) = 3/4 * h(J.s) * Nu_21cm(Hz) * A_21cm(s^-1) * Msol(kg)/m_proton(kg)
+//                         / (4 Pi) / (Dlum(Mpc) * 3.0857e22(m/Mpc))^2
+//          = 2.0e-28 * m / Dlum^2
+//-------------
 {
-  return  m / (2.35e+5 * d*d )
-          * Jansky2Watt_cst
-          * (Fr_HyperFin_Par*1e+9)/SpeedOfLight_Cst;
+  return  2.0e-28 * m / (d*d);
 }
-