60 fGlauberEnergy = 91.*
GeV;
62 fSAIDHighEnergyLimit = 2.6*
GeV;
66 for (
G4int i = 0; i < 93; i++) {
67 theGlauberFac[i] = 0.0;
68 theCoulombFac[i] = 0.0;
78 isInitialized =
false;
107 return (1 == Z && 2 >= A);
124 if(Z > 92) { Z = 92; }
126 if(ekin <= fLowEnergy) {
127 cross = theCoulombFac[
Z];
128 }
else if(ekin > fGlauberEnergy) {
135 G4cout <<
"G4BGGPionElasticXS::GetElementCrossSection for "
138 <<
" in nucleus Z= " << Z <<
" A= " << theA[
Z]
139 <<
" XS(b)= " << cross/
barn
157 if(ekin <= fSAIDHighEnergyLimit) {
179 G4cout <<
"G4BGGPionElasticXS::GetIsoCrossSection for "
182 <<
" in nucleus Z= " << Z <<
" A= " << A
183 <<
" XS(b)= " << cross/
barn
196 G4cout <<
"### G4BGGPionElasticXS WARNING: is not applicable to "
200 "G4BGGPionElasticXS::BuildPhysicsTable is used for wrong particle");
204 if(isInitialized) {
return; }
205 isInitialized =
true;
226 G4cout <<
"### G4BGGPionElasticXS::Initialise for "
229 for(
G4int iz=2; iz<93; iz++) {
237 theGlauberFac[iz] = csdn/csup;
239 G4cout <<
"Z= " << iz <<
" A= " << A
240 <<
" factor= " << theGlauberFac[iz] <<
G4endl;
255 for(
G4int iz=2; iz<93; iz++) {
258 G4cout <<
"Z= " << iz <<
" A= " << A
259 <<
" factor= " << theCoulombFac[iz] <<
G4endl;
267 outFile <<
"The Barashenkov-Glauber-Gribov cross section handles elastic\n"
268 <<
"scattering of pions from nuclei at all energies. The\n"
269 <<
"Barashenkov parameterization is used below 91 GeV and the\n"
270 <<
"Glauber-Gribov parameterization is used above 91 GeV.\n";
G4double GetElasticHadronNucleonXsc()
virtual ~G4BGGPionElasticXS()
G4double GetElasticCrossSection(const G4DynamicParticle *aParticle, G4int Z, G4int A)
G4double GetKineticEnergy() const
virtual G4double GetElasticIsotopeCrossSection(const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
G4ParticleDefinition * GetDefinition() const
static G4NistManager * Instance()
const G4String & GetParticleName() const
static constexpr double TeV
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
void SetMinKinEnergy(G4double value)
G4BGGPionElasticXS(const G4ParticleDefinition *)
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
static G4Proton * Proton()
static G4PionPlus * PionPlus()
void SetKineticEnergy(G4double aEnergy)
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
G4double GetElasticGlauberGribov(const G4DynamicParticle *, G4int Z, G4int A)
void SetMaxKinEnergy(G4double value)
virtual G4bool IsElementApplicable(const G4DynamicParticle *, G4int Z, const G4Material *)
static G4PionMinus * PionMinus()
G4double GetAtomicMassAmu(const G4String &symb) const
static constexpr double GeV
G4double GetHadronNucleonXscPDG(const G4DynamicParticle *, const G4ParticleDefinition *)
static constexpr double MeV
virtual G4bool IsIsoApplicable(const G4DynamicParticle *, G4int Z, G4int A, const G4Element *elm=0, const G4Material *mat=0)
virtual void BuildPhysicsTable(const G4ParticleDefinition &)
static constexpr double barn
virtual void CrossSectionDescription(std::ostream &) const
void BuildPhysicsTable(const G4ParticleDefinition &)