I’m doing the migration from Geant4.9.4 to Geant4.10.4.
With Geant4.9.4, we construct EM physics for all other charged particles except gamma, e+, e-, mu+,mu- like this:
G4hMultipleScattering* aMultipleScattering = new G4hMultipleScattering();
G4hLowEnergyIonisation* ahadronLowEIon = new G4hLowEnergyIonisation();
pmanager->AddProcess(aMultipleScattering,-1,1,1);
pmanager->AddProcess(ahadronLowEIon,-1,2,2);
ahadronLowEIon->SetFluorescence(false);
As I know, G4hLowEnergyIonisation has been removed since Geant4.9.5.
With Geant4.10.4, we just use G4EmLivermorePhysics. But during the validation between 9.4 and 10.4,
we found obvious differences for the energy loss of alpha and proton(including number of steps, dE/dx, energy deposit). Which one is more correct? What process is the sustitute of G4hLowEnergyIonisation in Geant4.10.4?
G4hIonisation is now capable both for high and low energy.
I would suggest you in 10.4 trying to use Geant4 components (physics constructors) available for 10.X seria.
See examples: $G4INSTALL/examples/extended/electromagnetic/TestEm7
For the reference Physics Lists $G4INSTALL/examples/extended/hadronic/Hadr00
Thank you for your reply.
I noticed that in G4EmLivermorePhysics,
For alpha:
G4ionIonisation* ionIoni = new G4ionIonisation();
ionIoni->SetStepFunction(0.1, 10*um);
For proton:
G4hIonisation* hIoni = new G4hIonisation();
hIoni->SetStepFunction(0.2, 50*um);
For different physics constructors, the parameters are different.
Could you please explain more about these two parameters of SetStepFunction?
Could these parameters affect the distribution of dE/dx, Step numbers, energy deposit?