| 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 //
27
28 #ifndef G4DNAScreenedRutherfordElasticModel_h
29 #define G4DNAScreenedRutherfordElasticModel_h 1
30
31 #include <CLHEP/Units/SystemOfUnits.h>
32
33 #include "G4VEmModel.hh"
34 #include "G4ParticleChangeForGamma.hh"
35 #include "G4ProductionCutsTable.hh"
36 #include "G4NistManager.hh"
37
38 class G4DNAScreenedRutherfordElasticModel : public G4VEmModel
39 {
40 public:
41 G4DNAScreenedRutherfordElasticModel(const G4ParticleDefinition* p = nullptr,
42 const G4String& nam = "DNAScreenedRutherfordElasticModel");
43
44 ~G4DNAScreenedRutherfordElasticModel() override;
45
46 G4DNAScreenedRutherfordElasticModel& operator=
47 (const G4DNAScreenedRutherfordElasticModel &right) = delete;
48 G4DNAScreenedRutherfordElasticModel(const G4DNAScreenedRutherfordElasticModel&) = delete;
49
50 void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
51
52 G4double CrossSectionPerVolume(const G4Material* material,
53 const G4ParticleDefinition* p,
54 G4double ekin,
55 G4double emin,
56 G4double emax) override;
57
58 void SampleSecondaries(std::vector<G4DynamicParticle*>*,
59 const G4MaterialCutsCouple*,
60 const G4DynamicParticle*,
61 G4double tmin,
62 G4double maxEnergy) override;
63
64 inline void SetKillBelowThreshold (G4double threshold);
65
66 inline void SelectFasterComputation(G4bool input);
67
68 protected:
69 G4ParticleChangeForGamma* fParticleChangeForGamma;
70
71 private:
72 G4double lowEnergyLimit;
73 G4double intermediateEnergyLimit;
74 G4double highEnergyLimit;
75
76 // -- Brenner & Zaider
77 std::vector<G4double> betaCoeff;
78 std::vector<G4double> deltaCoeff;
79 std::vector<G4double> gamma035_10Coeff;
80 std::vector<G4double> gamma10_100Coeff;
81 std::vector<G4double> gamma100_200Coeff;
82
83 // -- Water density table
84 const std::vector<G4double>* fpWaterDensity;
85
86 G4int verboseLevel;
87
88 G4bool isInitialised{false};
89 G4bool fasterCode;
90
91 // -- Cross section
92 G4double RutherfordCrossSection(G4double energy, G4double z);
93 G4double ScreeningFactor(G4double energy, G4double z);
94
95 // -- Final state according to Brenner & Zaider
96 G4double BrennerZaiderRandomizeCosTheta(G4double k);
97 G4double CalculatePolynomial(G4double k, std::vector<G4double>& vec);
98
99
100 // -- Final state according to Screened Rutherford
101 G4double ScreenedRutherfordRandomizeCosTheta(G4double k, G4double z);
102
103 };
104
105 inline void
106 G4DNAScreenedRutherfordElasticModel::SetKillBelowThreshold (G4double)
107 {
108 G4ExceptionDescription errMsg;
109 errMsg << "The method G4DNAScreenedRutherfordElasticModel::"
110 "SetKillBelowThreshold is deprecated";
111
112 G4Exception("G4DNAScreenedRutherfordElasticModel::SetKillBelowThreshold",
113 "deprecated",
114 JustWarning,
115 errMsg);
116 }
117
118 inline void
119 G4DNAScreenedRutherfordElasticModel::SelectFasterComputation(G4bool input)
120 {
121 fasterCode = input;
122 }
123
124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
125
126 #endif
127