Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/cross_sections/src/G4BGGPionElasticXS.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/cross_sections/src/G4BGGPionElasticXS.cc (Version 11.3.0) and /processes/hadronic/cross_sections/src/G4BGGPionElasticXS.cc (Version 10.7.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 // -------------------------------------------     27 // -------------------------------------------------------------------
 28 //                                                 28 //
 29 // GEANT4 Class file                               29 // GEANT4 Class file
 30 //                                                 30 //
 31 //                                                 31 //
 32 // File name:     G4BGGPionElasticXS               32 // File name:     G4BGGPionElasticXS
 33 //                                                 33 //
 34 // Author:        Vladimir Ivanchenko              34 // Author:        Vladimir Ivanchenko
 35 //                                                 35 //
 36 // Creation date: 01.10.2003                       36 // Creation date: 01.10.2003
 37 // Modifications:                                  37 // Modifications:
 38 //                                                 38 //
 39 // -------------------------------------------     39 // -------------------------------------------------------------------
 40 //                                                 40 //
 41                                                    41 
 42 #include "G4BGGPionElasticXS.hh"                   42 #include "G4BGGPionElasticXS.hh"
 43 #include "G4SystemOfUnits.hh"                      43 #include "G4SystemOfUnits.hh"
 44 #include "G4ComponentGGHadronNucleusXsc.hh"        44 #include "G4ComponentGGHadronNucleusXsc.hh"
 45 #include "G4UPiNuclearCrossSection.hh"             45 #include "G4UPiNuclearCrossSection.hh"
 46 #include "G4HadronNucleonXsc.hh"                   46 #include "G4HadronNucleonXsc.hh"
 47 #include "G4NuclearRadii.hh"                       47 #include "G4NuclearRadii.hh"
 48                                                    48 
 49 #include "G4Proton.hh"                             49 #include "G4Proton.hh"
 50 #include "G4PionPlus.hh"                           50 #include "G4PionPlus.hh"
 51 #include "G4PionMinus.hh"                          51 #include "G4PionMinus.hh"
 52 #include "G4NistManager.hh"                        52 #include "G4NistManager.hh"
 53 #include "G4HadronicParameters.hh"                 53 #include "G4HadronicParameters.hh"
 54 #include "G4Pow.hh"                                54 #include "G4Pow.hh"
 55                                                    55 
 56 G4double G4BGGPionElasticXS::theGlauberFacPiPl     56 G4double G4BGGPionElasticXS::theGlauberFacPiPlus[93] = {0.0};
 57 G4double G4BGGPionElasticXS::theCoulombFacPiPl     57 G4double G4BGGPionElasticXS::theCoulombFacPiPlus[93] = {0.0};
 58 G4double G4BGGPionElasticXS::theGlauberFacPiMi     58 G4double G4BGGPionElasticXS::theGlauberFacPiMinus[93] = {0.0};
 59 G4double G4BGGPionElasticXS::theCoulombFacPiMi     59 G4double G4BGGPionElasticXS::theCoulombFacPiMinus[93] = {0.0};
 60 G4int G4BGGPionElasticXS::theA[93] = {0};      <<  60 G4int    G4BGGPionElasticXS::theA[93] = {0};
                                                   >>  61 
                                                   >>  62 #ifdef G4MULTITHREADED
                                                   >>  63 G4Mutex G4BGGPionElasticXS::pionElasticXSMutex = G4MUTEX_INITIALIZER;
                                                   >>  64 #endif
 61                                                    65 
 62 //....oooOO0OOooo........oooOO0OOooo........oo     66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 63                                                    67 
 64 G4BGGPionElasticXS::G4BGGPionElasticXS(const G     68 G4BGGPionElasticXS::G4BGGPionElasticXS(const G4ParticleDefinition* p) 
 65  : G4VCrossSectionDataSet("BarashenkovGlauberG     69  : G4VCrossSectionDataSet("BarashenkovGlauberGribov") 
 66 {                                                  70 {
 67   verboseLevel = 0;                                71   verboseLevel = 0;
 68   fGlauberEnergy = 91.*CLHEP::GeV;             <<  72   fGlauberEnergy = 91.*GeV;
 69   fLowEnergy = 20.*CLHEP::MeV;                 <<  73   fLowEnergy = 20.*MeV;
 70   fLowestEnergy = 1.*CLHEP::MeV;               <<  74   fLowestEnergy = 1.*MeV;
 71   SetMinKinEnergy(0.0);                            75   SetMinKinEnergy(0.0);
 72   SetMaxKinEnergy( G4HadronicParameters::Insta     76   SetMaxKinEnergy( G4HadronicParameters::Instance()->GetMaxEnergy() );
 73                                                    77 
 74   fPion    = new G4UPiNuclearCrossSection();   <<  78   fPion = nullptr;
 75   fGlauber = new G4ComponentGGHadronNucleusXsc <<  79   fGlauber = nullptr;
 76   fHadron  = new G4HadronNucleonXsc();         <<  80   fHadron  = nullptr;
 77                                                    81 
 78   fG4pow = G4Pow::GetInstance();               <<  82   fG4pow   = G4Pow::GetInstance();
 79                                                    83 
 80   theProton = G4Proton::Proton();              <<  84   theProton= G4Proton::Proton();
 81   thePiPlus = G4PionPlus::PionPlus();          <<  85   thePiPlus= G4PionPlus::PionPlus();
 82   isPiplus = (p == thePiPlus);                     86   isPiplus = (p == thePiPlus);
                                                   >>  87   isMaster = false;
 83   SetForAllAtomsAndEnergies(true);                 88   SetForAllAtomsAndEnergies(true);
 84                                                << 
 85   if (0 == theA[0]) { Initialise(); }          << 
 86 }                                                  89 }
 87                                                    90 
 88 //....oooOO0OOooo........oooOO0OOooo........oo     91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 89                                                    92 
 90 G4BGGPionElasticXS::~G4BGGPionElasticXS()          93 G4BGGPionElasticXS::~G4BGGPionElasticXS()
 91 {                                                  94 {
 92   delete fHadron;                                  95   delete fHadron;
 93 }                                                  96 }
 94                                                    97 
 95 //....oooOO0OOooo........oooOO0OOooo........oo     98 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 96                                                    99 
 97 G4bool                                            100 G4bool 
 98 G4BGGPionElasticXS::IsElementApplicable(const     101 G4BGGPionElasticXS::IsElementApplicable(const G4DynamicParticle*, G4int,
 99                                         const     102                                         const G4Material*)
100 {                                                 103 {
101   return true;                                    104   return true;
102 }                                                 105 }
103                                                   106 
104 //....oooOO0OOooo........oooOO0OOooo........oo    107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
105                                                   108 
106 G4bool G4BGGPionElasticXS::IsIsoApplicable(con    109 G4bool G4BGGPionElasticXS::IsIsoApplicable(const G4DynamicParticle*, 
107                                            G4i    110                                            G4int Z, G4int, 
108                                            con    111                                            const G4Element*, const G4Material*)
109 {                                                 112 {
110   return (1 == Z);                                113   return (1 == Z);
111 }                                                 114 }
112                                                   115 
113 //....oooOO0OOooo........oooOO0OOooo........oo    116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
114                                                   117 
115 G4double                                          118 G4double
116 G4BGGPionElasticXS::GetElementCrossSection(con    119 G4BGGPionElasticXS::GetElementCrossSection(const G4DynamicParticle* dp,
117                                            G4i    120                                            G4int ZZ, const G4Material*)
118 {                                                 121 {
119   // this method should be called only for Z >    122   // this method should be called only for Z > 1
120   G4double cross = 0.0;                           123   G4double cross = 0.0;
121   G4double ekin = std::max(dp->GetKineticEnerg    124   G4double ekin = std::max(dp->GetKineticEnergy(), fLowestEnergy);
122   G4int Z = std::min(ZZ, 92);                     125   G4int Z = std::min(ZZ, 92);
123   if(1 == Z) {                                    126   if(1 == Z) {
124     cross = 1.0115*GetIsoCrossSection(dp,1,1);    127     cross = 1.0115*GetIsoCrossSection(dp,1,1);
125   } else {                                        128   } else {
126     if(ekin <= fLowEnergy) {                      129     if(ekin <= fLowEnergy) {
127       cross = (isPiplus) ? theCoulombFacPiPlus    130       cross = (isPiplus) ? theCoulombFacPiPlus[Z]*CoulombFactorPiPlus(ekin, Z)
128         : theCoulombFacPiMinus[Z]*FactorPiMinu    131         : theCoulombFacPiMinus[Z]*FactorPiMinus(ekin);
129     } else if(ekin > fGlauberEnergy) {            132     } else if(ekin > fGlauberEnergy) {
130       cross = (isPiplus) ? theGlauberFacPiPlus    133       cross = (isPiplus) ? theGlauberFacPiPlus[Z] : theGlauberFacPiMinus[Z];
131       cross *= fGlauber->GetElasticGlauberGrib    134       cross *= fGlauber->GetElasticGlauberGribov(dp, Z, theA[Z]);
132     } else {                                      135     } else {
133       cross = fPion->GetElasticCrossSection(dp    136       cross = fPion->GetElasticCrossSection(dp, Z, theA[Z]);
134     }                                             137     }
135   }                                               138   }
136 #ifdef G4VERBOSE                               << 
137   if(verboseLevel > 1) {                          139   if(verboseLevel > 1) {
138     G4cout << "G4BGGPionElasticXS::GetElementC    140     G4cout << "G4BGGPionElasticXS::GetElementCrossSection  for "
139            << dp->GetDefinition()->GetParticle    141            << dp->GetDefinition()->GetParticleName()
140            << "  Ekin(GeV)= " << dp->GetKineti    142            << "  Ekin(GeV)= " << dp->GetKineticEnergy()
141            << " in nucleus Z= " << Z << "  A=     143            << " in nucleus Z= " << Z << "  A= " << theA[Z]
142            << " XS(b)= " << cross/barn            144            << " XS(b)= " << cross/barn 
143            << G4endl;                             145            << G4endl;
144   }                                               146   }
145 #endif                                         << 
146   return cross;                                   147   return cross;
147 }                                                 148 }
148                                                   149 
149 //....oooOO0OOooo........oooOO0OOooo........oo << 
150                                                << 
151 G4double                                          150 G4double
152 G4BGGPionElasticXS::GetIsoCrossSection(const G    151 G4BGGPionElasticXS::GetIsoCrossSection(const G4DynamicParticle* dp, 
153                                        G4int,  << 152                                        G4int Z, G4int A, 
154                                        const G    153                                        const G4Isotope*,
155                                        const G    154                                        const G4Element*,
156                                        const G    155                                        const G4Material*)
157 {                                                 156 {
158   // this method should be called only for Z =    157   // this method should be called only for Z = 1
159   fHadron->HadronNucleonXscNS(dp->GetDefinitio    158   fHadron->HadronNucleonXscNS(dp->GetDefinition(), theProton, 
160                               dp->GetKineticEn    159                               dp->GetKineticEnergy());
161   G4double cross = A*fHadron->GetElasticHadron    160   G4double cross = A*fHadron->GetElasticHadronNucleonXsc();
162                                                   161 
163 #ifdef G4VERBOSE                               << 
164   if(verboseLevel > 1) {                          162   if(verboseLevel > 1) {
165     G4cout << "G4BGGPionElasticXS::GetIsoCross    163     G4cout << "G4BGGPionElasticXS::GetIsoCrossSection  for "
166            << dp->GetDefinition()->GetParticle    164            << dp->GetDefinition()->GetParticleName()
167            << "  Ekin(GeV)= " << dp->GetKineti    165            << "  Ekin(GeV)= " << dp->GetKineticEnergy()
168            << " in nucleus Z=1  A=" << A       << 166            << " in nucleus Z= " << Z << "  A= " << A
169            << " XS(b)= " << cross/barn            167            << " XS(b)= " << cross/barn 
170            << G4endl;                             168            << G4endl;
171   }                                               169   }
172 #endif                                         << 
173   return cross;                                   170   return cross;
174 }                                                 171 }
175                                                   172 
176 //....oooOO0OOooo........oooOO0OOooo........oo    173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
177                                                   174 
178 void G4BGGPionElasticXS::BuildPhysicsTable(con    175 void G4BGGPionElasticXS::BuildPhysicsTable(const G4ParticleDefinition& p)
179 {                                                 176 {
                                                   >> 177   if(fPion) { return; }
180   if(verboseLevel > 1) {                          178   if(verboseLevel > 1) {
181     G4cout << "G4BGGPionElasticXS::BuildPhysic    179     G4cout << "G4BGGPionElasticXS::BuildPhysicsTable for " 
182            << p.GetParticleName() << G4endl;      180            << p.GetParticleName() << G4endl;
183   }                                               181   } 
184   if(&p == G4PionPlus::PionPlus() || &p == G4P    182   if(&p == G4PionPlus::PionPlus() || &p == G4PionMinus::PionMinus()) {
185     isPiplus = (&p == G4PionPlus::PionPlus());    183     isPiplus = (&p == G4PionPlus::PionPlus());
186   } else {                                        184   } else {
187     G4ExceptionDescription ed;                    185     G4ExceptionDescription ed;
188     ed << "This BGG cross section is applicabl    186     ed << "This BGG cross section is applicable only to pions and not to " 
189        << p.GetParticleName() << G4endl;          187        << p.GetParticleName() << G4endl; 
190     G4Exception("G4BGGPionElasticXS::BuildPhys    188     G4Exception("G4BGGPionElasticXS::BuildPhysicsTable", "had001", 
191                 FatalException, ed);              189                 FatalException, ed);
                                                   >> 190     return;
192   }                                               191   }
193 }                                              << 
194                                                   192 
195 //....oooOO0OOooo........oooOO0OOooo........oo << 193   fPion    = new G4UPiNuclearCrossSection();
                                                   >> 194   fGlauber = new G4ComponentGGHadronNucleusXsc();
                                                   >> 195   fHadron  = new G4HadronNucleonXsc();
196                                                   196 
197 void G4BGGPionElasticXS::Initialise()          << 197   fPion->BuildPhysicsTable(p);
198 {                                              << 
199   theA[0] = theA[1] = 1;                       << 
200   G4ThreeVector mom(0.0,0.0,1.0);              << 
201   G4DynamicParticle dp(thePiPlus, mom, fGlaube << 
202                                                   198 
203   G4NistManager* nist = G4NistManager::Instanc << 199   if(0 == theA[0]) { 
                                                   >> 200 #ifdef G4MULTITHREADED
                                                   >> 201     G4MUTEXLOCK(&pionElasticXSMutex);
                                                   >> 202     if(0 == theA[0]) { 
                                                   >> 203 #endif
                                                   >> 204       isMaster = true;
                                                   >> 205 #ifdef G4MULTITHREADED
                                                   >> 206     }
                                                   >> 207     G4MUTEXUNLOCK(&pionElasticXSMutex);
                                                   >> 208 #endif
                                                   >> 209   } else {
                                                   >> 210     return;
                                                   >> 211   }
204                                                   212 
205   G4double csup, csdn;                         << 213   if(isMaster && 0 == theA[0]) {
206   for (G4int iz=2; iz<93; ++iz) {              << 
207                                                   214 
208     G4int A = G4lrint(nist->GetAtomicMassAmu(i << 215     theA[0] = theA[1] = 1;
209     theA[iz] = A;                              << 216     G4ThreeVector mom(0.0,0.0,1.0);
                                                   >> 217     G4DynamicParticle dp(thePiPlus, mom, fGlauberEnergy);
210                                                   218 
211     csup = fGlauber->GetElasticGlauberGribov(& << 219     G4NistManager* nist = G4NistManager::Instance();
212     csdn = fPion->GetElasticCrossSection(&dp,  << 
213     theGlauberFacPiPlus[iz] = csdn/csup;       << 
214   }                                            << 
215                                                   220 
216   dp.SetDefinition(G4PionMinus::PionMinus());  << 221     G4double csup, csdn;
217   for (G4int iz=2; iz<93; ++iz) {              << 222     for(G4int iz=2; iz<93; ++iz) {
218     csup = fGlauber->GetElasticGlauberGribov(& << 223 
219     csdn = fPion->GetElasticCrossSection(&dp,  << 224       G4int A = G4lrint(nist->GetAtomicMassAmu(iz));
220     theGlauberFacPiMinus[iz] = csdn/csup;      << 225       theA[iz] = A;
221                                                << 226 
222     if (verboseLevel > 1) {                    << 227       csup = fGlauber->GetElasticGlauberGribov(&dp, iz, A);
223       G4cout << "Z= " << iz <<  "  A= " << the << 228       csdn = fPion->GetElasticCrossSection(&dp, iz, A);
224        << " factorPiPlus= " << theGlauberFacPi << 229       theGlauberFacPiPlus[iz] = csdn/csup;
225        << " factorPiMinus= " << theGlauberFacP << 
226        << G4endl;                              << 
227     }                                             230     }
228   }                                            << 231 
229   theCoulombFacPiPlus[1] = 1.0;                << 232     dp.SetDefinition(G4PionMinus::PionMinus());
230   theCoulombFacPiMinus[1]= 1.0;                << 233     for(G4int iz=2; iz<93; ++iz) {
231   dp.SetKineticEnergy(fLowEnergy);             << 234       csup = fGlauber->GetElasticGlauberGribov(&dp, iz, theA[iz]);
232   dp.SetDefinition(thePiPlus);                 << 235       csdn = fPion->GetElasticCrossSection(&dp, iz, theA[iz]);
233   for (G4int iz=2; iz<93; ++iz) {              << 236       theGlauberFacPiMinus[iz] = csdn/csup;
234     theCoulombFacPiPlus[iz] = fPion->GetElasti << 237 
235       /CoulombFactorPiPlus(fLowEnergy, iz);    << 238       if(verboseLevel > 0) {
236   }                                            << 239         G4cout << "Z= " << iz <<  "  A= " << theA[iz] 
237   dp.SetDefinition(G4PionMinus::PionMinus());  << 240                << " factorPiPlus= " << theGlauberFacPiPlus[iz] 
238   for(G4int iz=2; iz<93; ++iz) {               << 241                << " factorPiMinus= " << theGlauberFacPiMinus[iz] 
239     theCoulombFacPiMinus[iz] = fPion->GetElast << 242                << G4endl; 
240       /FactorPiMinus(fLowEnergy);              << 243       }
241                                                << 244     }
242     if(verboseLevel > 1) {                     << 245     theCoulombFacPiPlus[1] = 1.0;
243       G4cout << "Z= " << iz <<  "  A= " << the << 246     theCoulombFacPiMinus[1]= 1.0;
244        << " CoulombFactorPiPlus= " << theCoulo << 247     dp.SetKineticEnergy(fLowEnergy);
245        << " CoulombFactorPiMinus= " << theCoul << 248     dp.SetDefinition(thePiPlus);
246        << G4endl;                              << 249     for(G4int iz=2; iz<93; ++iz) {
                                                   >> 250       theCoulombFacPiPlus[iz] = fPion->GetElasticCrossSection(&dp, iz, theA[iz])
                                                   >> 251         /CoulombFactorPiPlus(fLowEnergy, iz);
                                                   >> 252     }
                                                   >> 253     dp.SetDefinition(G4PionMinus::PionMinus());
                                                   >> 254     for(G4int iz=2; iz<93; ++iz) {
                                                   >> 255       theCoulombFacPiMinus[iz] = fPion->GetElasticCrossSection(&dp, iz, theA[iz])
                                                   >> 256         /FactorPiMinus(fLowEnergy);
                                                   >> 257 
                                                   >> 258       if(verboseLevel > 0) {
                                                   >> 259         G4cout << "Z= " << iz <<  "  A= " << theA[iz] 
                                                   >> 260                << " CoulombFactorPiPlus= " << theCoulombFacPiPlus[iz] 
                                                   >> 261                << " CoulombFactorPiMinus= " << theCoulombFacPiMinus[iz] 
                                                   >> 262                << G4endl; 
                                                   >> 263       }
247     }                                             264     }
248   }                                               265   }
249 }                                                 266 }
250                                                   267 
251 //....oooOO0OOooo........oooOO0OOooo........oo    268 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
252                                                   269 
253 G4double G4BGGPionElasticXS::CoulombFactorPiPl    270 G4double G4BGGPionElasticXS::CoulombFactorPiPlus(G4double kinEnergy, G4int Z)
254 {                                                 271 {
255   return (kinEnergy > 0.0) ?                      272   return (kinEnergy > 0.0) ? 
256     G4NuclearRadii::CoulombFactor(Z, theA[Z],     273     G4NuclearRadii::CoulombFactor(Z, theA[Z], thePiPlus, kinEnergy) : 0.0;
257 }                                                 274 }
258                                                   275 
259 //....oooOO0OOooo........oooOO0OOooo........oo    276 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
260                                                   277 
261 G4double G4BGGPionElasticXS::FactorPiMinus(G4d    278 G4double G4BGGPionElasticXS::FactorPiMinus(G4double kinEnergy)
262 {                                                 279 {
263   return 1.0/std::sqrt(kinEnergy);                280   return 1.0/std::sqrt(kinEnergy);
264 }                                                 281 }
265                                                   282 
266 //....oooOO0OOooo........oooOO0OOooo........oo    283 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
267                                                   284 
268 void                                              285 void
269 G4BGGPionElasticXS::CrossSectionDescription(st    286 G4BGGPionElasticXS::CrossSectionDescription(std::ostream& outFile) const 
270 {                                                 287 {
271   outFile << "The Barashenkov-Glauber-Gribov c    288   outFile << "The Barashenkov-Glauber-Gribov cross section handles elastic\n"
272           << "scattering of pions from nuclei     289           << "scattering of pions from nuclei at all energies. The\n"
273           << "Barashenkov parameterization is     290           << "Barashenkov parameterization is used below 91 GeV and the\n"
274           << "Glauber-Gribov parameterization     291           << "Glauber-Gribov parameterization is used above 91 GeV.\n";
275 }                                                 292 }
276                                                   293 
277 //....oooOO0OOooo........oooOO0OOooo........oo    294 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
278                                                   295