ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNAIonElasticModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4DNAIonElasticModel.hh
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: H. N. Tran (Ton Duc Thang University)
27 // p, H, He, He+ and He++ models are assumed identical
28 // NIMB 343, 132-137 (2015)
29 //
30 // The Geant4-DNA web site is available at http://geant4-dna.org
31 //
32 
33 #ifndef G4DNAIonElasticModel_h
34 #define G4DNAIonElasticModel_h 1
35 
36 #include <map>
38 #include "G4VEmModel.hh"
39 #include "G4Proton.hh"
42 #include "G4LogLogInterpolation.hh"
43 #include "G4ProductionCutsTable.hh"
44 #include "G4NistManager.hh"
45 
47 {
48 
49 public:
50 
52  const G4String& nam ="DNAIonElasticModel");
53 
54  virtual
56 
57  virtual void
58  Initialise (const G4ParticleDefinition* particuleDefinition,
59  const G4DataVector&);
60 
61  virtual G4double
63  const G4ParticleDefinition* p, G4double ekin,
65 
66  virtual void
67  SampleSecondaries (std::vector<G4DynamicParticle*>*,
69  G4double tmin, G4double maxEnergy);
70 
71  void
72  SetKillBelowThreshold (G4double threshold);
73 
74  G4double
76  {
77  return killBelowEnergy;
78  }
79 
80  inline void SelectStationary(G4bool input);
81 
82 protected:
83 
85 
86 private:
87 
89 
90  // Water density table
91  const std::vector<G4double>* fpMolWaterDensity;
92 
98 
100 
101  // Cross section
103 
104  // Final state
105 
106  G4double
107  Theta (G4ParticleDefinition * aParticleDefinition, G4double k,
108  G4double integrDiff);
109 
110  G4double
112  G4double xs2);
113 
114  G4double
116  G4double xs2);
117 
118  G4double
119  QuadInterpolator (G4double e11, G4double e12, G4double e21, G4double e22,
120  G4double x11, G4double x12, G4double x21, G4double x22,
122 
123  G4double
125  G4double xs2);
126 
127  typedef std::map<G4double, std::map<G4double, G4double> > TriDimensionMap;
129 
130  std::vector<G4double> eTdummyVec;
131 
132  typedef std::map<G4double, std::vector<G4double> > VecMap;
134 
135  G4double
136  RandomizeThetaCM (G4double k, G4ParticleDefinition * aParticleDefinition);
137 
138  //
139 
143 
144 };
145 
146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
147 
149 {
150  statCode = input;
151 }
152 
153 #endif
154