In Geant4.11.2 a new hadronic option lQMD was made available in the physics list shielding. It can be selected using an additional boolean parameter Shielding::Shielding(G4int verbose,const G4String&n_model, const G4String& HadrPhysVariant, G4bool useLightIonQMD) In Shielding.cc 175 // Ion Physics 176 RegisterPhysics( new G4IonElasticPhysics(verbose) ); 177 if (useLightIonQMD){ 178 RegisterPhysics( new G4LightIonQMDPhysics(verbose) ); 179 } else { 180 RegisterPhysics( new G4IonQMDPhysics(verbose) ); 181 } Unfortunately, in the physics list factory (G4PhysListFactory.cc), there is not yet a possibility to use it, requiring workarounds such as not using the recommended physics list factory. Perhaps one could introduce Shielding_LQMD or similar in the physics list factory? A suitable place could be next to the other Shielding options In G4PhysListFactory.cc) 182 else if(had_name == "Shielding") {p = new Shielding(verbose);} 183 else if(had_name == "ShieldingLEND") {p = new ShieldingLEND(verbose);} 184 else if(had_name == "ShieldingM") {p = new Shielding(verbose,"HP","M");}
Thanks for reporting this. For the coming beta-release, 11.3.beta, scheduled for June 28th, and for the next patch 11.2.p02 (not yet scheduled, but likely to appear sometimes this summer), there will be the physics list ShieldingLIQMD (and its variant, ShieldingLIQMD_HPT, with a special treatment of thermal neutrons) available from the physics list factory. As the name suggests, this physics list is similar to Shielding, except that it uses the new G4LightIonQMDReaction model for ions with kinetic energies between 30 MeV/u and 500 MeV/u.
