| 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 | } |
git://git.uio.no / u / mrichter / AliRoot.git / blame