Changeset 1361 for trunk/WuerzburgSoft/Thomas/mphys

Timestamp:

06/13/02 08:58:22 (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

@@ -116 +116 @@
 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 +128 @@
 
     Double_t lo = 1e4; 
-    Double_t hi = 1e10;
+    Double_t hi = 1e11;
     Double_t alpha = -2;
 
@@ -155 +155 @@
     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 +321 @@
 
                     // 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 +339 @@
 
                     if (fabs(e->GetTheta()*3437)<1 &&  // < 1min
-                        e->GetEnergy()>lo*2)
+                        e->GetEnergy()>lo          &&
+                        e->GetEnergy()>3e4)
                         continue;