Changeset 1507 for trunk/WuerzburgSoft/Thomas/mphys/MCascade.cc

Timestamp:

08/19/02 08:58:09 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/MCascade.cc (modified) (7 diffs)

trunk/WuerzburgSoft/Thomas/mphys/MCascade.cc

@@ -50 +50 @@
 }
 
-Double_t RandomTheta(Double_t Eg, Double_t Ep)
+Double_t RandomThetaG(Double_t Eg, Double_t Ep)
 {
     Double_t E0 = 511e-6; // [GeV]
@@ -59 +59 @@
         return 0;
 
-    static TF1 func("RndTheta", Sbar_sigmas, 0, 0, 0);
+    static TF1 func("RndThetaG", Sbar_sigmas, 0, 0, 0);
 
     func.SetRange(0, log10(f));
@@ -139 +139 @@
         }
 
-        // WRONG!
-        MPhoton *p = e.DoInvCompton(0);
+        MPhoton *p = e.DoInvCompton();
 
         fBranchElectrons->GetTree()->Fill();
@@ -200 +199 @@
 
         Double_t Ep = pow(10, pe);
-        Double_t theta = RandomTheta(Eg, Ep);
+        Double_t theta = RandomThetaG(Eg, Ep);
         if (theta==0)
         {
@@ -316 +315 @@
 MCascade::MCascade()
 {
+    if (gRandom)
+        delete gRandom;
+
     TRandom r(0);
-    delete gRandom;
     gRandom = new TRandom3(r.GetSeed());
 
@@ -332 +333 @@
 {
     delete gRandom;
+    gRandom = 0;
 }
 
@@ -373 +375 @@
 
     // ------------------------------
+
+    cout << endl;
 
     cout << "R = " << fSrcR << "kpc" << endl;