Index: /trunk/MagicSoft/Mars/Changelog
===================================================================
--- /trunk/MagicSoft/Mars/Changelog	(revision 1360)
+++ /trunk/MagicSoft/Mars/Changelog	(revision 1361)
@@ -1,3 +1,10 @@
                                                                   -*-*- END -*-*-
+ 2002/06/13: Thomas Bretz
+
+   * mdata/MData.cc:
+     - changed ostream.h to ofstream.h to make it compile on alpha
+
+
+
  2002/06/10: Thomas Bretz
 
Index: /trunk/MagicSoft/Mars/mdata/MData.cc
===================================================================
--- /trunk/MagicSoft/Mars/mdata/MData.cc	(revision 1360)
+++ /trunk/MagicSoft/Mars/mdata/MData.cc	(revision 1361)
@@ -48,5 +48,5 @@
 #include "MData.h"
 
-#include <ostream.h>
+#include <fstream.h>
 
 ClassImp(MData);
Index: /trunk/WuerzburgSoft/Thomas/mphys/phys.C
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1360)
+++ /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1361)
@@ -116,8 +116,8 @@
 void DoIt()
 {
-    Double_t R = 100; // [kpc]
-    Double_t startz = MParticle::ZofR(&R);
-
-    cout << "R = " << R << "kpc" << endl;
+    //    Double_t R = 1798; // [kpc]
+    Double_t startz = 0.003; //MParticle::ZofR(&R);
+
+    //    cout << "R = " << R << "kpc" << endl;
     cout << "Z = " << startz << endl;
 
@@ -128,5 +128,5 @@
 
     Double_t lo = 1e4; 
-    Double_t hi = 1e10;
+    Double_t hi = 1e11;
     Double_t alpha = -2;
 
@@ -155,5 +155,5 @@
     MBinning binspoly;
     binspolx.SetEdges(16, -180, 180);
-    binspoly.SetEdgesLog(20, 1e-13, 1e-4);
+    binspoly.SetEdgesLog(20, 1e-10, 1e-3);
     MH::SetBinning(&position, &binspolx, &binspoly);
     MH::SetBinning(&angle,    &binspolx, &binspoly);
@@ -321,5 +321,5 @@
 
                     // WRONG!
-                    Double_t theta = rand.Uniform(TMath::Pi()*2);
+                    Double_t theta = rand.Uniform(TMath::Pi()/2)+TMath::Pi()*3/4;
                     MPhoton *p = e->DoInvCompton(theta);
 
@@ -339,5 +339,6 @@
 
                     if (fabs(e->GetTheta()*3437)<1 &&  // < 1min
-                        e->GetEnergy()>lo*2)
+                        e->GetEnergy()>lo          &&
+                        e->GetEnergy()>3e4)
                         continue;