Message: Re: Charge, momentum and mass are not enough to compute the force Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Charge, momentum and mass are not enough to compute the force 

Forum: Fields: Magnetic and Otherwise
Re: Warning Charge, momentum and mass are not enough to compute the force (Bertalan Juhasz)
Date: 03 Mar, 2005
From: Gumplinger Peter <Gumplinger Peter>

Hello Bertalam,

Thank you for your interesting contribution to this forum.

Have you come across how a particle's spin is propagated in Geant4? I think the magnetic moment can be propagated in the same way.

Please, take a look at:

source/geometry/magneticfield/include/G4Mag_SpinEqRhs.hh source/geometry/magneticfield/include/G4Mag_SpinEqRhs.cc

and sorry, I should have written documentation about 'spin tracking' by now (since 7.0).

In DetectorConstruction:

#include "G4Mag_SpinEqRhs.hh"

G4Mag_EqRhs* fEquation = new G4Mag_SpinEqRhs(magField);

G4MagIntegratorStepper* pStepper = new G4ClassicalRK4(fEquation,12);

(notice the 12!)

In PrimaryGenerator:

particleGun->SetParticlePolarization(....)

where you set the initial mag. moment direction.

In your equivalent to G4Mag_SpinEqRhs the meaning of 'spin' becomes mag. moment and in general, the meaning in G4 of 'Polarization' becomes mag. moment.

I think the machinery is already there, all you got to do is write the G4Mag_MagMomentEqRhs (which sounds like you have already done) and set the 'polarization' initially.

Let me know how you are doing implementing propagation of mag. moments in this way.

Peter

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Charge, momentum and mass are not enough to compute the force   (Bertalan Juhasz - 04 Mar, 2005)
(_ Feedback: Re: Charge, momentum and mass are not enough to compute the force   (Peter Gumplinger - 04 Mar, 2005)
 Add Message Add Message
to: "Re: Charge, momentum and mass are not enough to compute the force"

 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 ]