Message: Re: Spin Tracking in B-field Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Spin Tracking in B-field 

Forum: Fields: Magnetic and Otherwise
Re: Question Spin Tracking in B-field (Victor)
Date: 29 Apr, 2014
From: Gumplinger Peter <Gumplinger Peter>

Hallo VK,

> Now, my really simple setup is Global Magnetic Field(in x-dir) with
> G4UniformMagneticField, plus I've included exactly the same code as in
> the Manual for Spin Tracking(or the post link just above). I've set
> Beam's polarization as (0, 0, 1) and it just hits a target. I use
> QGSP_BERT physics list. Upon generation of the primary(I shoot 100GeV
> e-) I can see that it has proper polarization vector. But when I catch
> it in the Sensetive Detector(World is SD for now...) and do
> GetPolarization it gives me (0,0,0) for the primary

How do you identify the 'primary' e-; with:

aTrack->GetParentID() == 0

and (&&) for good measure also:

aTrack->GetTrackID() == 1

You can also follow the progression of the spin with G4SteppingVerbose by setting the verbosity to 4.

/tracking/verbose 4

> and for the
> secondaries....

Secondaries won't have spin/polarization unless the process creating them also specifies their polarization. The only processes that do that are in:

/processes/electromagnetic/polarisation

and G4eplusAnnihilation, G4eeToTwoGammaModel, G4SynchrotronRadiation(InMat), G4Scintillation anhd G4Cerenkov - but the latter are all secondary photons.

Peter

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

1 Question: Re: Spin Tracking in B-field   (Victor - 29 Apr, 2014)
1 Question: Re: Spin Tracking in B-field   (Victor - 29 Apr, 2014)
1 Feedback: Re: Spin Tracking in B-field   (Gumplinger Peter - 01 May, 2014)
... 2 Message(s)
2 More: Re: Spin Tracking in B-field   (Gumplinger Peter - 01 May, 2014)
 Add Message Add Message
to: "Re: Spin Tracking in B-field"

 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 ]