In our application we use G4PrimaryParticle and SetDaughter method in order to force a specific decay channel. We noticed that the polarization which is set to the daughter by SetPolarization function is not used in the subsequent decay of the daughter particle (muon decay with spin). After enabling verbose tracking output we found that the polarization is indeed not set for the daughter particle. We found a possible bug in the G4PrimaryTransformer where in the transformation from Primary to Dynamic particle, the polarization is set for the mother particle but not for the decay products in the SetDecayProducts function. A simple fix would be to assign the polarization of the daughter primary particle to the newly created dynamic particle, like it is done for momentum, charge, etc.
Thank you for reporting the issue. We will fix it very soon, but the fix won't go with version 10.3 scheduled on this coming Friday. It will go with the next patch release, and we will also make a retroactive patch for 10.2. Kind regards, Makoto Asai
Created attachment 434 [details] G4PrimaryTransformer.cc with the fix (source code to be included in version 10.3)
The fix will be included in version 10.3. It may also be included in a retroactive patch for 10.2 as well. For your convenience, the fixed code is attached. Thank you for reporting this issue. Makoto