BmnRoot 26.10.0-18-g4e3414f9
Software framework of the BM@N experiment
Loading...
Searching...
No Matches
convertFieldCSV2ROOT_1900.C
Go to the documentation of this file.
1/*
2Macro takes files with magnetic field measurements (in V) and
3creates magnetic field map in root-format (in kGs)
4Author: S. Merts (JINR), S. Lebedeva (SPSU)
5*/
6
7R__ADD_INCLUDE_PATH($VMCWORKDIR)
8
9typedef struct
10{
11 Int_t N;
12 Float_t min;
13 Float_t max;
14 Float_t step;
15} Coordinate_info_t;
16
17// convertion from voltage to kGauss for NEW and OLD hall sensors
18Double_t ConvertVoltageToFeild(Double_t U);
19Double_t ConvertVoltageToFeild_X(Double_t U);
20Double_t ConvertVoltageToFeild_Y(Double_t U);
21Double_t ConvertVoltageToFeild_Z(Double_t U);
22
23void convertFieldCSV2ROOT_1900()
24{
25
26 const Int_t nZ = 47; // vertical planes
27 vector<Int_t> zCoord = {1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31,
28 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63,
29 65, 67, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91, 93};
30
31 vector<TString> files[nZ];
32
33 //==== Z=1 ====//
34 files[0].push_back("/home/merz/data/magfield/1900/1_new/2022x05x30_file_3_5250_CW.csv");
35 files[0].push_back("/home/merz/data/magfield/1900/1_new/2022x05x30_file_2_2000_CW.csv");
36 files[0].push_back("/home/merz/data/magfield/1900/1_new/2022x05x30_file_1_0_CW.csv");
37
38 //==== Z=3 ====//
39 // files[1].push_back("/home/merz/data/magfield/1900/3/2022x05x25_file_3_5775_CW.csv");
40 // files[1].push_back("/home/merz/data/magfield/1900/3/2022x05x25_file_2_2450_CW.csv");
41 // files[1].push_back("/home/merz/data/magfield/1900/3/2022x05x25_file_1_0_CW.csv");
42 files[1].push_back("/home/merz/data/magfield/1900/3/file_3_.csv");
43
44 //==== Z=5 ====//
45 files[2].push_back("/home/merz/data/magfield/1900/5_new/2022x05x30_file_3_4500_CW.csv");
46 files[2].push_back("/home/merz/data/magfield/1900/5_new/2022x05x30_file_2_3350_CW.csv");
47 files[2].push_back("/home/merz/data/magfield/1900/5_new/2022x05x30_file_1_0_CW.csv");
48
49 //==== Z=7 ====//
50 // files[3].push_back("/home/merz/data/magfield/1900/7/2022x05x25_file_2_3150_CW.csv");
51 // files[3].push_back("/home/merz/data/magfield/1900/7/2022x05x25_file_1_0_CW.csv");
52 files[3].push_back("/home/merz/data/magfield/1900/7/file_7_.csv");
53
54 //==== Z=9 ====//
55 files[4].push_back("/home/merz/data/magfield/1900/9_new/2022x05x30_file_4_4450_CW.csv");
56 files[4].push_back("/home/merz/data/magfield/1900/9_new/2022x05x29_file_3_1775_CW.csv");
57 files[4].push_back("/home/merz/data/magfield/1900/9_new/2022x05x29_file_2_0_CW.csv");
58 files[4].push_back("/home/merz/data/magfield/1900/9_new/2022x05x29_file_1_0_CW.csv");
59
60 //==== Z=11 ====//
61 // files[5].push_back("/home/merz/data/magfield/1900/11/2022x05x25_file_3_4150_CW.csv");
62 // files[5].push_back("/home/merz/data/magfield/1900/11/2022x05x25_file_2_3375_CW.csv");
63 // files[5].push_back("/home/merz/data/magfield/1900/11/2022x05x25_file_1_0_CW.csv");
64 files[5].push_back("/home/merz/data/magfield/1900/11/file_11_.csv");
65
66 //==== Z=13 ====//
67 files[6].push_back("/home/merz/data/magfield/1900/13_new/2022x05x29_file_5_5300_CW.csv");
68 files[6].push_back("/home/merz/data/magfield/1900/13_new/2022x05x29_file_4_2675_CW.csv");
69 files[6].push_back("/home/merz/data/magfield/1900/13_new/2022x05x29_file_3_2675_CW.csv");
70 files[6].push_back("/home/merz/data/magfield/1900/13_new/2022x05x29_file_2_25_CW.csv");
71 files[6].push_back("/home/merz/data/magfield/1900/13_new/2022x05x29_file_1_0_CW.csv");
72
73 //==== Z=15 ====//
74 // files[7].push_back("/home/merz/data/magfield/1900/15/2022x05x25_file_3_4500_CW.csv");
75 // files[7].push_back("/home/merz/data/magfield/1900/15/2022x05x25_file_2_3200_CW.csv");
76 // files[7].push_back("/home/merz/data/magfield/1900/15/2022x05x25_file_1_0_CW.csv");
77 files[7].push_back("/home/merz/data/magfield/1900/15/file_15_.csv");
78
79 //==== Z=17 ====//
80 files[8].push_back("/home/merz/data/magfield/1900/17_new/2022x05x29_file_3_3500_CW.csv");
81 files[8].push_back("/home/merz/data/magfield/1900/17_new/2022x05x29_file_2_3325_CW.csv");
82 files[8].push_back("/home/merz/data/magfield/1900/17_new/2022x05x29_file_1_0_CW.csv");
83
84 //==== Z=19 ====//
85 files[9].push_back("/home/merz/data/magfield/1900/19_new/2022x05x30_file_3_3275_CW.csv");
86 files[9].push_back("/home/merz/data/magfield/1900/19_new/2022x05x30_file_2_2000_CW.csv");
87 files[9].push_back("/home/merz/data/magfield/1900/19_new/2022x05x30_file_1_0_CW.csv");
88
89 //==== Z=21 ====//
90 files[10].push_back("/home/merz/data/magfield/1900/21_new/2022x05x29_file_2_3275_CW.csv");
91 files[10].push_back("/home/merz/data/magfield/1900/21_new/2022x05x29_file_1_0_CW.csv");
92
93 //==== Z=23 ====//
94 // files[11].push_back("/home/merz/data/magfield/1900/23/2022x05x25_file_2_3475_CW.csv");
95 // files[11].push_back("/home/merz/data/magfield/1900/23/2022x05x25_file_1_0_CW.csv");
96 files[11].push_back("/home/merz/data/magfield/1900/23/file_23_.csv");
97
98 //==== Z=25 ====//
99 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_10_5725_CW.csv");
100 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_9_3900_CW.csv");
101 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_8_2100_CW.csv");
102 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_7_1675_CW.csv");
103 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_6_1450_CW.csv");
104 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_5_1225_CW.csv");
105 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_4_1000_CW.csv");
106 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_3_550_CW.csv");
107 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_2_325_CW.csv");
108 files[12].push_back("/home/merz/data/magfield/1900/25/2022x05x26_file_1_0_CW.csv");
109
110 //==== Z=27 ====//
111 // files[13].push_back("/home/merz/data/magfield/1900/27/2022x05x25_file_2_3825_CW.csv");
112 // files[13].push_back("/home/merz/data/magfield/1900/27/2022x05x25_file_1_0_CW.csv");
113 files[13].push_back("/home/merz/data/magfield/1900/27/file_27_.csv");
114
115 //==== Z=29 ====//
116 files[14].push_back("/home/merz/data/magfield/1900/29/2022x05x26_file_4_6000_CW.csv");
117 files[14].push_back("/home/merz/data/magfield/1900/29/2022x05x26_file_3_5125_CW.csv");
118 files[14].push_back("/home/merz/data/magfield/1900/29/2022x05x26_file_2_3500_CW.csv");
119 files[14].push_back("/home/merz/data/magfield/1900/29/2022x05x26_file_1_0_CW.csv");
120
121 //==== Z=31 ====//
122 files[15].push_back("/home/merz/data/magfield/1900/31_new/2022x05x29_file_3_4725_CW.csv");
123 files[15].push_back("/home/merz/data/magfield/1900/31_new/2022x05x29_file_2_1800_CW.csv");
124 files[15].push_back("/home/merz/data/magfield/1900/31_new/2022x05x29_file_1_0_CW.csv");
125
126 //==== Z=33 ====//
127 files[16].push_back("/home/merz/data/magfield/1900/33/2022x05x26_file_4_4625_CW.csv");
128 files[16].push_back("/home/merz/data/magfield/1900/33/2022x05x26_file_3_4275_CW.csv");
129 files[16].push_back("/home/merz/data/magfield/1900/33/2022x05x26_file_2_2250_CW.csv");
130 files[16].push_back("/home/merz/data/magfield/1900/33/2022x05x26_file_1_0_CW.csv");
131
132 //==== Z=35 ====//
133 // files[17].push_back("/home/merz/data/magfield/1900/35/2022x05x24_file_3_4950_CW.csv");
134 // files[17].push_back("/home/merz/data/magfield/1900/35/2022x05x24_file_2_975_CW.csv");
135 // files[17].push_back("/home/merz/data/magfield/1900/35/2022x05x24_file_1_0_CW.csv");
136 files[17].push_back("/home/merz/data/magfield/1900/35/file_35_.csv");
137
138 //==== Z=37 ====//
139 files[18].push_back("/home/merz/data/magfield/1900/37_new/2022x05x29_file_3_5175_CW.csv");
140 files[18].push_back("/home/merz/data/magfield/1900/37_new/2022x05x29_file_2_2150_CW.csv");
141 files[18].push_back("/home/merz/data/magfield/1900/37_new/2022x05x29_file_1_0_CW.csv");
142
143 //==== Z=39 ====//
144 // files[19].push_back("/home/merz/data/magfield/1900/39_best/2022x05x24_file_3_3350_CW.csv");
145 // files[19].push_back("/home/merz/data/magfield/1900/39_best/2022x05x24_file_1_875_CW.csv");
146 // files[19].push_back("/home/merz/data/magfield/1900/39_best/2022x05x24_file_1_0_CW.csv");
147 files[19].push_back("/home/merz/data/magfield/1900/39/file_39_.csv");
148
149 //==== Z=41 ====//
150 files[20].push_back("/home/merz/data/magfield/1900/41_new/2022x05x29_file_3_5650_CW.csv");
151 files[20].push_back("/home/merz/data/magfield/1900/41_new/2022x05x29_file_2_2725_CW.csv");
152 files[20].push_back("/home/merz/data/magfield/1900/41_new/2022x05x29_file_1_0_CW.csv");
153
154 //==== Z=43 ====//
155 files[21].push_back("/home/merz/data/magfield/1900/43_new_updown/2022x05x29_file_5_5975_CW.csv");
156 files[21].push_back("/home/merz/data/magfield/1900/43_new_updown/2022x05x29_file_4_4350_CW.csv");
157 files[21].push_back("/home/merz/data/magfield/1900/43_new_updown/2022x05x29_file_3_3000_CW.csv");
158 files[21].push_back("/home/merz/data/magfield/1900/43_new_updown/2022x05x29_file_2_1450_CW.csv");
159 files[21].push_back("/home/merz/data/magfield/1900/43_new_updown/2022x05x29_file_1_0_CW.csv");
160
161 //==== Z=45 ====//
162 files[22].push_back("/home/merz/data/magfield/1900/45/2022x05x27_file_2_3825_CW.csv");
163 files[22].push_back("/home/merz/data/magfield/1900/45/2022x05x27_file_1_0_CW.csv");
164
165 //==== Z=47 ====//
166 files[23].push_back("/home/merz/data/magfield/1900/47_new/2022x05x29_file_2_3075_CW.csv");
167 files[23].push_back("/home/merz/data/magfield/1900/47_new/2022x05x29_file_1_0_CW.csv");
168
169 //==== Z=49 ====//
170 files[24].push_back("/home/merz/data/magfield/1900/49/2022x05x27_file_2_4200_CW.csv");
171 files[24].push_back("/home/merz/data/magfield/1900/49/2022x05x27_file_1_0_CW.csv");
172
173 //==== Z=51 ====//
174 files[25].push_back("/home/merz/data/magfield/1900/51_new/2022x05x29_file_2_3475_CW.csv");
175 files[25].push_back("/home/merz/data/magfield/1900/51_new/2022x05x28_file_1_0_CW.csv");
176
177 //==== Z=53 ====//
178 // files[26].push_back("/home/merz/data/magfield/1900/53/2022x05x27_file_3_5200_CW.csv");
179 // files[26].push_back("/home/merz/data/magfield/1900/53/2022x05x27_file_2_2600_CW.csv");
180 // files[26].push_back("/home/merz/data/magfield/1900/53/2022x05x27_file_1_0_CW.csv");
181 // files[26].push_back("/home/merz/data/magfield/1900/53/2022x05x27_file_1_1_CW.csv");
182 // files[26].push_back("/home/merz/data/magfield/1900/53/2022x05x27_file_1_2_CW.csv");
183 files[26].push_back("/home/merz/data/magfield/1900/53/file_53_.csv");
184
185 //==== Z=55 ====//
186 files[27].push_back("/home/merz/data/magfield/1900/55_new/2022x05x28_file_2_3100_CW.csv");
187 files[27].push_back("/home/merz/data/magfield/1900/55_new/2022x05x28_file_1_0_CW.csv");
188
189 //==== Z=57 ====//
190 files[28].push_back("/home/merz/data/magfield/1900/57/2022x05x27_file_3_4475_CW.csv");
191 files[28].push_back("/home/merz/data/magfield/1900/57/2022x05x27_file_2_3200_CW.csv");
192 files[28].push_back("/home/merz/data/magfield/1900/57/2022x05x27_file_1_0_CW.csv");
193
194 //==== Z=59 ====//
195 files[29].push_back("/home/merz/data/magfield/1900/59_new/2022x05x28_file_2_3275_CW.csv");
196 files[29].push_back("/home/merz/data/magfield/1900/59_new/2022x05x28_file_1_0_CW.csv");
197
198 //==== Z=61 ====//
199 files[30].push_back("/home/merz/data/magfield/1900/61/2022x05x27_file_2_4225_CW.csv");
200 files[30].push_back("/home/merz/data/magfield/1900/61/2022x05x27_file_1_0_CW.csv");
201
202 //==== Z=63 ====//
203 files[31].push_back("/home/merz/data/magfield/1900/63_new/2022x05x28_file_2_3325_CW.csv");
204 files[31].push_back("/home/merz/data/magfield/1900/63_new/2022x05x28_file_1_0_CW.csv");
205
206 //==== Z=65 ====//
207 files[32].push_back("/home/merz/data/magfield/1900/65/2022x05x27_file_2_3950_CW.csv");
208 files[32].push_back("/home/merz/data/magfield/1900/65/2022x05x27_file_1_0_CW.csv");
209
210 //==== Z=67 ====//
211 files[33].push_back("/home/merz/data/magfield/1900/67_new/2022x05x28_file_3_4825_CW.csv");
212 files[33].push_back("/home/merz/data/magfield/1900/67_new/2022x05x28_file_2_1825_CW.csv");
213 files[33].push_back("/home/merz/data/magfield/1900/67_new/2022x05x28_file_1_0_CW.csv");
214
215 //==== Z=69 ====//
216 files[34].push_back("/home/merz/data/magfield/1900/69/2022x05x27_file_2_3725_CW.csv");
217 files[34].push_back("/home/merz/data/magfield/1900/69/2022x05x27_file_1_0_CW.csv");
218
219 //==== Z=71 ====//
220 files[35].push_back("/home/merz/data/magfield/1900/71_new/2022x05x28_file_3_5150_CW.csv");
221 files[35].push_back("/home/merz/data/magfield/1900/71_new/2022x05x28_file_2_2850_CW.csv");
222 files[35].push_back("/home/merz/data/magfield/1900/71_new/2022x05x28_file_1_0_CW.csv");
223
224 //==== Z=73 ====//
225 files[36].push_back("/home/merz/data/magfield/1900/73/2022x05x27_file_3_3475_CW.csv");
226 files[36].push_back("/home/merz/data/magfield/1900/73/2022x05x27_file_2_500_CW.csv");
227 files[36].push_back("/home/merz/data/magfield/1900/73/2022x05x27_file_1_0_CW.csv");
228
229 //==== Z=75 ====//
230 // files[37].push_back("/home/merz/data/magfield/1900/75/2022x05x22_file_2_5075_CW.csv");
231 // files[37].push_back("/home/merz/data/magfield/1900/75/2022x05x22_file_1_0_CW.csv");
232 // files[37].push_back("/home/merz/data/magfield/1900/75/2022x05x22_1900A_hall_75_CW.csv");
233 // files[37].push_back("/home/merz/data/magfield/1900/75/2022x05x22_1900A_hall_75_part2_CW.csv");
234 files[37].push_back("/home/merz/data/magfield/1900/75/file_75_.csv");
235
236 //==== Z=77 ====//
237 files[38].push_back("/home/merz/data/magfield/1900/77/2022x05x27_file_2_3175_CW.csv");
238 files[38].push_back("/home/merz/data/magfield/1900/77/2022x05x27_file_1_0_CW.csv");
239
240 //==== Z=79 ====//
241 // files[39].push_back("/home/merz/data/magfield/1900/79/2022x05x22_file_2_4625_CW.csv");
242 // files[39].push_back("/home/merz/data/magfield/1900/79/2022x05x22_file_1_0_CW.csv");
243 // files[39].push_back("/home/merz/data/magfield/1900/79/2022x05x22_1900A_hall_79_CW.csv");
244 // files[39].push_back("/home/merz/data/magfield/1900/79/2022x05x22_1900A_hall_79_part2_CW.csv");
245 files[39].push_back("/home/merz/data/magfield/1900/79/file_79_.csv");
246
247 //==== Z=81 ====//
248 files[40].push_back("/home/merz/data/magfield/1900/81/2022x05x28_file_3_5275_CW.csv");
249 files[40].push_back("/home/merz/data/magfield/1900/81/2022x05x28_file_2_2225_CW.csv");
250 files[40].push_back("/home/merz/data/magfield/1900/81/2022x05x28_file_1_0_CW.csv");
251
252 //==== Z=83 ====//
253 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_file_4_6200_CW.csv");
254 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_file_3_2650_CW.csv");
255 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_file_2_175_CW.csv");
256 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_file_1_0_CW.csv");
257 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_1900A_hall_83_CW.csv");
258 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_1900A_hall_83_1_CW.csv");
259 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_1900A_hall_83_2_CW.csv");
260 // files[41].push_back("/home/merz/data/magfield/1900/83/2022x05x22_1900A_hall_83_3_CW.csv");
261 files[41].push_back("/home/merz/data/magfield/1900/83/file_83_.csv");
262
263 //==== Z=85 ====//
264 files[42].push_back("/home/merz/data/magfield/1900/85_repeat/2022x05x28_file_3_4975_CW.csv");
265 files[42].push_back("/home/merz/data/magfield/1900/85_repeat/2022x05x28_file_2_2000_CW.csv");
266 files[42].push_back("/home/merz/data/magfield/1900/85_repeat/2022x05x28_file_1_0_CW.csv");
267
268 //==== Z=87 ====//
269 files[43].push_back("/home/merz/data/magfield/1900/87_new/2022x05x29_file_2_3325_CW.csv");
270 files[43].push_back("/home/merz/data/magfield/1900/87_new/2022x05x29_file_1_0_CW.csv");
271
272 //==== Z=89 ====//
273 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_6_5175_CW.csv");
274 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_5_4975_CW.csv");
275 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_4_4775_CW.csv");
276 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_3_4575_CW.csv");
277 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_2_3175_CW.csv");
278 files[44].push_back("/home/merz/data/magfield/1900/89_new/2022x05x28_file_1_0_CW.csv");
279
280 //==== Z=91 ====//
281 files[45].push_back("/home/merz/data/magfield/1900/91_new/2022x05x28_file_3_4400_CW.csv");
282 files[45].push_back("/home/merz/data/magfield/1900/91_new/2022x05x28_file_2_3525_CW.csv");
283 files[45].push_back("/home/merz/data/magfield/1900/91_new/2022x05x28_file_1_0_CW.csv");
284
285 //==== Z=93 ====//
286 files[46].push_back("/home/merz/data/magfield/1900/93_new/2022x05x28_file_3_5525_CW.csv");
287 files[46].push_back("/home/merz/data/magfield/1900/93_new/2022x05x28_file_2_2700_CW.csv");
288 files[46].push_back("/home/merz/data/magfield/1900/93_new/2022x05x28_file_1_0_CW.csv");
289
290 BmnFieldParser* parser = new BmnFieldParser();
291
292 vector<Int_t> zOfOldSensor = {3, 7, 11, 15, 23, 27, 35, 39, 53, 75, 77, 79, 83};
293
294 TString out = "FieldMap_1900.root";
295 TFile* f = new TFile(out.Data(), "RECREATE");
296 TTree* t = new TTree("Field_map", "data from ASCII");
297
298 BmnFieldPoint pointToWrite;
299 t->Branch("field_points", &pointToWrite, "X/D:Y/D:Z/D:Bx/D:By/D:Bz/D");
300
301 for (Int_t i = 0; i < nZ; ++i) {
302 // TString out = Form("field_1900_%d.root", zCoord[i]);
303 // TFile* f = new TFile(out.Data(), "RECREATE");
304 // TTree* t = new TTree("Field_map", "data from ASCII");
305
306 // BmnFieldPoint pointToWrite;
307 // t->Branch("field_points", &pointToWrite, "X/D:Y/D:Z/D:Bx/D:By/D:Bz/D");
308 map<pair<Int_t, Int_t>, BmnFieldPoint> all_field_points;
309 // Bool_t isNewSensor = kFALSE;
310 for (auto file : files[i]) {
311 // if (find(zOfOldSensor.begin(), zOfOldSensor.end(), zCoord[i]) == zOfOldSensor.end()) isNewSensor = kTRUE;
312 vector<BmnFieldPoint> fieldPoints = parser->ParseData(file);
313 for (auto pnt : fieldPoints) {
314 all_field_points.insert(pair<pair<Int_t, Int_t>, BmnFieldPoint>(
315 pair<Int_t, Int_t>(Int_t(pnt.GetX()), Int_t(pnt.GetY())), pnt));
316 }
317 }
318
319 const Double_t By_ped = -0.158; // pedestal for By
320
321 Int_t nPoints = 0;
322 for (auto pnt : all_field_points) {
323
324 Double_t xFrame = pnt.second.GetX();
325 Double_t yFrame = pnt.second.GetY();
326 Double_t zFrame = pnt.second.GetZ();
327
328 if (xFrame < 0 || xFrame > 3000 || yFrame < 0 || yFrame > 6000) {
329 continue;
330 }
331
332 Double_t BxFrame = pnt.second.GetBx();
333 Double_t ByFrame = pnt.second.GetBy();
334 Double_t BzFrame = pnt.second.GetBz();
335
336 Double_t xBmn = (1.002 * (-xFrame) + 1448.218) * 0.1;
337 Double_t yBmn = (zFrame * 10 - 386.73) * 0.1;
338 Double_t zBmn = (1.002 * yFrame - 1621.615) * 0.1;
339
340 // Convert voltage from channels of measurement frame coordinates to magnetic field components of hall
341 // sensor coordinates
342 Double_t ByHall = ConvertVoltageToFeild_X(+BxFrame);
343 Double_t BxHall = ConvertVoltageToFeild_Y(+ByFrame);
344 Double_t BzHall = ConvertVoltageToFeild_Z(+BzFrame);
345
346 // Convert field components of hall sensor coordinates to field components of BM@N coordinates
347 Double_t BxBmn = -BxHall;
348 Double_t ByBmn = -BzHall;
349 Double_t BzBmn = -ByHall;
350
351 BmnFieldPoint tmpPnt(xBmn, yBmn, zBmn, BxBmn, ByBmn - By_ped, BzBmn);
352 pointToWrite = tmpPnt;
353 nPoints++;
354 t->Fill();
355 }
356
357 // f->Write();
358 }
359 f->Write();
360}
361
362Double_t ConvertVoltageToFeild(Double_t U)
363{
364 // U - voltage in V
365 // return field in kGs
366 Double_t k = 136.5;
367 Double_t b = 23.16;
368 Double_t scale = 1000;
369
370 return (k * U * scale + b) / 1000;
371}
372
373Double_t ConvertVoltageToFeild_X(Double_t U)
374{
375 // U - voltage in V
376 // return field in kGs
377 Double_t k = 136.652;
378 Double_t b = -11.782;
379 Double_t scale = 1000;
380 return (k * U * scale + b) / 1000;
381}
382
383Double_t ConvertVoltageToFeild_Y(Double_t U)
384{
385 // U - voltage in V
386 // return field in kGs
387 Double_t k = 139.666;
388 Double_t b = -2.07;
389 Double_t scale = 1000;
390 return (k * U * scale + b) / 1000;
391}
392
393Double_t ConvertVoltageToFeild_Z(Double_t U)
394{
395 // U - voltage in V
396 // return field in kGs
397 Double_t k = 144.824;
398 Double_t b = 152.728; // 50.0; //изначальное значеие из калибровок 152.728, но я изменил его для лучшей
399 // совместимости со старой картой
400 Double_t scale = 1000;
401 return (k * U * scale + b) / 1000;
402}
i
int i
Definition P4_F32vec4.h:22
min
friend F32vec4 min(const F32vec4 &a, const F32vec4 &b)
Definition P4_F32vec4.h:30
max
friend F32vec4 max(const F32vec4 &a, const F32vec4 &b)
Definition P4_F32vec4.h:31
Int_t
int Int_t
Float_t
float Float_t
Double_t
double Double_t
R__ADD_INCLUDE_PATH
#define R__ADD_INCLUDE_PATH(PATH)
BmnFieldParser::ParseData
vector< BmnFieldPoint > ParseData(TString file_name, Bool_t shift=kTRUE)
Definition BmnFieldParser.cxx:9
TString::Data
const char * Data() const
TTree::Fill
virtual Int_t Fill()
TTree::Branch
virtual Int_t Branch(const char *folder, Int_t bufsize=32000, Int_t splitlevel=99)
ConvertVoltageToFeild
Double_t ConvertVoltageToFeild(Double_t U)
Definition convertFieldCSV2ROOT_1900.C:362
ConvertVoltageToFeild_X
Double_t ConvertVoltageToFeild_X(Double_t U)
Definition convertFieldCSV2ROOT_1900.C:373
ConvertVoltageToFeild_Y
Double_t ConvertVoltageToFeild_Y(Double_t U)
Definition convertFieldCSV2ROOT_1900.C:383
ConvertVoltageToFeild_Z
Double_t ConvertVoltageToFeild_Z(Double_t U)
Definition convertFieldCSV2ROOT_1900.C:393
convertFieldCSV2ROOT_1900
void convertFieldCSV2ROOT_1900()
Definition convertFieldCSV2ROOT_1900.C:23
Coordinate_info_t
Definition convertFieldCSV2ROOT_1900.C:10