What's the substitute of G4hLowEnergyIonisation in Geant4 10.4

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?

Thanks in advance.

Hello,

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

More details on PhysicsLIsts : http://geant4-userdoc.web.cern.ch/geant4-userdoc/UsersGuides/PhysicsListGuide/html/index.html
and for EM physics: http://geant4.web.cern.ch/content/physics-lists/em-physics-lists

VI

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?

Thanks a lot.

see :
http://geant4-userdoc.web.cern.ch/geant4-userdoc/UsersGuides/PhysicsReferenceManual/html/electromagnetic/energy_loss/enloss.html

http://lappweb.in2p3.fr/~maire/tutorials/energyRange.pdf

You can take any example TestEm1, TestEm5, TestEm11, TestEm18, insert in you macro the command
/process/eLoss/StepFunction
and play with it