Changeset 1366

Timestamp:

06/14/02 09:14:49 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/phys.C (modified) (5 diffs)

trunk/WuerzburgSoft/Thomas/mphys/phys.C

@@ -35 +35 @@
     Double_t res = MPhoton::Int2(&Ep, k);
 
-    return res*1e55;
+    return res*1e55; //65/k[0];
+    // return MPhoton::Planck(&Ep, &k[1]);
 }
 
@@ -115 +116 @@
 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 +125 @@
     MPairProduction pair;
 
-    Double_t runtime = 7*60*60; // [s]
+    Double_t runtime = 15*60; // [s]
 
     Double_t lo = 1e4; 
@@ -155 +156 @@
     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 +344 @@
 
                     if (fabs(e->GetTheta()*3437)<1 &&  // < 1min
-                        e->GetEnergy()>lo          &&
-                        e->GetEnergy()>3e4)
+                        e->GetEnergy()>lo)
                         continue;