File indexing completed on 2025-01-18 09:57:58
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047
0048 #ifndef G4BGGNucleonElasticXS_h
0049 #define G4BGGNucleonElasticXS_h 1
0050
0051 #include "globals.hh"
0052 #include "G4VCrossSectionDataSet.hh"
0053 #include "G4ParticleDefinition.hh"
0054 #include "G4Threading.hh"
0055
0056 class G4ComponentGGHadronNucleusXsc;
0057 class G4NucleonNuclearCrossSection;
0058 class G4HadronNucleonXsc;
0059 class G4Material;
0060 class G4Element;
0061 class G4Isotope;
0062
0063 class G4BGGNucleonElasticXS final : public G4VCrossSectionDataSet
0064 {
0065 public:
0066
0067 explicit G4BGGNucleonElasticXS (const G4ParticleDefinition*);
0068
0069 ~G4BGGNucleonElasticXS() final;
0070
0071 G4bool IsElementApplicable(const G4DynamicParticle*, G4int Z,
0072 const G4Material* mat) final;
0073
0074 G4bool IsIsoApplicable(const G4DynamicParticle*, G4int Z, G4int A,
0075 const G4Element* elm,
0076 const G4Material* mat) final;
0077
0078 G4double GetElementCrossSection(const G4DynamicParticle*, G4int Z,
0079 const G4Material* mat) final;
0080
0081 G4double GetIsoCrossSection(const G4DynamicParticle*, G4int Z, G4int A,
0082 const G4Isotope* iso=nullptr,
0083 const G4Element* elm=nullptr,
0084 const G4Material* mat=nullptr) final;
0085
0086 void BuildPhysicsTable(const G4ParticleDefinition&) final;
0087
0088 void CrossSectionDescription(std::ostream&) const final;
0089
0090 private:
0091
0092 G4double CoulombFactor(G4double kinEnergy, G4int Z);
0093
0094 G4BGGNucleonElasticXS & operator=(const G4BGGNucleonElasticXS &right);
0095 G4BGGNucleonElasticXS(const G4BGGNucleonElasticXS&);
0096
0097 G4double fGlauberEnergy;
0098 G4double fLowEnergy;
0099
0100 static G4double theGlauberFacP[93];
0101 static G4double theCoulombFacP[93];
0102 static G4double theGlauberFacN[93];
0103 static G4double theCoulombFacN[93];
0104 static G4int theA[93];
0105
0106 const G4ParticleDefinition* theProton;
0107 G4ComponentGGHadronNucleusXsc* fGlauber;
0108 G4NucleonNuclearCrossSection* fNucleon;
0109 G4HadronNucleonXsc* fHadron;
0110 G4bool isProton;
0111 G4bool isMaster;
0112
0113 #ifdef G4MULTITHREADED
0114 static G4Mutex nucleonElasticXSMutex;
0115 #endif
0116 };
0117
0118 #endif