| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // Author: Alexei Sytov <<
27 26
28 #ifndef G4ChannelingFastSimInterpolation_h 27 #ifndef G4ChannelingFastSimInterpolation_h
29 #define G4ChannelingFastSimInterpolation_h 28 #define G4ChannelingFastSimInterpolation_h
30 29
31 #include "G4PhysicsVector.hh" 30 #include "G4PhysicsVector.hh"
32 #include "G4Physics2DVector.hh" 31 #include "G4Physics2DVector.hh"
33 #include "G4ThreeVector.hh" 32 #include "G4ThreeVector.hh"
34 #include "G4PhysicsLinearVector.hh" 33 #include "G4PhysicsLinearVector.hh"
35 34
36 /** \file G4ChannelingFastSimInterpolation.hh 35 /** \file G4ChannelingFastSimInterpolation.hh
37 * \brief Definition of the G4ChannelingFastSim 36 * \brief Definition of the G4ChannelingFastSimInterpolation class
38 * The class includes spline interpolation coef 37 * The class includes spline interpolation coefficients for the important functions
39 * needed in Channeling FastSimulation model, i 38 * needed in Channeling FastSimulation model, i.e. electric fields, nuclear and
40 * electron densities and minimum energy of ion 39 * electron densities and minimum energy of ionization. All the functions have
41 * transverse coordinates as arguments (x for c 40 * transverse coordinates as arguments (x for crystal planes, x and y for crystal axes).
42 * All the functions are calculated in the co-r 41 * All the functions are calculated in the co-rotating reference system (along crystal
43 * planes/axes). 42 * planes/axes).
44 */ 43 */
45 44
46 class G4ChannelingFastSimInterpolation { 45 class G4ChannelingFastSimInterpolation {
47 46
48 public: 47 public:
49 G4ChannelingFastSimInterpolation(G4double 48 G4ChannelingFastSimInterpolation(G4double dx0,
50 G4double dy 49 G4double dy0,
51 G4int nPoin 50 G4int nPointsx0,
52 G4int nPoin 51 G4int nPointsy0,
53 G4int iMode 52 G4int iModel0);
54 ~G4ChannelingFastSimInterpolation() = defa 53 ~G4ChannelingFastSimInterpolation() = default;
55 54
56 ///Get Spline Function 55 ///Get Spline Function
57 G4double GetIF(G4double xx, G4double yy); 56 G4double GetIF(G4double xx, G4double yy);
58 57
59 ///Set spline coefficients 58 ///Set spline coefficients
60 void SetCoefficients1D(G4double AI0, 59 void SetCoefficients1D(G4double AI0,
61 G4double BI0, 60 G4double BI0,
62 G4double CI0, 61 G4double CI0,
63 G4double DI0, 62 G4double DI0,
64 G4int i); 63 G4int i);
65 void SetCoefficients2D(G4double AI3D0, 64 void SetCoefficients2D(G4double AI3D0,
66 G4double BI3D0, 65 G4double BI3D0,
67 G4double CI3D0, 66 G4double CI3D0,
68 G4int i, 67 G4int i,
69 G4int j, 68 G4int j,
70 G4int k); 69 G4int k);
71 70
72 private: 71 private:
73 G4double Spline1D(G4double xx); 72 G4double Spline1D(G4double xx);
74 G4double Spline2D(G4double xx, G4double yy 73 G4double Spline2D(G4double xx, G4double yy);// cubic spline of 2-variable function
75 74
76 G4double fDx=0., fDy=0.; //channel width a 75 G4double fDx=0., fDy=0.; //channel width and height
77 G4double fStepi=0., fStepj=0.; //interpola 76 G4double fStepi=0., fStepj=0.; //interpolation steps in x and y, respectively
78 G4double fStepi2=0.; //=fStepi*fStepi 77 G4double fStepi2=0.; //=fStepi*fStepi
79 G4int nPointsx=0, nPointsy=0; //number of 78 G4int nPointsx=0, nPointsy=0; //number of interpolation nodes in x and y, respectively
80 79
81 std::vector <G4double> fAI; 80 std::vector <G4double> fAI;
82 std::vector <G4double> fBI; 81 std::vector <G4double> fBI;
83 std::vector <G4double> fCI; 82 std::vector <G4double> fCI;
84 std::vector <G4double> fDI; 83 std::vector <G4double> fDI;
85 84
86 std::vector<std::vector<G4double>> fAI3D; 85 std::vector<std::vector<G4double>> fAI3D;
87 std::vector<std::vector<G4double>> fBI3D; 86 std::vector<std::vector<G4double>> fBI3D;
88 std::vector<std::vector<G4double>> fCI3D; 87 std::vector<std::vector<G4double>> fCI3D;
89 std::vector<std::vector<G4double>> fAI3D3; 88 std::vector<std::vector<G4double>> fAI3D3;
90 std::vector<std::vector<G4double>> fBI3D3; 89 std::vector<std::vector<G4double>> fBI3D3;
91 std::vector<std::vector<G4double>> fCI3D3; 90 std::vector<std::vector<G4double>> fCI3D3;
92 91
93 G4int iModel=1; 92 G4int iModel=1;
94 93
95 }; 94 };
96 95
97 #endif 96 #endif
98 97