Index: /trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc	(revision 1357)
+++ /trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc	(revision 1358)
@@ -288,6 +288,4 @@
 }
 
-#include <iostream.h>
-
 MPhoton *MElectron::DoInvCompton(Double_t theta)
 {
Index: /trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc	(revision 1357)
+++ /trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc	(revision 1358)
@@ -120,5 +120,5 @@
     const Double_t f = sqrt(sqrbetah*sqrbetae)*salpha;
 
-    cout << " {" << f << "," << E << "," << atan(tan1) << "," << atan(-tan2) << "} " << flush;
+    // cout << " {" << f << "," << E << "," << atan(tan1) << "," << atan(-tan2) << "} " << flush;
 
     MElectron &p0 = *new MElectron(E*(1.-f), gamma->GetZ());
Index: /trunk/WuerzburgSoft/Thomas/mphys/phys.C
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1357)
+++ /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1358)
@@ -212,15 +212,9 @@
                 if (phot.Integral(-log10(Eg)-8, -10.5)==0)
                 {
-                    // ???????????????????
-                    hist.Fill(Eg, Eg*Eg*weight);
-                    angle.Fill(p->GetPsi()*kRad2Deg, p->GetTheta()*kRad2Deg, weight);
-                    hista.Fill(p->GetTheta()*kRad2Deg, weight);
-                    delete listg.Remove(p);
                     cout << "z" << flush;
                     continue;
                 }
 
-                Double_t Ep = pow(10, phot.GetRandom());
-
+                Double_t Ep    = pow(10, phot.GetRandom());
                 Double_t theta = RandomTheta(Eg, Ep);
                 if (theta==0)