]>
Commit | Line | Data |
---|---|---|
7d5a9359 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | /* $Id$ */ | |
17 | ||
7d5a9359 | 18 | #include <TCanvas.h> |
9edefa04 | 19 | #include <TFile.h> |
7d5a9359 | 20 | #include <TGraph.h> |
a1e17193 | 21 | #include <TH1F.h> |
9edefa04 | 22 | #include <TMath.h> |
7d5a9359 | 23 | #include <TMultiGraph.h> |
24 | #include <TNtuple.h> | |
25 | #include <TObjString.h> | |
26 | #include <TString.h> | |
27 | #include <TStyle.h> | |
9edefa04 | 28 | |
29 | #include "AliFieldReader.h" | |
30 | #include "AliMagFMaps.h" | |
31 | ||
7d5a9359 | 32 | ClassImp(AliFieldReader) |
33 | ||
34 | ||
6772b8d8 | 35 | |
7d5a9359 | 36 | //_______________________________________________________________________ |
6772b8d8 | 37 | AliFieldReader::AliFieldReader(): |
38 | fField(0), | |
39 | fMap(0), | |
40 | fCatalogue(0), | |
41 | fHtmlMain(0), | |
42 | fStepSize(0.), | |
43 | fZStart(1383.), | |
44 | fDd(0.08), | |
45 | fDz(0.064), | |
46 | fPolarity(1.), | |
47 | fCatalogueName("goodfiles.list") | |
7d5a9359 | 48 | { |
14b22234 | 49 | // |
50 | // Constructor | |
51 | // | |
6772b8d8 | 52 | } |
53 | ||
54 | AliFieldReader::AliFieldReader(const AliFieldReader& reader): | |
55 | TObject(reader), | |
56 | fField(0), | |
57 | fMap(0), | |
58 | fCatalogue(0), | |
59 | fHtmlMain(0), | |
60 | fStepSize(0.), | |
61 | fZStart(0.), | |
62 | fDd(0.), | |
63 | fDz(0.), | |
64 | fPolarity(0.), | |
65 | fCatalogueName(0) | |
66 | { | |
67 | // | |
68 | // Constructor | |
69 | // | |
70 | reader.Copy(*this); | |
7d5a9359 | 71 | } |
72 | ||
73 | //_______________________________________________________________________ | |
74 | AliFieldReader::~AliFieldReader() | |
75 | { | |
76 | // | |
77 | // Destructor | |
78 | // | |
79 | } | |
80 | ||
81 | //_______________________________________________________________________ | |
82 | void AliFieldReader::Init() | |
83 | { | |
14b22234 | 84 | // |
85 | // Initialize the reader | |
86 | // | |
7d5a9359 | 87 | // Calculated map |
88 | fField = new AliMagFMaps("Maps","Maps", 2, 1., 10., 2); | |
89 | // Catalogue | |
90 | fCatalogue = fopen(fCatalogueName, "r"); | |
91 | ||
92 | // HTML | |
93 | fHtmlMain = fopen("bmap.html", "w"); | |
94 | MakeHtmlHeaderMain(fHtmlMain); | |
95 | ||
96 | } | |
97 | ||
98 | ||
99 | void AliFieldReader::ReadMap() | |
100 | { | |
14b22234 | 101 | // |
102 | // Read the measured dipole field map | |
103 | // | |
7d5a9359 | 104 | |
105 | Float_t zA[450], bxzA[200], byzA[200], bzzA[200], bxzcA[200], byzcA[200], bzzcA[200]; | |
106 | Float_t yA[450], bxyA[200], byyA[200], bzyA[200], bxycA[200], byycA[200], bzycA[200]; | |
107 | Float_t xA[450], bxxA[200], byxA[200], bzxA[200], bxxcA[200], byxcA[200], bzxcA[200]; | |
108 | ||
109 | Char_t sLine[255]; | |
110 | Char_t fLine[255]; | |
111 | ||
112 | Float_t xpos, ypos, zpos; | |
113 | Int_t iboxpos; | |
114 | ||
115 | Float_t h[3], x[3], b[3]; | |
116 | Float_t temp; | |
117 | Float_t xnt[17]; | |
118 | Int_t ires; | |
119 | ||
120 | Int_t iret, ipos; | |
121 | Int_t ic115 = 0; | |
122 | Float_t dx = 0.; | |
123 | Float_t dy = 0.; | |
124 | ||
125 | Init(); | |
126 | // Read register | |
127 | ReadRegisterMap(); | |
128 | // n-tuple | |
129 | fMap = new TNtuple("Field Map", "Map", | |
130 | "x:y:z:ix:iy:iz:bx:by:bz:bxc:byc:bzc:ifile:ireg:irot:temp:cal", 4000); | |
131 | ||
132 | // | |
133 | // Loop over files | |
134 | // | |
135 | Int_t ifiles = 0; | |
136 | // File catalogue | |
137 | while ((fgets(fLine, 255, fCatalogue)) != NULL && ifiles <= 2000) { | |
138 | ||
139 | if (strncmp(fLine,"#HEADER", 7) == 0) { | |
140 | iret = sscanf(&fLine[7],"%f", &fZStart); | |
141 | continue; | |
142 | } | |
143 | ifiles++; | |
144 | ||
145 | // if (ifiles != 87) continue; | |
146 | char fileName[32]; | |
147 | iret = sscanf(fLine, "%s", fileName); | |
148 | printf("Reading File %s\n", fileName); | |
149 | ||
150 | // Get run name | |
151 | TString* tsFile = new TString(fileName); | |
152 | TObjArray * tokens = tsFile->Tokenize("."); | |
153 | const char* runName = (((TObjString*) tokens->At(0))->GetString()).Data(); | |
154 | FILE* file = fopen(fileName, "r"); | |
155 | ||
156 | ||
157 | Float_t bdl = 0.; | |
158 | Float_t bdlc = 0.; | |
159 | Int_t iz = 0; | |
160 | Int_t izz = 0; | |
161 | Int_t iy = 0; | |
162 | Int_t ix = 0; | |
163 | Int_t iheader = 0; | |
164 | Float_t stepsz = 0.; | |
165 | // Graphs go here | |
166 | TMultiGraph* bxzmg = new TMultiGraph("bxzmg", "B_{x}"); | |
167 | TMultiGraph* byzmg = new TMultiGraph("byzmg", "B_{y}"); | |
168 | TMultiGraph* bzzmg = new TMultiGraph("bzzmg", "B_{z}"); | |
169 | while ((fgets(sLine, 255, file)) != NULL) { | |
170 | // read z-position | |
171 | ||
172 | if (strncmp(sLine," Current", 8) == 0) { | |
173 | Float_t current; | |
174 | iret = sscanf(&sLine[9],"%f", ¤t); | |
175 | printf("Current %f \n", current); | |
176 | ||
177 | fPolarity = current / 6000.; | |
178 | } | |
179 | if (strncmp(sLine," z", 2) == 0) { | |
180 | Float_t zmin, zmax; | |
181 | Int_t nsteps; | |
182 | iret = sscanf(&sLine[3],"%f %f %d", &zmin, &zmax, &nsteps); | |
183 | printf("zmin zmax %13.3f %13.3f %13.3f\n", | |
184 | zmin, zmax, TMath::Abs(zmax - zmin)/Float_t(nsteps)); | |
185 | zmax = TMath::Max(zmin, zmax); | |
186 | stepsz = TMath::Abs(zmax - zmin)/Float_t(nsteps); | |
187 | } | |
188 | ||
189 | if (strncmp(sLine," X position", 11) == 0) | |
190 | { | |
191 | TString string; | |
192 | TString* tsLine = new TString(sLine); | |
193 | TObjArray * tokens = tsLine->Tokenize("="); | |
194 | string = ((TObjString*) tokens->At(1))->GetString(); | |
195 | iret = sscanf(string.Data(), "%f", &xpos); | |
196 | string = ((TObjString*) tokens->At(2))->GetString(); | |
197 | iret = sscanf(string.Data(), "%f", &ypos); | |
198 | string = ((TObjString*) tokens->At(3))->GetString(); | |
199 | iret = sscanf(string.Data(), "%d", &iboxpos); | |
200 | ||
201 | ||
202 | printf("This file is for x = %13.3f y = %13.3f Box Position %1d\n", xpos, ypos, iboxpos); | |
203 | iheader++; | |
204 | ||
205 | } | |
206 | ||
207 | if (strncmp(sLine," R Z-pos",8) == 0) | |
208 | { | |
209 | iret = sscanf(&sLine[8],"%e", &zpos); | |
210 | ipos = 0; | |
211 | ic115 = 0; | |
212 | izz++; | |
213 | } | |
214 | ||
215 | ||
216 | if (strncmp(sLine,"C",1) == 0) { | |
217 | ipos++; | |
218 | Int_t ireg; | |
219 | iret = sscanf(&sLine[2],"%d %f %f %f %f %d", &ireg, &h[0], &h[1], &h[2], &temp, &ires); | |
220 | // | |
221 | // fix for address 115 | |
222 | // | |
223 | if (ireg == 115) ic115++; | |
224 | if (ic115 == 2 && ireg == 115) ireg = 119; | |
225 | ||
226 | if (iheader == 1) { | |
227 | ||
228 | Float_t bx = 0., by = 0., bz = 0.; | |
229 | Int_t jx = fRegMap[ireg][0]; | |
230 | Int_t jy = fRegMap[ireg][1]; | |
231 | Int_t jz = fRegMap[ireg][2]; | |
232 | ||
233 | switch (iboxpos) { | |
234 | case 0: | |
235 | bx = h[1]; | |
236 | by = h[2]; | |
237 | dx = -0.36 + jx * fDd; | |
238 | dy = jy * fDd; | |
239 | break; | |
240 | case 1: | |
241 | bx = -h[2]; | |
242 | by = h[1]; | |
243 | dx = 0.36 + jy * fDd; | |
244 | dy = -(-0.36 + jx * fDd); | |
245 | break; | |
246 | case 2: | |
247 | bx = -h[1]; | |
248 | by = -h[2]; | |
249 | dx = 0.36 - jx * fDd; | |
250 | dy = -jy * fDd; | |
251 | break; | |
252 | case 3: | |
253 | bx = h[2]; | |
254 | by = -h[1]; | |
255 | dx = -0.36 - jy * fDd; | |
256 | dy = -(0.36 - jx * fDd); | |
257 | break; | |
258 | } // switch | |
259 | ||
260 | bz = h[0]; | |
261 | if (jz == 0) { | |
262 | bz = -bz; | |
263 | if (iboxpos == 1 || iboxpos == 3) { | |
264 | by = -by; | |
265 | } else { | |
266 | bx = -bx; | |
267 | } | |
268 | } | |
269 | ||
270 | ||
271 | Float_t dz = (jz == 0)? fDz/2. : - fDz/2.; | |
272 | dz *= 100.; | |
273 | ||
274 | ||
275 | Float_t xc = xpos + dx; | |
276 | Float_t yc = ypos + dy; | |
277 | ||
278 | ||
279 | x[0] = - xc * 100.; | |
280 | x[1] = + yc * 100.; | |
281 | x[2] = - (-zpos * 100. + fZStart + dz); | |
282 | ||
283 | fField->Field(x, b); | |
284 | b[0] *= fPolarity; | |
285 | b[1] *= fPolarity; | |
286 | b[2] *= fPolarity; | |
287 | ||
288 | xnt[ 0] = xc; | |
289 | xnt[ 1] = yc; | |
290 | xnt[ 2] = -x[2] / 100.; | |
291 | xnt[ 3] = Float_t (jx); | |
292 | xnt[ 4] = Float_t (jy); | |
293 | xnt[ 5] = Float_t (jz); | |
294 | xnt[ 6] = bx; | |
295 | xnt[ 7] = by; | |
296 | xnt[ 8] = bz; | |
297 | xnt[ 9] = b[0]/10.; | |
298 | xnt[10] = b[1]/10.; | |
299 | xnt[11] = b[2]/10.; | |
300 | xnt[12] = Float_t (ifiles); | |
301 | xnt[13] = Float_t (ireg); | |
302 | xnt[14] = Float_t(iboxpos); | |
303 | xnt[15] = temp; | |
304 | xnt[16] = Float_t(ires); | |
305 | ||
306 | fMap->Fill(xnt); | |
307 | ||
308 | ||
309 | // | |
310 | // Calculated field | |
311 | ||
312 | if (jy != -1 && jz == 1 && jy == 0 && izz == 40){ | |
313 | x[1] = x[1] + 10.; | |
314 | fField->Field(x, b); | |
315 | yA [jx] = yc; | |
316 | bxyA[jx] = bx; | |
317 | byyA[jx] = by; | |
318 | bzyA[jx] = bz; | |
319 | bxycA[jx] = b[0] / 10.; | |
320 | byycA[jx] = b[1] / 10.; | |
321 | bzycA[jx] = b[2] / 10.; | |
322 | iy++; | |
323 | } // if | |
324 | ||
325 | if (jy != -1 && jz == 1 && jy == 0 && izz == 1){ | |
326 | x[1] = x[1] + 10.; | |
327 | fField->Field(x, b); | |
328 | xA [jx] = xc; | |
329 | bxxA[jx] = bx; | |
330 | byxA[jx] = by; | |
331 | bzxA[jx] = bz; | |
332 | bxxcA[jx] = b[0] / 10.; | |
333 | byxcA[jx] = b[1] / 10.; | |
334 | bzxcA[jx] = b[2] / 10.; | |
335 | ix++; | |
336 | } // if | |
337 | ||
338 | ||
339 | if (ireg == 181) { | |
340 | fField->Field(x, b); | |
341 | // printf("Field %f %f %f %f %f %f \n", x[0], x[1], x[2], b[0], b[1] , b[2]); | |
342 | ||
343 | bdl += stepsz * bx; | |
344 | bdlc += stepsz * b[0]; | |
345 | ||
346 | zA [iz] = x[2]; | |
347 | bxzA[iz] = bx; | |
348 | byzA[iz] = by; | |
349 | bzzA[iz] = bz; | |
350 | bxzcA[iz] = b[0] / 10.; | |
351 | byzcA[iz] = b[1] / 10.; | |
352 | bzzcA[iz] = b[2] / 10.; | |
353 | iz++; | |
354 | } | |
355 | } // if 1st header | |
356 | } // if C | |
357 | } // next line | |
358 | gStyle->SetOptStat(0); | |
359 | char title[128]; | |
360 | sprintf(title, "File#: %5d, X = %13.2f m , Y = %13.2f m, Box Orientation: %2d", ifiles, xpos, ypos, iboxpos); | |
361 | TCanvas* c2 = new TCanvas("c2", title, 1200, 800); | |
362 | c2->Divide(2,2); | |
363 | c2->cd(1); | |
364 | TGraph* bxg = new TGraph(iz, zA, bxzA); | |
365 | TGraph* bxcg = new TGraph(iz, zA, bxzcA); | |
366 | bxcg->SetLineColor(2); | |
367 | bxzmg->Add(bxg); | |
368 | bxzmg->Add(bxcg); | |
369 | bxzmg->Draw("lA"); | |
370 | bxzmg->GetHistogram()->SetXTitle("z[m]"); | |
371 | bxzmg->GetHistogram()->SetYTitle("B_{x} [T]"); | |
372 | bxzmg->Draw("lA"); | |
373 | ||
374 | c2->cd(2); | |
375 | TGraph* byg = new TGraph(iz, zA, byzA); | |
376 | TGraph* bycg = new TGraph(iz, zA, byzcA); | |
377 | bycg->SetLineColor(2); | |
378 | byzmg->Add(byg); | |
379 | byzmg->Add(bycg); | |
380 | byzmg->Draw("lA"); | |
381 | byzmg->GetHistogram()->SetXTitle("z[m]"); | |
382 | byzmg->GetHistogram()->SetYTitle("B_{y} [T]"); | |
383 | byzmg->Draw("lA"); | |
384 | ||
385 | c2->cd(3); | |
386 | TGraph* bzg = new TGraph(iz, zA, bzzA); | |
387 | TGraph* bzcg = new TGraph(iz, zA, bzzcA); | |
388 | bzcg->SetLineColor(2); | |
389 | bzzmg->Add(bzg); | |
390 | bzzmg->Add(bzcg); | |
391 | bzzmg->Draw("lA"); | |
392 | bzzmg->GetHistogram()->SetXTitle("z[m]"); | |
393 | bzzmg->GetHistogram()->SetYTitle("B_{z} [T]"); | |
394 | bzzmg->Draw("lA"); | |
395 | ||
396 | c2->Update(); | |
397 | char pictFile[64]; | |
398 | sprintf(pictFile, "%s.gif", runName); | |
399 | c2->SaveAs(pictFile); | |
400 | // | |
401 | // Html generation | |
402 | // | |
403 | char htmlFile[64]; | |
404 | sprintf(htmlFile, "%s.html", runName); | |
405 | FILE* chtml = fopen(htmlFile, "w"); | |
406 | MakeHtmlHeaderPict(chtml); | |
407 | MakeHtmlPict(chtml, pictFile); | |
408 | MakeHtmlTableEntry(fHtmlMain, fileName, htmlFile, xpos, ypos, iboxpos, bdl, ifiles); | |
409 | ||
410 | // | |
411 | // | |
412 | printf("Bdl [Tm] %f %f \n", 2. * bdl, 2 * bdlc / 10.); | |
413 | } // files | |
414 | MakeHtmlTrailor(fHtmlMain); | |
415 | TFile* out = new TFile("fmap.root", "recreate"); | |
416 | fMap->Write(); | |
417 | out->Close(); | |
418 | } | |
419 | ||
420 | void AliFieldReader::ReadMapSolenoid(){ | |
14b22234 | 421 | // |
422 | // Read map for solenoid measurement | |
423 | // | |
7d5a9359 | 424 | Float_t phiA[450], bzPhiA[200], brPhiA[200], btPhiA[200], bbPhiA[200]; |
425 | Float_t bzcPhiA[200], brcPhiA[200], btcPhiA[200], bbcPhiA[200]; | |
426 | Char_t sLine[255]; | |
427 | Char_t fLine[255]; | |
428 | ||
429 | Float_t zpos, phipos, skewing, temp; | |
430 | Int_t ical; | |
431 | Float_t zmin, zmax; | |
432 | Float_t h[3], x[3], b[3]; | |
433 | Int_t iret; | |
434 | Int_t ipos = 0; | |
435 | ||
436 | Init(); | |
437 | ReadRegisterMapSolenoid(); | |
438 | fMap = new TNtuple("Field Map", "Map", | |
439 | "r:phi:z:br:bt:bz:brc:btc:bzc:ifile:ireg:temp:cal:arm", 4000); | |
440 | ||
441 | // | |
442 | // Loop over files | |
443 | // | |
444 | Int_t ifiles = 0; | |
445 | // File catalogue | |
446 | while ((fgets(fLine, 255, fCatalogue)) != NULL && ifiles <= 2000) { | |
447 | ||
448 | if (strncmp(fLine,"#HEADER", 7) == 0) { | |
449 | iret = sscanf(&fLine[7],"%f", &fZStart); | |
450 | continue; | |
451 | } | |
452 | ifiles++; | |
453 | ||
454 | char fileName[32]; | |
455 | iret = sscanf(fLine, "%s", fileName); | |
456 | printf("Reading File %s\n", fileName); | |
457 | ||
458 | // Get run name | |
459 | TString* tsFile = new TString(fileName); | |
460 | TObjArray * tokens = tsFile->Tokenize("."); | |
461 | Int_t n = tokens->GetEntries(); | |
462 | char* runName = new char[256]; | |
463 | sprintf(runName, "%s", (((TObjString*) tokens->At(0))->GetString()).Data()); | |
464 | if (n > 2) { | |
465 | for (Int_t i = 1; i < n-1; i++) | |
466 | { | |
467 | sprintf(runName, "%s.%s", | |
468 | runName, (((TObjString*) tokens->At(i))->GetString()).Data()); | |
469 | } | |
470 | } | |
471 | FILE* file = fopen(fileName, "r"); | |
472 | ||
473 | ||
474 | Float_t bdl = 0.; | |
475 | Float_t bdlc = 0.; | |
476 | Int_t iA = 0; | |
477 | Int_t izz = 0; | |
478 | Int_t iphi = 0; | |
479 | // Graphs go here | |
480 | TMultiGraph* bxzmg = new TMultiGraph("bxzmg", "B_{z}"); | |
481 | TMultiGraph* byzmg = new TMultiGraph("byzmg", "B_{r}"); | |
482 | TMultiGraph* bzzmg = new TMultiGraph("bzzmg", "B_{t}"); | |
483 | TMultiGraph* bbmg = new TMultiGraph("bbmg", "|B|"); | |
484 | ||
485 | while ((fgets(sLine, 255, file)) != NULL) { | |
486 | if (strncmp(sLine," z", 2) == 0) { | |
487 | Int_t nsteps; | |
488 | iret = sscanf(&sLine[3],"%f %f %d", &zmin, &zmax, &nsteps); | |
489 | printf("zmin zmax %13.3f %13.3f %13.3f\n", | |
490 | zmin, zmax, TMath::Abs(zmax - zmin)/Float_t(nsteps)); | |
491 | } | |
492 | if (strncmp(sLine," R\tPOSITION NUMBER", 18) == 0) | |
493 | { | |
494 | // | |
495 | // Current z-position | |
496 | TString string; | |
497 | TString* tsLine = new TString(sLine); | |
498 | TObjArray * tokens = tsLine->Tokenize("="); | |
499 | string = ((TObjString*) tokens->At(1))->GetString(); | |
500 | iret = sscanf(string.Data(), "%f", &zpos); | |
501 | printf("POSITION NUMBER Z: %f\n", zpos); | |
502 | izz ++; | |
503 | delete tsLine; | |
504 | } | |
505 | ||
506 | ||
507 | if (strncmp(sLine," SKEWING ON Z:", 14) == 0) | |
508 | { | |
509 | // | |
510 | // Skew in z | |
511 | iret = sscanf(&sLine[14],"%f", &skewing); | |
512 | printf("SKEWING ON Z: %f\n", skewing); | |
513 | } | |
514 | ||
515 | ||
516 | if (strncmp(sLine,"Phi", 3) == 0) | |
517 | { | |
518 | Float_t phiStart, phiStop; | |
519 | ||
520 | iret = sscanf(&sLine[3],"%f %f", &phiStart, &phiStop); | |
521 | ||
522 | printf("phiStart phiStop %f %f\n", phiStart, phiStop); | |
523 | ||
524 | } | |
525 | ||
526 | ||
527 | if (strncmp(sLine," R\tPhi-Angle",12) == 0) | |
528 | { | |
529 | iret = sscanf(&sLine[12],"%e", &phipos); | |
530 | ipos = 0; | |
531 | iphi++; | |
532 | } | |
533 | ||
534 | ||
535 | if (strncmp(sLine,"C",1) == 0) { | |
536 | ||
537 | Int_t ireg; | |
538 | iret = sscanf(&sLine[2],"%d %f %f %f %f %d", &ireg, &h[0], &h[1], &h[2], &temp, &ical); | |
539 | ipos++; | |
540 | ||
541 | Int_t ir = fRegMap[ireg][0]; | |
542 | Int_t ia = fRegMap[ireg][1]; | |
543 | Float_t rpos = 0.; | |
544 | ||
545 | ||
546 | if (ia == 0) { | |
547 | rpos = 0.2295 + ir * 0.16; | |
548 | } else { | |
549 | if (ireg == 81) rpos = 0.2295; | |
550 | if (ireg == 59) rpos = 1.0295; | |
551 | if (ireg == 142) rpos = 2.1495; | |
552 | if (ireg == 180) rpos = 3.1095; | |
553 | if (ireg == 69) rpos = 4.2295; | |
554 | ||
555 | // if (ireg == 55) rpos = 0.2295; | |
556 | //if (ireg == 195) rpos = 1.0295; | |
557 | //if (ireg == 129) rpos = 2.1495; | |
558 | //if (ireg == 167) rpos = 3.1095; | |
559 | //if (ireg == 142) rpos = 4.2295; | |
560 | } | |
561 | ||
562 | Float_t phi = phipos; | |
563 | if (ia == 1) { | |
564 | phi += 180.; | |
565 | if (phi > 360.) phi -= 360.; | |
566 | } | |
567 | ||
568 | phi = - phi * TMath::Pi() / 180.; | |
569 | Float_t xpos = rpos * TMath::Cos(phi); | |
570 | Float_t ypos = rpos * TMath::Sin(phi); | |
571 | x[0] = - xpos * 100.; | |
572 | x[1] = ypos * 100.; | |
573 | x[2] = -400. + zpos * 100.; | |
574 | ||
575 | fField->Field(x, b); | |
576 | Float_t phi0 = TMath::Pi() - phi; | |
577 | ||
578 | Float_t brc = b[0] * TMath::Cos(phi0) + b[1] * TMath::Sin(phi0); | |
579 | Float_t btc = - b[0] * TMath::Sin(phi0) + b[1] * TMath::Cos(phi0); | |
580 | Float_t bzc = b[2]; | |
581 | ||
582 | brc /= 10.; | |
583 | btc /= 10.; | |
584 | bzc /= 10.; | |
585 | ||
586 | fMap->Fill(rpos, -phi, -(615.5 - fZStart) / 100. + zpos, h[2], -h[1], h[0], brc, btc, bzc, ifiles, ireg, temp, Float_t(ical), Int_t(ia)); | |
587 | ||
588 | if (ireg == 174) { | |
589 | printf("Field (Bx, By, Bz) at position %d: %13.3f %13.3f %13.3f %13.3f %13.3f %13.3f %5d\n", | |
590 | ipos, zpos, rpos, phipos, h[0], h[1], h[2], ia); | |
591 | if (izz == 1) { | |
592 | phiA[iA] = phi; | |
593 | bzPhiA[iA] = -h[0]; | |
594 | brPhiA[iA] = h[2]; | |
595 | btPhiA[iA] = -h[1]; | |
596 | bbPhiA[iA] = TMath::Sqrt(h[0] * h[0] + h[1] * h[1] + h[2] * h[2]); | |
597 | bzcPhiA[iA] = bzc; | |
598 | brcPhiA[iA] = brc; | |
599 | btcPhiA[iA] = btc; | |
600 | bbcPhiA[iA] = TMath::Sqrt(brc * brc + bzc * bzc + btc * btc); | |
601 | iA++; | |
602 | } | |
603 | } | |
604 | } // if R | |
605 | } // next line | |
606 | ||
607 | gStyle->SetOptStat(0); | |
608 | char title[128]; | |
609 | sprintf(title, "Z = %13.2f m , Phi = %13.2f m", zpos, phipos); | |
610 | ||
611 | ||
612 | TCanvas* c2 = new TCanvas("c2", title, 1200, 800); | |
613 | ||
614 | c2->Divide(2,2); | |
615 | c2->cd(1); | |
616 | TGraph* bzg = new TGraph(iA, phiA, bzPhiA); | |
617 | TGraph* bzcg = new TGraph(iA, phiA, bzcPhiA); | |
618 | bzcg->SetLineColor(2); | |
619 | bxzmg->Add(bzg); | |
620 | bxzmg->Add(bzcg); | |
621 | bxzmg->Draw("lA"); | |
622 | bxzmg->GetHistogram()->SetXTitle("#phi[rad]"); | |
623 | bxzmg->GetHistogram()->SetYTitle("B_{z} [T]"); | |
624 | bxzmg->Draw("lA"); | |
625 | ||
626 | c2->cd(2); | |
627 | TGraph* brg = new TGraph(iA, phiA, brPhiA); | |
628 | TGraph* brcg = new TGraph(iA, phiA, brcPhiA); | |
629 | brcg->SetLineColor(2); | |
630 | byzmg->Add(brcg); | |
631 | byzmg->Add(brg); | |
632 | ||
633 | byzmg->SetMaximum(0.03); | |
634 | byzmg->SetMinimum(-0.03); | |
635 | byzmg->Draw("lA"); | |
636 | byzmg->GetHistogram()->SetXTitle("#phi[rad]"); | |
637 | byzmg->GetHistogram()->SetYTitle("B_{r} [T]"); | |
638 | byzmg->Draw("lA"); | |
639 | ||
640 | c2->cd(3); | |
641 | TGraph* btg = new TGraph(iA, phiA, btPhiA); | |
642 | TGraph* btcg = new TGraph(iA, phiA, btcPhiA); | |
643 | btcg->SetLineColor(2); | |
644 | bzzmg->Add(btcg); | |
645 | bzzmg->Add(btg); | |
646 | bzzmg->Draw("lA"); | |
647 | bzzmg->SetMaximum(0.03); | |
648 | bzzmg->SetMinimum(-0.03); | |
649 | bzzmg->GetHistogram()->SetXTitle("#phi[rad]"); | |
650 | bzzmg->GetHistogram()->SetYTitle("B_{t} [T]"); | |
651 | bzzmg->Draw("lA"); | |
652 | ||
653 | ||
654 | c2->cd(4); | |
655 | TGraph* bg = new TGraph(iA, phiA, bbPhiA); | |
656 | TGraph* bcg = new TGraph(iA, phiA, bbcPhiA); | |
657 | bcg->SetLineColor(2); | |
658 | bbmg->Add(bg); | |
659 | bbmg->Add(bcg); | |
660 | bbmg->Draw("lA"); | |
661 | bbmg->GetHistogram()->SetXTitle("#phi[rad]"); | |
662 | bbmg->GetHistogram()->SetYTitle("|B| [T]"); | |
663 | bbmg->Draw("lA"); | |
664 | ||
665 | ||
666 | ||
667 | ||
668 | char pictFile[64]; | |
669 | sprintf(pictFile, "%s.gif", runName); | |
670 | c2->SaveAs(pictFile); | |
671 | ||
672 | // | |
673 | // Html generation | |
674 | // | |
675 | char htmlFile[64]; | |
676 | sprintf(htmlFile, "%s.html", runName); | |
677 | FILE* chtml = fopen(htmlFile, "w"); | |
678 | MakeHtmlHeaderPict(chtml); | |
679 | MakeHtmlPict(chtml, pictFile); | |
680 | MakeHtmlTableEntry(fHtmlMain, fileName, htmlFile, zmin, zmax, ifiles, 0., 0); | |
681 | ||
682 | // | |
683 | // | |
684 | printf("Bdl [Tm] %f %f \n", 2. * bdl, 2 * bdlc / 10.); | |
685 | } // files | |
686 | MakeHtmlTrailor(fHtmlMain); | |
687 | TFile* out = new TFile("fmap.root", "recreate"); | |
688 | fMap->Write(); | |
689 | out->Close(); | |
690 | } | |
691 | ||
692 | ||
693 | ||
694 | void AliFieldReader::MakeHtmlHeaderMain(FILE* file) | |
695 | { | |
14b22234 | 696 | // |
697 | // Write the header of the heml output | |
698 | // | |
7d5a9359 | 699 | fprintf(file,"<!DOCTYPE html PUBLIC \"-//W3C//DTD HTML 4.01 Transitional//EN\">\n"); |
700 | fprintf(file, "<html>\n"); | |
701 | fprintf(file, "<head>\n"); | |
702 | fprintf(file, "<meta http-equiv=\"content-type\"\n"); | |
703 | fprintf(file, "content=\"text/html; charset=ISO-8859-1\"\n"); | |
704 | fprintf(file, "<title>main.html</title>\n"); | |
705 | fprintf(file, "<\\head>\n"); | |
706 | fprintf(file, "<body>\n"); | |
707 | fprintf(file, "<table cellpadding=\"1\" cellspacing=\"1\" border=\"1\"\n"); | |
708 | fprintf(file, "style=\"text-align: left; width: 80;\">\n"); | |
709 | fprintf(file, "<tbody>\n"); | |
710 | fprintf(file, "<td style=\"vertical-align: top;\">File# <br></td>\n"); | |
711 | fprintf(file, "<td style=\"vertical-align: top;\">File Name <br></td>\n"); | |
712 | fprintf(file, "<td style=\"vertical-align: top;\">X-Position <br></td>\n"); | |
713 | fprintf(file, "<td style=\"vertical-align: top;\">Y-Position <br></td>\n"); | |
714 | fprintf(file, "<td style=\"vertical-align: top;\">Box Orientation <br></td>\n"); | |
715 | fprintf(file, "<td style=\"vertical-align: top;\">B.dl <br></td>\n"); | |
716 | fprintf(file, "<td style=\"vertical-align: top;\">Link to plots <br></td>\n"); | |
717 | ||
718 | } | |
719 | ||
720 | void AliFieldReader::MakeHtmlHeaderPict(FILE* file) | |
721 | { | |
14b22234 | 722 | // |
723 | // Write header for picture | |
724 | // | |
7d5a9359 | 725 | fprintf(file,"<!DOCTYPE html PUBLIC \"-//W3C//DTD HTML 4.01 Transitional//EN\">\n"); |
726 | fprintf(file, "<html>\n"); | |
727 | fprintf(file, "<head>\n"); | |
728 | fprintf(file, "<meta http-equiv=\"content-type\"\n"); | |
729 | fprintf(file, "content=\"text/html; charset=ISO-8859-1\"\n"); | |
730 | fprintf(file, "<title>main.html</title>\n"); | |
731 | fprintf(file, "</head>\n"); | |
732 | fprintf(file, "<body>\n"); | |
733 | } | |
734 | ||
735 | void AliFieldReader:: MakeHtmlPict(FILE* chtml, char* pictFile) | |
736 | { | |
14b22234 | 737 | // |
738 | // Write html for including picture | |
739 | // | |
7d5a9359 | 740 | fprintf(chtml, "<img src=\"./%s\" alt=\"%s\" style=\"width: 1196px; height: 772px;\">\n", |
741 | pictFile, pictFile); | |
742 | ||
743 | fprintf(chtml, "<br> <br> <br>\n"); | |
744 | fprintf(chtml, "<a href=\"./bmap.html\">Back to main page</a><br>\n"); | |
745 | ||
746 | fprintf(chtml, "</body>\n"); | |
747 | fprintf(chtml, "</header>\n"); | |
748 | fclose(chtml); | |
749 | } | |
750 | ||
751 | void AliFieldReader::MakeHtmlTableEntry(FILE* htmlmain, char* fileName, char* htmlFile, Float_t x, Float_t y, Int_t i, Float_t bdl, Int_t ifile) | |
752 | { | |
753 | fprintf(htmlmain, "<tr>\n"); | |
754 | // | |
755 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%5d</td>\n", ifile); | |
756 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%s</td>\n", fileName); | |
757 | ||
758 | // | |
759 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%13.2f</td>\n",x); | |
760 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%13.2f</td>\n",y); | |
761 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%3d </td>\n",i); | |
762 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">%13.3f</td>\n",bdl); | |
763 | // | |
764 | fprintf(htmlmain, "<td style=\"vertical-align: top;\">\n"); | |
765 | fprintf(htmlmain, "<span style=\"text-decoration: underline;\"><a href=\"./%s\">gif</a></span></td>\n", htmlFile); | |
766 | fprintf(htmlmain, "</tr>\n"); | |
767 | } | |
768 | ||
769 | ||
770 | void AliFieldReader::MakeHtmlTrailor(FILE* htmlmain) | |
771 | { | |
14b22234 | 772 | // |
773 | // Write the html trailor | |
774 | // | |
7d5a9359 | 775 | fprintf(htmlmain, "</tbody>\n"); |
776 | fprintf(htmlmain, "</table>\n"); | |
777 | fprintf(htmlmain, "</body>\n"); | |
778 | fprintf(htmlmain, "</header>\n"); | |
779 | fclose(htmlmain); | |
780 | } | |
781 | ||
782 | void AliFieldReader::ReadRegisterMap() | |
783 | { | |
14b22234 | 784 | // |
785 | // Read the register map | |
786 | // | |
7d5a9359 | 787 | FILE* regmap = fopen("register.map", "r"); |
788 | Int_t ireg; | |
789 | for (ireg = 0; ireg < 200; ireg++) { | |
790 | fRegMap[ireg][0] = -1; | |
791 | fRegMap[ireg][1] = -1; | |
792 | fRegMap[ireg][2] = -1; | |
793 | } | |
794 | for (Int_t iz = 0; iz < 2; iz++) { | |
795 | for (Int_t iy = 2; iy >= 0; iy--) { | |
796 | for (Int_t ix = 0; ix < 10; ix++) { | |
797 | fscanf(regmap, "%d\n", &ireg); | |
798 | printf("Address %5d %5d %5d %5d \n", iz, iy, ix, ireg); | |
799 | fRegMap[ireg][1] = iy; | |
800 | fRegMap[ireg][2] = iz; | |
801 | if (iz == 1) { | |
802 | fRegMap[ireg][0] = ix; | |
803 | } else { | |
804 | fRegMap[ireg][0] = 9 - ix; | |
805 | } | |
806 | } // ix | |
807 | } // iy | |
808 | } // iz | |
809 | fclose(regmap); | |
810 | printf("-> ReadRegisterMap()\n\n"); | |
811 | } | |
812 | ||
813 | ||
814 | void AliFieldReader::ReadRegisterMapSolenoid() | |
815 | { | |
14b22234 | 816 | // |
817 | // Read the register map | |
818 | // | |
7d5a9359 | 819 | FILE* regmap = fopen("register.map", "r"); |
820 | Int_t ireg; | |
821 | ||
822 | // Initialize | |
823 | for (ireg = 0; ireg < 200; ireg++) { | |
824 | fRegMap[ireg][0] = -1; | |
825 | fRegMap[ireg][1] = -1; | |
826 | fRegMap[ireg][2] = -1; | |
827 | } | |
828 | ||
829 | // Main arm | |
830 | for (Int_t ir = 0; ir < 33; ir++) { | |
831 | fscanf(regmap, "%d\n", &ireg); | |
832 | fRegMap[ireg][0] = ir; | |
833 | fRegMap[ireg][1] = 0; | |
834 | } | |
835 | // Opposite arm | |
836 | for (Int_t ir = 0; ir < 5; ir++) { | |
837 | fscanf(regmap, "%d\n", &ireg); | |
838 | fRegMap[ireg][0] = ir; | |
839 | fRegMap[ireg][1] = 1; | |
840 | } | |
841 | ||
842 | fclose(regmap); | |
843 | ||
844 | } | |
6772b8d8 | 845 | |
846 | ||
847 | AliFieldReader& AliFieldReader::operator=(const AliFieldReader& rhs) | |
848 | { | |
849 | // Assignment operator | |
850 | rhs.Copy(*this); | |
851 | return *this; | |
852 | } | |
853 | ||
854 | ||
855 | void AliFieldReader::Copy( TObject&) const | |
856 | { | |
857 | // | |
858 | // Copy | |
859 | // | |
860 | Fatal("Copy","Not implemented!\n"); | |
861 | } |