Index: /trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc	(revision 1364)
+++ /trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc	(revision 1365)
@@ -12,18 +12,9 @@
 ClassImp(MParticle);
 
-/***********************
+/**************************************************
  *
- * calc r from z:
+ * H0 = 50./3.0857e19; // [km / Mpc s] --> [s^-1]
  *
- *        R = 2*c/H0 *(1+z - sqrt(1+z))
- *
- *        z = -0.5 * (3 + R' +/- sqrt(R'+1))   R' = R*H0/c
- *
- *        H0 = 50./3.0857e19; // [km / Mpc s] --> [s^-1]
- *
- ************************
- */
-#include <iostream.h>
-#include <iomanip.h>
+ **************************************************/
 
 Double_t MParticle::ZofR(Double_t *x, Double_t *k)
@@ -50,5 +41,4 @@
     const Double_t R = 1./(1-r*H0/c/2);   // [1]
 
-
     return R*R - 1;
 }
@@ -67,5 +57,5 @@
      const Double_t R = c/H0 * 2 * (z1 - sqrt(z1)); // [m]
 
-     return  R * pc/ly/1e3;                   // [kpc]
+     return  R * pc/ly/1e3;                         // [kpc]
      */
      Double_t z1 = x[0] + 1;
@@ -79,5 +69,5 @@
      const Double_t R = c/H0 * 2 * (1 - 1./sqrt(z1)); // [m]
 
-     return  R * pc/ly/1e3;                           // [kpc]
+     return R * pc/ly/1e3;                            // [kpc]
 }
 
@@ -107,5 +97,4 @@
 
     return E*E / (exp(EkT)-1.); // [GeV^2]
-
 }
 
@@ -115,15 +104,5 @@
     //
     //  default constructor
-    //  creates an a list of histograms for all pixels and both gain channels
     //
-
-    //
-    //   set the name and title of this object
-    //
-
-    /*
-     fName  = name  ? name  : "MParticle";
-     fTitle = title ? title : "Container for a particle";
-     */
 }
 
@@ -159,6 +138,8 @@
     fPsi   = atan2(r(1), r(0));
 
-    if (fTheta*2 > TMath::Pi())
+    /*
+     if (fTheta*2 > TMath::Pi())
         fTheta = fabs(fTheta-TMath::Pi());
+     */
 }