Changeset 1361

Timestamp:

06/13/02 08:58:22 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk

Files:

• MagicSoft/Mars/Changelog (modified) (1 diff)
• MagicSoft/Mars/mdata/MData.cc (modified) (1 diff)
• WuerzburgSoft/Thomas/mphys/phys.C (modified) (5 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                                   -*-*- 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
 

trunk/MagicSoft/Mars/mdata/MData.cc

@@ -48 +48 @@
 #include "MData.h"
 
-#include <ostream.h>
+#include <fstream.h>
 
 ClassImp(MData);

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;