Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNAIonElasticModel.hh
Go to the documentation of this file.
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,
64  G4double emin, G4double emax);
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 
88  G4bool statCode;
89 
90  // Water density table
91  const std::vector<G4double>* fpMolWaterDensity;
92 
93  G4double killBelowEnergy;
94  G4double lowEnergyLimit;
95  G4double highEnergyLimit;
96  G4bool isInitialised;
97  G4int verboseLevel;
98 
99  G4double fParticle_Mass;
100 
101  // Cross section
102  G4DNACrossSectionDataSet* fpTableData;
103 
104  // Final state
105 
106  G4double
107  Theta (G4ParticleDefinition * aParticleDefinition, G4double k,
108  G4double integrDiff);
109 
110  G4double
111  LinLogInterpolate (G4double e1, G4double e2, G4double e, G4double xs1,
112  G4double xs2);
113 
114  G4double
115  LogLogInterpolate (G4double e1, G4double e2, G4double e, G4double xs1,
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,
121  G4double t1, G4double t2, G4double t, G4double e);
122 
123  G4double
124  LinLinInterpolate (G4double e1, G4double e2, G4double e, G4double xs1,
125  G4double xs2);
126 
127  typedef std::map<double, std::map<double, double> > TriDimensionMap;
128  TriDimensionMap fDiffCrossSectionData;
129 
130  std::vector<double> eTdummyVec;
131 
132  typedef std::map<double, std::vector<double> > VecMap;
133  VecMap eVecm;
134 
135  G4double
136  RandomizeThetaCM (G4double k, G4ParticleDefinition * aParticleDefinition);
137 
138  //
139 
141  operator= (const G4DNAIonElasticModel &right);
143 
144 };
145 
146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
147 
149 {
150  statCode = input;
151 }
152 
153 #endif
154 
const char * p
Definition: xmltok.h:285
int G4int
Definition: G4Types.hh:78
void SetKillBelowThreshold(G4double threshold)
string material
Definition: eplot.py:19
G4ParticleChangeForGamma * fParticleChangeForGamma
bool G4bool
Definition: G4Types.hh:79
virtual void Initialise(const G4ParticleDefinition *particuleDefinition, const G4DataVector &)
static const G4double emax
tuple t1
Definition: plottest35.py:33
virtual G4double CrossSectionPerVolume(const G4Material *material, const G4ParticleDefinition *p, G4double ekin, G4double emin, G4double emax)
double G4double
Definition: G4Types.hh:76
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
void SelectStationary(G4bool input)
G4DNAIonElasticModel(const G4ParticleDefinition *p=0, const G4String &nam="DNAIonElasticModel")