/// Callback to construct processes (uses the G4 particle table)
void Geant4PhysicsList::constructProcesses(G4VUserPhysicsList* physics_pointer) {
debug("constructProcesses %p", physics_pointer);
for (PhysicsProcesses::const_iterator i = m_discreteProcesses.begin(); i != m_discreteProcesses.end(); ++i) {
const string& part_name = (*i).first;
const ParticleProcesses& procs = (*i).second;
vector<G4ParticleDefinition*> defs(Geant4ParticleHandle::g4DefinitionsRegEx(part_name));
if (defs.empty()) {
except("Particle:%s Cannot find the corresponding entry in the particle table.", part_name.c_str());
}
for (vector<G4ParticleDefinition*>::const_iterator id = defs.begin(); id != defs.end(); ++id) {
G4ParticleDefinition* particle = *id;
G4ProcessManager* mgr = particle->GetProcessManager();
for (ParticleProcesses::const_iterator ip = procs.begin(); ip != procs.end(); ++ip) {
const Process& p = (*ip);
G4VProcess* g4 = PluginService::Create<G4VProcess*>(p.name);
if (!g4) { // Error no factory for this process
except("Particle:%s -> [%s] Cannot create discrete physics process %s",
part_name.c_str(), particle->GetParticleName().c_str(), p.name.c_str());
}
mgr->AddDiscreteProcess(g4);
info("Particle:%s -> [%s] added discrete process %s",
part_name.c_str(), particle->GetParticleName().c_str(), p.name.c_str());
}
}
}
for (PhysicsProcesses::const_iterator i = m_processes.begin(); i != m_processes.end(); ++i) {
const string& part_name = (*i).first;
const ParticleProcesses& procs = (*i).second;
vector<G4ParticleDefinition*> defs(Geant4ParticleHandle::g4DefinitionsRegEx(part_name));
if (defs.empty()) {
except("Particle:%s Cannot find the corresponding entry in the particle table.", part_name.c_str());
}
for (vector<G4ParticleDefinition*>::const_iterator id = defs.begin(); id != defs.end(); ++id) {
G4ParticleDefinition* particle = *id;
G4ProcessManager* mgr = particle->GetProcessManager();
for (ParticleProcesses::const_iterator ip = procs.begin(); ip != procs.end(); ++ip) {
const Process& p = (*ip);
G4VProcess* g4 = PluginService::Create<G4VProcess*>(p.name);
if (!g4) { // Error no factory for this process
except("Particle:%s -> [%s] Cannot create physics process %s",
part_name.c_str(), particle->GetParticleName().c_str(), p.name.c_str());
}
mgr->AddProcess(g4, p.ordAtRestDoIt, p.ordAlongSteptDoIt, p.ordPostStepDoIt);
info("Particle:%s -> [%s] added process %s with flags (%d,%d,%d)",
part_name.c_str(), particle->GetParticleName().c_str(), p.name.c_str(),
p.ordAtRestDoIt, p.ordAlongSteptDoIt, p.ordPostStepDoIt);
}
}
}
}
void PhysListEmStandardSS::ConstructProcess()
{
// Add standard EM Processes
aParticleIterator->reset();
while( (*aParticleIterator)() ){
G4ParticleDefinition* particle = aParticleIterator->value();
G4ProcessManager* pmanager = particle->GetProcessManager();
G4String particleName = particle->GetParticleName();
if (particleName == "gamma") {
// gamma
pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
pmanager->AddDiscreteProcess(new G4ComptonScattering);
pmanager->AddDiscreteProcess(new G4GammaConversion);
} else if (particleName == "e-") {
//electron
pmanager->AddProcess(new G4eIonisation, -1, 1, 1);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 2, 2);
pmanager->AddDiscreteProcess(new G4CoulombScattering);
} else if (particleName == "e+") {
//positron
pmanager->AddProcess(new G4eIonisation, -1, 1, 1);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 2, 2);
pmanager->AddProcess(new G4eplusAnnihilation, 0,-1, 3);
pmanager->AddDiscreteProcess(new G4CoulombScattering);
} else if (particleName == "mu+" ||
particleName == "mu-" ) {
//muon
pmanager->AddProcess(new G4MuIonisation, -1, 1, 1);
pmanager->AddProcess(new G4MuBremsstrahlung, -1, 2, 2);
pmanager->AddProcess(new G4MuPairProduction, -1, 3, 3);
pmanager->AddDiscreteProcess(new G4CoulombScattering);
} else if (particleName == "alpha" || particleName == "He3") {
pmanager->AddProcess(new G4ionIonisation, -1, 1, 1);
G4CoulombScattering* cs = new G4CoulombScattering();
//cs->AddEmModel(0, new G4IonCoulombScatteringModel());
cs->SetBuildTableFlag(false);
pmanager->AddDiscreteProcess(cs);
} else if (particleName == "GenericIon" ) {
pmanager->AddProcess(new G4ionIonisation, -1, 1, 1);
G4CoulombScattering* cs = new G4CoulombScattering();
cs->AddEmModel(0, new G4IonCoulombScatteringModel());
cs->SetBuildTableFlag(false);
pmanager->AddDiscreteProcess(cs);
} else if ((!particle->IsShortLived()) &&
(particle->GetPDGCharge() != 0.0) &&
(particle->GetParticleName() != "chargedgeantino")) {
//all others charged particles except geantino
pmanager->AddDiscreteProcess(new G4CoulombScattering);
pmanager->AddProcess(new G4hIonisation, -1, 1, 1);
}
}
// Em options
//
// Main options and setting parameters are shown here.
// Several of them have default values.
//
G4EmProcessOptions emOptions;
//physics tables
//
emOptions.SetMinEnergy(100*eV); //default
emOptions.SetMaxEnergy(100*TeV); //default
emOptions.SetDEDXBinning(12*20); //default=12*7
emOptions.SetLambdaBinning(12*20); //default=12*7
emOptions.SetSplineFlag(true); //default
//energy loss
//
emOptions.SetStepFunction(0.2, 100*um); //default=(0.2, 1*mm)
emOptions.SetLinearLossLimit(1.e-2); //default
//ionization
//
emOptions.SetSubCutoff(false); //default
// scattering
emOptions.SetPolarAngleLimit(0.0);
}
void Physics::ConstructEM() {
theParticleIterator->reset();
while( ( *theParticleIterator )() ) {
G4ParticleDefinition * particle = theParticleIterator->value();
G4ProcessManager * pmanager = particle->GetProcessManager();
G4String particleName = particle->GetParticleName();
if (particleName == "gamma") {
pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
pmanager->AddDiscreteProcess(new G4ComptonScattering);
pmanager->AddDiscreteProcess(new G4GammaConversion);
pmanager->AddDiscreteProcess(new G4CoulombScattering);
pmanager->AddDiscreteProcess(new G4RayleighScattering);
}
else if( particleName == "e-" ) {
if( Run_Type == "Recon" ) {
pmanager->AddProcess(new G4eIonisation, -1, 1, 1 );
}
else {
pmanager->AddProcess(new G4eMultipleScattering, -1, 1, 1 );
pmanager->AddProcess(new G4eIonisation, -1, 2, 2 );
if (fullBrem && supBrem)
{
std::cout<<"\n\nCANNOT USE FULL AND SUPPRESSED PHOTON BREMSSTRAHLUNG SIMULTANEOUSLY!\n\n";
std::cout<<"USING FULL BREMSSTRAHLUNG FOR ELECTRONS!\n\n";
supBrem = false;
}
else if (fullBrem)
{
pmanager->AddProcess(new G4eBremsstrahlung, -1, -1, 3 );
}
else if (supBrem)
{
G4eBremsstrahlung * eBremP = new G4eBremsstrahlung;
bremPhotonSuppression * bremSup = new bremPhotonSuppression;
bremSup->RegisterProcess(eBremP);
pmanager->AddProcess(bremSup, -1, -1, 3 );
}
else
{
std::cout<<"\n\nNO ELECTRON BREMSSTRAHLUNG INCLUDED!\n\n";
}
}
}
else if( particleName == "e+" ) {
if( Run_Type == "Recon" ) {
pmanager->AddProcess(new G4eIonisation, -1, 1, 1 );
}
else {
pmanager->AddProcess(new G4eMultipleScattering, -1, 1, 1 );
pmanager->AddProcess(new G4eIonisation, -1, 2, 2 );
if (fullBrem && supBrem)
{
std::cout<<"\n\nCANNOT USE FULL AND SUPPRESSED PHOTON BREMSSTRAHLUNG SIMULTANEOUSLY!\n\n";
std::cout<<"USING FULL BREMSSTRAHLUNG FOR POSITRONS!\n\n";
supBrem = false;
}
else if (fullBrem)
{
pmanager->AddProcess(new G4eBremsstrahlung, -1, -1, 3 );
}
else if (supBrem)
{
G4eBremsstrahlung * eBremP = new G4eBremsstrahlung;
bremPhotonSuppression * bremSup = new bremPhotonSuppression;
bremSup->RegisterProcess(eBremP);
pmanager->AddProcess(bremSup, -1, -1, 3 );
}
else
{
std::cout<<"\n\nNO POSITRON BREMSSTRAHLUNG INCLUDED!\n\n";
}
pmanager->AddProcess(new G4eplusAnnihilation, 0, -1, 4 );
}
}
else if( particleName == "mu+" ||
particleName == "mu-" ) {
pmanager->AddProcess(new G4MuMultipleScattering, -1, 1, 1 );
pmanager->AddProcess(new G4MuIonisation, -1, 2, 2 );
if (fullBrem && supBrem)
{
std::cout<<"\n\nCANNOT USE FULL AND SUPPRESSED PHOTON BREMSSTRAHLUNG SIMULTANEOUSLY!\n\n";
std::cout<<"USING FULL BREMSSTRAHLUNG FOR MUONS!\n\n";
supBrem = false;
}
else if (fullBrem)
{
pmanager->AddProcess(new G4MuBremsstrahlung, -1, -1, 3 );
}
else if (supBrem)
{
G4MuBremsstrahlung * eBremP = new G4MuBremsstrahlung;
bremPhotonSuppression * bremSup = new bremPhotonSuppression;
bremSup->RegisterProcess(eBremP);
pmanager->AddProcess(bremSup, -1, -1, 3 );
}
else
{
std::cout<<"\n\nNO MUON BREMSSTRAHLUNG INCLUDED!\n\n";
}
pmanager->AddProcess(new G4MuPairProduction, -1, -1, 4 );
}
// All others charged particles except chargedgeantino.
else if( !( particle->IsShortLived() ) &&
particle->GetPDGCharge() != 0.0 &&
particle->GetParticleName() != "chargedgeantino" ) {
if( Run_Type == "Recon" ) {
pmanager->AddProcess(new G4hIonisation, -1, 1, 1 );
}
else {
pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1 );
pmanager->AddProcess(new G4hIonisation, -1, 2, 2 );
if (fullBrem && supBrem)
{
std::cout<<"\n\nCANNOT USE FULL AND SUPPRESSED PHOTON BREMSSTRAHLUNG SIMULTANEOUSLY!\n\n";
std::cout<<"USING FULL BREMSSTRAHLUNG FOR HADRONS!\n\n";
supBrem = false;
}
else if (fullBrem)
{
pmanager->AddProcess(new G4hBremsstrahlung, -1, -1, 3 );
}
else if (supBrem)
{
G4hBremsstrahlung * eBremP = new G4hBremsstrahlung;
bremPhotonSuppression * bremSup = new bremPhotonSuppression;
bremSup->RegisterProcess(eBremP);
pmanager->AddProcess(bremSup, -1, -1, 3 );
}
else
{
std::cout<<"\n\nNO HADRON BREMSSTRAHLUNG INCLUDED!\n\n";
}
pmanager->AddProcess(new G4hPairProduction, -1, -1, 4 );
}
// Step limit.
//pmanager->AddProcess(new G4StepLimiter, -1, -1, 3 );
//pmanager->AddProcess(new G4UserSpecialCuts, -1, -1, 4 );
}
}
}
void PhysListEmStandardGS::ConstructProcess()
{
// Add standard EM Processes
//
aParticleIterator->reset();
while( (*aParticleIterator)() ){
G4ParticleDefinition* particle = aParticleIterator->value();
G4ProcessManager* pmanager = particle->GetProcessManager();
G4String particleName = particle->GetParticleName();
if (particleName == "gamma") {
// gamma
pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
pmanager->AddDiscreteProcess(new G4ComptonScattering);
pmanager->AddDiscreteProcess(new G4GammaConversion);
} else if (particleName == "e-") {
//electron
G4eMultipleScattering* msc = new G4eMultipleScattering();
msc->AddEmModel(0, new G4GoudsmitSaundersonMscModel());
pmanager->AddProcess(msc, -1, 1, 1);
pmanager->AddProcess(new G4eIonisation, -1, 2, 2);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 3, 3);
} else if (particleName == "e+") {
//positron
G4eMultipleScattering* msc = new G4eMultipleScattering();
msc->AddEmModel(0, new G4GoudsmitSaundersonMscModel());
pmanager->AddProcess(msc, -1, 1, 1);
pmanager->AddProcess(new G4eIonisation, -1, 2, 2);
pmanager->AddProcess(new G4eBremsstrahlung, -1, 3, 3);
pmanager->AddProcess(new G4eplusAnnihilation, 0,-1, 4);
} else if (particleName == "mu+" ||
particleName == "mu-" ) {
//muon
pmanager->AddProcess(new G4MuMultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4MuIonisation, -1, 2, 2);
pmanager->AddProcess(new G4MuBremsstrahlung, -1, 3, 3);
pmanager->AddProcess(new G4MuPairProduction, -1, 4, 4);
} else if( particleName == "proton" ||
particleName == "pi-" ||
particleName == "pi+" ) {
//proton
pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4hIonisation, -1, 2, 2);
pmanager->AddProcess(new G4hBremsstrahlung, -1, 3, 3);
pmanager->AddProcess(new G4hPairProduction, -1, 4, 4);
} else if( particleName == "alpha" ||
particleName == "He3" ) {
//alpha
pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4ionIonisation, -1, 2, 2);
pmanager->AddProcess(new G4NuclearStopping, -1, 3,-1);
} else if( particleName == "GenericIon" ) {
//Ions
pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
G4ionIonisation* ionIoni = new G4ionIonisation();
ionIoni->SetEmModel(new G4IonParametrisedLossModel());
pmanager->AddProcess(ionIoni, -1, 2, 2);
pmanager->AddProcess(new G4NuclearStopping, -1, 3,-1);
} else if ((!particle->IsShortLived()) &&
(particle->GetPDGCharge() != 0.0) &&
(particle->GetParticleName() != "chargedgeantino")) {
//all others charged particles except geantino
pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4hIonisation, -1, 2, 2);
}
}
// Em options
//
// Main options and setting parameters are shown here.
// Several of them have default values.
//
G4EmProcessOptions emOptions;
//physics tables
//
emOptions.SetMinEnergy(100*eV); //default
emOptions.SetMaxEnergy(100*TeV); //default
emOptions.SetDEDXBinning(12*20); //default=12*7
emOptions.SetLambdaBinning(12*20); //default=12*7
emOptions.SetSplineFlag(true); //default
//multiple coulomb scattering
//
emOptions.SetMscStepLimitation(fUseDistanceToBoundary); //default=fUseSafety
emOptions.SetMscRangeFactor(0.04); //default
emOptions.SetMscGeomFactor (2.5); //default
emOptions.SetSkin(3.); //default
//energy loss
//
emOptions.SetStepFunction(0.2, 100*um); //default=(0.2, 1*mm)
emOptions.SetLinearLossLimit(1.e-2); //default
//ionization
//
emOptions.SetSubCutoff(false); //default
}