Changeset 1363 for trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc

Timestamp:

06/13/02 16:36:19 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc (modified) (3 diffs)

trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc

@@ -79 +79 @@
     //  theta: interaction angle [rad]
     //
-
-
     const Double_t E0     = 511e-6;              // [GeV]
     const Double_t Eg     = gamma->GetEnergy();  // [GeV]
@@ -98 +96 @@
     const Double_t GammaH = (Eg+Ep)/sqrt(s);
 
-    Double_t psi = stheta/(GammaH*(ctheta-sqrt(sqrbetah)));
+    const Double_t psi = stheta/(GammaH*(ctheta-sqrt(sqrbetah)));
 
     fAngle->SetParameter(0, sqrt(sqrbetae));
@@ -104 +102 @@
     const Double_t alpha = psi-acos(fAngle->GetRandom());
 
-    Double_t salpha = sin(alpha);
-    Double_t calpha = cos(alpha);
+    const Double_t salpha = sin(alpha);
+    const Double_t calpha = cos(alpha);
 
     const Double_t tphi = stheta/(Eg/Ep+ctheta); // tan(phi)
 
-    Double_t bb = sqrt(sqrbetah/sqrbetae);
+    const Double_t bb = sqrt(sqrbetah/sqrbetae);
 
-    Double_t s1 = calpha/GammaH;
-    Double_t s2 = tphi*s1 - salpha - bb;
+    const Double_t s1 = calpha/GammaH;
+    const Double_t s2 = tphi*s1 - salpha - bb;
 
-    Double_t tan1 = ((salpha+bb)*tphi+s1)/s2;
-    Double_t tan2 = ((salpha-bb)*tphi+s1)/s2;
+    const Double_t tan1 = ((salpha+bb)*tphi+s1)/s2;
+    const Double_t tan2 = ((salpha-bb)*tphi+s1)/s2;
 
     const Double_t E = (Eg+Ep)/2;;