Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4INCL::CrossSectionsTruncatedMultiPions Class Reference

Truncated multipion cross sections. More...

#include <G4INCLCrossSectionsTruncatedMultiPions.hh>

Inheritance diagram for G4INCL::CrossSectionsTruncatedMultiPions:
Collaboration diagram for G4INCL::CrossSectionsTruncatedMultiPions:

Public Member Functions

 CrossSectionsTruncatedMultiPions (const G4int nPi=std::numeric_limits< G4int >::max())
 
virtual G4double elastic (Particle const *const p1, Particle const *const p2)
 Elastic particle-particle cross section. More...
 
virtual G4double piNToDelta (Particle const *const p1, Particle const *const p2)
 Cross section for Delta production - piN Channel. More...
 
virtual G4double piNToxPiN (const G4int xpi, Particle const *const p1, Particle const *const p2)
 Cross section for X pion production - piN Channel. More...
 
virtual G4double NNToxPiNN (const G4int xpi, Particle const *const p1, Particle const *const p2)
 Cross section for X pion production - NN Channel. More...
 
- Public Member Functions inherited from G4INCL::CrossSectionsMultiPions
 CrossSectionsMultiPions ()
 
virtual G4double total (Particle const *const p1, Particle const *const p2)
 Total (elastic+inelastic) particle-particle cross section. More...
 
virtual G4double NDeltaToNN (Particle const *const p1, Particle const *const p2)
 Cross section for NDelta->NN. More...
 
virtual G4double NNToNDelta (Particle const *const p1, Particle const *const p2)
 Cross section for Delta production - NN Channel. More...
 
virtual G4double calculateNNAngularSlope (G4double energyCM, G4int iso)
 Calculate the slope of the NN DDXS. More...
 
virtual G4double piNToEtaN (Particle const *const p1, Particle const *const p2)
 Cross sections for mesonic resonance production - piN Channel. More...
 
virtual G4double piNToOmegaN (Particle const *const p1, Particle const *const p2)
 Cross section for PiN->OmegaN. More...
 
virtual G4double piNToEtaPrimeN (Particle const *const p1, Particle const *const p2)
 Cross section for PiN->EtaPrimeN. More...
 
virtual G4double etaNToPiPiN (Particle const *const p1, Particle const *const p2)
 Cross sections for mesonic resonance absorption on nucleon - pipiN Channel. More...
 
virtual G4double omegaNToPiPiN (Particle const *const p1, Particle const *const p2)
 Cross section for OmegaN->PiPiN. More...
 
virtual G4double etaNToPiN (Particle const *const p1, Particle const *const p2)
 Cross sections for mesonic resonance absorption on nucleon - piN Channel. More...
 
virtual G4double omegaNToPiN (Particle const *const p1, Particle const *const p2)
 Cross section for OmegaN->PiN. More...
 
virtual G4double etaPrimeNToPiN (Particle const *const p1, Particle const *const p2)
 Cross section for EtaPrimeN->PiN. More...
 
virtual G4double NNToNNEta (Particle const *const particle1, Particle const *const particle2)
 Cross section for Eta production - NN entrance channel. More...
 
virtual G4double NNToNNEtaExclu (Particle const *const particle1, Particle const *const particle2)
 Cross section for Eta production (exclusive) - NN entrance channel. More...
 
virtual G4double NNToNNEtaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2)
 Cross section for X pion production - NNEta Channel. More...
 
virtual G4double NNToNDeltaEta (Particle const *const p1, Particle const *const p2)
 Cross section for N-Delta-Eta production - NNEta Channel. More...
 
virtual G4double NNToNNOmega (Particle const *const particle1, Particle const *const particle2)
 Cross section for Eta production - NN entrance channel. More...
 
virtual G4double NNToNNOmegaExclu (Particle const *const particle1, Particle const *const particle2)
 Cross section for Eta production (exclusive) - NN entrance channel. More...
 
virtual G4double NNToNNOmegaxPi (const G4int xpi, Particle const *const p1, Particle const *const p2)
 Cross section for X pion production - NNEta Channel. More...
 
virtual G4double NNToNDeltaOmega (Particle const *const p1, Particle const *const p2)
 Cross section for N-Delta-Eta production - NNEta Channel. More...
 
- Public Member Functions inherited from G4INCL::ICrossSections
 ICrossSections ()
 
virtual ~ICrossSections ()
 

Protected Attributes

const G4int nMaxPi
 Maximum number of pions produced in TruncatedMultiPion collisions. More...
 
- Protected Attributes inherited from G4INCL::CrossSectionsMultiPions
const HornerC7 s11pzHC
 Horner coefficients for s11pz. More...
 
const HornerC8 s01ppHC
 Horner coefficients for s01pp. More...
 
const HornerC4 s01pzHC
 Horner coefficients for s01pz. More...
 
const HornerC4 s11pmHC
 Horner coefficients for s11pm. More...
 
const HornerC5 s12pmHC
 Horner coefficients for s12pm. More...
 
const HornerC3 s12ppHC
 Horner coefficients for s12pp. More...
 
const HornerC4 s12zzHC
 Horner coefficients for s12zz. More...
 
const HornerC4 s02pzHC
 Horner coefficients for s02pz. More...
 
const HornerC6 s02pmHC
 Horner coefficients for s02pm. More...
 
const HornerC4 s12mzHC
 Horner coefficients for s12mz. More...
 

Additional Inherited Members

- 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. More...
 
G4double NNElasticFixed (const G4double s, const G4int i)
 Internal implementation of the NN elastic cross section with fixed isospin. More...
 
G4double NNTot (Particle const *const part1, Particle const *const part2)
 Internal implementation of the NN total cross section. More...
 
G4double NNTotFixed (const G4double s, const G4int i)
 Internal implementation of the NN total cross section with fixed isospin. More...
 
G4double NNInelasticIso (const G4double ener, const G4int iso)
 Internal implementation of the isospin dependent NN reaction cross section. More...
 
virtual G4double NNOnePiOrDelta (const G4double ener, const G4int iso, const G4double xsiso)
 Cross section for direct 1-pion production + delta production - NN entrance channel. More...
 
virtual G4double NNTwoPi (const G4double ener, const G4int iso, const G4double xsiso)
 Cross section for direct 2-pion production - NN entrance channel. More...
 
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. More...
 
virtual G4double NNOnePi (Particle const *const part1, Particle const *const part2)
 Cross section for direct 1-pion production - NN entrance channel. More...
 
virtual G4double NNOnePiOrDelta (Particle const *const part1, Particle const *const part2)
 Cross section for direct 1-pion production - NN entrance channel. More...
 
virtual G4double NNTwoPi (Particle const *const part1, Particle const *const part2)
 Cross section for direct 2-pion production - NN entrance channel. More...
 
virtual G4double NNThreePi (Particle const *const part1, Particle const *const part2)
 Cross section for direct 3-pion production - NN entrance channel. More...
 
virtual G4double NNFourPi (Particle const *const part1, Particle const *const part2)
 Cross section for direct 4-pion production - NN entrance channel. More...
 
G4double spnPiPlusPHE (const G4double x)
 Internal function for pion cross sections. More...
 
G4double spnPiMinusPHE (const G4double x)
 Internal function for pion cross sections. More...
 
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. More...
 
virtual G4double piNTwoPi (Particle const *const p1, Particle const *const p2)
 Cross section for Two (more) pion production - piN entrance channel. More...
 
- Static Protected Attributes inherited from G4INCL::CrossSectionsMultiPions
static const G4int nMaxPiNN = 4
 Maximum number of outgoing pions in NN collisions. More...
 
static const G4int nMaxPiPiN = 4
 Maximum number of outgoing pions in piN collisions. More...
 
static const G4double s11pzOOT = 0.0035761542037692665889
 One over threshold for s11pz. More...
 
static const G4double s01ppOOT = 0.003421025623481919853
 One over threshold for s01pp. More...
 
static const G4double s01pzOOT = 0.0035739814152966403123
 One over threshold for s01pz. More...
 
static const G4double s11pmOOT = 0.0034855350296270480281
 One over threshold for s11pm. More...
 
static const G4double s12pmOOT = 0.0016672224074691565119
 One over threshold for s12pm. More...
 
static const G4double s12ppOOT = 0.0016507643038726931312
 One over threshold for s12pp. More...
 
static const G4double s12zzOOT = 0.0011111111111111111111
 One over threshold for s12zz. More...
 
static const G4double s02pzOOT = 0.00125
 One over threshold for s02pz. More...
 
static const G4double s02pmOOT = 0.0016661112962345883443
 One over threshold for s02pm. More...
 
static const G4double s12mzOOT = 0.0017047391749062392793
 One over threshold for s12mz. More...
 

Detailed Description

Truncated multipion cross sections.

Definition at line 54 of file G4INCLCrossSectionsTruncatedMultiPions.hh.

Constructor & Destructor Documentation

G4INCL::CrossSectionsTruncatedMultiPions::CrossSectionsTruncatedMultiPions ( const G4int  nPi = std::numeric_limits<G4int>::max())

Definition at line 50 of file G4INCLCrossSectionsTruncatedMultiPions.cc.

50  :
51  nMaxPi(nPi)
52  {}
const G4int nMaxPi
Maximum number of pions produced in TruncatedMultiPion collisions.

Member Function Documentation

G4double G4INCL::CrossSectionsTruncatedMultiPions::elastic ( Particle const *const  p1,
Particle const *const  p2 
)
virtual

Elastic particle-particle cross section.

The elastic pi-N cross section is calculated by difference: elastic = total - inelastic(pion production) - delta, but we need to use the unmodified piNToDelta cross section, or else we will violate the cross-section sum.

This modification is only necessary for piN collisions.

Reimplemented from G4INCL::CrossSectionsMultiPions.

Definition at line 54 of file G4INCLCrossSectionsTruncatedMultiPions.cc.

54  {
55 // if(!p1->isPion() && !p2->isPion()) {
56  if((p1->isNucleon()||p1->isDelta()) && (p2->isNucleon()||p2->isDelta())){
57  return NNElastic(p1, p2);
58  }
59 // else if (p1->isNucleon() || p2->isNucleon()) {
60  else if ((p1->isNucleon() && p2->isPion()) || (p2->isNucleon() && p1->isPion())){
61  G4double pielas = piNTot(p1,p2) - piNIne(p1,p2) - CrossSectionsMultiPions::piNToDelta(p1,p2);
62  if (pielas < 0.)
63  pielas = 0.;
64  return pielas;
65  } else
66  return 0.0;
67  }
G4double piNIne(Particle const *const p1, Particle const *const p2)
virtual G4double piNToDelta(Particle const *const p1, Particle const *const p2)
Cross section for Delta production - piN Channel.
G4double NNElastic(Particle const *const part1, Particle const *const part2)
Internal implementation of the NN elastic cross section.
G4double piNTot(Particle const *const p1, Particle const *const p2)
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double G4INCL::CrossSectionsTruncatedMultiPions::NNToxPiNN ( const G4int  xpi,
Particle const *const  p1,
Particle const *const  p2 
)
virtual

Cross section for X pion production - NN Channel.

Reimplemented from G4INCL::CrossSectionsMultiPions.

Definition at line 78 of file G4INCLCrossSectionsTruncatedMultiPions.cc.

78  {
79  if(xpi<nMaxPi)
80  return CrossSectionsMultiPions::NNToxPiNN(xpi, particle1, particle2);
81  else if(xpi==nMaxPi) {
82  G4double sum = 0.;
83  for(G4int i=xpi; i<=nMaxPiNN; ++i)
84  sum += CrossSectionsMultiPions::NNToxPiNN(i, particle1, particle2);
85  return sum;
86  } else
87  return 0.;
88  }
virtual G4double NNToxPiNN(const G4int xpi, Particle const *const p1, Particle const *const p2)
Cross section for X pion production - NN Channel.
const G4int nMaxPi
Maximum number of pions produced in TruncatedMultiPion collisions.
int G4int
Definition: G4Types.hh:78
static const G4int nMaxPiNN
Maximum number of outgoing pions in NN collisions.
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double G4INCL::CrossSectionsTruncatedMultiPions::piNToDelta ( Particle const *const  p1,
Particle const *const  p2 
)
virtual

Cross section for Delta production - piN Channel.

Reimplemented from G4INCL::CrossSectionsMultiPions.

Definition at line 69 of file G4INCLCrossSectionsTruncatedMultiPions.cc.

69  {
70  G4double sum = CrossSectionsMultiPions::piNToDelta(particle1, particle2);
71  if(nMaxPi<=1) {
72  for(G4int i=nMaxPi+1; i<=nMaxPiPiN; ++i)
73  sum += CrossSectionsMultiPions::piNToxPiN(i, particle1, particle2);
74  }
75  return sum;
76  }
virtual G4double piNToDelta(Particle const *const p1, Particle const *const p2)
Cross section for Delta production - piN Channel.
const G4int nMaxPi
Maximum number of pions produced in TruncatedMultiPion collisions.
virtual G4double piNToxPiN(const G4int xpi, Particle const *const p1, Particle const *const p2)
Cross section for X pion production - piN Channel.
int G4int
Definition: G4Types.hh:78
double G4double
Definition: G4Types.hh:76
static const G4int nMaxPiPiN
Maximum number of outgoing pions in piN collisions.

Here is the call graph for this function:

G4double G4INCL::CrossSectionsTruncatedMultiPions::piNToxPiN ( const G4int  xpi,
Particle const *const  p1,
Particle const *const  p2 
)
virtual

Cross section for X pion production - piN Channel.

Reimplemented from G4INCL::CrossSectionsMultiPions.

Definition at line 90 of file G4INCLCrossSectionsTruncatedMultiPions.cc.

90  {
91  if(xpi<nMaxPi)
92  return CrossSectionsMultiPions::piNToxPiN(xpi, particle1, particle2);
93  else if(xpi==nMaxPi) {
94  G4double sum = 0.;
95  for(G4int i=xpi; i<=nMaxPiPiN; ++i)
96  sum += CrossSectionsMultiPions::piNToxPiN(i, particle1, particle2);
97  return sum;
98  } else
99  return 0.;
100  }
const G4int nMaxPi
Maximum number of pions produced in TruncatedMultiPion collisions.
virtual G4double piNToxPiN(const G4int xpi, Particle const *const p1, Particle const *const p2)
Cross section for X pion production - piN Channel.
int G4int
Definition: G4Types.hh:78
double G4double
Definition: G4Types.hh:76
static const G4int nMaxPiPiN
Maximum number of outgoing pions in piN collisions.

Here is the call graph for this function:

Member Data Documentation

const G4int G4INCL::CrossSectionsTruncatedMultiPions::nMaxPi
protected

Maximum number of pions produced in TruncatedMultiPion collisions.

Definition at line 81 of file G4INCLCrossSectionsTruncatedMultiPions.hh.


The documentation for this class was generated from the following files: