File indexing completed on 2025-01-18 09:58:24
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
0049 #ifndef G4hNuclearStoppingModel_h
0050 #define G4hNuclearStoppingModel_h 1
0051
0052 #include "globals.hh"
0053 #include "G4VLowEnergyModel.hh"
0054 #include "G4VhNuclearStoppingPower.hh"
0055
0056 class G4ParticleDefinition;
0057 class G4Material;
0058 class G4DynamicParticle;
0059
0060 class G4hNuclearStoppingModel : public G4VLowEnergyModel
0061 {
0062 public:
0063
0064 explicit G4hNuclearStoppingModel(const G4String& name);
0065
0066 ~G4hNuclearStoppingModel() ;
0067
0068 G4double TheValue(const G4DynamicParticle* particle,
0069 const G4Material* material) override;
0070
0071 G4double TheValue(const G4ParticleDefinition* aParticle,
0072 const G4Material* material,
0073 G4double kineticEnergy) override;
0074
0075 G4double HighEnergyLimit(const G4ParticleDefinition* ,
0076 const G4Material* ) const override
0077 {return highEnergyLimit;};
0078
0079 G4double LowEnergyLimit(const G4ParticleDefinition* ,
0080 const G4Material* ) const override
0081 {return lowEnergyLimit;};
0082
0083 G4double HighEnergyLimit(const G4ParticleDefinition* ) const override
0084 {return highEnergyLimit;};
0085
0086 G4double LowEnergyLimit(const G4ParticleDefinition* ) const override
0087 {return lowEnergyLimit;};
0088
0089 G4bool IsInCharge(const G4DynamicParticle* ,
0090 const G4Material* ) const override
0091 {return true;};
0092
0093 G4bool IsInCharge(const G4ParticleDefinition* ,
0094 const G4Material* ) const override
0095 {return true;};
0096
0097 void SetNuclearStoppingFluctuationsOn()
0098 {nStopingPowerTable->SetNuclearStoppingFluctuationsOn();};
0099
0100 void SetNuclearStoppingFluctuationsOff()
0101 {nStopingPowerTable->SetNuclearStoppingFluctuationsOff();};
0102
0103 G4hNuclearStoppingModel(G4hNuclearStoppingModel &) = delete;
0104 G4hNuclearStoppingModel & operator=(const G4hNuclearStoppingModel &right) = delete;
0105
0106 private:
0107 void InitializeMe();
0108
0109 G4double StoppingPower(const G4Material* material,
0110 G4double kineticEnergy,
0111 G4double z1,
0112 G4double m1) const;
0113
0114
0115 G4VhNuclearStoppingPower* nStopingPowerTable;
0116 G4String modelName;
0117
0118 G4double factorPDG2AMU;
0119
0120
0121 G4double theZieglerFactor;
0122
0123
0124 G4double lowEnergyLimit;
0125 G4double highEnergyLimit;
0126 };
0127
0128 #endif