TPC/AliTPCMonitorMappingHandler.cxx

   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 /*
  17   $Log$
  18   Revision 1.2  2007/10/12 13:36:27  cvetan
  19   Coding convention fixes from Stefan
  20
  21   Revision 1.1  2007/09/17 10:23:31  cvetan
  22   New TPC monitoring package from Stefan Kniege. The monitoring package can be started by running TPCMonitor.C macro located in macros folder.
  23
  24 */
  25
  26 ////////////////////////////////////////////////////////////////////////
  27 ////
  28 //// AliTPCMonitorMappingHandler class
  29 ////
  30 //// Class for handling mapping information TPC
  31 ////
  32 //// The mapping information for the TPC front end electornics (pads, front end cards)
  33 //// are handled by this class.
  34 //// The information from the design mapping and from the hardware address can be
  35 //// cross checked in the TPCMonitor.C.
  36 //// Malfunctioning front end cards can be identified by looking at single channels
  37 //// displayed with  the TPCMonitor.
  38 ////
  39 ////
  40 //// Authors: Roland Bramm,
  41 ////          Stefan Kniege, IKF, Frankfurt
  42 ////
  43 /////////////////////////////////////////////////////////////////////////
  44
  45
  46 #include <cstdlib>
  47 #include "AliTPCMonitorMappingHandler.h"
  48 #include "TH1.h"
  49 #include "TLegend.h"
  50 #include "AliLog.h"
  51 #include <Riostream.h>
  52
  53 ClassImp(AliTPCMonitorMappingHandler)
  54
  55 //_____________________________________________________________________________________________
  56 AliTPCMonitorMappingHandler::AliTPCMonitorMappingHandler(Char_t* name, Char_t* title):
  57   TNamed(name,title),
  58   fnumofChannels(0),
  59   fmaxHWAdress(0),
  60   fsizeofArray(0),
  61   fmapping(new Short_t*[24000]),
  62   fmappingChannelinRow(new Int_t*[160]),
  63   fu2ftestmapping(new Short_t*[7000]),
  64   fMapHwFECglobal(new Int_t*[24000]),
  65   fecGainMap(new Float_t*[7000])
  66 {
  67   // Constructor : Initialize mapping arrays
  68
  69   for(Int_t in = 0; in<160; in++)
  70     {
  71       Int_t* hold = new Int_t[150];
  72       for(Int_t jn = 0; jn<150;jn++)  hold[jn]=0;
  73       fmappingChannelinRow[in]= hold;
  74     }
  75
  76   for(Int_t i = 0; i<7000; i++)
  77     {
  78       Short_t* hold = new Short_t[8];
  79       for(Int_t j = 0; j<8;j++)  hold[j]=0;
  80       fu2ftestmapping[i]= hold;
  81     }
  82
  83   for(Int_t i = 0; i<24000; i++)
  84     {
  85       Int_t* hold = new Int_t[2];
  86       for(Int_t j = 0; j<2;j++)  hold[j]=0;
  87       fMapHwFECglobal[i]= hold;
  88     }
  89
  90
  91   for(Int_t i = 0; i<7000; i++)
  92     {
  93       Float_t* hold = new Float_t[128];
  94       for(Int_t j = 0; j<128;j++)  hold[j]=0;
  95       fecGainMap[i]= hold;
  96     }
  97 }
  98
  99 //____________________________________________________________________________
 100 AliTPCMonitorMappingHandler::AliTPCMonitorMappingHandler(const  AliTPCMonitorMappingHandler &maphand) :
 101   TNamed(maphand.GetName(),maphand.GetTitle()),
 102   fnumofChannels(maphand.fnumofChannels),
 103   fmaxHWAdress(maphand.fmaxHWAdress),
 104   fsizeofArray(maphand.fsizeofArray),
 105   fmapping(new Short_t*[24000]),
 106   fmappingChannelinRow(new Int_t*[160]),
 107   fu2ftestmapping(new Short_t*[7000]),
 108   fMapHwFECglobal(new Int_t*[24000]),
 109   fecGainMap(new Float_t*[7000])
 110 {
 111   // copy constructor
 112
 113
 114   for(Int_t in = 0; in<160; in++)
 115     {
 116       Int_t* hold = new Int_t[150];
 117       for(Int_t jn = 0; jn<150;jn++)  hold[jn]=maphand.fmappingChannelinRow[in][jn];
 118       fmappingChannelinRow[in]= hold;
 119     }
 120
 121   for(Int_t i = 0; i<7000; i++)
 122     {
 123       Short_t* hold = new Short_t[8];
 124       for(Int_t j = 0; j<8;j++)  hold[j]=maphand.fu2ftestmapping[i][j];
 125       fu2ftestmapping[i]= hold;
 126     }
 127
 128   for(Int_t i = 0; i<24000; i++)
 129     {
 130       Int_t* hold = new Int_t[2];
 131       for(Int_t j = 0; j<2;j++)  hold[j]=maphand.fMapHwFECglobal[i][j];
 132       fMapHwFECglobal[i]= hold;
 133     }
 134
 135   for(Int_t i = 0; i<7000; i++)
 136     {
 137       Float_t* hold = new Float_t[128];
 138       for(Int_t j = 0; j<128;j++)  hold[j]=maphand.fecGainMap[i][j];
 139       fecGainMap[i]= hold;
 140     }
 141
 142 }
 143
 144
 145 //____________________________________________________________________________
 146 AliTPCMonitorMappingHandler &AliTPCMonitorMappingHandler:: operator= (const AliTPCMonitorMappingHandler& maphand)
 147 {
 148   // assignment operator
 149   if(this!=&maphand)
 150     {
 151       fnumofChannels=maphand.fnumofChannels;
 152       fmaxHWAdress=maphand.fmaxHWAdress;
 153       fsizeofArray=maphand.fsizeofArray;
 154
 155       fmapping = new Short_t*[24000]; // empty
 156
 157
 158       fmappingChannelinRow = new Int_t*[160];
 159       for(Int_t in = 0; in<160; in++)
 160         {
 161           Int_t* hold = new Int_t[150];
 162           for(Int_t jn = 0; jn<150;jn++)  hold[jn]=maphand.fmappingChannelinRow[in][jn];
 163           fmappingChannelinRow[in]= hold;
 164         }
 165
 166       fu2ftestmapping      = new Short_t*[7000];
 167       for(Int_t i = 0; i<7000; i++)
 168         {
 169           Short_t* hold = new Short_t[8];
 170           for(Int_t j = 0; j<8;j++)  hold[j]=maphand.fu2ftestmapping[i][j];
 171           fu2ftestmapping[i]= hold;
 172         }
 173
 174       fMapHwFECglobal = new Int_t*[24000];
 175       for(Int_t i = 0; i<24000; i++)
 176         {
 177           Int_t* hold = new Int_t[2];
 178           for(Int_t j = 0; j<2;j++)  hold[j]=maphand.fMapHwFECglobal[i][j];
 179           fMapHwFECglobal[i]= hold;
 180         }
 181
 182       fecGainMap = new Float_t*[7000];
 183       for(Int_t i = 0; i<7000; i++)
 184         {
 185           Float_t* hold = new Float_t[128];
 186           for(Int_t j = 0; j<128;j++)  hold[j]=maphand.fecGainMap[i][j];
 187           fecGainMap[i]= hold;
 188         }
 189
 190
 191     }
 192   return *this;
 193 }
 194
 195
 196 //_____________________________________________________________________________________________
 197 AliTPCMonitorMappingHandler::~AliTPCMonitorMappingHandler()
 198 {
 199   // Destructor
 200   Short_t* temp;
 201   for(Int_t i = 0; i < 24000 ; i++) {
 202     temp = (Short_t *) fmapping[i];
 203     if(temp[1] != -1)
 204       delete temp;
 205   }
 206   temp = (Short_t *) fmapping;
 207   delete temp;
 208
 209
 210   Int_t* temp1;
 211   for(Int_t i = 0; i < 24000 ; i++) {
 212     temp1 = (Int_t *) fMapHwFECglobal[i];
 213     delete temp1;
 214   }
 215   temp1 = (Int_t *) fMapHwFECglobal;
 216   delete temp1;
 217
 218   for(Int_t i = 0; i < 7000 ; i++) {
 219     temp = (Short_t *) fu2ftestmapping[i];
 220     delete temp;
 221   }
 222
 223   temp = (Short_t *) fu2ftestmapping;
 224   delete temp;
 225
 226   for(Int_t i = 0; i < 150 ; i++) {
 227     temp1= (Int_t *) fmappingChannelinRow[i];
 228     delete temp1;
 229   }
 230   temp1 = (Int_t *) fmappingChannelinRow;
 231   delete temp1;
 232
 233   Float_t* temp2;
 234   for(Int_t i = 0; i<7000; i++)
 235     {
 236       temp2= (Float_t *) fecGainMap[i];
 237       delete temp2;
 238     }
 239   temp2 = (Float_t *) fecGainMap;
 240   delete temp2;
 241
 242
 243 }
 244
 245 //_____________________________________________________________________________________________
 246 void AliTPCMonitorMappingHandler::ReadMapping(char* mapfile)
 247 {
 248   // Read global Mapping file
 249   // Format of data in mapping file:
 250   // column 0:  hadrware address
 251   // column 1:  readout index in IROC/OROC
 252   // column 2:  global pad row number (0-158)
 253   // column 3   pad number in row
 254   // column 4:  connector
 255   // column 5:  pin
 256   // column 6:  fec number in IROC/OROC
 257   // column 7:  fec channel
 258   // column 8:  fec connector
 259   // column 9:  altro channel
 260   // column 10: altro chip
 261
 262   // Mapping data for a given hardware address are placed at the
 263   // index corresponding to the value of the hardware address in
 264   // the fmapping array..
 265   // The mapping information for the hardware address 'hwaddr'
 266   // can hence be found in fmapping[hwaddr]
 267
 268
 269   Short_t*  mappingRow;
 270   char readcarry[255];
 271   Int_t version = -1;
 272   Int_t actPos  = 0;
 273   Int_t oldPos  = 0;
 274
 275   ifstream *in = new ifstream();
 276   in->open(mapfile);
 277   *in >> readcarry;
 278   version = atoi(readcarry);
 279
 280   *in >> readcarry;
 281   fnumofChannels = atoi(readcarry);
 282   *in >> readcarry;
 283   fmaxHWAdress = atoi(readcarry);
 284   fsizeofArray = fmaxHWAdress;
 285
 286   Short_t *fmappingEmptyRow = new Short_t[11];
 287   for(Int_t i = 0; i < 11 ; i++) {
 288     fmappingEmptyRow[i] = 0;
 289   }
 290   fmappingEmptyRow[1] = -1;
 291   for(Int_t i = 0; i < fnumofChannels ; i++) {
 292     mappingRow = new Short_t[11];
 293     for(Int_t j = 0 ; j < 11 ; j++) {
 294       *in >> readcarry;
 295       mappingRow[j] = atoi(readcarry);
 296     }
 297     actPos = mappingRow[0];
 298     fmapping[actPos] = mappingRow;
 299     if( (actPos - oldPos) > 1) {
 300       for(Int_t j = (oldPos+1); j < actPos; j++) {
 301         fmapping[j] = fmappingEmptyRow;
 302       }
 303     }
 304     oldPos = actPos;
 305     fmapping[actPos][0] = i;
 306   }
 307
 308   in->close();
 309   delete in;
 310 }
 311
 312 //_____________________________________________________________________________________________
 313 Int_t  AliTPCMonitorMappingHandler::ReadFECMapping(char* u2ftestfile)
 314 {
 315   // Read in Mapping of global FEC numbers to branches, rcu patches etc.
 316   // Format in Mapping file
 317   // column 0: global card number (set to 0 if card number does not exist)
 318   // column 1: side of the TPC
 319   // column 2: sector on the side   (0-17)
 320   // column 3: rcu in the sector    (0-5)
 321   // column 4: fec number in rcu    (0-24)
 322   // column 5: fec number in branch (0-12)
 323   // column 6: branch number        (0,1)
 324   // column 7: local hardware address of first channel (not used)
 325
 326   // Order of data is kept in fu2ftestmapping
 327
 328   ifstream datin(u2ftestfile);
 329
 330   Int_t carry     = 0;
 331   Int_t ncards    = 0;
 332
 333   for(Int_t ind = 0; ind<7000; ind++)
 334     {
 335       for(Int_t entr = 0; entr<8; entr++)
 336         {
 337           datin >> carry ;
 338           fu2ftestmapping[ind][entr] = carry  ;
 339
 340           if(entr==0 && carry!=0) ncards++;
 341         }
 342     }
 343   return ncards ;
 344 }
 345
 346
 347
 348 //_____________________________________________________________________________________________
 349 void AliTPCMonitorMappingHandler::ReadfecHwMap(Int_t sector)
 350 {
 351   // Create mapping of global FEC numbers to hardware addresses for a given sector
 352
 353   Int_t  fside           =0;
 354   Int_t  fsector        = 0;
 355   Int_t  fec            = 0;
 356   Int_t  branch         = 0;
 357   Int_t  rcupatch      = 0;
 358   Int_t  altrchann      = 0;
 359   Int_t  altrchip       = 0;
 360   Int_t  nextHwAddress  = 0;
 361   Int_t  nfecs          = 0;
 362
 363   if(sector/18==0) fside =65;
 364   else             fside =67;
 365
 366   if(sector>18)    fsector= sector-18;
 367   else             fsector= sector   ;
 368
 369   for(Int_t ind = 0; ind<7000; ind++)
 370     {
 371       if((Int_t)U2fGetSide(ind)==fside && U2fGetSector(ind)==fsector)
 372         {
 373           nfecs++;
 374           fec            = U2fGetFECinBranch(ind);
 375           branch         = U2fGetBranch(ind);
 376           rcupatch      = U2fGetRCU(ind);
 377
 378           for(Int_t ch = 0; ch<128; ch++)
 379             {
 380               altrchann      = ch%16;
 381               altrchip       = ch/16;
 382
 383               nextHwAddress  = (   ((branch&1)<<11) + (fec<<7) + (altrchip<<4) + (altrchann)  + ((rcupatch-1)<<12) );
 384
 385               fMapHwFECglobal[nextHwAddress][0] = ind;
 386               fMapHwFECglobal[nextHwAddress][1] = ch ;
 387             }
 388         }
 389     }
 390 }
 391 //_____________________________________________________________________________________________
 392 void AliTPCMonitorMappingHandler::ReadfecGainMap(char* fecgainmap)
 393 {
 394   // Read global gain calibration pap
 395   // Format in file :
 396   // colummn 0      : FEC number
 397   // colums  1-128  : gain calibration factors
 398   ifstream datin(fecgainmap);
 399
 400   Int_t   fecnr  = 0;
 401   Float_t val    = 0.0 ;
 402
 403   while(!datin.eof())
 404     {
 405       datin >> fecnr ;
 406       for(Int_t in = 0; in<128; in++)
 407         {
 408           datin >> val ;
 409           fecGainMap[fecnr][in] = val;
 410         }
 411     }
 412 }
 413
 414 //_____________________________________________________________________________________________
 415 void  AliTPCMonitorMappingHandler::ReadRowMappingGlob(char* fpathtoMappingRowfile)
 416 {
 417   // Read mapping of hardware addresses in rows
 418   // Format of file:
 419   // column 0:        global row number (0-158)
 420   // column 1:        number of pads in this row (npads)
 421   // column 2-npads:  hardware addresses for these pads
 422
 423   char readcarry[256];
 424   char readcarry2[256];
 425   ifstream *in = new ifstream();
 426   in->open(fpathtoMappingRowfile);
 427
 428   for(Int_t i = 0; i < 159 ; i++) {
 429     *in >> readcarry;   // row number
 430     *in >> readcarry2;  // numof pads
 431     fmappingChannelinRow[i][0] = atoi(readcarry2);
 432     fmappingChannelinRow[i][1] = atoi(readcarry);
 433
 434     for(Int_t j = 2 ; j < fmappingChannelinRow[i][0]+2 ; j++) {
 435       *in >> readcarry;
 436       fmappingChannelinRow[i][j] = atoi(readcarry);
 437     }
 438   }
 439   in->close();
 440 }
 441
 442
 443
 444 //_____________________________________________________________________________________________
 445 Int_t AliTPCMonitorMappingHandler::GetNumOfChannels() const
 446 {
 447   // Return number of channels
 448   return fnumofChannels;
 449 }
 450
 451 //_____________________________________________________________________________________________
 452 Int_t AliTPCMonitorMappingHandler::GetSizeofArray() const
 453 {
 454   // Return sise of global mapping fmapping array.
 455   // Value orresponds to max value of hardware addresses
 456   return fsizeofArray;
 457 }
 458
 459
 460 //_____________________________________________________________________________________________
 461 Short_t* AliTPCMonitorMappingHandler::GetLine(Int_t hwaddr)const
 462 {
 463   // Return pointer to mapping array for the hardware address hwaddr
 464   Short_t* retval;
 465   if(hwaddr <= fsizeofArray)
 466     retval = fmapping[hwaddr];
 467   else
 468     retval = 0;
 469   return retval;
 470 }
 471
 472 //_____________________________________________________________________________________________
 473 Int_t AliTPCMonitorMappingHandler::GetIndex(Int_t hwaddr) const
 474 {
 475   // Return readout index for the hardware address hwaddr
 476   Int_t retval;
 477   if(hwaddr <= fsizeofArray)
 478     retval = fmapping[hwaddr][1];
 479   else
 480     retval = 0;
 481   return retval;
 482 }
 483
 484 //_____________________________________________________________________________________________
 485 Int_t AliTPCMonitorMappingHandler::GetPadRow(Int_t hwaddr) const
 486 {
 487   // Return global pad row  for the hardware address hwaddr
 488   Int_t retval;
 489   if(hwaddr <= fsizeofArray)
 490     retval = fmapping[hwaddr][2];
 491   else
 492     retval = 0;
 493   return retval;
 494 }
 495
 496 //_____________________________________________________________________________________________
 497 Int_t AliTPCMonitorMappingHandler::GetPad(Int_t hwaddr) const
 498 {
 499   // Return pad number in row for the hardware address hwaddr
 500   Int_t retval;
 501   if(hwaddr < fsizeofArray)
 502     retval = fmapping[hwaddr][3];
 503   else
 504     retval = 0;
 505   return retval;
 506 }
 507
 508 //_____________________________________________________________________________________________
 509 Int_t AliTPCMonitorMappingHandler::GetConnector(Int_t hwaddr) const
 510 {
 511   // Return connector for the hardware address hwaddr
 512   Int_t retval;
 513   if(hwaddr <= fsizeofArray)
 514     retval = fmapping[hwaddr][4];
 515   else
 516     retval = 0;
 517   return retval;
 518 }
 519
 520 //_____________________________________________________________________________________________
 521 Int_t AliTPCMonitorMappingHandler::GetPin(Int_t hwaddr) const
 522 {
 523   // Return pin for the hardware address hwaddr
 524   Int_t retval;
 525   if(hwaddr <= fsizeofArray)
 526     retval = fmapping[hwaddr][5];
 527   else
 528     retval = 0;
 529   return retval;
 530 }
 531
 532 //_____________________________________________________________________________________________
 533 Int_t AliTPCMonitorMappingHandler::GetFEC(Int_t hwaddr) const
 534 {
 535   // Return fec number in IROC/OROC  for the hardware address hwaddr
 536   Int_t retval;
 537   if(hwaddr <= fsizeofArray)
 538     retval = fmapping[hwaddr][6];
 539   else
 540     retval = 0;
 541   return retval;
 542 }
 543
 544 //_____________________________________________________________________________________________
 545 Int_t AliTPCMonitorMappingHandler::GetFECchannel(Int_t hwaddr) const
 546 {
 547   // Return FEC channel for the hardware address hwaddr
 548   Int_t retval;
 549   if(hwaddr < fsizeofArray)
 550     retval = fmapping[hwaddr][7];
 551   else
 552     retval = 0;
 553   return retval;
 554 }
 555
 556 //_____________________________________________________________________________________________
 557 Int_t AliTPCMonitorMappingHandler::GetFECconnector(Int_t hwaddr) const
 558 {
 559   // Return FEC connector for the hardware address hwaddr
 560   Int_t retval;
 561   if(hwaddr <= fsizeofArray)
 562     retval = fmapping[hwaddr][8];
 563   else
 564     retval = 0;
 565   return retval;
 566 }
 567
 568 //_____________________________________________________________________________________________
 569 Int_t AliTPCMonitorMappingHandler::GetAltroChannel(Int_t hwaddr)  const
 570 {
 571   // Return Altro channel for the hardware address hwaddr
 572   Int_t retval;
 573   if(hwaddr <= fsizeofArray)
 574     retval = fmapping[hwaddr][9];
 575   else
 576     retval = 0;
 577   return retval;
 578 }
 579
 580 //_____________________________________________________________________________________________
 581 Int_t AliTPCMonitorMappingHandler::GetAltro(Int_t hwaddr)  const
 582 {
 583   // Return Altro chip number in FEC for the hardware address hwaddr
 584   Int_t retval;
 585   if(hwaddr <= fsizeofArray)
 586     retval = fmapping[hwaddr][10];
 587   else
 588     retval = 0;
 589   return retval;
 590 }
 591
 592
 593 //_____________________________________________________________________________________________
 594 Int_t AliTPCMonitorMappingHandler::GetNumofPads(Int_t row)
 595 {
 596   // Return number of pads in row
 597   if(row<159)
 598     return fmappingChannelinRow[row][0];
 599   else  {
 600     AliError("Row number to high");
 601     return 0;
 602   }
 603 }
 604
 605 //_____________________________________________________________________________________________
 606 Int_t AliTPCMonitorMappingHandler::GetPadAddInRow(Int_t row,Int_t pad )
 607 {
 608   // Return hardware address for given pad in row
 609   if(row<159)
 610     return fmappingChannelinRow[row][pad+2];
 611   else {
 612     AliError("Row number to high");
 613     return 0;
 614   }
 615 }
 616
 617
 618
 619 //_____________________________________________________________________________________________
 620 Int_t AliTPCMonitorMappingHandler::U2fGetFECnr(Int_t index) const
 621 {
 622   // Return FEC number for index  (FEC number should be equal to index)
 623   return  fu2ftestmapping[index][0];
 624 }
 625
 626 //_____________________________________________________________________________________________
 627 Int_t AliTPCMonitorMappingHandler::U2fGetSide(Int_t fecnr) const
 628 {
 629   // Return side on which FEC is installed
 630   return  fu2ftestmapping[fecnr][1];
 631 }
 632
 633 //_____________________________________________________________________________________________
 634 Int_t AliTPCMonitorMappingHandler::U2fGetSector(Int_t fecnr) const
 635 {
 636   // Return sector in which FEC is installed
 637   return  fu2ftestmapping[fecnr][2];
 638 }
 639
 640 //_____________________________________________________________________________________________
 641 Int_t AliTPCMonitorMappingHandler::U2fGetRCU(Int_t fecnr) const
 642 {
 643   // Rerurn rcu patch in which FEC is installed
 644   return  fu2ftestmapping[fecnr][3];
 645 }
 646
 647 //_____________________________________________________________________________________________
 648 Int_t AliTPCMonitorMappingHandler::U2fGetFECinRCU(Int_t fecnr) const
 649 {
 650   // Return index of FEC in RCU (0-25)
 651   return  fu2ftestmapping[fecnr][4];
 652 }
 653
 654 //_____________________________________________________________________________________________
 655 Int_t AliTPCMonitorMappingHandler::U2fGetFECinBranch(Int_t fecnr) const
 656 {
 657   // Return index of FEC in branch (0-12)
 658   return  fu2ftestmapping[fecnr][5];
 659 }
 660
 661 //_____________________________________________________________________________________________
 662 Int_t AliTPCMonitorMappingHandler::U2fGetBranch(Int_t fecnr) const
 663 {
 664   // Return branch in which FEC is installed (0,1)
 665   return  fu2ftestmapping[fecnr][6];
 666
 667 }