Loading...
Searching...
No Matches
Public Member Functions |
Protected Member Functions |
Protected Attributes |
Static Protected Attributes |
List of all members
G4INCL::CrossSectionsMultiPionsAndResonances Class Reference
Multipion and mesonic Resonances cross sections. More...
#include <G4INCLCrossSectionsMultiPionsAndResonances.hh>
Inheritance diagram for G4INCL::CrossSectionsMultiPionsAndResonances:
Public Member Functions | |
| CrossSectionsMultiPionsAndResonances () | |
| virtual G4double | elastic (Particle const *const p1, Particle const *const p2) |
| new elastic particle-particle cross section | |
| virtual G4double | total (Particle const *const p1, Particle const *const p2) |
| new total particle-particle cross section | |
| virtual G4double | piNToxPiN (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - piN Channel (modified due to the mesonic resonances) | |
| virtual G4double | piNToEtaN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance production - piN Channel. | |
| virtual G4double | piNToOmegaN (Particle const *const p1, Particle const *const p2) |
| Cross section for PiN->OmegaN. | |
| virtual G4double | piNToEtaPrimeN (Particle const *const p1, Particle const *const p2) |
| Cross section for PiN->EtaPrimeN. | |
| virtual G4double | etaNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - piN Channel. | |
| virtual G4double | omegaNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross section for OmegaN->PiN. | |
| virtual G4double | etaPrimeNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross section for EtaPrimeN->PiN. | |
| virtual G4double | etaNToPiPiN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - pipiN Channel. | |
| virtual G4double | NNToNNEta (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production (inclusive) - NN entrance channel. | |
| virtual G4double | NNToNNEtaExclu (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToNNOmega (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Omega production (inclusive) - NN entrance channel. | |
| virtual G4double | NNToNNOmegaExclu (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Omega production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToxPiNN (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NN Channel. | |
| virtual G4double | NNToNNEtaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NNEta Channel. | |
| virtual G4double | NNToNDeltaEta (Particle const *const p1, Particle const *const p2) |
| Cross section for N-Delta-Eta production - NNEta Channel. | |
| virtual G4double | NNToNNOmegaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NNOmega Channel. | |
| virtual G4double | NNToNDeltaOmega (Particle const *const p1, Particle const *const p2) |
| Cross section for N-Delta-Eta production - NNOmega Channel. | |
| Public Member Functions inherited from G4INCL::CrossSectionsMultiPions | |
| CrossSectionsMultiPions () | |
| virtual G4double | elastic (Particle const *const p1, Particle const *const p2) |
| Elastic particle-particle cross section. | |
| virtual G4double | total (Particle const *const p1, Particle const *const p2) |
| Total (elastic+inelastic) particle-particle cross section. | |
| virtual G4double | NDeltaToNN (Particle const *const p1, Particle const *const p2) |
| Cross section for NDelta->NN. | |
| virtual G4double | NNToNDelta (Particle const *const p1, Particle const *const p2) |
| Cross section for Delta production - NN Channel. | |
| virtual G4double | piNToDelta (Particle const *const p1, Particle const *const p2) |
| Cross section for Delta production - piN Channel. | |
| virtual G4double | piNToxPiN (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - piN Channel. | |
| virtual G4double | NNToxPiNN (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NN Channel. | |
| virtual G4double | calculateNNAngularSlope (G4double energyCM, G4int iso) |
| Calculate the slope of the NN DDXS. | |
| virtual G4double | piNToEtaN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance production - piN Channel. | |
| virtual G4double | piNToOmegaN (Particle const *const p1, Particle const *const p2) |
| Cross section for PiN->OmegaN. | |
| virtual G4double | piNToEtaPrimeN (Particle const *const p1, Particle const *const p2) |
| Cross section for PiN->EtaPrimeN. | |
| virtual G4double | etaNToPiPiN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - pipiN Channel. | |
| virtual G4double | omegaNToPiPiN (Particle const *const p1, Particle const *const p2) |
| Cross section for OmegaN->PiPiN. | |
| virtual G4double | etaNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - piN Channel. | |
| virtual G4double | omegaNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross section for OmegaN->PiN. | |
| virtual G4double | etaPrimeNToPiN (Particle const *const p1, Particle const *const p2) |
| Cross section for EtaPrimeN->PiN. | |
| virtual G4double | NNToNNEta (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production - NN entrance channel. | |
| virtual G4double | NNToNNEtaExclu (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToNNEtaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NNEta Channel. | |
| virtual G4double | NNToNDeltaEta (Particle const *const p1, Particle const *const p2) |
| Cross section for N-Delta-Eta production - NNEta Channel. | |
| virtual G4double | NNToNNOmega (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production - NN entrance channel. | |
| virtual G4double | NNToNNOmegaExclu (Particle const *const particle1, Particle const *const particle2) |
| Cross section for Eta production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToNNOmegaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2) |
| Cross section for X pion production - NNEta Channel. | |
| virtual G4double | NNToNDeltaOmega (Particle const *const p1, Particle const *const p2) |
| Cross section for N-Delta-Eta production - NNEta Channel. | |
| virtual G4double | NYelastic (Particle const *const p1, Particle const *const p2) |
| elastic scattering for Nucleon-Strange Particles cross sections | |
| virtual G4double | NKbelastic (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKelastic (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNLK (Particle const *const p1, Particle const *const p2) |
| Nucleon-Nucleon to Stange particles cross sections. | |
| virtual G4double | NNToNSK (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNLKpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNSKpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNLK2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNSK2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToNNKKb (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NNToMissingStrangeness (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NDeltaToNLK (Particle const *const p1, Particle const *const p2) |
| Nucleon-Delta to Stange particles cross sections. | |
| virtual G4double | NDeltaToNSK (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NDeltaToDeltaLK (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NDeltaToDeltaSK (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NDeltaToNNKKb (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToLK (Particle const *const p1, Particle const *const p2) |
| Nucleon-Pion to Stange particles cross sections. | |
| virtual G4double | NpiToSK (Particle const *const p1, Particle const *const p2) |
| virtual G4double | p_pimToSzKz (Particle const *const p1, Particle const *const p2) |
| virtual G4double | p_pimToSmKp (Particle const *const p1, Particle const *const p2) |
| virtual G4double | p_pizToSzKp (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToLKpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToSKpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToLK2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToSK2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToNKKb (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NpiToMissingStrangeness (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NLToNS (Particle const *const p1, Particle const *const p2) |
| Nucleon-Hyperon cross sections. | |
| virtual G4double | NSToNL (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NSToNS (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKToNK (Particle const *const p1, Particle const *const p2) |
| Nucleon-Kaon quasi-elastic and inelastic cross sections. | |
| virtual G4double | NKToNKpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKToNK2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToNKb (Particle const *const p1, Particle const *const p2) |
| Nucleon-antiKaon quasi-elastic and inelastic cross sections. | |
| virtual G4double | NKbToSpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToLpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToS2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToL2pi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToNKbpi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | NKbToNKb2pi (Particle const *const p1, Particle const *const p2) |
| Public Member Functions inherited from G4INCL::ICrossSections | |
| ICrossSections () | |
| virtual | ~ICrossSections () |
| virtual G4double | elastic (Particle const *const p1, Particle const *const p2)=0 |
| Elastic particle-particle cross section. | |
| virtual G4double | total (Particle const *const p1, Particle const *const p2)=0 |
| Total (elastic+inelastic) particle-particle cross section. | |
| virtual G4double | NDeltaToNN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NDelta->NN. | |
| virtual G4double | NNToNDelta (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NDelta. | |
| virtual G4double | NNToxPiNN (const G4int xpi, Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->xpiN. | |
| virtual G4double | piNToDelta (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for piN->NDelta. | |
| virtual G4double | piNToxPiN (const G4int xpi, Particle const *const p1, Particle const *const p2)=0 |
| Cross section for piN->piNpi. | |
| virtual G4double | piNToEtaN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for PiN->EtaN. | |
| virtual G4double | piNToOmegaN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for PiN->OmegaN. | |
| virtual G4double | piNToEtaPrimeN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for PiN->EtaPrimeN. | |
| virtual G4double | etaNToPiN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for EtaN->PiN. | |
| virtual G4double | etaNToPiPiN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for EtaN->PiPiN. | |
| virtual G4double | omegaNToPiN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for OmegaN->PiN. | |
| virtual G4double | omegaNToPiPiN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for OmegaN->PiPiN. | |
| virtual G4double | etaPrimeNToPiN (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for EtaPrimeN->PiN. | |
| virtual G4double | NNToNNEta (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEta (inclusive) | |
| virtual G4double | NNToNNEtaExclu (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEta (exclusive) | |
| virtual G4double | NNToNNEtaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEtaxPi. | |
| virtual G4double | NNToNDeltaEta (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for N-Delta-Eta production - NNEta Channel. | |
| virtual G4double | NNToNNOmega (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEta (inclusive) | |
| virtual G4double | NNToNNOmegaExclu (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEta (exclusive) | |
| virtual G4double | NNToNNOmegaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2)=0 |
| Cross section for NN->NNEtaxPi. | |
| virtual G4double | NNToNDeltaOmega (Particle const *const p1, Particle const *const p2)=0 |
| Cross section for N-Delta-Eta production - NNEta Channel. | |
| virtual G4double | NYelastic (Particle const *const p1, Particle const *const p2)=0 |
| elastic scattering for Nucleon-Strange Particles cross sections | |
| virtual G4double | NKbelastic (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKelastic (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNLK (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-Nucleon to Stange particles cross sections. | |
| virtual G4double | NNToNSK (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNLKpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNSKpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNLK2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNSK2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToNNKKb (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NNToMissingStrangeness (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NDeltaToNLK (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-Delta to Stange particles cross sections. | |
| virtual G4double | NDeltaToNSK (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NDeltaToDeltaLK (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NDeltaToDeltaSK (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NDeltaToNNKKb (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToLK (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-Pion to Stange particles cross sections. | |
| virtual G4double | NpiToSK (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | p_pimToSzKz (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | p_pimToSmKp (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | p_pizToSzKp (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToLKpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToSKpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToLK2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToSK2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToNKKb (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NpiToMissingStrangeness (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NLToNS (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-Hyperon cross sections. | |
| virtual G4double | NSToNL (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NSToNS (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKToNK (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-Kaon inelastic cross sections. | |
| virtual G4double | NKToNKpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKToNK2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToNKb (Particle const *const p1, Particle const *const p2)=0 |
| Nucleon-antiKaon inelastic cross sections. | |
| virtual G4double | NKbToSpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToLpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToS2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToL2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToNKbpi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | NKbToNKb2pi (Particle const *const p1, Particle const *const p2)=0 |
| virtual G4double | calculateNNAngularSlope (G4double energyCM, G4int iso)=0 |
| Calculate the slope of the NN DDXS. | |
Protected Member Functions | |
| G4double | piMinuspToEtaN (Particle const *const p1, Particle const *const p2) |
| Internal function for pion cross sections. | |
| G4double | piMinuspToEtaN (const G4double ECM) |
| G4double | piMinuspToOmegaN (Particle const *const p1, Particle const *const p2) |
| G4double | piMinuspToOmegaN (const G4double ECM) |
| virtual G4double | NNToNNEtaIso (const G4double ener, const G4int iso) |
| Cross section for One (more) pion production - piN entrance channel. | |
| virtual G4double | NNToNNEtaExcluIso (const G4double ener, const G4int iso) |
| Isotopic Cross section for Eta production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToNNEtaOnePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NNEta channel. | |
| virtual G4double | NNToNNEtaOnePiOrDelta (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NNEta channel. | |
| virtual G4double | NNToNNEtaTwoPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 2-pion production - NNEta channel. | |
| virtual G4double | NNToNNEtaThreePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 3-pion production - NNEta channel. | |
| virtual G4double | NNToNNEtaFourPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 4-pion production - NNEta channel. | |
| virtual G4double | NNToNNOmegaIso (const G4double ener, const G4int iso) |
| Isotopic Cross section for Omega production (inclusive) - NN entrance channel. | |
| virtual G4double | NNToNNOmegaExcluIso (const G4double ener, const G4int iso) |
| Isotopic Cross section for Omega production (exclusive) - NN entrance channel. | |
| virtual G4double | NNToNNOmegaOnePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NNOmega channel. | |
| virtual G4double | NNToNNOmegaOnePiOrDelta (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NNOmega channel. | |
| virtual G4double | NNToNNOmegaTwoPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 2-pion production - NNOmega channel. | |
| virtual G4double | NNToNNOmegaThreePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 3-pion production - NNOmega channel. | |
| virtual G4double | NNToNNOmegaFourPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 4-pion production - NNOmega channel. | |
| virtual G4double | etaNElastic (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - elastic Channel. | |
| virtual G4double | omegaNElastic (Particle const *const p1, Particle const *const p2) |
| virtual G4double | omegaNInelastic (Particle const *const p1, Particle const *const p2) |
| Cross sections for mesonic resonance absorption on nucleon - inelastic Channel. | |
| virtual G4double | omegaNToPiPiN (Particle const *const p1, Particle const *const p2) |
| Cross sections for omega-induced 2Pi emission on nucleon. | |
| Protected Member Functions inherited from G4INCL::CrossSectionsMultiPions | |
| G4double | NNElastic (Particle const *const part1, Particle const *const part2) |
| Internal implementation of the NN elastic cross section. | |
| G4double | NNElasticFixed (const G4double s, const G4int i) |
| Internal implementation of the NN elastic cross section with fixed isospin. | |
| G4double | NNTot (Particle const *const part1, Particle const *const part2) |
| Internal implementation of the NN total cross section. | |
| G4double | NNTotFixed (const G4double s, const G4int i) |
| Internal implementation of the NN total cross section with fixed isospin. | |
| G4double | NNInelasticIso (const G4double ener, const G4int iso) |
| Internal implementation of the isospin dependent NN reaction cross section. | |
| virtual G4double | NNOnePiOrDelta (const G4double ener, const G4int iso, const G4double xsiso) |
| Cross section for direct 1-pion production + delta production - NN entrance channel. | |
| virtual G4double | NNTwoPi (const G4double ener, const G4int iso, const G4double xsiso) |
| Cross section for direct 2-pion production - NN entrance channel. | |
| virtual G4double | NNThreePi (const G4double ener, const G4int iso, const G4double xsiso, const G4double xs1pi, const G4double xs2pi) |
| Cross section for direct 3-pion production - NN entrance channel. | |
| virtual G4double | NNOnePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NN entrance channel. | |
| virtual G4double | NNOnePiOrDelta (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 1-pion production - NN entrance channel. | |
| virtual G4double | NNTwoPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 2-pion production - NN entrance channel. | |
| virtual G4double | NNThreePi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 3-pion production - NN entrance channel. | |
| virtual G4double | NNFourPi (Particle const *const part1, Particle const *const part2) |
| Cross section for direct 4-pion production - NN entrance channel. | |
| G4double | spnPiPlusPHE (const G4double x) |
| Internal function for pion cross sections. | |
| G4double | spnPiMinusPHE (const G4double x) |
| Internal function for pion cross sections. | |
| G4double | piNIne (Particle const *const p1, Particle const *const p2) |
| G4double | piNTot (Particle const *const p1, Particle const *const p2) |
| G4double | piNTopiN (Particle const *const p1, Particle const *const p2) |
| G4double | piPluspIne (Particle const *const p1, Particle const *const p2) |
| G4double | piMinuspIne (Particle const *const p1, Particle const *const p2) |
| G4double | piPluspOnePi (Particle const *const p1, Particle const *const p2) |
| G4double | piMinuspOnePi (Particle const *const p1, Particle const *const p2) |
| G4double | piPluspTwoPi (Particle const *const p1, Particle const *const p2) |
| G4double | piMinuspTwoPi (Particle const *const p1, Particle const *const p2) |
| virtual G4double | piNOnePi (Particle const *const p1, Particle const *const p2) |
| Cross section for One (more) pion production - piN entrance channel. | |
| virtual G4double | piNTwoPi (Particle const *const p1, Particle const *const p2) |
| Cross section for Two (more) pion production - piN entrance channel. | |
Protected Attributes | |
| const HornerC7 | s11pzHC |
| Horner coefficients for s11pz. | |
| const HornerC8 | s01ppHC |
| Horner coefficients for s01pp. | |
| const HornerC4 | s01pzHC |
| Horner coefficients for s01pz. | |
| const HornerC4 | s11pmHC |
| Horner coefficients for s11pm. | |
| const HornerC5 | s12pmHC |
| Horner coefficients for s12pm. | |
| const HornerC3 | s12ppHC |
| Horner coefficients for s12pp. | |
| const HornerC4 | s12zzHC |
| Horner coefficients for s12zz. | |
| const HornerC4 | s02pzHC |
| Horner coefficients for s02pz. | |
| const HornerC6 | s02pmHC |
| Horner coefficients for s02pm. | |
| const HornerC4 | s12mzHC |
| Horner coefficients for s12mz. | |
| Protected Attributes inherited from G4INCL::CrossSectionsMultiPions | |
| const HornerC7 | s11pzHC |
| Horner coefficients for s11pz. | |
| const HornerC8 | s01ppHC |
| Horner coefficients for s01pp. | |
| const HornerC4 | s01pzHC |
| Horner coefficients for s01pz. | |
| const HornerC4 | s11pmHC |
| Horner coefficients for s11pm. | |
| const HornerC5 | s12pmHC |
| Horner coefficients for s12pm. | |
| const HornerC3 | s12ppHC |
| Horner coefficients for s12pp. | |
| const HornerC4 | s12zzHC |
| Horner coefficients for s12zz. | |
| const HornerC4 | s02pzHC |
| Horner coefficients for s02pz. | |
| const HornerC6 | s02pmHC |
| Horner coefficients for s02pm. | |
| const HornerC4 | s12mzHC |
| Horner coefficients for s12mz. | |
Static Protected Attributes | |
| static const G4int | nMaxPiNN = 4 |
| Maximum number of outgoing pions in NN collisions. | |
| static const G4int | nMaxPiPiN = 4 |
| Maximum number of outgoing pions in piN collisions. | |
| static const G4double | s11pzOOT = 0.0035761542037692665889 |
| One over threshold for s11pz. | |
| static const G4double | s01ppOOT = 0.003421025623481919853 |
| One over threshold for s01pp. | |
| static const G4double | s01pzOOT = 0.0035739814152966403123 |
| One over threshold for s01pz. | |
| static const G4double | s11pmOOT = 0.0034855350296270480281 |
| One over threshold for s11pm. | |
| static const G4double | s12pmOOT = 0.0016672224074691565119 |
| One over threshold for s12pm. | |
| static const G4double | s12ppOOT = 0.0016507643038726931312 |
| One over threshold for s12pp. | |
| static const G4double | s12zzOOT = 0.0011111111111111111111 |
| One over threshold for s12zz. | |
| static const G4double | s02pzOOT = 0.00125 |
| One over threshold for s02pz. | |
| static const G4double | s02pmOOT = 0.0016661112962345883443 |
| One over threshold for s02pm. | |
| static const G4double | s12mzOOT = 0.0017047391749062392793 |
| One over threshold for s12mz. | |
| Static Protected Attributes inherited from G4INCL::CrossSectionsMultiPions | |
| static const G4int | nMaxPiNN = 4 |
| Maximum number of outgoing pions in NN collisions. | |
| static const G4int | nMaxPiPiN = 4 |
| Maximum number of outgoing pions in piN collisions. | |
| static const G4double | s11pzOOT = 0.0035761542037692665889 |
| One over threshold for s11pz. | |
| static const G4double | s01ppOOT = 0.003421025623481919853 |
| One over threshold for s01pp. | |
| static const G4double | s01pzOOT = 0.0035739814152966403123 |
| One over threshold for s01pz. | |
| static const G4double | s11pmOOT = 0.0034855350296270480281 |
| One over threshold for s11pm. | |
| static const G4double | s12pmOOT = 0.0016672224074691565119 |
| One over threshold for s12pm. | |
| static const G4double | s12ppOOT = 0.0016507643038726931312 |
| One over threshold for s12pp. | |
| static const G4double | s12zzOOT = 0.0011111111111111111111 |
| One over threshold for s12zz. | |
| static const G4double | s02pzOOT = 0.00125 |
| One over threshold for s02pz. | |
| static const G4double | s02pmOOT = 0.0016661112962345883443 |
| One over threshold for s02pm. | |
| static const G4double | s12mzOOT = 0.0017047391749062392793 |
| One over threshold for s12mz. | |
Detailed Description
Multipion and mesonic Resonances cross sections.
Definition at line 54 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
Constructor & Destructor Documentation
◆ CrossSectionsMultiPionsAndResonances()
| G4INCL::CrossSectionsMultiPionsAndResonances::CrossSectionsMultiPionsAndResonances | ( | ) |
Definition at line 77 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
77 :
78 s11pzHC(-2.228000000000294018,8.7560000000005723725,-0.61000000000023239325,-5.4139999999999780324,3.3338333333333348023,-0.75835000000000022049,0.060623611111111114688),
79 s01ppHC(2.0570000000126518344,-6.029000000012135826,36.768500000002462784,-45.275666666666553533,25.112666666666611953,-7.2174166666666639187,1.0478875000000000275,-0.060804365079365080846),
80 s01pzHC(0.18030000000000441851,7.8700999999999953598,-4.0548999999999990425,0.555199999999999959),
81 s11pmHC(0.20590000000000031866,3.3450999999999993936,-1.4401999999999997825,0.17076666666666664973),
82 s12pmHC(-0.77235999999999901328,4.2626599999999991117,-1.9008899999999997323,0.30192266666666663379,-0.012270833333333331986),
83 s12ppHC(-0.75724999999999975664,2.0934399999999998565,-0.3803099999999999814),
84 s12zzHC(-0.89599999999996965072,7.882999999999978632,-7.1049999999999961928,1.884333333333333089),
85 s02pzHC(-1.0579999999999967036,11.113999999999994089,-8.5259999999999990196,2.0051666666666666525),
86 s02pmHC(2.4009000000012553286,-7.7680000000013376183,20.619000000000433505,-16.429666666666723928,5.2525708333333363472,-0.58969166666666670206),
87 s12mzHC(-0.21858699999999976269,1.9148999999999999722,-0.31727500000000001065,-0.027695000000000000486)
88 {
89 }
const HornerC4 s01pzHC
Horner coefficients for s01pz.
const HornerC7 s11pzHC
Horner coefficients for s11pz.
const HornerC4 s02pzHC
Horner coefficients for s02pz.
const HornerC5 s12pmHC
Horner coefficients for s12pm.
const HornerC4 s11pmHC
Horner coefficients for s11pm.
const HornerC8 s01ppHC
Horner coefficients for s01pp.
const HornerC4 s12zzHC
Horner coefficients for s12zz.
const HornerC3 s12ppHC
Horner coefficients for s12pp.
const HornerC4 s12mzHC
Horner coefficients for s12mz.
const HornerC6 s02pmHC
Horner coefficients for s02pm.
Member Function Documentation
◆ elastic()
|
virtual |
new elastic particle-particle cross section
Reimplemented from G4INCL::CrossSectionsMultiPions.
Reimplemented in G4INCL::CrossSectionsStrangeness.
Definition at line 118 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
118 {
119 if((p1->isNucleon()||p1->isDelta()) && (p2->isNucleon()||p2->isDelta())){
121 }
122 else if ((p1->isNucleon() && p2->isPion()) || (p2->isNucleon() && p1->isPion())){
124 }
125 else if ((p1->isNucleon() && p2->isEta()) || (p2->isNucleon() && p1->isEta())){
127 }
128 else if ((p1->isNucleon() && p2->isOmega()) || (p2->isNucleon() && p1->isOmega())){
130 }
131 else {
132 return 0.0;
133 }
134 }
virtual G4double etaNElastic(Particle const *const p1, Particle const *const p2)
Cross sections for mesonic resonance absorption on nucleon - elastic Channel.
virtual G4double omegaNElastic(Particle const *const p1, Particle const *const p2)
virtual G4double elastic(Particle const *const p1, Particle const *const p2)
Elastic particle-particle cross section.
Definition: G4INCLCrossSectionsMultiPions.cc:807
Referenced by total().
◆ etaNElastic()
|
protectedvirtual |
Cross sections for mesonic resonance absorption on nucleon - elastic Channel.
Definition at line 340 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
340 {
341 //
342 // Eta-Nucleon elastic cross sections
343 //
344// assert((particle1->isNucleon() && particle2->isEta()) || (particle1->isEta() && particle2->isNucleon()));
345
346 G4double sigma=0.;
347
348 const Particle *eta;
350
351 if (particle1->isEta()) {
352 eta = particle1;
353 nucleon = particle2;
354 }
355 else {
356 eta = particle2;
357 nucleon = particle1;
358 }
359
361
362 if (pLab < 700.)
363 sigma = 3.6838e-15*std::pow(pLab,6) - 9.7815e-12*std::pow(pLab,5) + 9.7914e-9*std::pow(pLab,4) -
364 4.3222e-06*std::pow(pLab,3) + 7.9188e-04*std::pow(pLab,2) - 1.8379e-01*pLab + 84.965;
365 else if (pLab < 1400.)
366 sigma = 3.562630e-16*std::pow(pLab,6) - 2.384766e-12*std::pow(pLab,5) + 6.601312e-9*std::pow(pLab,4) -
367 9.667078e-06*std::pow(pLab,3) + 7.894845e-03*std::pow(pLab,2) - 3.409200*pLab + 609.8501;
368 else if (pLab < 2025.)
369 sigma = -1.041950E-03*pLab + 2.110529E+00;
370 else
371 sigma=0.;
372
373 if (sigma < 0.) sigma = 0.;
374 return sigma; // Parameterization from the ANL-Osaka DCC model [PRC88(2013)035209]
375 }
G4double momentumInLab(Particle const *const p1, Particle const *const p2)
gives the momentum in the lab frame of two particles.
Definition: G4INCLKinematicsUtils.cc:135
Referenced by elastic(), and G4INCL::CrossSectionsStrangeness::elastic().
◆ etaNToPiN()
|
virtual |
Cross sections for mesonic resonance absorption on nucleon - piN Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 256 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
256 {
257 //
258 // Eta-Nucleon producing Pion cross sections
259 //
260// assert((particle1->isNucleon() && particle2->isEta()) || (particle1->isEta() && particle2->isNucleon()));
261
262 const Particle *eta;
264
265 if (particle1->isEta()) {
266 eta = particle1;
267 nucleon = particle2;
268 }
269 else {
270 eta = particle2;
271 nucleon = particle1;
272 }
273
275 G4double sigma=0.;
276
277 if (pLab <= 574.)
278 sigma= 1.511147E-13*std::pow(pLab,6)- 3.603636E-10*std::pow(pLab,5)+ 3.443487E-07*std::pow(pLab,4)- 1.681980E-04*std::pow(pLab,3)+ 4.437913E-02*std::pow(pLab,2)- 6.172108E+00*pLab+ 4.031449E+02;
279 else if (pLab <= 850.)
280 sigma= -8.00018E-14*std::pow(pLab,6)+ 3.50041E-10*std::pow(pLab,5)- 6.33891E-07*std::pow(pLab,4)+ 6.07658E-04*std::pow(pLab,3)- 3.24936E-01*std::pow(pLab,2)+ 9.18098E+01*pLab- 1.06943E+04;
281 else if (pLab <= 1300.)
282 sigma= 6.56364E-09*std::pow(pLab,3)- 2.07653E-05*std::pow(pLab,2)+ 1.84148E-02*pLab- 1.70427E+00;
283 else {
288 G4double masseta;
289 G4double massnucleon;
290 G4double pCM_eta;
291 masseta=eta->getMass();
292 massnucleon=nucleon->getMass();
293 pCM_eta=KinematicsUtils::momentumInCM(ECM, masseta, massnucleon);
295 G4double pCM_PiMinus=KinematicsUtils::momentumInCM(ECM, massPiMinus, massProton); // = pCM_PiPlus (because massPiMinus = massPiPlus)
296 sigma = (piMinuspToEtaN(ECM)/2.) * std::pow((pCM_PiZero/pCM_eta), 2) + piMinuspToEtaN(ECM) * std::pow((pCM_PiMinus/pCM_eta), 2);
297 }
298 if (sigma < 0.) sigma=0.;
299
300 return sigma;
301 }
G4double piMinuspToEtaN(Particle const *const p1, Particle const *const p2)
Internal function for pion cross sections.
G4double totalEnergyInCM(Particle const *const p1, Particle const *const p2)
Definition: G4INCLKinematicsUtils.cc:94
G4double momentumInCM(Particle const *const p1, Particle const *const p2)
gives the momentum in the CM frame of two particles.
Definition: G4INCLKinematicsUtils.cc:107
G4double getINCLMass(const G4int A, const G4int Z, const G4int S)
Get INCL nuclear mass (in MeV/c^2)
Definition: G4INCLParticleTable.cc:794
Referenced by etaNToPiPiN(), total(), and G4INCL::CrossSectionsStrangeness::total().
◆ etaNToPiPiN()
|
virtual |
Cross sections for mesonic resonance absorption on nucleon - pipiN Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 303 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
303 {
304 //
305 // Eta-Nucleon producing Two Pions cross sections
306 //
307// assert((particle1->isNucleon() && particle2->isEta()) || (particle1->isEta() && particle2->isNucleon()));
308
309 G4double sigma=0.;
310
311 const Particle *eta;
313
314 if (particle1->isEta()) {
315 eta = particle1;
316 nucleon = particle2;
317 }
318 else {
319 eta = particle2;
320 nucleon = particle1;
321 }
322
324
325 if (pLab < 450.)
326 sigma = 2.01854221E-13*std::pow(pLab,6) - 3.49750459E-10*std::pow(pLab,5) + 2.46011585E-07*std::pow(pLab,4) - 9.01422901E-05*std::pow(pLab,3) + 1.83382964E-02*std::pow(pLab,2) - 2.03113098E+00*pLab + 1.10358550E+02;
327 else if (pLab < 600.)
328 sigma = 2.01854221E-13*std::pow(450.,6) - 3.49750459E-10*std::pow(450.,5) + 2.46011585E-07*std::pow(450.,4) - 9.01422901E-05*std::pow(450.,3) + 1.83382964E-02*std::pow(450.,2) - 2.03113098E+00*450. + 1.10358550E+02;
329 else if (pLab <= 1300.)
330 sigma = -6.32793049e-16*std::pow(pLab,6) + 3.95985900e-12*std::pow(pLab,5) - 1.01727714e-8*std::pow(pLab,4) +
331 1.37055547e-05*std::pow(pLab,3) - 1.01830486e-02*std::pow(pLab,2) + 3.93492126*pLab - 609.447145;
332 else
333 sigma = etaNToPiN(particle1,particle2);
334
335 if (sigma < 0.) sigma = 0.;
336 return sigma; // Parameterization from the ANL-Osaka DCC model [PRC88(2013)035209] - eta p --> "pi+pi0 n" + "pi0 pi0 p" total XS
337 }
virtual G4double etaNToPiN(Particle const *const p1, Particle const *const p2)
Cross sections for mesonic resonance absorption on nucleon - piN Channel.
Referenced by total(), and G4INCL::CrossSectionsStrangeness::total().
◆ etaPrimeNToPiN()
|
virtual |
Cross section for EtaPrimeN->PiN.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 493 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
493 {
494#else
495 G4double CrossSectionsMultiPionsAndResonances::etaPrimeNToPiN(Particle const * const particle1, Particle const * const particle2) {
496#endif
497 //
498 // EtaPrime-Nucleon producing Pion cross sections
499 //
500// assert((particle1->isNucleon() && particle2->isEtaPrime()) || (particle1->isEtaPrime() && particle2->isNucleon()));
501
502 return 0.;
503 }
virtual G4double etaPrimeNToPiN(Particle const *const p1, Particle const *const p2)
Cross section for EtaPrimeN->PiN.
Referenced by etaPrimeNToPiN(), total(), and G4INCL::CrossSectionsStrangeness::total().
◆ NNToNDeltaEta()
|
virtual |
Cross section for N-Delta-Eta production - NNEta Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 1092 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1092 {
1093// assert(p1->isNucleon() && p2->isNucleon());
1094 const G4int i = ParticleTable::getIsospin(p1->getType()) + ParticleTable::getIsospin(p2->getType());
1095 const G4double ener=KinematicsUtils::totalEnergyInCM(p1, p2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
1096 if (ener < 2018.563) return 0.;
1097 G4double xsinelas;
1098 if (i!=0)
1099 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1100 else
1101 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1102 if (xsinelas <= 1.e-9) return 0.;
1105 if(i==0)
1106 sigma *= 0.5;
1107 return sigma;
1108 }
virtual G4double NNToNNEta(Particle const *const particle1, Particle const *const particle2)
Cross section for Eta production (inclusive) - NN entrance channel.
virtual G4double NNToNNEtaExclu(Particle const *const particle1, Particle const *const particle2)
Cross section for Eta production (exclusive) - NN entrance channel.
virtual G4double NNToNNEtaOnePiOrDelta(Particle const *const part1, Particle const *const part2)
Cross section for direct 1-pion production - NNEta channel.
G4double NNInelasticIso(const G4double ener, const G4int iso)
Internal implementation of the isospin dependent NN reaction cross section.
Definition: G4INCLCrossSectionsMultiPions.cc:224
G4int getIsospin(const ParticleType t)
Get the isospin of a particle.
Definition: G4INCLParticleTable.cc:478
◆ NNToNDeltaOmega()
|
virtual |
Cross section for N-Delta-Eta production - NNOmega Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 1241 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1241 {
1242// assert(p1->isNucleon() && p2->isNucleon());
1243//jcd to be removed
1244// return 0.;
1245//jcd
1246 const G4int i = ParticleTable::getIsospin(p1->getType()) + ParticleTable::getIsospin(p2->getType());
1247 const G4double ener=KinematicsUtils::totalEnergyInCM(p1, p2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1248 if (ener < 2018.563) return 0.;
1249 G4double xsinelas;
1250 if (i!=0)
1251 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1252 else
1253 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1254 if (xsinelas <= 1.e-9) return 0.;
1257 if(i==0)
1258 sigma *= 0.5;
1259 return sigma;
1260 }
virtual G4double NNToNNOmegaExclu(Particle const *const particle1, Particle const *const particle2)
Cross section for Omega production (exclusive) - NN entrance channel.
virtual G4double NNToNNOmegaOnePiOrDelta(Particle const *const part1, Particle const *const part2)
Cross section for direct 1-pion production - NNOmega channel.
virtual G4double NNToNNOmega(Particle const *const particle1, Particle const *const particle2)
Cross section for Omega production (inclusive) - NN entrance channel.
◆ NNToNNEta()
|
virtual |
Cross section for Eta production (inclusive) - NN entrance channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 691 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
691 {
692
694 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
695
696 if (iso != 0) {
698 }
699 else {
701 }
702 }
virtual G4double NNToNNEtaIso(const G4double ener, const G4int iso)
Cross section for One (more) pion production - piN entrance channel.
Referenced by NNToNDeltaEta(), NNToNNEtaFourPi(), NNToNNEtaxPi(), NNToxPiNN(), and G4INCL::CrossSectionsStrangeness::NNToxPiNN().
◆ NNToNNEtaExclu()
|
virtual |
Cross section for Eta production (exclusive) - NN entrance channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 767 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
767 {
768
770 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
771
772 if (iso != 0) {
774 }
775 else {
777 }
778 }
virtual G4double NNToNNEtaExcluIso(const G4double ener, const G4int iso)
Isotopic Cross section for Eta production (exclusive) - NN entrance channel.
Referenced by NNToNDeltaEta(), NNToNNEtaFourPi(), and NNToNNEtaxPi().
◆ NNToNNEtaExcluIso()
|
protectedvirtual |
Isotopic Cross section for Eta production (exclusive) - NN entrance channel.
Definition at line 704 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
704 {
705
710
711 if(Ecm>=3.875) { // By hand (JCD)
712 sNNEta = -13.008*Ecm*Ecm + 84.531*Ecm + 36.234;
713 }
714 else if(Ecm>=2.725) { // By hand (JCD)
715 sNNEta = -913.2809*std::pow(Ecm,5) + 15564.27*std::pow(Ecm,4) - 105054.9*std::pow(Ecm,3) + 351294.2*std::pow(Ecm,2) - 582413.9*Ecm + 383474.7;
716 }
717 else if(Ecm>=2.575) { // By hand (JCD)
718 sNNEta = -2640.3*Ecm*Ecm + 14692*Ecm - 20225;
719 }
720 else {
721 sNNEta = -147043.497285*std::pow(Ecm,4) + 1487222.5438123*std::pow(Ecm,3) - 5634399.900744*std::pow(Ecm,2) + 9477290.199378*Ecm - 5972174.353438;
722 }
723
727 G4double Thr0=0.;
728 if (iso > 0) {
729 Thr0=2.*Mp+Meta;
730 }
731 else if (iso < 0) {
732 Thr0=2.*Mn+Meta;
733 }
734 else {
735 Thr0=Mn+Mp+Meta;
736 }
737
738 if (sNNEta < 1.e-9 || Ecm < Thr0) sNNEta = 0.; // Thr0: Ecm threshold
739
740 if (iso != 0) {
741 return sNNEta/1000.; // parameterization in microbarn (not millibarn)!
742 }
743
744 if(Ecm>=3.9) {
745 sNNEta1 = sNNEta;
746 }
747 else if (Ecm >= 3.5) {
748 sNNEta1 = -1916.2*Ecm*Ecm*Ecm + 21556.0*Ecm*Ecm - 80828.0*Ecm + 101200.0;
749 }
750 else if (Ecm >= 2.525) {
751 sNNEta1 = -4433.586*Ecm*Ecm*Ecm*Ecm + 56581.54*Ecm*Ecm*Ecm - 270212.6*Ecm*Ecm + 571650.6*Ecm - 451091.6;
752 }
753 else {
754 sNNEta1 = 17570.217219*Ecm*Ecm - 84910.985402*Ecm + 102585.55847;
755 }
756
757 sNNEta2 = -10220.89518466*Ecm*Ecm+51227.30841724*Ecm-64097.96025731;
758 if (sNNEta2 < 0.) sNNEta2=0.;
759
760 sNNEta = 2*(sNNEta1+sNNEta2)-sNNEta;
761 if (sNNEta < 1.e-9 || Ecm < Thr0) sNNEta = 0.; // Thr0: Ecm threshold
762
763 return sNNEta/1000.; // parameterization in microbarn (not millibarn)!
764
765 }
G4double getRealMass(const G4INCL::ParticleType t)
Get particle mass (in MeV/c^2)
Definition: G4INCLParticleTable.cc:710
Referenced by NNToNNEtaExclu(), and NNToNNEtaIso().
◆ NNToNNEtaFourPi()
|
protectedvirtual |
Cross section for direct 4-pion production - NNEta channel.
Definition at line 1039 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1039 {
1040 //
1041 // Nucleon-Nucleon producing one eta and four pions
1042 //
1043
1044 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
1045 if (ener < 2018.563) return 0.;
1047 const G4int i = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1048 G4double xsinelas;
1049 if (i!=0)
1050 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1051 else
1052 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1053 if (xsinelas <= 1.e-9) return 0.;
1054 G4double ratio=(NNToNNEta(particle1, particle2)-NNToNNEtaExclu(particle1, particle2))/xsinelas;
1055 if(s<6.25E6)
1056 return 0.;
1057 const G4double sigma = NNToNNEta(particle1, particle2) - NNToNNEtaExclu(particle1, particle2) - ratio*(NNToNNEtaOnePiOrDelta(particle1, particle2) + NNToNNEtaTwoPi(particle1, particle2) + NNToNNEtaThreePi(particle1, particle2));
1058 return ((sigma>1.e-9) ? sigma : 0.);
1059 }
virtual G4double NNToNNEtaThreePi(Particle const *const part1, Particle const *const part2)
Cross section for direct 3-pion production - NNEta channel.
virtual G4double NNToNNEtaTwoPi(Particle const *const part1, Particle const *const part2)
Cross section for direct 2-pion production - NNEta channel.
Referenced by NNToNNEtaxPi().
◆ NNToNNEtaIso()
|
protectedvirtual |
Cross section for One (more) pion production - piN entrance channel.
Cross section for Two (more) pion production - piN entrance channel
Cross section for Three (more) pion production - piN entrance channel virtual G4double piNThreePi(Particle const * const p1, Particle const * const p2);
Isotopic Cross section for Eta production (inclusive) - NN entrance channel
Definition at line 638 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
638 {
639
644 G4double x=Ecm*Ecm/5.88;
645
646//jcd
647 if (Ecm >= 3.05) {
648 sNNEta = 2.5*std::pow((x-1.),1.47)*std::pow(x,-1.25)*1000.;
649 }
650 else if(Ecm >= 2.6) {
651 sNNEta = -327.29*Ecm*Ecm*Ecm + 2870.*Ecm*Ecm - 7229.3*Ecm + 5273.3;
653 }
654 else {
655 sNNEta = NNToNNEtaExcluIso(ener, 2)*1000.;
656 }
657//jcd
658 if (sNNEta < 1.e-9) sNNEta = 0.;
659
660 if (iso != 0) {
661 return sNNEta/1000.; // parameterization in microbarn (not millibarn)!
662 }
663
664 if(Ecm >= 6.25) {
665 sNNEta1 = sNNEta;
666 }
667 else if (Ecm >= 2.6) {
668 sNNEta1 = sNNEta*std::exp(-(-5.53151576/Ecm+0.8850425));
669 }
670 else if (Ecm >= 2.525) { // = exclusive pn
671 sNNEta1 = -4433.586*Ecm*Ecm*Ecm*Ecm + 56581.54*Ecm*Ecm*Ecm - 270212.6*Ecm*Ecm + 571650.6*Ecm - 451091.6;
672 }
673 else { // = exclusive pn
674 sNNEta1 = 17570.217219*Ecm*Ecm - 84910.985402*Ecm + 102585.55847;
675 }
676
677 sNNEta2 = -10220.89518466*Ecm*Ecm+51227.30841724*Ecm-64097.96025731;
678 if (sNNEta2 < 0.) sNNEta2=0.;
679
680 sNNEta = 2*(sNNEta1+sNNEta2)-sNNEta;
681
685 if (sNNEta < 1.e-9 || Ecm < Mn+Mp+Meta) sNNEta = 0.;
686
687 return sNNEta/1000.; // parameterization in microbarn (not millibarn)!
688 }
Referenced by NNToNNEta().
◆ NNToNNEtaOnePi()
|
protectedvirtual |
Cross section for direct 1-pion production - NNEta channel.
Definition at line 967 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
967 {
968 // Cross section for nucleon-nucleon producing one eta and one pion
969
970 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
971 if (iso!=0)
972 return 0.;
973
974 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta (= 2705.55 - 2018.563; 4074595.287720512986=2018.563*2018.563)
975 if (ener < 2018.563) return 0.;
976
979
980 return 0.25*(CrossSectionsMultiPions::NNOnePiOrDelta(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNOnePiOrDelta(ener, 2, xsiso2));
981 }
virtual G4double NNOnePiOrDelta(const G4double ener, const G4int iso, const G4double xsiso)
Cross section for direct 1-pion production + delta production - NN entrance channel.
Definition: G4INCLCrossSectionsMultiPions.cc:248
Referenced by NNToNNEtaxPi().
◆ NNToNNEtaOnePiOrDelta()
|
protectedvirtual |
Cross section for direct 1-pion production - NNEta channel.
Definition at line 983 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
983 {
984 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
985 if (ener < 2018.563) return 0.;
986 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
987
989 if (iso != 0)
991 else {
993 return 0.5*(CrossSectionsMultiPions::NNOnePiOrDelta(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNOnePiOrDelta(ener, 2, xsiso2));
994 }
995 }
Referenced by NNToNDeltaEta(), and NNToNNEtaFourPi().
◆ NNToNNEtaThreePi()
|
protectedvirtual |
Cross section for direct 3-pion production - NNEta channel.
Definition at line 1016 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1016 {
1017 //
1018 // Nucleon-Nucleon producing one eta and three pions
1019 //
1020
1021 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
1022 if (ener < 2018.563) return 0.;
1023 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1024
1025
1029 if (iso != 0)
1031 else {
1035 return 0.5*(CrossSectionsMultiPions::NNThreePi(ener, 0, xsiso0, xs1pi0, xs2pi0)+ CrossSectionsMultiPions::NNThreePi(ener, 2, xsiso2, xs1pi2, xs2pi2));
1036 }
1037 }
virtual G4double NNTwoPi(const G4double ener, const G4int iso, const G4double xsiso)
Cross section for direct 2-pion production - NN entrance channel.
Definition: G4INCLCrossSectionsMultiPions.cc:331
virtual G4double NNThreePi(const G4double ener, const G4int iso, const G4double xsiso, const G4double xs1pi, const G4double xs2pi)
Cross section for direct 3-pion production - NN entrance channel.
Definition: G4INCLCrossSectionsMultiPions.cc:423
Referenced by NNToNNEtaFourPi(), and NNToNNEtaxPi().
◆ NNToNNEtaTwoPi()
|
protectedvirtual |
Cross section for direct 2-pion production - NNEta channel.
Definition at line 997 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
997 {
998 //
999 // Nucleon-Nucleon producing one eta and two pions
1000 //
1001 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
1002 if (ener < 2018.563) return 0.;
1003 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1004
1005
1007 if (iso != 0) {
1009 }
1010 else {
1012 return 0.5*(CrossSectionsMultiPions::NNTwoPi(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNTwoPi(ener, 2, xsiso2));
1013 }
1014 }
Referenced by NNToNNEtaFourPi(), and NNToNNEtaxPi().
◆ NNToNNEtaxPi()
|
virtual |
Cross section for X pion production - NNEta Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 1061 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1061 {
1062 //
1063 // Nucleon-Nucleon producing one eta and xpi pions
1064 //
1065// assert(xpi>0 && xpi<=nMaxPiNN);
1066// assert(particle1->isNucleon() && particle2->isNucleon());
1067
1068 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 581.437; // 581.437 MeV translation to open pion production in NNEta
1069 if (ener < 2018.563) return 0.;
1070 const G4int i = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1071 G4double xsinelas;
1072 if (i!=0)
1073 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1074 else
1075 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1076 if (xsinelas <= 1.e-9) return 0.;
1077 G4double ratio=(NNToNNEta(particle1, particle2)-NNToNNEtaExclu(particle1, particle2))/xsinelas;
1078
1079 if (xpi == 1)
1081 else if (xpi == 2)
1083 else if (xpi == 3)
1085 else if (xpi == 4)
1087 else // should never reach this point
1088 return 0.;
1089 }
virtual G4double NNToNNEtaFourPi(Particle const *const part1, Particle const *const part2)
Cross section for direct 4-pion production - NNEta channel.
virtual G4double NNToNNEtaOnePi(Particle const *const part1, Particle const *const part2)
Cross section for direct 1-pion production - NNEta channel.
◆ NNToNNOmega()
|
virtual |
Cross section for Omega production (inclusive) - NN entrance channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 815 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
815 {
816
818 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
819//jcd to be removed
820// return 0.;
821//jcd
822 if (iso != 0) {
824 }
825 else {
827 }
828 }
virtual G4double NNToNNOmegaIso(const G4double ener, const G4int iso)
Isotopic Cross section for Omega production (inclusive) - NN entrance channel.
Referenced by NNToNDeltaOmega(), NNToNNOmegaFourPi(), NNToNNOmegaxPi(), NNToxPiNN(), and G4INCL::CrossSectionsStrangeness::NNToxPiNN().
◆ NNToNNOmegaExclu()
|
virtual |
Cross section for Omega production (exclusive) - NN entrance channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 875 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
875 {
876//jcd to be removed
877// return 0.;
878//jcd
879
881 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
882
883 if (iso != 0) {
885 }
886 else {
888 }
889 }
virtual G4double NNToNNOmegaExcluIso(const G4double ener, const G4int iso)
Isotopic Cross section for Omega production (exclusive) - NN entrance channel.
Referenced by NNToNDeltaOmega(), NNToNNOmegaFourPi(), and NNToNNOmegaxPi().
◆ NNToNNOmegaExcluIso()
|
protectedvirtual |
Isotopic Cross section for Omega production (exclusive) - NN entrance channel.
Definition at line 830 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
830 {
831
835 G4double sthroot=std::sqrt(7.06);
836
837 if(Ecm>=3.0744) { // By hand (JCD)
838 sNNOmega = 330.*(Ecm-sthroot)/(1.05+std::pow((Ecm-sthroot),2));
839 }
840 else if(Ecm>=2.65854){
841 sNNOmega = -1208.09757*std::pow(Ecm,3) + 10773.3322*std::pow(Ecm,2) - 31661.0223*Ecm + 30728.7241 ;
842 }
843 else {
844 sNNOmega = 0. ;
845 }
846
850 G4double Thr0=0.;
851 if (iso > 0) {
852 Thr0=2.*Mp+Momega;
853 }
854 else if (iso < 0) {
855 Thr0=2.*Mn+Momega;
856 }
857 else {
858 Thr0=Mn+Mp+Momega;
859 }
860
861 if (sNNOmega < 1.e-9 || Ecm < Thr0) sNNOmega = 0.; // Thr0: Ecm threshold
862
863 if (iso != 0) {
864 return sNNOmega/1000.; // parameterization in microbarn (not millibarn)!
865 }
866
867 sNNOmega1 = 3*sNNOmega; // 3.0: ratio pn/pp
868
869 sNNOmega = 2*sNNOmega1-sNNOmega;
870 if (sNNOmega < 1.e-9 || Ecm < Thr0) sNNOmega = 0.;
871
872 return sNNOmega/1000.; // parameterization in microbarn (not millibarn)!
873 }
Referenced by NNToNNOmegaExclu(), and NNToNNOmegaIso().
◆ NNToNNOmegaFourPi()
|
protectedvirtual |
Cross section for direct 4-pion production - NNOmega channel.
Definition at line 1182 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1182 {
1183 //
1184 // Nucleon-Nucleon producing one omega and four pions
1185 //
1186//jcd to be removed
1187// return 0.;
1188//jcd
1189
1190 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1191 if (ener < 2018.563) return 0.;
1193 const G4int i = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1194 G4double xsinelas;
1195 if (i!=0)
1196 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1197 else
1198 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1199 if (xsinelas <= 1.e-9) return 0.;
1200 G4double ratio=(NNToNNOmega(particle1, particle2)-NNToNNOmegaExclu(particle1, particle2))/xsinelas;
1201 if(s<6.25E6)
1202 return 0.;
1203 const G4double sigma = NNToNNOmega(particle1, particle2) - NNToNNOmegaExclu(particle1, particle2) - ratio*(NNToNNOmegaOnePiOrDelta(particle1, particle2) + NNToNNOmegaTwoPi(particle1, particle2) + NNToNNOmegaThreePi(particle1, particle2));
1204 return ((sigma>1.e-9) ? sigma : 0.);
1205 }
virtual G4double NNToNNOmegaTwoPi(Particle const *const part1, Particle const *const part2)
Cross section for direct 2-pion production - NNOmega channel.
virtual G4double NNToNNOmegaThreePi(Particle const *const part1, Particle const *const part2)
Cross section for direct 3-pion production - NNOmega channel.
Referenced by NNToNNOmegaxPi().
◆ NNToNNOmegaIso()
|
protectedvirtual |
Isotopic Cross section for Omega production (inclusive) - NN entrance channel.
Definition at line 781 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
781 {
782
786 G4double x=Ecm*Ecm/7.06;
787
788 if(Ecm>4.0) {
789 sNNOmega = 2.5*std::pow(x-1, 1.47)*std::pow(x, -1.11) ;
790 }
791 else if(Ecm>2.802) { // 2802 MeV -> threshold to open inclusive (based on multipion threshold and omega mass)
792 sNNOmega = (568.5254*Ecm*Ecm - 2694.045*Ecm + 3106.247)/1000.;
794 }
795 else {
796 sNNOmega = NNToNNOmegaExcluIso(ener, 2);
797 }
798
799 if (sNNOmega < 1.e-9) sNNOmega = 0.;
800
801 if (iso != 0) {
802 return sNNOmega;
803 }
804
805 sNNOmega1 = 3.*sNNOmega; // 3.0: ratio pn/pp (5 from theory; 2 from experiments)
806
807 sNNOmega = 2.*sNNOmega1-sNNOmega;
808
809 if (sNNOmega < 1.e-9) sNNOmega = 0.;
810
811 return sNNOmega;
812 }
Referenced by NNToNNOmega().
◆ NNToNNOmegaOnePi()
|
protectedvirtual |
Cross section for direct 1-pion production - NNOmega channel.
Definition at line 1111 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1111 {
1112 // Cross section for nucleon-nucleon producing one omega and one pion
1113
1114 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1115 if (iso!=0)
1116 return 0.;
1117
1118 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega (= 2802. - 2018.563; 4074595.287720512986=2018.563*2018.563)
1119 if (ener < 2018.563) return 0.;
1120
1123
1124 return 0.25*(CrossSectionsMultiPions::NNOnePiOrDelta(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNOnePiOrDelta(ener, 2, xsiso2));
1125 }
Referenced by NNToNNOmegaxPi().
◆ NNToNNOmegaOnePiOrDelta()
|
protectedvirtual |
Cross section for direct 1-pion production - NNOmega channel.
Definition at line 1127 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1127 {
1128 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1129 if (ener < 2018.563) return 0.;
1130 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1131
1133 if (iso != 0)
1135 else {
1137 return 0.5*(CrossSectionsMultiPions::NNOnePiOrDelta(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNOnePiOrDelta(ener, 2, xsiso2));
1138 }
1139 }
Referenced by NNToNDeltaOmega(), and NNToNNOmegaFourPi().
◆ NNToNNOmegaThreePi()
|
protectedvirtual |
Cross section for direct 3-pion production - NNOmega channel.
Definition at line 1159 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1159 {
1160 //
1161 // Nucleon-Nucleon producing one omega and three pions
1162 //
1163
1164 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1165 if (ener < 2018.563) return 0.;
1166 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1167
1168
1172 if (iso != 0)
1174 else {
1178 return 0.5*(CrossSectionsMultiPions::NNThreePi(ener, 0, xsiso0, xs1pi0, xs2pi0)+ CrossSectionsMultiPions::NNThreePi(ener, 2, xsiso2, xs1pi2, xs2pi2));
1179 }
1180 }
Referenced by NNToNNOmegaFourPi(), and NNToNNOmegaxPi().
◆ NNToNNOmegaTwoPi()
|
protectedvirtual |
Cross section for direct 2-pion production - NNOmega channel.
Definition at line 1141 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1141 {
1142 //
1143 // Nucleon-Nucleon producing one omega and two pions
1144 //
1145 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1146 if (ener < 2018.563) return 0.;
1147 const G4int iso=ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1148
1150 if (iso != 0) {
1152 }
1153 else {
1155 return 0.5*(CrossSectionsMultiPions::NNTwoPi(ener, 0, xsiso0)+ CrossSectionsMultiPions::NNTwoPi(ener, 2, xsiso2));
1156 }
1157 }
Referenced by NNToNNOmegaFourPi(), and NNToNNOmegaxPi().
◆ NNToNNOmegaxPi()
|
virtual |
Cross section for X pion production - NNOmega Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 1207 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
1207 {
1208 //
1209 // Nucleon-Nucleon producing one omega and xpi pions
1210 //
1211// assert(xpi>0 && xpi<=nMaxPiNN);
1212// assert(particle1->isNucleon() && particle2->isNucleon());
1213//jcd to be removed
1214// return 0.;
1215//jcd
1216
1217 const G4double ener=KinematicsUtils::totalEnergyInCM(particle1, particle2) - 783.437; // 783.437 MeV translation to open pion production in NNOmega
1218 if (ener < 2018.563) return 0.;
1219 const G4int i = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
1220 G4double xsinelas;
1221 if (i!=0)
1222 xsinelas=CrossSectionsMultiPions::NNInelasticIso(ener, 2);
1223 else
1224 xsinelas=0.5*(CrossSectionsMultiPions::NNInelasticIso(ener, 0)+CrossSectionsMultiPions::NNInelasticIso(ener, 2));
1225 if (xsinelas <= 1.e-9) return 0.;
1226 G4double ratio=(NNToNNOmega(particle1, particle2)-NNToNNOmegaExclu(particle1, particle2))/xsinelas;
1227
1228 if (xpi == 1)
1230 else if (xpi == 2)
1232 else if (xpi == 3)
1234 else if (xpi == 4)
1236 else // should never reach this point
1237 return 0.;
1238 }
virtual G4double NNToNNOmegaFourPi(Particle const *const part1, Particle const *const part2)
Cross section for direct 4-pion production - NNOmega channel.
virtual G4double NNToNNOmegaOnePi(Particle const *const part1, Particle const *const part2)
Cross section for direct 1-pion production - NNOmega channel.
◆ NNToxPiNN()
|
virtual |
Cross section for X pion production - NN Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Reimplemented in G4INCL::CrossSectionsStrangeness.
Definition at line 892 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
892 {
893 //
894 // Nucleon-Nucleon producing xpi pions cross sections
895 //
896// assert(xpi>0 && xpi<=nMaxPiNN);
897// assert(particle1->isNucleon() && particle2->isNucleon());
898
904 G4double newXS1Pi=0.;
905 G4double newXS2Pi=0.;
906 G4double newXS3Pi=0.;
907 G4double newXS4Pi=0.;
908
909 if (xpi == 1) {
910 if (oldXS4Pi != 0. || oldXS3Pi != 0.)
911 newXS1Pi=oldXS1Pi;
912 else if (oldXS2Pi != 0.) {
913 newXS2Pi=oldXS2Pi-xsEtaOmega;
914 if (newXS2Pi < 0.)
915 newXS1Pi=oldXS1Pi-(xsEtaOmega-oldXS2Pi);
916 else
917 newXS1Pi=oldXS1Pi;
918 }
919 else
920 newXS1Pi=oldXS1Pi-xsEtaOmega;
921 return newXS1Pi;
922 }
923 else if (xpi == 2) {
924 if (oldXS4Pi != 0.)
925 newXS2Pi=oldXS2Pi;
926 else if (oldXS3Pi != 0.) {
927 newXS3Pi=oldXS3Pi-xsEtaOmega;
928 if (newXS3Pi < 0.)
929 newXS2Pi=oldXS2Pi-(xsEtaOmega-oldXS3Pi);
930 else
931 newXS2Pi=oldXS2Pi;
932 }
933 else {
934 newXS2Pi=oldXS2Pi-xsEtaOmega;
935 if (newXS2Pi < 0.)
936 newXS2Pi=0.;
937 }
938 return newXS2Pi;
939 }
940 else if (xpi == 3) {
941 if (oldXS4Pi != 0.) {
942 newXS4Pi=oldXS4Pi-xsEtaOmega;
943 if (newXS4Pi < 0.)
944 newXS3Pi=oldXS3Pi-(xsEtaOmega-oldXS4Pi);
945 else
946 newXS3Pi=oldXS3Pi;
947 }
948 else {
949 newXS3Pi=oldXS3Pi-xsEtaOmega;
950 if (newXS3Pi < 0.)
951 newXS3Pi=0.;
952 }
953 return newXS3Pi;
954 }
955 else if (xpi == 4) {
956 newXS4Pi=oldXS4Pi-xsEtaOmega;
957 if (newXS4Pi < 0.)
958 newXS4Pi=0.;
959 return newXS4Pi;
960 }
961
962 else // should never reach this point
963 return 0.;
964 }
virtual G4double NNToxPiNN(const G4int xpi, Particle const *const p1, Particle const *const p2)
Cross section for X pion production - NN Channel.
Definition: G4INCLCrossSectionsMultiPions.cc:532
◆ omegaNElastic()
|
protectedvirtual |
Definition at line 405 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
405 {
406 //
407 // Omega-Nucleon elastic cross sections
408 //
409// assert((particle1->isNucleon() && particle2->isOmega()) || (particle1->isOmega() && particle2->isNucleon()));
410
411 G4double sigma=0.;
412
413 const Particle *omega;
415
416 if (particle1->isOmega()) {
417 omega = particle1;
418 nucleon = particle2;
419 }
420 else {
421 omega = particle2;
422 nucleon = particle1;
423 }
424
426
427 sigma = 5.4 + 10.*std::exp(-0.6*pLab); // Eq.(21) in G.I. Lykasov et al., EPJA 6, 71-81 (1999)
428
429 return sigma;
430 }
Referenced by elastic().
◆ omegaNInelastic()
|
protectedvirtual |
Cross sections for mesonic resonance absorption on nucleon - inelastic Channel.
Definition at line 377 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
377 {
378 //
379 // Omega-Nucleon inelastic cross sections
380 //
381// assert((particle1->isNucleon() && particle2->isOmega()) || (particle1->isOmega() && particle2->isNucleon()));
382
383 G4double sigma=0.;
384
385 const Particle *omega;
387
388 if (particle1->isOmega()) {
389 omega = particle1;
390 nucleon = particle2;
391 }
392 else {
393 omega = particle2;
394 nucleon = particle1;
395 }
396
398
399 sigma = 20. + 4.0/pLab; // Eq.(24) in G.I. Lykasov et al., EPJA 6, 71-81 (1999)
400
401 return sigma;
402 }
Referenced by omegaNToPiN(), omegaNToPiPiN(), total(), and G4INCL::CrossSectionsStrangeness::total().
◆ omegaNToPiN()
|
virtual |
Cross section for OmegaN->PiN.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 433 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
433 {
434 //
435 // Omega-Nucleon producing Pion cross sections
436 //
437// assert((particle1->isNucleon() && particle2->isOmega()) || (particle1->isOmega() && particle2->isNucleon()));
438
440
444
445 G4double massomega;
446 G4double massnucleon;
447 G4double pCM_omega;
448 G4double pLab_omega;
449
450 G4double sigma=0.;
451
452 if (particle1->isOmega()) {
453 massomega=particle1->getMass();
454 massnucleon=particle2->getMass();
455 }
456 else {
457 massomega=particle2->getMass();
458 massnucleon=particle1->getMass();
459 }
460 pCM_omega=KinematicsUtils::momentumInCM(ECM, massomega, massnucleon);
461 pLab_omega=KinematicsUtils::momentumInLab(ECM*ECM, massomega, massnucleon);
462
464 G4double pCM_PiMinus=KinematicsUtils::momentumInCM(ECM, massPiMinus, massProton); // = pCM_PiPlus (because massPiMinus = massPiPlus)
465
466 sigma = (piMinuspToOmegaN(ECM)/2.) * std::pow((pCM_PiZero/pCM_omega), 2)
467 + piMinuspToOmegaN(ECM) * std::pow((pCM_PiMinus/pCM_omega), 2);
468
470// if (sigma > omegaNInelastic(particle1, particle2)) {
471 sigma = omegaNInelastic(particle1, particle2);
472 }
473
474 return sigma;
475 }
virtual G4double omegaNInelastic(Particle const *const p1, Particle const *const p2)
Cross sections for mesonic resonance absorption on nucleon - inelastic Channel.
G4double piMinuspToOmegaN(Particle const *const p1, Particle const *const p2)
Referenced by omegaNToPiPiN().
◆ omegaNToPiPiN()
|
protectedvirtual |
Cross sections for omega-induced 2Pi emission on nucleon.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 478 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
478 {
479 //
480 // Omega-Nucleon producing 2 PionS cross sections
481 //
482// assert((particle1->isNucleon() && particle2->isOmega()) || (particle1->isOmega() && particle2->isNucleon()));
483
484 G4double sigma=0.;
485
487
488 return sigma;
489 }
virtual G4double omegaNToPiN(Particle const *const p1, Particle const *const p2)
Cross section for OmegaN->PiN.
◆ piMinuspToEtaN() [1/2]
|
protected |
Definition at line 550 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
550 {
551 //
552 // Pion-Nucleon producing Eta cross sections
553 //
554
557
559
560 G4double sigma;
561
562// new parameterization (JCD) - end of february 2016
563 if (ECM < 1486.5) sigma=0.;
564 else
565 {
566 if (ECM < 1535.)
567 {
568 sigma = -0.0000003689197974814*std::pow(ECM,4) + 0.002260193900097*std::pow(ECM,3) - 5.193105877187*std::pow(ECM,2) + 5303.505273919*ECM - 2031265.900648;
569 }
570 else if (ECM < 1670.)
571 {
572 sigma = -0.0000000337986446*std::pow(ECM,4) + 0.000218279989*std::pow(ECM,3) - 0.528276144*std::pow(ECM,2) + 567.828367*ECM - 228709.42;
573 }
574 else if (ECM < 1714.)
575 {
576 sigma = 0.000003737765*std::pow(ECM,2) - 0.005664062*ECM;
577 }
578 else sigma=1.47*std::pow(plab, -1.68);
579 }
580
581 return sigma;
582 }
◆ piMinuspToEtaN() [2/2]
|
protected |
Internal function for pion cross sections.
Definition at line 505 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
505 {
506 //
507 // Pion-Nucleon producing Eta cross sections
508 //
509// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
510
511 G4double masspion;
512 G4double massnucleon;
513 if (particle1->isPion()) {
514 masspion=particle1->getMass();
515 massnucleon=particle2->getMass();
516 }
517 else {
518 masspion=particle2->getMass();
519 massnucleon=particle1->getMass();
520 }
521
524
525 G4double sigma;
526
527// new parameterization (JCD) - end of february 2016
528 if (ECM < 1486.5) sigma=0.;
529 else
530 {
531 if (ECM < 1535.)
532 {
533 sigma = -0.0000003689197974814*std::pow(ECM,4) + 0.002260193900097*std::pow(ECM,3) - 5.193105877187*std::pow(ECM,2) + 5303.505273919*ECM - 2031265.900648;
534 }
535 else if (ECM < 1670.)
536 {
537 sigma = -0.0000000337986446*std::pow(ECM,4) + 0.000218279989*std::pow(ECM,3) - 0.528276144*std::pow(ECM,2) + 567.828367*ECM - 228709.42;
538 }
539 else if (ECM < 1714.)
540 {
541 sigma = 0.000003737765*std::pow(ECM,2) - 0.005664062*ECM;
542 }
543 else sigma=1.47*std::pow(plab, -1.68);
544 }
545//
546
547 return sigma;
548 }
Referenced by etaNToPiN(), and piNToEtaN().
◆ piMinuspToOmegaN() [1/2]
|
protected |
Definition at line 615 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
615 {
616 //
617 // Pion-Nucleon producing Omega cross sections
618 //
619//jcd to be removed
620// return 0.;
621//jcd
622
623// G4double param=1.095 ; // Deneye (Thesis)
625
628
630
631 G4double sigma;
632 if (plab < param) sigma=0.;
633 else sigma=13.76*(plab-param)/(std::pow(plab, 3.33)-1.07);
634
635 return sigma;
636 }
◆ piMinuspToOmegaN() [2/2]
|
protected |
Definition at line 584 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
584 {
585 //
586 // Pion-Nucleon producing Omega cross sections
587 //
588// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
589//jcd to be removed
590// return 0.;
591//jcd
592
593// G4double param=1.095 ; // Deneye (Thesis)
595
596 G4double masspion;
597 G4double massnucleon;
598 if (particle1->isPion()) {
599 masspion=particle1->getMass();
600 massnucleon=particle2->getMass();
601 }
602 else {
603 masspion=particle2->getMass();
604 massnucleon=particle1->getMass();
605 }
608
609 G4double sigma;
610 if (plab < param) sigma=0.;
611 else sigma=13.76*(plab-param)/(std::pow(plab, 3.33) - 1.07); // Phys. Rev. C 41, 1701–1718 (1990)
612
613 return sigma;
614}
Referenced by omegaNToPiN(), and piNToOmegaN().
◆ piNToEtaN()
|
virtual |
Cross sections for mesonic resonance production - piN Channel.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 195 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
195 {
196 //
197 // Pion-Nucleon producing Eta cross sections
198 //
199// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
200
201 G4double sigma;
202 sigma=piMinuspToEtaN(particle1,particle2);
203
204 const G4int isoin = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
205
206 if (isoin == -1) {
208 else return 0.5 * sigma;
209 }
210 else if (isoin == 1) {
212 else return 0.5 * sigma;
213 }
214 else return 0. ; // should never return 0. (?) // pi+ p and pi- n return 0.
215
216// return sigma;
217 }
Referenced by piNToxPiN(), and G4INCL::CrossSectionsStrangeness::piNToxPiN().
◆ piNToEtaPrimeN()
|
virtual |
Cross section for PiN->EtaPrimeN.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 244 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
244 {
245#else
246 G4double CrossSectionsMultiPionsAndResonances::piNToEtaPrimeN(Particle const * const particle1, Particle const * const particle2) {
247#endif
248 //
249 // Pion-Nucleon producing EtaPrime cross sections
250 //
251// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
252
253 return 0.;
254 }
virtual G4double piNToEtaPrimeN(Particle const *const p1, Particle const *const p2)
Cross section for PiN->EtaPrimeN.
Referenced by piNToEtaPrimeN().
◆ piNToOmegaN()
|
virtual |
Cross section for PiN->OmegaN.
Reimplemented from G4INCL::CrossSectionsMultiPions.
Definition at line 219 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
219 {
220 //
221 // Pion-Nucleon producing Omega cross sections
222 //
223// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
224
225 G4double sigma;
226 sigma=piMinuspToOmegaN(particle1,particle2);
227
228 const G4int isoin = ParticleTable::getIsospin(particle1->getType()) + ParticleTable::getIsospin(particle2->getType());
229
230 if (isoin == -1) {
232 else return 0.5 * sigma;
233 }
234 else if (isoin == 1) {
236 else return 0.5 * sigma;
237 }
238 else return 0. ; // should never return 0. (?) // pi+ p and pi- n return 0.
239
240// return sigma;
241 }
Referenced by piNToxPiN(), and G4INCL::CrossSectionsStrangeness::piNToxPiN().
◆ piNToxPiN()
|
virtual |
Cross section for X pion production - piN Channel (modified due to the mesonic resonances)
Reimplemented from G4INCL::CrossSectionsMultiPions.
Reimplemented in G4INCL::CrossSectionsStrangeness.
Definition at line 137 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
137 {
138 //
139 // pion-Nucleon producing xpi pions cross sections (corrected due to eta and omega)
140 //
141// assert(xpi>1 && xpi<=nMaxPiPiN);
142// assert((particle1->isNucleon() && particle2->isPion()) || (particle1->isPion() && particle2->isNucleon()));
143
149 G4double newXS2Pi=0.;
150 G4double newXS3Pi=0.;
151 G4double newXS4Pi=0.;
152
153 if (xpi == 2) {
154 if (oldXS4Pi != 0.)
155 newXS2Pi=oldXS2Pi;
156 else if (oldXS3Pi != 0.) {
157 newXS3Pi=oldXS3Pi-xsEta-xsOmega;
158 if (newXS3Pi < 1.e-09)
159 newXS2Pi=oldXS2Pi-(xsEta+xsOmega-oldXS3Pi);
160 else
161 newXS2Pi=oldXS2Pi;
162 }
163 else {
164 newXS2Pi=oldXS2Pi-xsEta-xsOmega;
165 if (newXS2Pi < 1.e-09)
166 newXS2Pi=0.;
167 }
168 return newXS2Pi;
169 }
170 else if (xpi == 3) {
171 if (oldXS4Pi != 0.) {
172 newXS4Pi=oldXS4Pi-xsEta-xsOmega;
173 if (newXS4Pi < 1.e-09)
174 newXS3Pi=oldXS3Pi-(xsEta+xsOmega-oldXS4Pi);
175 else
176 newXS3Pi=oldXS3Pi;
177 }
178 else {
179 newXS3Pi=oldXS3Pi-xsEta-xsOmega;
180 if (newXS3Pi < 1.e-09)
181 newXS3Pi=0.;
182 }
183 return newXS3Pi;
184 }
185 else if (xpi == 4) {
186 newXS4Pi=oldXS4Pi-xsEta-xsOmega;
187 if (newXS4Pi < 1.e-09)
188 newXS4Pi=0.;
189 return newXS4Pi;
190 }
191 else // should never reach this point
192 return 0.;
193 }
virtual G4double piNToEtaN(Particle const *const p1, Particle const *const p2)
Cross sections for mesonic resonance production - piN Channel.
virtual G4double piNToOmegaN(Particle const *const p1, Particle const *const p2)
Cross section for PiN->OmegaN.
virtual G4double piNToxPiN(const G4int xpi, Particle const *const p1, Particle const *const p2)
Cross section for X pion production - piN Channel.
Definition: G4INCLCrossSectionsMultiPions.cc:849
◆ total()
|
virtual |
new total particle-particle cross section
Reimplemented from G4INCL::CrossSectionsMultiPions.
Reimplemented in G4INCL::CrossSectionsStrangeness.
Definition at line 91 of file G4INCLCrossSectionsMultiPionsAndResonances.cc.
91 {
92 G4double inelastic;
93 if(p1->isNucleon() && p2->isNucleon()) {
95 } else if((p1->isNucleon() && p2->isDelta()) ||
96 (p1->isDelta() && p2->isNucleon())) {
97 inelastic = CrossSectionsMultiPions::NDeltaToNN(p1, p2);
98 } else if((p1->isNucleon() && p2->isPion()) ||
99 (p1->isPion() && p2->isNucleon())) {
101 } else if((p1->isNucleon() && p2->isEta()) ||
102 (p1->isEta() && p2->isNucleon())) {
104 } else if((p1->isNucleon() && p2->isOmega()) ||
105 (p1->isOmega() && p2->isNucleon())) {
106 inelastic = omegaNInelastic(p1,p2);
107 } else if((p1->isNucleon() && p2->isEtaPrime()) ||
108 (p1->isEtaPrime() && p2->isNucleon())) {
109 inelastic = etaPrimeNToPiN(p1,p2);
110 } else {
111 inelastic = 0.;
112 }
113
115 }
virtual G4double etaNToPiPiN(Particle const *const p1, Particle const *const p2)
Cross sections for mesonic resonance absorption on nucleon - pipiN Channel.
virtual G4double elastic(Particle const *const p1, Particle const *const p2)
new elastic particle-particle cross section
G4double piNTot(Particle const *const p1, Particle const *const p2)
Definition: G4INCLCrossSectionsMultiPions.cc:712
G4double NNTot(Particle const *const part1, Particle const *const part2)
Internal implementation of the NN total cross section.
Definition: G4INCLCrossSectionsMultiPions.cc:164
virtual G4double NDeltaToNN(Particle const *const p1, Particle const *const p2)
Cross section for NDelta->NN.
Definition: G4INCLCrossSectionsMultiPions.cc:749
Member Data Documentation
◆ nMaxPiNN
|
staticprotected |
Maximum number of outgoing pions in NN collisions.
Definition at line 110 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ nMaxPiPiN
|
staticprotected |
Maximum number of outgoing pions in piN collisions.
Definition at line 113 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s01ppHC
|
protected |
Horner coefficients for s01pp.
Definition at line 118 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s01ppOOT
|
staticprotected |
One over threshold for s01pp.
Definition at line 139 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s01pzHC
|
protected |
Horner coefficients for s01pz.
Definition at line 120 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s01pzOOT
|
staticprotected |
One over threshold for s01pz.
Definition at line 141 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s02pmHC
|
protected |
Horner coefficients for s02pm.
Definition at line 132 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s02pmOOT
|
staticprotected |
One over threshold for s02pm.
Definition at line 153 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s02pzHC
|
protected |
Horner coefficients for s02pz.
Definition at line 130 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s02pzOOT
|
staticprotected |
One over threshold for s02pz.
Definition at line 151 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s11pmHC
|
protected |
Horner coefficients for s11pm.
Definition at line 122 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s11pmOOT
|
staticprotected |
One over threshold for s11pm.
Definition at line 143 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s11pzHC
|
protected |
Horner coefficients for s11pz.
Definition at line 116 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s11pzOOT
|
staticprotected |
One over threshold for s11pz.
Definition at line 137 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12mzHC
|
protected |
Horner coefficients for s12mz.
Definition at line 134 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12mzOOT
|
staticprotected |
One over threshold for s12mz.
Definition at line 155 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12pmHC
|
protected |
Horner coefficients for s12pm.
Definition at line 124 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12pmOOT
|
staticprotected |
One over threshold for s12pm.
Definition at line 145 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12ppHC
|
protected |
Horner coefficients for s12pp.
Definition at line 126 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12ppOOT
|
staticprotected |
One over threshold for s12pp.
Definition at line 147 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12zzHC
|
protected |
Horner coefficients for s12zz.
Definition at line 128 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
◆ s12zzOOT
|
staticprotected |
One over threshold for s12zz.
Definition at line 149 of file G4INCLCrossSectionsMultiPionsAndResonances.hh.
The documentation for this class was generated from the following files:
