Changeset 1369 for trunk/WuerzburgSoft/Thomas/mphys

Timestamp:

06/17/02 13:55:42 (23 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/WuerzburgSoft/Thomas/mphys

Files:

• Changelog (modified) (1 diff)
• MElectron.cc (modified) (1 diff)
• MParticle.cc (modified) (1 diff)

trunk/WuerzburgSoft/Thomas/mphys/Changelog

@@ -1 +1 @@
                                                                   -*-*- END -*-*-
+ 2002/06/17: Thomas Bretz
+
+   * MElectron.cc:
+     - Added SetNewDirection for the electron itself
+
+   * MParticle.cc:
+     - changed sqrt() to asin(sqrt()) - no effect (but correct)
+
+
+
  2002/06/14: Thomas Bretz
 

trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc

@@ -319 +319 @@
     fEnergy -= e;
 
+    const Double_t phi = rand.Uniform(TMath::Pi()*2);
+
     MPhoton &p = *new MPhoton(e, fZ);
     p = *this;
-    p.SetNewDirection(theta1, rand.Uniform(TMath::Pi()*2));
+    p.SetNewDirection(theta1, phi);
+
+    const Double_t beta2  = sqrt(1.-E0/fEnergy*E0/fEnergy);
+    const Double_t theta2 = asin((epsilon*sin(t)-e*sin(theta1))/fEnergy/beta2);
+
+    SetNewDirection(theta2, phi);
 
     return &p;

trunk/WuerzburgSoft/Thomas/mphys/MParticle.cc

@@ -135 +135 @@
     r *= B;
 
-    fTheta = sqrt(r(0)*r(0)+r(1)*r(1)); // Numerically bad: acos(r(2));
+    fTheta = asin(sqrt(r(0)*r(0)+r(1)*r(1))); // Numerically bad: acos(r(2));
     fPsi   = atan2(r(1), r(0));