1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
23 #include <TMultiGraph.h>
25 #include <TObjString.h>
29 #include "AliFieldReader.h"
30 #include "AliMagFMaps.h"
32 ClassImp(AliFieldReader)
36 //_______________________________________________________________________
37 AliFieldReader::AliFieldReader():
47 fCatalogueName("goodfiles.list")
54 AliFieldReader::AliFieldReader(const AliFieldReader& reader):
73 //_______________________________________________________________________
74 AliFieldReader::~AliFieldReader()
81 //_______________________________________________________________________
82 void AliFieldReader::Init()
85 // Initialize the reader
88 fField = new AliMagFMaps("Maps","Maps", 2, 1., 10., 2);
90 fCatalogue = fopen(fCatalogueName, "r");
93 fHtmlMain = fopen("bmap.html", "w");
94 MakeHtmlHeaderMain(fHtmlMain);
99 void AliFieldReader::ReadMap()
102 // Read the measured dipole field map
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];
112 Float_t xpos, ypos, zpos;
115 Float_t h[3], x[3], b[3];
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);
137 while ((fgets(fLine, 255, fCatalogue)) != NULL && ifiles <= 2000) {
139 if (strncmp(fLine,"#HEADER", 7) == 0) {
140 iret = sscanf(&fLine[7],"%f", &fZStart);
145 // if (ifiles != 87) continue;
147 iret = sscanf(fLine, "%s", fileName);
148 printf("Reading File %s\n", fileName);
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");
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) {
172 if (strncmp(sLine," Current", 8) == 0) {
174 iret = sscanf(&sLine[9],"%f", ¤t);
175 printf("Current %f \n", current);
177 fPolarity = current / 6000.;
179 if (strncmp(sLine," z", 2) == 0) {
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);
189 if (strncmp(sLine," X position", 11) == 0)
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);
202 printf("This file is for x = %13.3f y = %13.3f Box Position %1d\n", xpos, ypos, iboxpos);
207 if (strncmp(sLine," R Z-pos",8) == 0)
209 iret = sscanf(&sLine[8],"%e", &zpos);
216 if (strncmp(sLine,"C",1) == 0) {
219 iret = sscanf(&sLine[2],"%d %f %f %f %f %d", &ireg, &h[0], &h[1], &h[2], &temp, &ires);
221 // fix for address 115
223 if (ireg == 115) ic115++;
224 if (ic115 == 2 && ireg == 115) ireg = 119;
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];
237 dx = -0.36 + jx * fDd;
243 dx = 0.36 + jy * fDd;
244 dy = -(-0.36 + jx * fDd);
249 dx = 0.36 - jx * fDd;
255 dx = -0.36 - jy * fDd;
256 dy = -(0.36 - jx * fDd);
263 if (iboxpos == 1 || iboxpos == 3) {
271 Float_t dz = (jz == 0)? fDz/2. : - fDz/2.;
275 Float_t xc = xpos + dx;
276 Float_t yc = ypos + dy;
281 x[2] = - (-zpos * 100. + fZStart + dz);
290 xnt[ 2] = -x[2] / 100.;
291 xnt[ 3] = Float_t (jx);
292 xnt[ 4] = Float_t (jy);
293 xnt[ 5] = Float_t (jz);
300 xnt[12] = Float_t (ifiles);
301 xnt[13] = Float_t (ireg);
302 xnt[14] = Float_t(iboxpos);
304 xnt[16] = Float_t(ires);
312 if (jy != -1 && jz == 1 && jy == 0 && izz == 40){
319 bxycA[jx] = b[0] / 10.;
320 byycA[jx] = b[1] / 10.;
321 bzycA[jx] = b[2] / 10.;
325 if (jy != -1 && jz == 1 && jy == 0 && izz == 1){
332 bxxcA[jx] = b[0] / 10.;
333 byxcA[jx] = b[1] / 10.;
334 bzxcA[jx] = b[2] / 10.;
341 // printf("Field %f %f %f %f %f %f \n", x[0], x[1], x[2], b[0], b[1] , b[2]);
344 bdlc += stepsz * b[0];
350 bxzcA[iz] = b[0] / 10.;
351 byzcA[iz] = b[1] / 10.;
352 bzzcA[iz] = b[2] / 10.;
358 gStyle->SetOptStat(0);
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);
364 TGraph* bxg = new TGraph(iz, zA, bxzA);
365 TGraph* bxcg = new TGraph(iz, zA, bxzcA);
366 bxcg->SetLineColor(2);
370 bxzmg->GetHistogram()->SetXTitle("z[m]");
371 bxzmg->GetHistogram()->SetYTitle("B_{x} [T]");
375 TGraph* byg = new TGraph(iz, zA, byzA);
376 TGraph* bycg = new TGraph(iz, zA, byzcA);
377 bycg->SetLineColor(2);
381 byzmg->GetHistogram()->SetXTitle("z[m]");
382 byzmg->GetHistogram()->SetYTitle("B_{y} [T]");
386 TGraph* bzg = new TGraph(iz, zA, bzzA);
387 TGraph* bzcg = new TGraph(iz, zA, bzzcA);
388 bzcg->SetLineColor(2);
392 bzzmg->GetHistogram()->SetXTitle("z[m]");
393 bzzmg->GetHistogram()->SetYTitle("B_{z} [T]");
398 sprintf(pictFile, "%s.gif", runName);
399 c2->SaveAs(pictFile);
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);
412 printf("Bdl [Tm] %f %f \n", 2. * bdl, 2 * bdlc / 10.);
414 MakeHtmlTrailor(fHtmlMain);
415 TFile* out = new TFile("fmap.root", "recreate");
420 void AliFieldReader::ReadMapSolenoid(){
422 // Read map for solenoid measurement
424 Float_t phiA[450], bzPhiA[200], brPhiA[200], btPhiA[200], bbPhiA[200];
425 Float_t bzcPhiA[200], brcPhiA[200], btcPhiA[200], bbcPhiA[200];
429 Float_t zpos, phipos, skewing, temp;
432 Float_t h[3], x[3], b[3];
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);
446 while ((fgets(fLine, 255, fCatalogue)) != NULL && ifiles <= 2000) {
448 if (strncmp(fLine,"#HEADER", 7) == 0) {
449 iret = sscanf(&fLine[7],"%f", &fZStart);
455 iret = sscanf(fLine, "%s", fileName);
456 printf("Reading File %s\n", fileName);
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());
465 for (Int_t i = 1; i < n-1; i++)
467 sprintf(runName, "%s.%s",
468 runName, (((TObjString*) tokens->At(i))->GetString()).Data());
471 FILE* file = fopen(fileName, "r");
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|");
485 while ((fgets(sLine, 255, file)) != NULL) {
486 if (strncmp(sLine," z", 2) == 0) {
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));
492 if (strncmp(sLine," R\tPOSITION NUMBER", 18) == 0)
495 // Current z-position
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);
507 if (strncmp(sLine," SKEWING ON Z:", 14) == 0)
511 iret = sscanf(&sLine[14],"%f", &skewing);
512 printf("SKEWING ON Z: %f\n", skewing);
516 if (strncmp(sLine,"Phi", 3) == 0)
518 Float_t phiStart, phiStop;
520 iret = sscanf(&sLine[3],"%f %f", &phiStart, &phiStop);
522 printf("phiStart phiStop %f %f\n", phiStart, phiStop);
527 if (strncmp(sLine," R\tPhi-Angle",12) == 0)
529 iret = sscanf(&sLine[12],"%e", &phipos);
535 if (strncmp(sLine,"C",1) == 0) {
538 iret = sscanf(&sLine[2],"%d %f %f %f %f %d", &ireg, &h[0], &h[1], &h[2], &temp, &ical);
541 Int_t ir = fRegMap[ireg][0];
542 Int_t ia = fRegMap[ireg][1];
547 rpos = 0.2295 + ir * 0.16;
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;
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;
562 Float_t phi = phipos;
565 if (phi > 360.) phi -= 360.;
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.;
573 x[2] = -400. + zpos * 100.;
576 Float_t phi0 = TMath::Pi() - phi;
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);
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));
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);
596 bbPhiA[iA] = TMath::Sqrt(h[0] * h[0] + h[1] * h[1] + h[2] * h[2]);
600 bbcPhiA[iA] = TMath::Sqrt(brc * brc + bzc * bzc + btc * btc);
607 gStyle->SetOptStat(0);
609 sprintf(title, "Z = %13.2f m , Phi = %13.2f m", zpos, phipos);
612 TCanvas* c2 = new TCanvas("c2", title, 1200, 800);
616 TGraph* bzg = new TGraph(iA, phiA, bzPhiA);
617 TGraph* bzcg = new TGraph(iA, phiA, bzcPhiA);
618 bzcg->SetLineColor(2);
622 bxzmg->GetHistogram()->SetXTitle("#phi[rad]");
623 bxzmg->GetHistogram()->SetYTitle("B_{z} [T]");
627 TGraph* brg = new TGraph(iA, phiA, brPhiA);
628 TGraph* brcg = new TGraph(iA, phiA, brcPhiA);
629 brcg->SetLineColor(2);
633 byzmg->SetMaximum(0.03);
634 byzmg->SetMinimum(-0.03);
636 byzmg->GetHistogram()->SetXTitle("#phi[rad]");
637 byzmg->GetHistogram()->SetYTitle("B_{r} [T]");
641 TGraph* btg = new TGraph(iA, phiA, btPhiA);
642 TGraph* btcg = new TGraph(iA, phiA, btcPhiA);
643 btcg->SetLineColor(2);
647 bzzmg->SetMaximum(0.03);
648 bzzmg->SetMinimum(-0.03);
649 bzzmg->GetHistogram()->SetXTitle("#phi[rad]");
650 bzzmg->GetHistogram()->SetYTitle("B_{t} [T]");
655 TGraph* bg = new TGraph(iA, phiA, bbPhiA);
656 TGraph* bcg = new TGraph(iA, phiA, bbcPhiA);
657 bcg->SetLineColor(2);
661 bbmg->GetHistogram()->SetXTitle("#phi[rad]");
662 bbmg->GetHistogram()->SetYTitle("|B| [T]");
669 sprintf(pictFile, "%s.gif", runName);
670 c2->SaveAs(pictFile);
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);
684 printf("Bdl [Tm] %f %f \n", 2. * bdl, 2 * bdlc / 10.);
686 MakeHtmlTrailor(fHtmlMain);
687 TFile* out = new TFile("fmap.root", "recreate");
694 void AliFieldReader::MakeHtmlHeaderMain(FILE* file)
697 // Write the header of the heml output
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");
720 void AliFieldReader::MakeHtmlHeaderPict(FILE* file)
723 // Write header for picture
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");
735 void AliFieldReader:: MakeHtmlPict(FILE* chtml, char* pictFile)
738 // Write html for including picture
740 fprintf(chtml, "<img src=\"./%s\" alt=\"%s\" style=\"width: 1196px; height: 772px;\">\n",
743 fprintf(chtml, "<br> <br> <br>\n");
744 fprintf(chtml, "<a href=\"./bmap.html\">Back to main page</a><br>\n");
746 fprintf(chtml, "</body>\n");
747 fprintf(chtml, "</header>\n");
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)
753 fprintf(htmlmain, "<tr>\n");
755 fprintf(htmlmain, "<td style=\"vertical-align: top;\">%5d</td>\n", ifile);
756 fprintf(htmlmain, "<td style=\"vertical-align: top;\">%s</td>\n", fileName);
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);
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");
770 void AliFieldReader::MakeHtmlTrailor(FILE* htmlmain)
773 // Write the html trailor
775 fprintf(htmlmain, "</tbody>\n");
776 fprintf(htmlmain, "</table>\n");
777 fprintf(htmlmain, "</body>\n");
778 fprintf(htmlmain, "</header>\n");
782 void AliFieldReader::ReadRegisterMap()
785 // Read the register map
787 FILE* regmap = fopen("register.map", "r");
789 for (ireg = 0; ireg < 200; ireg++) {
790 fRegMap[ireg][0] = -1;
791 fRegMap[ireg][1] = -1;
792 fRegMap[ireg][2] = -1;
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;
802 fRegMap[ireg][0] = ix;
804 fRegMap[ireg][0] = 9 - ix;
810 printf("-> ReadRegisterMap()\n\n");
814 void AliFieldReader::ReadRegisterMapSolenoid()
817 // Read the register map
819 FILE* regmap = fopen("register.map", "r");
823 for (ireg = 0; ireg < 200; ireg++) {
824 fRegMap[ireg][0] = -1;
825 fRegMap[ireg][1] = -1;
826 fRegMap[ireg][2] = -1;
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;
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;
847 AliFieldReader& AliFieldReader::operator=(const AliFieldReader& rhs)
849 // Assignment operator
855 void AliFieldReader::Copy( TObject&) const
860 Fatal("Copy","Not implemented!\n");