// // ******************************************************************** // * License and Disclaimer * // * * // * The Geant4 software is copyright of the Copyright Holders of * // * the Geant4 Collaboration. It is provided under the terms and * // * conditions of the Geant4 Software License, included in the file * // * LICENSE and available at http://cern.ch/geant4/license . These * // * include a list of copyright holders. * // * * // * Neither the authors of this software system, nor their employing * // * institutes,nor the agencies providing financial support for this * // * work make any representation or warranty, express or implied, * // * regarding this software system or assume any liability for its * // * use. Please see the license in the file LICENSE and URL above * // * for the full disclaimer and the limitation of liability. * // * * // * This code implementation is the result of the scientific and * // * technical work of the GEANT4 collaboration. * // * By using, copying, modifying or distributing the software (or * // * any work based on the software) you agree to acknowledge its * // * use in resulting scientific publications, and indicate your * // * acceptance of all terms of the Geant4 Software license. * // ******************************************************************** // // // $Id: PlhPhysicsList.cc,v 1.2 2008/12/18 13:02:26 gunter Exp $ // GEANT4 tag $Name: geant4-09-03-cand-01 $ // #include "globals.hh" #include "PlhPhysicsList.hh" #include "G4ParticleDefinition.hh" #include "G4ParticleWithCuts.hh" #include "G4ProcessManager.hh" #include "G4ProcessVector.hh" #include "G4ParticleTypes.hh" #include "G4ParticleTable.hh" #include "G4Material.hh" #include "G4MaterialTable.hh" #include "G4ios.hh" #include PlhPhysicsList::PlhPhysicsList(): G4VUserPhysicsList() { SetVerboseLevel(1); } PlhPhysicsList::~PlhPhysicsList() { } void PlhPhysicsList::ConstructParticle() { // In this method, static member functions should be called // for all particles which you want to use. // This ensures that objects of these particle types will be // created in the program. ConstructBosons(); ConstructLeptons(); ConstructMesons(); ConstructBarions(); ConstructIons(); } void PlhPhysicsList::ConstructBosons() { // pseudo-particles G4Geantino::GeantinoDefinition(); G4ChargedGeantino::ChargedGeantinoDefinition(); // gamma G4Gamma::GammaDefinition(); // optical photon G4OpticalPhoton::OpticalPhotonDefinition(); } void PlhPhysicsList::ConstructLeptons() { // leptons G4MuonPlus::MuonPlusDefinition(); G4MuonMinus::MuonMinusDefinition(); G4TauMinus::TauMinusDefinition(); G4TauPlus::TauPlusDefinition(); G4Electron::ElectronDefinition(); G4Positron::PositronDefinition(); G4NeutrinoTau::NeutrinoTauDefinition(); G4AntiNeutrinoTau::AntiNeutrinoTauDefinition(); G4NeutrinoMu::NeutrinoMuDefinition(); G4AntiNeutrinoMu::AntiNeutrinoMuDefinition(); G4NeutrinoE::NeutrinoEDefinition(); G4AntiNeutrinoE::AntiNeutrinoEDefinition(); } void PlhPhysicsList::ConstructMesons() { // mesons G4PionPlus::PionPlusDefinition(); G4PionMinus::PionMinusDefinition(); G4PionZero::PionZeroDefinition(); G4Eta::EtaDefinition(); G4EtaPrime::EtaPrimeDefinition(); G4RhoZero::RhoZeroDefinition(); G4KaonPlus::KaonPlusDefinition(); G4KaonMinus::KaonMinusDefinition(); G4KaonZero::KaonZeroDefinition(); G4AntiKaonZero::AntiKaonZeroDefinition(); G4KaonZeroLong::KaonZeroLongDefinition(); G4KaonZeroShort::KaonZeroShortDefinition(); G4DMesonPlus::DMesonPlusDefinition(); G4DMesonMinus::DMesonMinusDefinition(); G4DMesonZero::DMesonZeroDefinition(); G4AntiDMesonZero::AntiDMesonZeroDefinition(); G4DsMesonPlus::DsMesonPlusDefinition(); G4DsMesonMinus::DsMesonMinusDefinition(); G4JPsi::JPsiDefinition(); G4BMesonPlus::BMesonPlusDefinition(); G4BMesonMinus::BMesonMinusDefinition(); G4BMesonZero::BMesonZeroDefinition(); G4AntiBMesonZero::AntiBMesonZeroDefinition(); G4BsMesonZero::BsMesonZeroDefinition(); G4AntiBsMesonZero::AntiBsMesonZeroDefinition(); } void PlhPhysicsList::ConstructBarions() { // barions G4Proton::ProtonDefinition(); G4AntiProton::AntiProtonDefinition(); G4Neutron::NeutronDefinition(); G4AntiNeutron::AntiNeutronDefinition(); G4Lambda::LambdaDefinition(); G4AntiLambda::AntiLambdaDefinition(); G4SigmaZero::SigmaZeroDefinition(); G4AntiSigmaZero::AntiSigmaZeroDefinition(); G4SigmaPlus::SigmaPlusDefinition(); G4AntiSigmaPlus::AntiSigmaPlusDefinition(); G4SigmaMinus::SigmaMinusDefinition(); G4AntiSigmaMinus::AntiSigmaMinusDefinition(); G4XiZero::XiZeroDefinition(); G4AntiXiZero::AntiXiZeroDefinition(); G4XiMinus::XiMinusDefinition(); G4AntiXiMinus::AntiXiMinusDefinition(); G4OmegaMinus::OmegaMinusDefinition(); G4AntiOmegaMinus::AntiOmegaMinusDefinition(); G4LambdacPlus::LambdacPlusDefinition(); G4SigmacPlusPlus::SigmacPlusPlusDefinition(); G4SigmacPlus::SigmacPlusDefinition(); G4SigmacZero::SigmacZeroDefinition(); G4XicPlus::XicPlusDefinition(); G4XicZero::XicZeroDefinition(); G4OmegacZero::OmegacZeroDefinition(); G4AntiLambdacPlus::AntiLambdacPlusDefinition(); G4AntiSigmacPlusPlus::AntiSigmacPlusPlusDefinition(); G4AntiSigmacPlus::AntiSigmacPlusDefinition(); G4AntiSigmacZero::AntiSigmacZeroDefinition(); G4AntiXicPlus::AntiXicPlusDefinition(); G4AntiXicZero::AntiXicZeroDefinition(); G4AntiOmegacZero::AntiOmegacZeroDefinition(); } void PlhPhysicsList::ConstructIons() { // nuclei G4Alpha::AlphaDefinition(); G4Deuteron::DeuteronDefinition(); G4Triton::TritonDefinition(); G4He3::He3Definition(); } void PlhPhysicsList::ConstructProcess() { AddTransportation(); ConstructEM(); ConstructLeptHad(); // ConstructHad(); // ConstructGeneral(); } #include "G4ComptonScattering.hh" #include "G4GammaConversion.hh" #include "G4PhotoElectricEffect.hh" #include "G4MultipleScattering.hh" #include "G4eIonisation.hh" #include "G4eBremsstrahlung.hh" #include "G4eplusAnnihilation.hh" #include "G4MuIonisation.hh" #include "G4MuBremsstrahlung.hh" #include "G4MuPairProduction.hh" #include "G4hIonisation.hh" void PlhPhysicsList::ConstructEM() { theParticleIterator->reset(); while( (*theParticleIterator)() ){ G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if (particleName == "gamma") { // gamma // Construct processes for gamma pmanager->AddDiscreteProcess(new G4GammaConversion()); pmanager->AddDiscreteProcess(new G4ComptonScattering()); pmanager->AddDiscreteProcess(new G4PhotoElectricEffect()); } else if (particleName == "e-") { //electron // Construct processes for electron pmanager->AddProcess(new G4MultipleScattering(),-1,1,1); pmanager->AddProcess(new G4eIonisation(),-1,2,2); pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3); } else if (particleName == "e+") { //positron // Construct processes for positron pmanager->AddProcess(new G4MultipleScattering(),-1,1,1); pmanager->AddProcess(new G4eIonisation(),-1,2,2); pmanager->AddProcess(new G4eBremsstrahlung(),-1,-1,3); pmanager->AddProcess(new G4eplusAnnihilation(),0,-1,4); } else if( particleName == "mu+" || particleName == "mu-" ) { //muon // Construct processes for muon+ pmanager->AddProcess(new G4MultipleScattering(),-1,1,1); pmanager->AddProcess(new G4MuIonisation(),-1,2,2); pmanager->AddProcess(new G4MuBremsstrahlung(),-1,-1,3); pmanager->AddProcess(new G4MuPairProduction(),-1,-1,4); } else { if ((particle->GetPDGCharge() != 0.0) && (particle->GetParticleName() != "chargedgeantino")) { // all others charged particles except geantino pmanager->AddProcess(new G4MultipleScattering(),-1,1,1); pmanager->AddProcess(new G4hIonisation(),-1,2,2); } } } } // Hadron Processes #include "G4HadronElasticProcess.hh" #include "G4HadronFissionProcess.hh" #include "G4HadronCaptureProcess.hh" #include "G4PionPlusInelasticProcess.hh" #include "G4PionMinusInelasticProcess.hh" #include "G4KaonPlusInelasticProcess.hh" #include "G4KaonZeroSInelasticProcess.hh" #include "G4KaonZeroLInelasticProcess.hh" #include "G4KaonMinusInelasticProcess.hh" #include "G4ProtonInelasticProcess.hh" #include "G4AntiProtonInelasticProcess.hh" #include "G4NeutronInelasticProcess.hh" #include "G4AntiNeutronInelasticProcess.hh" #include "G4LambdaInelasticProcess.hh" #include "G4AntiLambdaInelasticProcess.hh" #include "G4SigmaPlusInelasticProcess.hh" #include "G4SigmaMinusInelasticProcess.hh" #include "G4AntiSigmaPlusInelasticProcess.hh" #include "G4AntiSigmaMinusInelasticProcess.hh" #include "G4XiZeroInelasticProcess.hh" #include "G4XiMinusInelasticProcess.hh" #include "G4AntiXiZeroInelasticProcess.hh" #include "G4AntiXiMinusInelasticProcess.hh" #include "G4DeuteronInelasticProcess.hh" #include "G4TritonInelasticProcess.hh" #include "G4AlphaInelasticProcess.hh" #include "G4OmegaMinusInelasticProcess.hh" #include "G4AntiOmegaMinusInelasticProcess.hh" // Low-energy Models #include "G4LElastic.hh" #include "G4LFission.hh" #include "G4LCapture.hh" #include "G4LEPionPlusInelastic.hh" #include "G4LEPionMinusInelastic.hh" #include "G4LEKaonPlusInelastic.hh" #include "G4LEKaonZeroSInelastic.hh" #include "G4LEKaonZeroLInelastic.hh" #include "G4LEKaonMinusInelastic.hh" #include "G4LEProtonInelastic.hh" #include "G4LEAntiProtonInelastic.hh" #include "G4LENeutronInelastic.hh" #include "G4LEAntiNeutronInelastic.hh" #include "G4LELambdaInelastic.hh" #include "G4LEAntiLambdaInelastic.hh" #include "G4LESigmaPlusInelastic.hh" #include "G4LESigmaMinusInelastic.hh" #include "G4LEAntiSigmaPlusInelastic.hh" #include "G4LEAntiSigmaMinusInelastic.hh" #include "G4LEXiZeroInelastic.hh" #include "G4LEXiMinusInelastic.hh" #include "G4LEAntiXiZeroInelastic.hh" #include "G4LEAntiXiMinusInelastic.hh" #include "G4LEDeuteronInelastic.hh" #include "G4LETritonInelastic.hh" #include "G4LEAlphaInelastic.hh" #include "G4LEOmegaMinusInelastic.hh" #include "G4LEAntiOmegaMinusInelastic.hh" // // ConstructHad() // // Makes discrete physics processes for the hadrons, at present limited // to those particles with GHEISHA interactions (INTRC > 0). // The processes are: Elastic scattering, Inelastic scattering, // Fission (for neutron only), and Capture (neutron). // // F.W.Jones 06-JUL-1998 // void PlhPhysicsList::ConstructHad() { G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess; G4LElastic* theElasticModel = new G4LElastic; theElasticProcess->RegisterMe(theElasticModel); theParticleIterator->reset(); while ((*theParticleIterator)()) { G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if (particleName == "pi+") { pmanager->AddDiscreteProcess(theElasticProcess); G4PionPlusInelasticProcess* theInelasticProcess = new G4PionPlusInelasticProcess("inelastic"); G4LEPionPlusInelastic* theInelasticModel = new G4LEPionPlusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "pi-") { pmanager->AddDiscreteProcess(theElasticProcess); G4PionMinusInelasticProcess* theInelasticProcess = new G4PionMinusInelasticProcess("inelastic"); G4LEPionMinusInelastic* theInelasticModel = new G4LEPionMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "kaon+") { pmanager->AddDiscreteProcess(theElasticProcess); G4KaonPlusInelasticProcess* theInelasticProcess = new G4KaonPlusInelasticProcess("inelastic"); G4LEKaonPlusInelastic* theInelasticModel = new G4LEKaonPlusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "kaon0S") { pmanager->AddDiscreteProcess(theElasticProcess); G4KaonZeroSInelasticProcess* theInelasticProcess = new G4KaonZeroSInelasticProcess("inelastic"); G4LEKaonZeroSInelastic* theInelasticModel = new G4LEKaonZeroSInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "kaon0L") { pmanager->AddDiscreteProcess(theElasticProcess); G4KaonZeroLInelasticProcess* theInelasticProcess = new G4KaonZeroLInelasticProcess("inelastic"); G4LEKaonZeroLInelastic* theInelasticModel = new G4LEKaonZeroLInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "kaon-") { pmanager->AddDiscreteProcess(theElasticProcess); G4KaonMinusInelasticProcess* theInelasticProcess = new G4KaonMinusInelasticProcess("inelastic"); G4LEKaonMinusInelastic* theInelasticModel = new G4LEKaonMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "proton") { pmanager->AddDiscreteProcess(theElasticProcess); G4ProtonInelasticProcess* theInelasticProcess = new G4ProtonInelasticProcess("inelastic"); G4LEProtonInelastic* theInelasticModel = new G4LEProtonInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_proton") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiProtonInelasticProcess* theInelasticProcess = new G4AntiProtonInelasticProcess("inelastic"); G4LEAntiProtonInelastic* theInelasticModel = new G4LEAntiProtonInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "neutron") { pmanager->AddDiscreteProcess(theElasticProcess); G4NeutronInelasticProcess* theInelasticProcess = new G4NeutronInelasticProcess("inelastic"); G4LENeutronInelastic* theInelasticModel = new G4LENeutronInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); G4HadronFissionProcess* theFissionProcess = new G4HadronFissionProcess; G4LFission* theFissionModel = new G4LFission; theFissionProcess->RegisterMe(theFissionModel); pmanager->AddDiscreteProcess(theFissionProcess); G4HadronCaptureProcess* theCaptureProcess = new G4HadronCaptureProcess; G4LCapture* theCaptureModel = new G4LCapture; theCaptureProcess->RegisterMe(theCaptureModel); pmanager->AddDiscreteProcess(theCaptureProcess); } else if (particleName == "anti_neutron") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiNeutronInelasticProcess* theInelasticProcess = new G4AntiNeutronInelasticProcess("inelastic"); G4LEAntiNeutronInelastic* theInelasticModel = new G4LEAntiNeutronInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "lambda") { pmanager->AddDiscreteProcess(theElasticProcess); G4LambdaInelasticProcess* theInelasticProcess = new G4LambdaInelasticProcess("inelastic"); G4LELambdaInelastic* theInelasticModel = new G4LELambdaInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_lambda") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiLambdaInelasticProcess* theInelasticProcess = new G4AntiLambdaInelasticProcess("inelastic"); G4LEAntiLambdaInelastic* theInelasticModel = new G4LEAntiLambdaInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "sigma+") { pmanager->AddDiscreteProcess(theElasticProcess); G4SigmaPlusInelasticProcess* theInelasticProcess = new G4SigmaPlusInelasticProcess("inelastic"); G4LESigmaPlusInelastic* theInelasticModel = new G4LESigmaPlusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "sigma-") { pmanager->AddDiscreteProcess(theElasticProcess); G4SigmaMinusInelasticProcess* theInelasticProcess = new G4SigmaMinusInelasticProcess("inelastic"); G4LESigmaMinusInelastic* theInelasticModel = new G4LESigmaMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_sigma+") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiSigmaPlusInelasticProcess* theInelasticProcess = new G4AntiSigmaPlusInelasticProcess("inelastic"); G4LEAntiSigmaPlusInelastic* theInelasticModel = new G4LEAntiSigmaPlusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_sigma-") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiSigmaMinusInelasticProcess* theInelasticProcess = new G4AntiSigmaMinusInelasticProcess("inelastic"); G4LEAntiSigmaMinusInelastic* theInelasticModel = new G4LEAntiSigmaMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "xi0") { pmanager->AddDiscreteProcess(theElasticProcess); G4XiZeroInelasticProcess* theInelasticProcess = new G4XiZeroInelasticProcess("inelastic"); G4LEXiZeroInelastic* theInelasticModel = new G4LEXiZeroInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "xi-") { pmanager->AddDiscreteProcess(theElasticProcess); G4XiMinusInelasticProcess* theInelasticProcess = new G4XiMinusInelasticProcess("inelastic"); G4LEXiMinusInelastic* theInelasticModel = new G4LEXiMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_xi0") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiXiZeroInelasticProcess* theInelasticProcess = new G4AntiXiZeroInelasticProcess("inelastic"); G4LEAntiXiZeroInelastic* theInelasticModel = new G4LEAntiXiZeroInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_xi-") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiXiMinusInelasticProcess* theInelasticProcess = new G4AntiXiMinusInelasticProcess("inelastic"); G4LEAntiXiMinusInelastic* theInelasticModel = new G4LEAntiXiMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "deuteron") { pmanager->AddDiscreteProcess(theElasticProcess); G4DeuteronInelasticProcess* theInelasticProcess = new G4DeuteronInelasticProcess("inelastic"); G4LEDeuteronInelastic* theInelasticModel = new G4LEDeuteronInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "triton") { pmanager->AddDiscreteProcess(theElasticProcess); G4TritonInelasticProcess* theInelasticProcess = new G4TritonInelasticProcess("inelastic"); G4LETritonInelastic* theInelasticModel = new G4LETritonInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "alpha") { pmanager->AddDiscreteProcess(theElasticProcess); G4AlphaInelasticProcess* theInelasticProcess = new G4AlphaInelasticProcess("inelastic"); G4LEAlphaInelastic* theInelasticModel = new G4LEAlphaInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "omega-") { pmanager->AddDiscreteProcess(theElasticProcess); G4OmegaMinusInelasticProcess* theInelasticProcess = new G4OmegaMinusInelasticProcess("inelastic"); G4LEOmegaMinusInelastic* theInelasticModel = new G4LEOmegaMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } else if (particleName == "anti_omega-") { pmanager->AddDiscreteProcess(theElasticProcess); G4AntiOmegaMinusInelasticProcess* theInelasticProcess = new G4AntiOmegaMinusInelasticProcess("inelastic"); G4LEAntiOmegaMinusInelastic* theInelasticModel = new G4LEAntiOmegaMinusInelastic; theInelasticProcess->RegisterMe(theInelasticModel); pmanager->AddDiscreteProcess(theInelasticProcess); } } } #include "G4MuonNucleusProcess.hh" void PlhPhysicsList::ConstructLeptHad() { G4MuonNucleusProcess* theMuonNucleusProcess = new G4MuonNucleusProcess(); theParticleIterator->reset(); while( (*theParticleIterator)() ){ G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); G4String particleName = particle->GetParticleName(); if( particleName == "mu+" || particleName == "mu-" ) { pmanager->AddDiscreteProcess(theMuonNucleusProcess); } } } #include "G4Decay.hh" void PlhPhysicsList::ConstructGeneral() { G4Decay* theDecayProcess = new G4Decay(); theParticleIterator->reset(); while( (*theParticleIterator)() ){ G4ParticleDefinition* particle = theParticleIterator->value(); G4ProcessManager* pmanager = particle->GetProcessManager(); if (theDecayProcess->IsApplicable(*particle)) { pmanager ->AddProcess(theDecayProcess); pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep); pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest); } } } void PlhPhysicsList::SetCuts(G4double cut) { if (verboseLevel >0){ G4cout << "PlhPhysicsList::SetCuts:"; G4cout << "CutLength : " << cut/mm << " (mm)" << G4endl; } // set cut values for gamma at first and for e- second and next for e+, // because some processes for e+/e- need cut values for gamma SetCutValue(cut, "gamma"); SetCutValue(cut, "e-"); SetCutValue(cut, "e+"); // set cut values for proton and anti_proton before all other hadrons // because some processes for hadrons need cut values for proton/anti_proton SetCutValue(cut, "proton"); SetCutValue(cut, "anti_proton"); SetCutValueForOthers(cut); if (verboseLevel>1) { DumpCutValuesTable(); } }