Dear all,
for what I understood, the WentzelVI model is a mixed scheme to simulate the single and multiple scattering, depending on a threshold scattering angle.
I have read that the threshold angle is 0.2 rad, for which collisions are considered as “single scattering” if theta>0.2 and multple scattering, if below.
Is there a way to change this default parameter to see how close I get to single scattering without increasing too much the computational time?
Thanks for helping!!
Hello,
please, note that this is a maximal limit angle. The real limit is dynamic - defined by momentum of a particle. This maximal limit is effective only at low momentum. It can be changed via UI command:
/process/msc/ThetaLimit 0.1 radian
The default value 0.2 is optimized using existing Physics Lists and number of test setups. It cannot be considered as the best value, extra study is needed to make a conclusion. The optimum may be different for different particle type - this is not known so far.
VI
Ok thank you very much.
I also read that together wiht the G4WentzelVIModel I need to have implemented the G4eSingleCoulombScatteringModel or the G4eCoulombScatteringModel. Correct?
What happens in case this is not implemented? Everything is treated as pure MSC?
Thanks again
In the reference EM physics configurations G4WentzelVIModel used in multiple scattering process, G4eCoulombScatteringModel - in G4CoulombScattering process. They share that same parameters, so angular limits are the same. EM constructors Opt0 and others configured correctly.
In custom physics all kind of mistakes are possible: double counting, holes in applicability of models…
VI