Changeset 1448 for trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc

Timestamp:

07/26/02 14:49:58 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc (modified) (5 diffs)

trunk/WuerzburgSoft/Thomas/mphys/MElectron.cc

@@ -371 +371 @@
 Bool_t MElectron::SetNewPositionB(Double_t B)
 {
+    if (B==0)
+        return SetNewPosition();
+
     static TRandom rand(0);
     Double_t x = rand.Exp(GetInteractionLength());
@@ -456 +459 @@
         p(0) = beta_p/beta*x;
         p(1) = GetPType()==kEElectron?r*sin(rho):-r*sin(rho);
-        p(2) = r*(1-cos(rho));
+        p(2) = r*(1.-cos(rho));
     }
     else
@@ -463 +466 @@
         p(1) = beta_o/beta*x;
         p(2) = 0;
+        cout << "------------- HEY! --------------" << endl;
     }
 
@@ -474 +478 @@
     p *= N2;
     p *= M2;
+
+    if (p(2)<x) // happens sometimes in case B==0
+    {
+        cout << "----- HA: " << B << " " << x << " " << p(2) << " " << x-p(2) << endl;
+        p(2)=x;
+    }
+
     s -= p;
 
@@ -479 +490 @@
     fPhi = atan2(s(1), s(0));
     fZ   = ZofR(&s(2));
+    fX  += x-p(2);
 
     // -----------------------------