Changeset 921 for trunk/source/geometry/magneticfield/src/G4EqEMFieldWithSpin.cc

Timestamp:

Feb 16, 2009, 10:14:30 AM (15 years ago)

Author:

garnier

Message:

en test de gl2ps. Problemes de libraries

File:

• trunk/source/geometry/magneticfield/src/G4EqEMFieldWithSpin.cc (modified) (4 diffs)

trunk/source/geometry/magneticfield/src/G4EqEMFieldWithSpin.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4EqEMFieldWithSpin.cc,v 1.2 2008/04/24 12:33:08 tnikitin Exp $
-// GEANT4 tag $Name: HEAD $
+// $Id: G4EqEMFieldWithSpin.cc,v 1.4 2008/11/21 21:17:03 gum Exp $
+// GEANT4 tag $Name: geant4-09-02-cand-01 $
 //
 //
@@ -40 +40 @@
 
 #include "G4EqEMFieldWithSpin.hh"
+#include "G4ElectroMagneticField.hh"
 #include "G4ThreeVector.hh"
 #include "globals.hh"
 
 G4EqEMFieldWithSpin::G4EqEMFieldWithSpin(G4ElectroMagneticField *emField )
-      : G4EquationOfMotion( emField ) { anomaly = 1.165923e-3; }
+      : G4EquationOfMotion( emField )
+{ 
+  anomaly = 0.0011659208;
+}
+
+G4EqEMFieldWithSpin::~G4EqEMFieldWithSpin()
+{
+} 
 
 void  
 G4EqEMFieldWithSpin::SetChargeMomentumMass(G4double particleCharge, // e+ units
-                                            G4double MomentumXc,
+                                            G4double MomentumXc,
                                             G4double particleMass)
 {
@@ -61 +69 @@
    beta  = MomentumXc/E;
    gamma = E/particleMass;
+
 }
-
-
 
 void
 G4EqEMFieldWithSpin::EvaluateRhsGivenB(const G4double y[],
-                                        const G4double Field[],
-                                              G4double dydx[] ) const
+                                       const G4double Field[],
+                                             G4double dydx[] ) const
 {
 
@@ -114 +121 @@
 
    G4ThreeVector Spin(y[9],y[10],y[11]);
+
+   if (Spin.mag() > 0.) Spin = Spin.unit();
+
    G4ThreeVector dSpin;