Changeset 1364 for trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc

Timestamp:

06/14/02 09:03:11 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc (modified) (1 diff)

trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc

@@ -71 +71 @@
      Double_t z1 = x[0] + 1;
 
-     const Double_t c  = 299792458;                 // [m/s]
-     const Double_t H0 = 50./3.0857e19;             // [km / Mpc s] --> [s^-1]
-
-     const Double_t ly = 3600.*24.*365.*c;          // [m/ly]
-     const Double_t pc = 1./3.258;                  // [pc/ly]
+     const Double_t c  = 299792458;                   // [m/s]
+     const Double_t H0 = 50./3.0857e19;               // [km / Mpc s] --> [s^-1]
+
+     const Double_t ly = 3600.*24.*365.*c;            // [m/ly]
+     const Double_t pc = 1./3.258;                    // [pc/ly]
 
      const Double_t R = c/H0 * 2 * (1 - 1./sqrt(z1)); // [m]
 
-     return  R * pc/ly/1e3;                   // [kpc]
+     return  R * pc/ly/1e3;                           // [kpc]
+}
+
+Double_t MParticle::Planck(Double_t *x, Double_t *k)
+{
+    //
+    // Planck, per unit volume, per unit energy
+    //
+    // constants (see below) moved out of function
+    //
+    const Double_t E   = x[0];                    // [GeV]
+    const Double_t z   = k ? k[0] : 0;
+
+    const Double_t T   = 2.96*(z+1);              // [K]
+    const Double_t e   = 1.602176462e-19;         // [C]
+    const Double_t kB  = 1e-9/e*1.3806503e-23;    // [GeV/K]
+
+    const Double_t EkT = E/kB/T;
+
+    /*
+     //Double_t c   = 299792458;             // [m/s]
+     //Double_t h   = 1e-9/e*6.62606876e-34; // [GeVs]
+     //Double_t hc  = h*c;                   // [GeVm]
+     Double_t konst = 4.*TMath::Pi() * 2. / (hc*hc*hc);
+     return konst * E*E / (exp(EkT)-1.); // [1 / GeV / m^3 ]
+     */
+
+    return E*E / (exp(EkT)-1.); // [GeV^2]
+
 }