Index: /trunk/WuerzburgSoft/Thomas/mphys/phys.C
===================================================================
--- /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1365)
+++ /trunk/WuerzburgSoft/Thomas/mphys/phys.C	(revision 1366)
@@ -35,5 +35,6 @@
     Double_t res = MPhoton::Int2(&Ep, k);
 
-    return res*1e55;
+    return res*1e55; //65/k[0];
+    // return MPhoton::Planck(&Ep, &k[1]);
 }
 
@@ -115,6 +116,6 @@
 void DoIt()
 {
-    Double_t R = 500; // [kpc]
-    Double_t startz = MParticle::ZofR(&R);
+    Double_t startz = 0.003;                   //MParticle::ZofR(&R);
+    Double_t R      = MParticle::RofZ(&startz); // [kpc]
 
     cout << "R = " << R << "kpc" << endl;
@@ -124,5 +125,5 @@
     MPairProduction pair;
 
-    Double_t runtime = 7*60*60; // [s]
+    Double_t runtime = 15*60; // [s]
 
     Double_t lo = 1e4; 
@@ -155,5 +156,5 @@
     MBinning binspoly2;
     binspolx.SetEdges(16, -180, 180);
-    binspoly1.SetEdges(20, 0, 2e-6);
+    binspoly1.SetEdges(20, 0, 5e-6);
     binspoly2.SetEdges(20, 0, 1e-5);
     MH::SetBinning(&angle,    &binspolx, &binspoly1);
@@ -343,6 +344,5 @@
 
                     if (fabs(e->GetTheta()*3437)<1 &&  // < 1min
-                        e->GetEnergy()>lo          &&
-                        e->GetEnergy()>3e4)
+                        e->GetEnergy()>lo)
                         continue;