Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // author: Vladimir.Grichine@cern.ch 27 // 28 // Implements data from: Barashenkov V.S., Nucleon-Nucleus Cross Section, 29 // Preprint JINR P2-89-770, p. 12, Dubna 1989 (scanned version from KEK) 30 // Based on G. Folger version of G4PiNuclearCrossSection class 31 // 32 // Modified: 33 // 05.03.2024 V.Ivanchenko removed obsolete methods and calls 34 // 35 36 #ifndef G4NucleonNuclearCrossSection_h 37 #define G4NucleonNuclearCrossSection_h 38 39 #include "G4VCrossSectionDataSet.hh" 40 #include "G4ParticleDefinition.hh" 41 #include "globals.hh" 42 43 class G4ComponentBarNucleonNucleusXsc; 44 class G4ParticleDefinition; 45 46 class G4NucleonNuclearCrossSection : public G4VCrossSectionDataSet 47 { 48 public: 49 50 G4NucleonNuclearCrossSection(); 51 ~G4NucleonNuclearCrossSection() override = default; 52 53 static const char* Default_Name() { return "BarashenkovNucleonXS"; } 54 55 G4bool IsElementApplicable(const G4DynamicParticle* aParticle, 56 G4int Z, const G4Material* mat) final; 57 58 // return inelastic x-section 59 G4double GetElementCrossSection(const G4DynamicParticle* aParticle, 60 G4int Z, const G4Material* mat=nullptr) final; 61 62 void CrossSectionDescription(std::ostream&) const final; 63 64 // return elastic x-section 65 inline G4double GetElasticCrossSection(const G4DynamicParticle* aParticle, 66 G4int Z); 67 68 // access methods should be called after ComputeCrossSection(...) 69 inline G4double GetTotalXsc() { return fTotalXsc; }; 70 inline G4double GetInelasticXsc() { return fInelasticXsc; }; 71 inline G4double GetElasticXsc() { return fElasticXsc; }; 72 73 G4NucleonNuclearCrossSection& operator= 74 (const G4NucleonNuclearCrossSection &right) = delete; 75 G4NucleonNuclearCrossSection(const G4NucleonNuclearCrossSection&) = delete; 76 77 private: 78 79 void ComputeCrossSections(const G4ParticleDefinition*, 80 G4double kinEnergy, G4int Z); 81 82 G4ComponentBarNucleonNucleusXsc* fBarash; 83 84 G4double fTotalXsc{0.0}; 85 G4double fInelasticXsc{0.0}; 86 G4double fElasticXsc{0.0}; 87 }; 88 89 inline 90 G4double G4NucleonNuclearCrossSection::GetElasticCrossSection( 91 const G4DynamicParticle* dp, G4int Z) 92 { 93 ComputeCrossSections(dp->GetDefinition(), dp->GetKineticEnergy(), Z); 94 return fElasticXsc; 95 } 96 97 #endif 98