// Here we declare what the objects of our class 
// MyPhysicsList actually do. In the header file .hh we 
// only declared it.

#include "physics.hh"

// We define the constructor of the MyPhysicsList class
// with the same name.

MyPhysicsList::MyPhysicsList()
{
// By calling RegisterPhysics with these modules we tell Geant to include
// these sets of physical processes in our simulation.  
  //  RegisterPhysics (new G4EmStandardPhysics()); // lowest electron energies 1keV
  //  RegisterPhysics (new G4OpticalPhysics());
  //  RegisterPhysics (new G4EmLivermorePhysics()); // electrons and photons [250 eV - 1 GeV]
  //  RegisterPhysics (new G4EmPenelopePhysics());  // electrons, positrons and photons [100 eV - 1 GeV]
    RegisterPhysics(new G4EmStandardPhysics_option4()); // best current advanced EM options standard + lowenergy
  //  RegisterPhysics(new G4EmStandardPhysicsSS());       // standard EM physics and single scattering model
  //  RegisterPhysics(new G4EmLowEPPhysics());            // low-energy EM physics implementing experimental low-energy models
  //  RegisterPhysics(new G4eIonisation());               //
  //  RegisterPhysics(new G4PenelopeIonisationModel());   //
}

MyPhysicsList::~MyPhysicsList()
{}