Line | |
---|
1 | #include <iostream> |
---|
2 | #include <cmath> |
---|
3 | |
---|
4 | #include "crystal.h" |
---|
5 | using namespace std; |
---|
6 | |
---|
7 | |
---|
8 | |
---|
9 | //==================================Constructeurs, destructeur==================================== |
---|
10 | |
---|
11 | Crystal::Crystal(int C_orient, int IS, double C_xmax, double C_ymax, double Cry_length, double Rcurv) |
---|
12 | : miscut(miscut), Rcurv(Rcurv), Cry_length(Cry_length), C_xmax(C_xmax), C_ymax(C_ymax), Alayer(Alayer), IS(IS), C_orient(C_orient), ZN(ZN), NAM(NAM) |
---|
13 | { |
---|
14 | |
---|
15 | //Some usefull parameters !! |
---|
16 | |
---|
17 | Alayer = 0; |
---|
18 | |
---|
19 | cry_bend = (Cry_length / Rcurv); //total bending for this crystal; |
---|
20 | s_length = Rcurv * (sin(cry_bend)); |
---|
21 | xmax = C_xmax ; |
---|
22 | L_chan = Cry_length; //Careful : the channeling length is the length of the particle trajectory, not the long. one |
---|
23 | tchan = cry_bend; |
---|
24 | ymin = -C_ymax / 2 ; |
---|
25 | ymax = C_ymax / 2; |
---|
26 | |
---|
27 | ZN = 1; |
---|
28 | NAM = 1; |
---|
29 | miscut = 0; |
---|
30 | |
---|
31 | DLAI[0] = 1.6; |
---|
32 | DLAI[1] = 0.57; |
---|
33 | DLAI[2] = 2.2; |
---|
34 | DLAI[3] = 1.0; |
---|
35 | SAI[0] = 42.0; |
---|
36 | SAI[1] = 140.0; |
---|
37 | SAI[2] = 42.0; |
---|
38 | SAI[3] = 50.0; |
---|
39 | DLRI[0] = 0.09336; |
---|
40 | DLRI[1] = 0.0035; |
---|
41 | DLRI[2] = 0.188; |
---|
42 | DLRI[3] = 0.023; |
---|
43 | DLYI[0] = 0.455; |
---|
44 | DLYI[1] = 0.096; |
---|
45 | DLYI[2] = 0.400; |
---|
46 | DLYI[3] = 0.162; |
---|
47 | AI[0] = 0.96e-7; |
---|
48 | AI[1] = 0.56e-7; |
---|
49 | AI[2] = 0.63e-7; |
---|
50 | AI[3] = 1.0e-7; |
---|
51 | DES[0] = 0.56; |
---|
52 | DES[1] = 3.0; |
---|
53 | DES[2] = 0.6; |
---|
54 | DES[3] = 1.0; |
---|
55 | eUm[0] = 21.34; |
---|
56 | eUm[1] = 21.0; |
---|
57 | eUm[2] = 21.0; |
---|
58 | eUm[3] = 21.0; |
---|
59 | |
---|
60 | } |
---|
61 | |
---|
62 | /* |
---|
63 | void Crystal::afficheCrystal(){ |
---|
64 | |
---|
65 | cout << "Here is the implicit parameters : " <<endl; |
---|
66 | cout << "L_chan:" << L_chan << endl; |
---|
67 | cout << "tchan:" << tchan << endl; |
---|
68 | |
---|
69 | } |
---|
70 | */ |
---|
71 | |
---|
Note: See
TracBrowser
for help on using the repository browser.