|Message: Re: Spin tracking of muons at rest||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Thanks Kevin, for both of your two insightful postings.
2) The G4Decay class itself (G4DecayWithSpin inherits from G4Decay) does not define decay modes of the particle. Geant4 provides two ways of doing this: using G4DecayChannel in G4DecayTable, and using thePreAssignedDecayProducts of G4DynamicParticle. The G4Decay class calculates the PhysicalInteractionLength and boosts decay products created by G4VDecayChannel or event generators.
3) Yes, G4DecayWithSpin retrieves or sets a particles polarization and evolves it through the remainder of a particle's lifetime. (this is crucial to your ultimate question).
> But if I attach G4DecayWithSpin to idxPostStep, do I simultaneously get BOTH
> G4Mag_EqRhs BMT evolution, AND the G4DecayWithSpin evolution?
The method G4DecayWithSpin::Spin_Precession takes the argument fRemainderLifeTime. For decays in flight, this will be zero, and consequently, the calculated remaining spin precession angle will also be zero.
> Should I attach G4Decay to idxPostStep instead of G4DecayWithSpin?
This will not work because only G4DecayWithSpin has the code:
parent_polarization = aParticle->GetPolarization(); decaychannel->SetPolarization(parent_polarization);
necessary to set the spin direction in the DecayChannelWithSpins.
> can you point out how the code in the distribution ensures
> the G4DecayWithSpin terms don't come in to play for moving muons?
rotationangle = fRemainderLifeTime * omega = 0. * omega
You can verify that fRemainderLifeTime has already been updated (reduced to zero) by the time the G4DecayWithSpin::DecayIt is called for a decay in flight.
|Inline Depth:||Outline Depth:||Add message:|