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

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

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

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

@@ -25 +25 @@
 //
 //
-// $Id: G4EqEMFieldWithSpin.cc,v 1.8 2009/11/06 22:31:35 gum Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4EqEMFieldWithSpin.cc,v 1.9 2010/07/14 10:00:36 gcosmo Exp $
+// GEANT4 tag $Name: field-V09-03-03 $
 //
 //
@@ -33 +33 @@
 //  30.08.2007 Chris Gong, Peter Gumplinger
 //  14.02.2009 Kevin Lynch
-//  06.11.2009 Hiromi Iinuma see:
-//  http://hypernews.slac.stanford.edu/HyperNews/geant4/get/emfields/161.html
+//  06.11.2009 Hiromi Iinuma
 //
 // -------------------------------------------------------------------
@@ -44 +43 @@
 
 G4EqEMFieldWithSpin::G4EqEMFieldWithSpin(G4ElectroMagneticField *emField )
-  : G4EquationOfMotion( emField )
+  : G4EquationOfMotion( emField ), fElectroMagCof(0.), fMassCof(0.),
+    omegac(0.), anomaly(0.0011659208), pcharge(0.), E(0.), gamma(0.), beta(0.)
 {
-  anomaly = 0.0011659208;
 }
 
@@ -63 +62 @@
    omegac = 0.105658387*GeV/particleMass * 2.837374841e-3*(rad/cm/kilogauss);
 
-   ParticleCharge = particleCharge;
+   pcharge = particleCharge;
 
    E = std::sqrt(sqr(MomentumXc)+sqr(particleMass));
@@ -135 +134 @@
 
    G4ThreeVector dSpin
-     = ParticleCharge*omegac*( ucb*(Spin.cross(BField))-udb*(Spin.cross(u))
+     = pcharge*omegac*( ucb*(Spin.cross(BField))-udb*(Spin.cross(u))
                                // from Jackson
                                // -uce*Spin.cross(u.cross(EField)) );
                                // but this form has one less operation
-                               - uce*(u*(Spin*EField) - EField*(Spin*u)) );
+                      - uce*(u*(Spin*EField) - EField*(Spin*u)) );
 
    dydx[ 9] = dSpin.x();