Probably I found out the solution:
So according to these posts Big differences in muon generation between 11.0.4 and 11.1.1 and Range of low-energy muons changed between 11.0.p0 and 11.2.p01 I downgrade geant4 to version 11.0.4 and apply next command before Initialize.
The main thing is I forced GammaToMuonsFactor to 1e16!!!
UImanager->ApplyCommand(“/physics_lists/em/GammaToMuons true”);
UImanager->ApplyCommand(“/physics_lists/em/GammaToMuonsFactor 1e16”);
After this my ratio electrons to muons now ~0.25 as in CRY and agrees to literature. Energy and angular distribution of muons and electrons are also good!
In my physics list I use:
// EM Physics
RegisterPhysics(new G4EmStandardPhysics_option4());
// Extra
RegisterPhysics( new G4EmExtraPhysics(ver) );
// Decays
RegisterPhysics( new G4DecayPhysics(ver) );
// Hadron Elastic scattering
RegisterPhysics( new G4HadronElasticPhysicsHP(ver) );
// Hadron Physics
RegisterPhysics( new G4HadronPhysicsINCLXX(ver));
// Stopping Physics
RegisterPhysics( new G4StoppingPhysics(ver) );
// Ion Physics
RegisterPhysics( new G4IonINCLXXPhysics(ver));
But actually now it also works with general FTFP_BERT physics list, and result almost the same.
This method doesn’t work in the last geant4 version, but I will test one more time and other versions.
So any comments?