Overview
LCPhys is a Geant4 physics list written by Dennis Wright for simulation of events in ILC detectors.
LCPhys HomepageLCPhys Source Code Documentation (doxygen)
...
This class creates the mesons using G4MesonConstructor, baryons with G4BaryonConstructor, and resonances and quarks with G4ShortLivedConstructor.
Hadron Models
- hadronic elastic process and model - same for all hadrons
- G4PionPlus, G4PionMinus
- G4MultipleScattering
- G4hIonisation
- G4PionPlus
- G4PionPlusInelasticProcess
- G4PionMinus
- G4KaonPlus
- G4KaonMinus
- G4KaonZeroLong
- G4KaonZeroShort
- G4Proton
- G4AntiProton
- G4Neutron
- G4AntiNeutron
- G4Lambda
- G4AntiLambda
- G4SigmaMinus
- G4AntiSigmaMinus
- G4SigmaPlus
- G4AntiSigmaPlus
- G4XiMinus
- G4AntiXiMinus
- G4XiZero
- G4AntiXiZero
- G4OmegaMinus
- G4AntiOmegaMinus
No Format |
---|
// Hadronic Elastic Process and Model (the same for all hadrons) G4HadronElasticProcess* elasticProcess = new G4HadronElasticProcess(); G4LElastic* elasticModel = new G4LElastic(); elasticProcess->RegisterMe(elasticModel); // Hadronic inelastic models // Bertini cascade model: use between 0 and 9.9 GeV // for |
- Bertini cascade model for p,n,pi+,pi- between 0 and 9.9 GeV
No Format |
---|
G4CascadeInterface* bertiniModel = new G4CascadeInterface(); bertiniModel->SetMaxEnergy(9.9*GeV); // Bertini cascade model: use between 0 and 13 GeV for // |
- Bertini cascade model for K+,K-,K0L,K0S,Lambda,Sigma+,Sigma-,Xi0,Xi- between 0 and 13 GeV
No Format |
---|
G4CascadeInterface* bertiniModelStrange = new G4CascadeInterface(); bertiniModelStrange->SetMaxEnergy(13*GeV); // Low energy parameterized models : use between |
- Low energy parameterized models between 9.5
...
- and
...
- 25
...
- GeV
- The LEP models are defined for each type of particle. (see below)
- Quark-Gluon String Model (QGSP) for p,n,pi+,pi-,K+,K-,K0L,K0S
...
- between
...
- 12
...
- GeV
...
- and
...
- 100 TeV
- G4TheoFSGenerator
- G4GeneratorPrecompoundInterface - transport
- G4PreCompoundModel - deexcitation
- G4QGSModel - high energy generator
- G4ExcitedStringDecay - fragmentation model
- G4QGSMFragmentation - Quark-Gluon String model fragmentation
- G4ExcitedStringDecay - fragmentation model
- G4GeneratorPrecompoundInterface - transport
- G4TheoFSGenerator
Defining the QGSP process and model in LCHadronPhysics.cc.
No Format |
---|
TeV G4TheoFSGenerator* QGSPModel = new G4TheoFSGenerator(); G4GeneratorPrecompoundInterface* theCascade = new G4GeneratorPrecompoundInterface(); G4ExcitationHandler* exHandler = new G4ExcitationHandler(); G4PreCompoundModel* preCompound = new G4PreCompoundModel(exHandler); theCascade->SetDeExcitation(preCompound); QGSPModel->SetTransport(theCascade); G4QGSMFragmentation* frag = new G4QGSMFragmentation(); G4ExcitedStringDecay* stringDecay = new G4ExcitedStringDecay(frag); G4QGSModel<G4QGSParticipants>* stringModel = new G4QGSModel<G4QGSParticipants>(); stringModel->SetFragmentationModel(stringDecay); QGSPModel->SetHighEnergyGenerator(stringModel); QGSPModel->SetMinEnergy(12*GeV); QGSPModel->SetMaxEnergy(100*TeV); // G4ProcessManager * pManager = 0; /////////////////// // // // pi+ physics // // // /////////////////// pManager = G4PionPlus::PionPlus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4PionPlusInelasticProcess* pipinelProc = new G4PionPlusInelasticProcess(); G4PiNuclearCrossSection* pion_XC = new G4PiNuclearCrossSection(); pipinelProc->AddDataSet(pion_XC); pipinelProc->RegisterMe(bertiniModel); G4LEPionPlusInelastic* LEPpipModel = new G4LEPionPlusInelastic(); LEPpipModel->SetMinEnergy(LEPpnpiLimit); LEPpipModel->SetMaxEnergy(LEPUpperLimit); pipinelProc->RegisterMe(LEPpipModel); pipinelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(pipinelProc); /////////////////// // // // pi- physics // // // /////////////////// pManager = G4PionMinus::PionMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4PionMinusInelasticProcess* piminelProc = new G4PionMinusInelasticProcess(); piminelProc->AddDataSet(pion_XC); piminelProc->RegisterMe(bertiniModel); G4LEPionMinusInelastic* LEPpimModel = new G4LEPionMinusInelastic(); LEPpimModel->SetMinEnergy(LEPpnpiLimit); LEPpimModel->SetMaxEnergy(LEPUpperLimit); piminelProc->RegisterMe(LEPpimModel); piminelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(piminelProc); // pi- absorption at rest G4PionMinusAbsorptionAtRest* pimAbsorb = new G4PionMinusAbsorptionAtRest(); pManager->AddRestProcess(pimAbsorb); /////////////////// // // // K+ physics // // // /////////////////// pManager = G4KaonPlus::KaonPlus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4KaonPlusInelasticProcess* kpinelProc = new G4KaonPlusInelasticProcess(); kpinelProc->RegisterMe(bertiniModelStrange); kpinelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(kpinelProc); /////////////////// // // // K- physics // // // /////////////////// pManager = G4KaonMinus::KaonMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4KaonMinusInelasticProcess* kminelProc = new G4KaonMinusInelasticProcess(); kminelProc->RegisterMe(bertiniModelStrange); kminelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(kminelProc); // K- absorption at rest G4KaonMinusAbsorption* kmAbsorb = new G4KaonMinusAbsorption(); pManager->AddRestProcess(kmAbsorb); /////////////////// // // // K0L physics // // // /////////////////// pManager = G4KaonZeroLong::KaonZeroLong()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4KaonZeroLInelasticProcess* k0LinelProc = new G4KaonZeroLInelasticProcess(); k0LinelProc->RegisterMe(bertiniModelStrange); k0LinelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(k0LinelProc); /////////////////// // // // K0S physics // // // /////////////////// pManager = G4KaonZeroShort::KaonZeroShort()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4KaonZeroSInelasticProcess* k0SinelProc = new G4KaonZeroSInelasticProcess(); k0SinelProc->RegisterMe(bertiniModelStrange); k0SinelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(k0SinelProc); /////////////////// // // // Proton // // // /////////////////// pManager = G4Proton::Proton()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4ProtonInelasticProcess* pinelProc = new G4ProtonInelasticProcess(); G4ProtonInelasticCrossSection* proton_XC = new G4ProtonInelasticCrossSection(); pinelProc->AddDataSet(proton_XC); pinelProc->RegisterMe(bertiniModel); G4LEProtonInelastic* LEPpModel = new G4LEProtonInelastic(); LEPpModel->SetMinEnergy(LEPpnpiLimit); LEPpModel->SetMaxEnergy(LEPUpperLimit); pinelProc->RegisterMe(LEPpModel); pinelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(pinelProc); /////////////////// // // // Anti-Proton // // // /////////////////// pManager = G4AntiProton::AntiProton()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiProtonInelasticProcess* apinelProc = new G4AntiProtonInelasticProcess(); G4LEAntiProtonInelastic* LEPapModel = new G4LEAntiProtonInelastic(); apinelProc->RegisterMe(LEPapModel); G4HEAntiProtonInelastic* HEPapModel = new G4HEAntiProtonInelastic(); apinelProc->RegisterMe(HEPapModel); pManager->AddDiscreteProcess(apinelProc); // anti-proton annihilation at rest G4AntiProtonAnnihilationAtRest* apAnnihil = new G4AntiProtonAnnihilationAtRest(); pManager->AddRestProcess(apAnnihil); /////////////////// // // // Neutron // // // /////////////////// pManager = G4Neutron::Neutron()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4NeutronInelasticProcess* ninelProc = new G4NeutronInelasticProcess(); G4NeutronInelasticCrossSection* neutron_XC = new G4NeutronInelasticCrossSection(); ninelProc->AddDataSet(neutron_XC); ninelProc->RegisterMe(bertiniModel); G4LENeutronInelastic* LEPnModel = new G4LENeutronInelastic(); LEPnModel->SetMinEnergy(LEPpnpiLimit); LEPnModel->SetMaxEnergy(LEPUpperLimit); ninelProc->RegisterMe(LEPnModel); ninelProc->RegisterMe(QGSPModel); pManager->AddDiscreteProcess(ninelProc); // neutron-induced fission G4HadronFissionProcess* neutronFission = new G4HadronFissionProcess(); G4LFission* neutronFissionModel = new G4LFission(); neutronFissionModel->SetMinEnergy(0.); neutronFissionModel->SetMaxEnergy(20*TeV); neutronFission->RegisterMe(neutronFissionModel); pManager->AddDiscreteProcess(neutronFission); // neutron capture G4HadronCaptureProcess* neutronCapture = new G4HadronCaptureProcess(); G4LCapture* neutronCaptureModel = new G4LCapture(); neutronCaptureModel->SetMinEnergy(0.); neutronCaptureModel->SetMaxEnergy(20*TeV); neutronCapture->RegisterMe(neutronCaptureModel); pManager->AddDiscreteProcess(neutronCapture); /////////////////// // // // Anti-Neutron // // // /////////////////// pManager = G4AntiNeutron::AntiNeutron()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiNeutronInelasticProcess* aninelProc = new G4AntiNeutronInelasticProcess(); G4LEAntiNeutronInelastic* LEPanModel = new G4LEAntiNeutronInelastic(); aninelProc->RegisterMe(LEPanModel); G4HEAntiNeutronInelastic* HEPanModel = new G4HEAntiNeutronInelastic(); aninelProc->RegisterMe(HEPanModel); pManager->AddDiscreteProcess(aninelProc); // anti-neutron annihilation at rest G4AntiNeutronAnnihilationAtRest* anAnnihil = new G4AntiNeutronAnnihilationAtRest(); pManager->AddRestProcess(anAnnihil); /////////////////// // // // Lambda // // // /////////////////// pManager = G4Lambda::Lambda()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4LambdaInelasticProcess* linelProc = new G4LambdaInelasticProcess(); linelProc->RegisterMe(bertiniModelStrange); G4LELambdaInelastic* LEPlModel = new G4LELambdaInelastic(); LEPlModel->SetMinEnergy(LEPLowerLimitForHyperons); linelProc->RegisterMe(LEPlModel); G4HELambdaInelastic* HEPlModel = new G4HELambdaInelastic(); linelProc->RegisterMe(HEPlModel); pManager->AddDiscreteProcess(linelProc); /////////////////// // // // Anti-Lambda // // // /////////////////// pManager = G4AntiLambda::AntiLambda()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiLambdaInelasticProcess* alinelProc = new G4AntiLambdaInelasticProcess(); G4LEAntiLambdaInelastic* LEPalModel = new G4LEAntiLambdaInelastic(); alinelProc->RegisterMe(LEPalModel); G4HEAntiLambdaInelastic* HEPalModel = new G4HEAntiLambdaInelastic(); alinelProc->RegisterMe(HEPalModel); pManager->AddDiscreteProcess(alinelProc); /////////////////// // // // Sigma- // // // /////////////////// pManager = G4SigmaMinus::SigmaMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4SigmaMinusInelasticProcess* sminelProc = new G4SigmaMinusInelasticProcess(); sminelProc->RegisterMe(bertiniModelStrange); G4LESigmaMinusInelastic* LEPsmModel = new G4LESigmaMinusInelastic(); LEPsmModel->SetMinEnergy(LEPLowerLimitForHyperons); sminelProc->RegisterMe(LEPsmModel); G4HESigmaMinusInelastic* HEPsmModel = new G4HESigmaMinusInelastic(); sminelProc->RegisterMe(HEPsmModel); pManager->AddDiscreteProcess(sminelProc); /////////////////// // // // Anti-Sigma- // // // /////////////////// pManager = G4AntiSigmaMinus::AntiSigmaMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiSigmaMinusInelasticProcess* asminelProc = new G4AntiSigmaMinusInelasticProcess(); G4LEAntiSigmaMinusInelastic* LEPasmModel = new G4LEAntiSigmaMinusInelastic(); asminelProc->RegisterMe(LEPasmModel); G4HEAntiSigmaMinusInelastic* HEPasmModel = new G4HEAntiSigmaMinusInelastic(); asminelProc->RegisterMe(HEPasmModel); pManager->AddDiscreteProcess(asminelProc); /////////////////// // // // Sigma+ // // // /////////////////// pManager = G4SigmaPlus::SigmaPlus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4SigmaPlusInelasticProcess* spinelProc = new G4SigmaPlusInelasticProcess(); spinelProc->RegisterMe(bertiniModelStrange); G4LESigmaPlusInelastic* LEPspModel = new G4LESigmaPlusInelastic(); LEPspModel->SetMinEnergy(LEPLowerLimitForHyperons); spinelProc->RegisterMe(LEPspModel); G4HESigmaPlusInelastic* HEPspModel = new G4HESigmaPlusInelastic(); spinelProc->RegisterMe(HEPspModel); pManager->AddDiscreteProcess(spinelProc); /////////////////// // // // Anti-Sigma+ // // // /////////////////// pManager = G4AntiSigmaPlus::AntiSigmaPlus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiSigmaPlusInelasticProcess* aspinelProc = new G4AntiSigmaPlusInelasticProcess(); G4LEAntiSigmaPlusInelastic* LEPaspModel = new G4LEAntiSigmaPlusInelastic(); aspinelProc->RegisterMe(LEPaspModel); G4HEAntiSigmaPlusInelastic* HEPaspModel = new G4HEAntiSigmaPlusInelastic(); aspinelProc->RegisterMe(HEPaspModel); pManager->AddDiscreteProcess(aspinelProc); /////////////////// // // // Xi- // // // /////////////////// pManager = G4XiMinus::XiMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4XiMinusInelasticProcess* xminelProc = new G4XiMinusInelasticProcess(); xminelProc->RegisterMe(bertiniModelStrange); G4LEXiMinusInelastic* LEPxmModel = new G4LEXiMinusInelastic(); LEPxmModel->SetMinEnergy(LEPLowerLimitForHyperons); xminelProc->RegisterMe(LEPxmModel); G4HEXiMinusInelastic* HEPxmModel = new G4HEXiMinusInelastic(); xminelProc->RegisterMe(HEPxmModel); pManager->AddDiscreteProcess(xminelProc); /////////////////// // // // Anti-Xi- // // // /////////////////// pManager = G4AntiXiMinus::AntiXiMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiXiMinusInelasticProcess* axminelProc = new G4AntiXiMinusInelasticProcess(); G4LEAntiXiMinusInelastic* LEPaxmModel = new G4LEAntiXiMinusInelastic(); axminelProc->RegisterMe(LEPaxmModel); G4HEAntiXiMinusInelastic* HEPaxmModel = new G4HEAntiXiMinusInelastic(); axminelProc->RegisterMe(HEPaxmModel); pManager->AddDiscreteProcess(axminelProc); /////////////////// // // // Xi0 // // // /////////////////// pManager = G4XiZero::XiZero()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4XiZeroInelasticProcess* x0inelProc = new G4XiZeroInelasticProcess(); x0inelProc->RegisterMe(bertiniModelStrange); G4LEXiZeroInelastic* LEPx0Model = new G4LEXiZeroInelastic(); LEPx0Model->SetMinEnergy(LEPLowerLimitForHyperons); x0inelProc->RegisterMe(LEPx0Model); G4HEXiZeroInelastic* HEPx0Model = new G4HEXiZeroInelastic(); x0inelProc->RegisterMe(HEPx0Model); pManager->AddDiscreteProcess(x0inelProc); /////////////////// // // // Anti-Xi0 // // // /////////////////// pManager = G4AntiXiZero::AntiXiZero()->GetProcessManager(); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiXiZeroInelasticProcess* ax0inelProc = new G4AntiXiZeroInelasticProcess(); G4LEAntiXiZeroInelastic* LEPax0Model = new G4LEAntiXiZeroInelastic(); ax0inelProc->RegisterMe(LEPax0Model); G4HEAntiXiZeroInelastic* HEPax0Model = new G4HEAntiXiZeroInelastic(); ax0inelProc->RegisterMe(HEPax0Model); pManager->AddDiscreteProcess(ax0inelProc); /////////////////// // // // Omega- // // // /////////////////// pManager = G4OmegaMinus::OmegaMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4OmegaMinusInelasticProcess* ominelProc = new G4OmegaMinusInelasticProcess(); G4LEOmegaMinusInelastic* LEPomModel = new G4LEOmegaMinusInelastic(); ominelProc->RegisterMe(LEPomModel); G4HEOmegaMinusInelastic* HEPomModel = new G4HEOmegaMinusInelastic(); ominelProc->RegisterMe(HEPomModel); pManager->AddDiscreteProcess(ominelProc); /////////////////// // // // Anti-Omega- // // // /////////////////// pManager = G4AntiOmegaMinus::AntiOmegaMinus()->GetProcessManager(); // EM processes pManager->AddProcess(new G4MultipleScattering(), -1, 1, 1); pManager->AddProcess(new G4hIonisation(), -1, 2, 2); // hadron elastic pManager->AddDiscreteProcess(elasticProcess); // hadron inelastic G4AntiOmegaMinusInelasticProcess* aominelProc = new G4AntiOmegaMinusInelasticProcess(); G4LEAntiOmegaMinusInelastic* LEPaomModel = new G4LEAntiOmegaMinusInelastic(); aominelProc->RegisterMe(LEPaomModel); G4HEAntiOmegaMinusInelastic* HEPaomModel = new G4HEAntiOmegaMinusInelastic(); aominelProc->RegisterMe(HEPaomModel); pManager->AddDiscreteProcess(aominelProc); |
Ion Physics
...
- G4PionPlus
- G4PionMinus
- G4KaonPlus
- G4KaonMinus
- G4KaonZeroLong
- G4KaonZeroShort
- G4Proton
- G4AntiProton
- G4Neutron
- G4AntiNeutron
- G4Lambda
- G4AntiLambda
- G4SigmaMinus
- G4AntiSigmaMinus
- G4SigmaPlus
- G4AntiSigmaPlus
- G4XiMinus
- G4AntiXiMinus
- G4XiZero
- G4AntiXiZero
- G4OmegaMinus
- G4AntiOmegaMinus
Ion Physics
...