Physics Lists to produce Secondary Particles

[Adwaya_Gupta]

Geant4 Version: geant4-11-02-patch-02
Operating System: MacOS
Compiler/Version: Apple clang version 15.0.0
CMake Version: 3.31.0

---

I currently have a simple simulation of a proton colliding with a cube of Technetium.
I’m interested in the secondary particles produced from this, however currently only electrons and positrons are being produced.
How do I need to modify my current physics list (shown below) to allow the production of other secondary particles (pions, kaons etc.)

#include "PhysicsList.hh"

PhysicsList::PhysicsList() {
    RegisterPhysics( new G4EmStandardPhysics() );
    RegisterPhysics( new G4OpticalPhysics() ); 
}

PhysicsList::~PhysicsList() {}

[maire1]

see examples/extended/hahronic/Hadr06/src/PhysicsList.cc
(or Hadr05 or Hadr07 etc …)

[Adwaya_Gupta]

Even after including the PhysicsLists used in the specified file (my updated PhysicsList is shown below), its still only electrons/positrons being produced.
Is there anywhere I can find further information about each of the physics lists?

#include "PhysicsList.hh"

PhysicsList::PhysicsList() {
    G4int verb = 1;
    SetVerboseLevel(verb);

    const G4double meanLife = 1*nanosecond, halfLife = meanLife*std::log(2);
    G4NuclideTable::GetInstance()->SetThresholdOfHalfLife(halfLife);

    RegisterPhysics( new G4HadronElasticPhysicsXS(verb) );
    RegisterPhysics( new G4HadronElasticPhysics(verb) );
    RegisterPhysics( new G4IonElasticPhysics(verb) );
    RegisterPhysics( new G4IonPhysicsXS(verb) );
    RegisterPhysics( new G4StoppingPhysics(verb) );
    RegisterPhysics( new G4EmStandardPhysics_option3() );
    RegisterPhysics( new G4DecayPhysics() );
    RegisterPhysics( new G4RadioactiveDecayPhysics() ); 
    RegisterPhysics( new G4OpticalPhysics() );  
}

PhysicsList::~PhysicsList() {}

[mkelsey]

Could you try just using one of the Geant4 reference physics lists (FTFP_BERT or Shielding are my two go-to’s). While the custom physics list you posted above ought to do everything you want, starting with a reference list and then comparing that to yours might be easier.

[jrellin]

How are you generating the protons? The default energy using gps, for example, is 1 MeV. Which is unlikely to overcome any possible coloumb barrier.