Message: Re: Problems with transportation with spin Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Problems with transportation with spin 

Forum: Fast Simulation, Transportation & Others
Re: None Problems with transportation with spin (Tom Roberts)
Date: 15 Mar, 2012
From: Tom Roberts <Tom Roberts>

On Thu, 15 Mar 2012 20:50:42 GMT, Tom Roberts wrote:
> 3. The constructor of G4Mag_SpinEqRhs takes a G4MagneticField, not a
>    G4ElectroMagneticField. As a result, the particles do not see any electric
>    field -- neither spin nor momentum are changed by an E field.

I just stumbled across G4EqEMFieldWithSpin. It appears to resolve this -- a quick test shows that a mu+ moves in a constant E field just as for the usual equation of motion without spin, and the spin precessed a little bit.

In my class derived from G4ElectroMagneticField, I simply replaced
    G4EquationOfMotion *e = new G4Mag_SpinEqRhs((G4MagneticField*)this);
    pStepper = new G4ClassicalRK4(e,12);
    G4EquationOfMotion *e = new G4EqEMFieldWithSpin(this);
    pStepper = new G4ClassicalRK4(e,12);

Tom Roberts

 Add Message Add Message
to: "Re: Problems with transportation with spin"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews

[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]