ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4DNACPA100ElasticModel.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4DNACPA100ElasticModel.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 // CPA100 elastic model class for electrons
27 //
28 // Based on the work of M. Terrissol and M. C. Bordage
29 //
30 // Users are requested to cite the following papers:
31 // - M. Terrissol, A. Baudre, Radiat. Prot. Dosim. 31 (1990) 175-177
32 // - M.C. Bordage, J. Bordes, S. Edel, M. Terrissol, X. Franceries,
33 // M. Bardies, N. Lampe, S. Incerti, Phys. Med. 32 (2016) 1833-1840
34 //
35 // Authors of this class:
36 // M.C. Bordage, M. Terrissol, S. Edel, J. Bordes, S. Incerti
37 //
38 // 15.01.2014: creation
39 //
40 
41 #ifndef G4DNACPA100ElasticModel_h
42 #define G4DNACPA100ElasticModel_h 1
43 
44 #include <map>
46 #include "G4VEmModel.hh"
47 #include "G4Electron.hh"
49 #include "G4LogLogInterpolation.hh"
50 //#include "G4DNACPA100LogLogInterpolation.hh"
51 #include "G4ProductionCutsTable.hh"
52 #include "G4NistManager.hh"
53 
55 {
56 
57 public:
58 
60  const G4String& nam = "DNACPA100ElasticModel");
61 
62  virtual ~G4DNACPA100ElasticModel();
63 
64  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
65 
67  const G4ParticleDefinition* p,
68  G4double ekin,
69  G4double emin,
70  G4double emax);
71 
72  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
73  const G4MaterialCutsCouple*,
74  const G4DynamicParticle*,
75  G4double tmin,
76  G4double maxEnergy);
77 
78  //---
79  // kept for backward compatibility
80  inline void SetKillBelowThreshold (G4double threshold);
82  //---
83 
84  inline void SelectStationary(G4bool input);
85 
86 protected:
87 
89 
90 private:
91 
93 
94  // Water density table
95  const std::vector<G4double>* fpMolWaterDensity;
96 
102 
103  // Cross section
104 
105  typedef std::map<G4String,G4String,std::less<G4String> > MapFile;
107 
108  typedef std::map<G4String,G4DNACrossSectionDataSet*,std::less<G4String> > MapData;
110 
111  // Final state
112 
113  //G4double DifferentialCrossSection(G4ParticleDefinition * aParticleDefinition, G4double k, G4double theta);
114 
115  G4double Theta(G4ParticleDefinition * aParticleDefinition, G4double k, G4double integrDiff);
116 
118 
120 
122 
124  G4double e12,
125  G4double e21,
126  G4double e22,
127  G4double x11,
128  G4double x12,
129  G4double x21,
130  G4double x22,
131  G4double t1,
132  G4double t2,
133  G4double t,
134  G4double e);
135 
136  typedef std::map<G4double, std::map<G4double, G4double> > TriDimensionMap;
137 
139  std::vector<G4double> eTdummyVec;
140 
141  typedef std::map<G4double, std::vector<G4double> > VecMap;
143 
145 
146  //
147 
150 
151 };
152 
153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
154 
156 {
157  statCode = input;
158 }
159 
160 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
161 
162 //---
163 // kept for backward compatibility
164 
165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
166 
168 {
169 
170  G4ExceptionDescription errMsg;
171  errMsg << "*** WARNING : "
172  << "G4DNACPA100ElasticModel::SetKillBelowThreshold"
173  << "is deprecated, the kill threshold won't be taken into account";
174 
175  G4Exception ("*** WARNING : G4DNACPA100ElasticModel::SetKillBelowThreshold DEPRECATED","",JustWarning,"") ;
176 
177 }
178 
179 #endif