HLT/src/AliL3Transform.cxx

   1 /* $Id$
   2 // Author: Anders Vestbo <mailto:vestbo@fi.uib.no>, Uli Frankenfeld <mailto:franken@fi.uib.no>
   3 // -- Copyright &copy AV
   4 // changes done by Constantin Loizides <mailto:loizides@ikf.physik.uni-frankfurt.de>
   5 */
   6
   7 #include "AliL3Logging.h"
   8 #include "AliL3Transform.h"
   9 #include <math.h>
  10
  11 //_____________________________________________________________
  12 // AliL3Transform
  13 //
  14 // Transformation class for ALICE TPC.
  15
  16 ClassImp(AliL3Transform);
  17
  18
  19 AliL3Transform::AliL3Transform()
  20 {
  21   //constructor
  22
  23   LOG(AliL3Log::kWarning,"AliL3Transform::AliL3Transform()","Init")
  24     <<"You are using the _old_ version of parameter inits. V3.06 only! Use the constructor AliL3Transform(char*) instead."<<ENDLOG;
  25   Init();
  26 }
  27
  28 AliL3Transform::AliL3Transform(const Char_t* pathname){
  29   //constructor
  30   Init(pathname);
  31 }
  32
  33 AliL3Transform::~AliL3Transform(){
  34 }
  35
  36 /// old init function used by Anders for AliRoot version 3.06
  37 void AliL3Transform::Init()
  38 {
  39   fVersion = 0;
  40
  41   //sector:
  42   fNTimeBins = 446;
  43   fNRowLow = 64;
  44   fNRowUp = 112;
  45   fNSectorLow = 36;
  46   fNSectorUp = 36;
  47   fNSector = 72;
  48   fPadPitchWidthLow = 0.400000;
  49   fPadPitchWidthUp = 0.600000;
  50   fZWidth = 0.56599998474121093750;
  51   fZSigma = 0.22880849748219134199;
  52
  53   //slices:
  54   fNSlice = 36;
  55   fNRow = 176;
  56   fPi = 3.141592653589793;
  57
  58   fNRotShift=0.5;
  59   for(Int_t i=0;i<36;i++){
  60     fCos[i] = cos( (2*fPi/18) * (i+fNRotShift) );
  61     fSin[i] = sin( (2*fPi/18) * (i+fNRotShift) );
  62   }
  63
  64   fX[0] = 84.570007324218750;
  65   fX[1] = 85.320007324218750;
  66   fX[2] = 86.070007324218750;
  67   fX[3] = 86.820007324218750;
  68   fX[4] = 87.570007324218750;
  69   fX[5] = 88.320007324218750;
  70   fX[6] = 89.070007324218750;
  71   fX[7] = 89.820007324218750;
  72   fX[8] = 90.570007324218750;
  73   fX[9] = 91.320007324218750;
  74   fX[10] = 92.070007324218750;
  75   fX[11] = 92.820007324218750;
  76   fX[12] = 93.570007324218750;
  77   fX[13] = 94.320007324218750;
  78   fX[14] = 95.070007324218750;
  79   fX[15] = 95.820007324218750;
  80   fX[16] = 96.570007324218750;
  81   fX[17] = 97.320007324218750;
  82   fX[18] = 98.070007324218750;
  83   fX[19] = 98.820007324218750;
  84   fX[20] = 99.570007324218750;
  85   fX[21] = 100.320007324218750;
  86   fX[22] = 101.070007324218750;
  87   fX[23] = 101.820007324218750;
  88   fX[24] = 102.570007324218750;
  89   fX[25] = 103.320007324218750;
  90   fX[26] = 104.070007324218750;
  91   fX[27] = 104.820007324218750;
  92   fX[28] = 105.570007324218750;
  93   fX[29] = 106.320007324218750;
  94   fX[30] = 107.070007324218750;
  95   fX[31] = 107.820007324218750;
  96   fX[32] = 108.570007324218750;
  97   fX[33] = 109.320007324218750;
  98   fX[34] = 110.070007324218750;
  99   fX[35] = 110.820007324218750;
 100   fX[36] = 111.570007324218750;
 101   fX[37] = 112.320007324218750;
 102   fX[38] = 113.070007324218750;
 103   fX[39] = 113.820007324218750;
 104   fX[40] = 114.570007324218750;
 105   fX[41] = 115.320007324218750;
 106   fX[42] = 116.070007324218750;
 107   fX[43] = 116.820007324218750;
 108   fX[44] = 117.570007324218750;
 109   fX[45] = 118.320007324218750;
 110   fX[46] = 119.070007324218750;
 111   fX[47] = 119.820007324218750;
 112   fX[48] = 120.570007324218750;
 113   fX[49] = 121.320007324218750;
 114   fX[50] = 122.070007324218750;
 115   fX[51] = 122.820007324218750;
 116   fX[52] = 123.570007324218750;
 117   fX[53] = 124.320007324218750;
 118   fX[54] = 125.070007324218750;
 119   fX[55] = 125.820007324218750;
 120   fX[56] = 126.570007324218750;
 121   fX[57] = 127.320007324218750;
 122   fX[58] = 128.070007324218750;
 123   fX[59] = 128.820007324218750;
 124   fX[60] = 129.570007324218750;
 125   fX[61] = 130.320007324218750;
 126   fX[62] = 131.070007324218750;
 127   fX[63] = 131.820007324218750;
 128   fX[64] = 135.054992675781250;
 129   fX[65] = 136.054992675781250;
 130   fX[66] = 137.054992675781250;
 131   fX[67] = 138.054992675781250;
 132   fX[68] = 139.054992675781250;
 133   fX[69] = 140.054992675781250;
 134   fX[70] = 141.054992675781250;
 135   fX[71] = 142.054992675781250;
 136   fX[72] = 143.054992675781250;
 137   fX[73] = 144.054992675781250;
 138   fX[74] = 145.054992675781250;
 139   fX[75] = 146.054992675781250;
 140   fX[76] = 147.054992675781250;
 141   fX[77] = 148.054992675781250;
 142   fX[78] = 149.054992675781250;
 143   fX[79] = 150.054992675781250;
 144   fX[80] = 151.054992675781250;
 145   fX[81] = 152.054992675781250;
 146   fX[82] = 153.054992675781250;
 147   fX[83] = 154.054992675781250;
 148   fX[84] = 155.054992675781250;
 149   fX[85] = 156.054992675781250;
 150   fX[86] = 157.054992675781250;
 151   fX[87] = 158.054992675781250;
 152   fX[88] = 159.054992675781250;
 153   fX[89] = 160.054992675781250;
 154   fX[90] = 161.054992675781250;
 155   fX[91] = 162.054992675781250;
 156   fX[92] = 163.054992675781250;
 157   fX[93] = 164.054992675781250;
 158   fX[94] = 165.054992675781250;
 159   fX[95] = 166.054992675781250;
 160   fX[96] = 167.054992675781250;
 161   fX[97] = 168.054992675781250;
 162   fX[98] = 169.054992675781250;
 163   fX[99] = 170.054992675781250;
 164   fX[100] = 171.054992675781250;
 165   fX[101] = 172.054992675781250;
 166   fX[102] = 173.054992675781250;
 167   fX[103] = 174.054992675781250;
 168   fX[104] = 175.054992675781250;
 169   fX[105] = 176.054992675781250;
 170   fX[106] = 177.054992675781250;
 171   fX[107] = 178.054992675781250;
 172   fX[108] = 179.054992675781250;
 173   fX[109] = 180.054992675781250;
 174   fX[110] = 181.054992675781250;
 175   fX[111] = 182.054992675781250;
 176   fX[112] = 183.054992675781250;
 177   fX[113] = 184.054992675781250;
 178   fX[114] = 185.054992675781250;
 179   fX[115] = 186.054992675781250;
 180   fX[116] = 187.054992675781250;
 181   fX[117] = 188.054992675781250;
 182   fX[118] = 189.054992675781250;
 183   fX[119] = 190.054992675781250;
 184   fX[120] = 191.054992675781250;
 185   fX[121] = 192.054992675781250;
 186   fX[122] = 193.054992675781250;
 187   fX[123] = 194.054992675781250;
 188   fX[124] = 195.054992675781250;
 189   fX[125] = 196.054992675781250;
 190   fX[126] = 197.054992675781250;
 191   fX[127] = 198.054992675781250;
 192   fX[128] = 199.054992675781250;
 193   fX[129] = 200.054992675781250;
 194   fX[130] = 201.054992675781250;
 195   fX[131] = 202.054992675781250;
 196   fX[132] = 203.054992675781250;
 197   fX[133] = 204.054992675781250;
 198   fX[134] = 205.054992675781250;
 199   fX[135] = 206.054992675781250;
 200   fX[136] = 207.054992675781250;
 201   fX[137] = 208.054992675781250;
 202   fX[138] = 209.054992675781250;
 203   fX[139] = 210.054992675781250;
 204   fX[140] = 211.054992675781250;
 205   fX[141] = 212.054992675781250;
 206   fX[142] = 213.054992675781250;
 207   fX[143] = 214.054992675781250;
 208   fX[144] = 215.054992675781250;
 209   fX[145] = 216.054992675781250;
 210   fX[146] = 217.054992675781250;
 211   fX[147] = 218.054992675781250;
 212   fX[148] = 219.054992675781250;
 213   fX[149] = 220.054992675781250;
 214   fX[150] = 221.054992675781250;
 215   fX[151] = 222.054992675781250;
 216   fX[152] = 223.054992675781250;
 217   fX[153] = 224.054992675781250;
 218   fX[154] = 225.054992675781250;
 219   fX[155] = 226.054992675781250;
 220   fX[156] = 227.054992675781250;
 221   fX[157] = 228.054992675781250;
 222   fX[158] = 229.054992675781250;
 223   fX[159] = 230.054992675781250;
 224   fX[160] = 231.054992675781250;
 225   fX[161] = 232.054992675781250;
 226   fX[162] = 233.054992675781250;
 227   fX[163] = 234.054992675781250;
 228   fX[164] = 235.054992675781250;
 229   fX[165] = 236.054992675781250;
 230   fX[166] = 237.054992675781250;
 231   fX[167] = 238.054992675781250;
 232   fX[168] = 239.054992675781250;
 233   fX[169] = 240.054992675781250;
 234   fX[170] = 241.054992675781250;
 235   fX[171] = 242.054992675781250;
 236   fX[172] = 243.054992675781250;
 237   fX[173] = 244.054992675781250;
 238   fX[174] = 245.054992675781250;
 239   fX[175] = 246.054992675781250;
 240   fNPads[0] = 67;
 241   fNPads[1] = 67;
 242   fNPads[2] = 67;
 243   fNPads[3] = 69;
 244   fNPads[4] = 69;
 245   fNPads[5] = 69;
 246   fNPads[6] = 71;
 247   fNPads[7] = 71;
 248   fNPads[8] = 71;
 249   fNPads[9] = 73;
 250   fNPads[10] = 73;
 251   fNPads[11] = 73;
 252   fNPads[12] = 75;
 253   fNPads[13] = 75;
 254   fNPads[14] = 75;
 255   fNPads[15] = 77;
 256   fNPads[16] = 77;
 257   fNPads[17] = 77;
 258   fNPads[18] = 79;
 259   fNPads[19] = 79;
 260   fNPads[20] = 79;
 261   fNPads[21] = 81;
 262   fNPads[22] = 81;
 263   fNPads[23] = 81;
 264   fNPads[24] = 83;
 265   fNPads[25] = 83;
 266   fNPads[26] = 83;
 267   fNPads[27] = 85;
 268   fNPads[28] = 85;
 269   fNPads[29] = 85;
 270   fNPads[30] = 87;
 271   fNPads[31] = 87;
 272   fNPads[32] = 87;
 273   fNPads[33] = 89;
 274   fNPads[34] = 89;
 275   fNPads[35] = 89;
 276   fNPads[36] = 91;
 277   fNPads[37] = 91;
 278   fNPads[38] = 91;
 279   fNPads[39] = 93;
 280   fNPads[40] = 93;
 281   fNPads[41] = 93;
 282   fNPads[42] = 95;
 283   fNPads[43] = 95;
 284   fNPads[44] = 95;
 285   fNPads[45] = 97;
 286   fNPads[46] = 97;
 287   fNPads[47] = 97;
 288   fNPads[48] = 99;
 289   fNPads[49] = 99;
 290   fNPads[50] = 99;
 291   fNPads[51] = 101;
 292   fNPads[52] = 101;
 293   fNPads[53] = 101;
 294   fNPads[54] = 103;
 295   fNPads[55] = 103;
 296   fNPads[56] = 103;
 297   fNPads[57] = 105;
 298   fNPads[58] = 105;
 299   fNPads[59] = 105;
 300   fNPads[60] = 107;
 301   fNPads[61] = 107;
 302   fNPads[62] = 107;
 303   fNPads[63] = 109;
 304   fNPads[64] = 73;
 305   fNPads[65] = 75;
 306   fNPads[66] = 75;
 307   fNPads[67] = 75;
 308   fNPads[68] = 77;
 309   fNPads[69] = 77;
 310   fNPads[70] = 77;
 311   fNPads[71] = 77;
 312   fNPads[72] = 79;
 313   fNPads[73] = 79;
 314   fNPads[74] = 79;
 315   fNPads[75] = 81;
 316   fNPads[76] = 81;
 317   fNPads[77] = 81;
 318   fNPads[78] = 83;
 319   fNPads[79] = 83;
 320   fNPads[80] = 83;
 321   fNPads[81] = 83;
 322   fNPads[82] = 85;
 323   fNPads[83] = 85;
 324   fNPads[84] = 85;
 325   fNPads[85] = 87;
 326   fNPads[86] = 87;
 327   fNPads[87] = 87;
 328   fNPads[88] = 87;
 329   fNPads[89] = 89;
 330   fNPads[90] = 89;
 331   fNPads[91] = 89;
 332   fNPads[92] = 91;
 333   fNPads[93] = 91;
 334   fNPads[94] = 91;
 335   fNPads[95] = 93;
 336   fNPads[96] = 93;
 337   fNPads[97] = 93;
 338   fNPads[98] = 93;
 339   fNPads[99] = 95;
 340   fNPads[100] = 95;
 341   fNPads[101] = 95;
 342   fNPads[102] = 97;
 343   fNPads[103] = 97;
 344   fNPads[104] = 97;
 345   fNPads[105] = 97;
 346   fNPads[106] = 99;
 347   fNPads[107] = 99;
 348   fNPads[108] = 99;
 349   fNPads[109] = 101;
 350   fNPads[110] = 101;
 351   fNPads[111] = 101;
 352   fNPads[112] = 103;
 353   fNPads[113] = 103;
 354   fNPads[114] = 103;
 355   fNPads[115] = 103;
 356   fNPads[116] = 105;
 357   fNPads[117] = 105;
 358   fNPads[118] = 105;
 359   fNPads[119] = 107;
 360   fNPads[120] = 107;
 361   fNPads[121] = 107;
 362   fNPads[122] = 107;
 363   fNPads[123] = 109;
 364   fNPads[124] = 109;
 365   fNPads[125] = 109;
 366   fNPads[126] = 111;
 367   fNPads[127] = 111;
 368   fNPads[128] = 111;
 369   fNPads[129] = 113;
 370   fNPads[130] = 113;
 371   fNPads[131] = 113;
 372   fNPads[132] = 113;
 373   fNPads[133] = 115;
 374   fNPads[134] = 115;
 375   fNPads[135] = 115;
 376   fNPads[136] = 117;
 377   fNPads[137] = 117;
 378   fNPads[138] = 117;
 379   fNPads[139] = 117;
 380   fNPads[140] = 119;
 381   fNPads[141] = 119;
 382   fNPads[142] = 119;
 383   fNPads[143] = 121;
 384   fNPads[144] = 121;
 385   fNPads[145] = 121;
 386   fNPads[146] = 123;
 387   fNPads[147] = 123;
 388   fNPads[148] = 123;
 389   fNPads[149] = 123;
 390   fNPads[150] = 125;
 391   fNPads[151] = 125;
 392   fNPads[152] = 125;
 393   fNPads[153] = 127;
 394   fNPads[154] = 127;
 395   fNPads[155] = 127;
 396   fNPads[156] = 127;
 397   fNPads[157] = 129;
 398   fNPads[158] = 129;
 399   fNPads[159] = 129;
 400   fNPads[160] = 131;
 401   fNPads[161] = 131;
 402   fNPads[162] = 131;
 403   fNPads[163] = 133;
 404   fNPads[164] = 133;
 405   fNPads[165] = 133;
 406   fNPads[166] = 133;
 407   fNPads[167] = 135;
 408   fNPads[168] = 135;
 409   fNPads[169] = 135;
 410   fNPads[170] = 137;
 411   fNPads[171] = 137;
 412   fNPads[172] = 137;
 413   fNPads[173] = 137;
 414   fNPads[174] = 139;
 415   fNPads[175] = 139;
 416 }
 417
 418 void AliL3Transform::Init(const Char_t* path){
 419   Init(); //call the old init function to have some initial values
 420
 421   Char_t *pathname=new Char_t[1024];
 422   strcpy(pathname,path);
 423   strcat(pathname,"/l3transform.config");
 424
 425   FILE *fptr=fopen(pathname,"r");
 426   if(!fptr){
 427     LOG(AliL3Log::kWarning,"AliL3Transform::Init","File Open")
 428       <<"Pointer to Config File \""<<pathname<<"\" 0x0!"<<ENDLOG;
 429     return;
 430   }
 431
 432   Char_t d1[250], d2[100], d3[100];
 433   Int_t dummy=0;
 434   Double_t ddummy=0.0;
 435
 436   while(!feof(fptr)) {
 437     fscanf(fptr,"%s",d1);
 438
 439     if(strcmp(d1,"fNTimeBins")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNTimeBins=(Int_t)dummy;}
 440     else if(strcmp(d1,"fNRowLow")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRowLow=(Int_t)dummy;}
 441     else if(strcmp(d1,"fNRowUp")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRowUp=(Int_t)dummy;}
 442     else if(strcmp(d1,"fNSectorLow")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSectorLow=(Int_t)dummy;}
 443     else if(strcmp(d1,"fNSectorUp")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSectorUp=(Int_t)dummy;}
 444     else if(strcmp(d1,"fNSector")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSector=(Int_t)dummy;}
 445     else if(strcmp(d1,"fPadPitchWidthLow")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPadPitchWidthLow=(Double_t)ddummy;}
 446     else if(strcmp(d1,"fPadPitchWidthUp")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPadPitchWidthUp=(Double_t)ddummy;}
 447     else if(strcmp(d1,"fZWidth")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fZWidth=(Double_t)ddummy;}
 448     else if(strcmp(d1,"fZSigma")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fZSigma=(Double_t)ddummy;}
 449     else if(strcmp(d1,"fNSlice")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSlice=(Int_t)dummy;}
 450     else if(strcmp(d1,"fNRow")==0){
 451       fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRow=(Int_t)dummy;
 452       if(fNRow!=176){
 453         LOG(AliL3Log::kError,"AliL3Transform::Init","Overflow")<<"Number of PadRows should be 176! Check and fgrep the code for 176 to see the consequences of this major change!"<<ENDLOG;
 454       }
 455     }
 456     else if(strcmp(d1,"fNRotShift")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fNRotShift=(Double_t)ddummy;}
 457     else if(strcmp(d1,"fPi")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPi=(Double_t)ddummy;}
 458     else if(strcmp(d1,"fX[0]")==0){
 459       fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fX[0]=(Double_t)ddummy;
 460       for(Int_t i=1;i<fNRow;i++){fscanf(fptr,"%s %s %lf %s",d1,d2,&ddummy,d3);fX[i]=(Double_t)dummy;}
 461     }
 462     else if(strcmp(d1,"fNPads[0]")==0){
 463       fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNPads[0]=(Int_t)dummy;
 464       for(Int_t i=1;i<fNRow;i++){fscanf(fptr,"%s %s %d %s",d1,d2,&dummy,d3);fNPads[i]=(Int_t)dummy;}
 465     }
 466   }
 467
 468   for(Int_t i=0;i<36;i++){
 469     fCos[i] = cos( (2*fPi/18) * (i+fNRotShift) );
 470     fSin[i] = sin( (2*fPi/18) * (i+fNRotShift) );
 471   }
 472
 473   fclose(fptr);
 474   delete pathname;
 475   fVersion=1; //new version
 476 }
 477
 478
 479 Double_t AliL3Transform::GetEta(Float_t *xyz)
 480 {
 481   Double_t r3 = sqrt(xyz[0]*xyz[0]+xyz[1]*xyz[1]+xyz[2]*xyz[2]);
 482   Double_t eta = 0.5 * log((r3+xyz[2])/(r3-xyz[2]));
 483   return eta;
 484 }
 485
 486 Double_t AliL3Transform::GetEta(Int_t padrow,Int_t pad,Int_t time)
 487 {
 488   Float_t xyz[3];
 489   Int_t sector,row;
 490   Slice2Sector(0,padrow,sector,row);
 491   Raw2Local(xyz,sector,row,pad,time);
 492
 493   return GetEta(xyz);
 494 }
 495
 496 Double_t AliL3Transform::GetPhi(Float_t *xyz)
 497 {
 498
 499   Double_t phi = atan2(xyz[1],xyz[0]);
 500   //if(phi<0) phi=phi+2*TMath::Pi();
 501   return phi;
 502 }
 503
 504 Bool_t AliL3Transform::Slice2Sector(Int_t slice, Int_t slicerow, Int_t & sector, Int_t &row) const{
 505   if(slicerow<0&&slicerow>=fNRow) return kFALSE;
 506   if(slice<0||slice>=fNSlice) return kFALSE;
 507
 508   if(slicerow<fNRowLow){
 509     sector = slice;
 510     row    = slicerow;
 511   }
 512   else {
 513     sector = slice+fNSlice;
 514     row    = slicerow-fNRowLow;
 515   }
 516   return kTRUE;
 517 }
 518
 519 Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t  sector) const{
 520   if(sector<0||sector>=fNSector) return kFALSE;
 521   if(sector<fNSectorLow) slice = sector;
 522   else          slice = sector - fNSectorLow;
 523   return kTRUE;
 524 }
 525
 526 Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t & slicerow,Int_t  sector, Int_t row) const{
 527   if(sector<0||sector>=fNSector||row<0) return kFALSE;
 528   if(sector<fNSectorLow){
 529     if(row>=fNRowLow) return kFALSE;
 530     slice = sector;
 531     slicerow = row;
 532   }
 533   else{
 534     if(row>=fNRowUp) return kFALSE;
 535     slice = sector - fNSectorLow;
 536     slicerow = row + fNRowLow;
 537   }
 538   return kTRUE;
 539 }
 540
 541 Double_t AliL3Transform::Row2X(Int_t slicerow){
 542   if(slicerow<0||slicerow>=fNRow) return 0;
 543   return fX[slicerow];
 544 }
 545
 546 void AliL3Transform::Local2Global(Float_t *xyz,Int_t slice)
 547 {
 548   //Transformation to global coordinate system
 549   Float_t x0 = xyz[0];
 550   Float_t y0 = xyz[1];
 551   Float_t cs,sn;
 552   cs = fCos[slice];
 553   sn = fSin[slice];
 554   xyz[0]=x0*cs-y0*sn;
 555   xyz[1]=x0*sn+y0*cs;
 556   xyz[2]=xyz[2];//global z=local z
 557 }
 558
 559 void AliL3Transform::Local2GlobalAngle(Float_t *angle,Int_t slice){
 560   angle[0] = fmod(angle[0]+(slice+0.5)*(2*fPi/18),2*fPi);
 561 }
 562
 563 void AliL3Transform::Global2LocalAngle(Float_t *angle,Int_t slice){
 564   angle[0] = angle[0]-(slice+0.5)*(2*fPi/18);
 565   if(angle[0]<0) angle[0]+=2*fPi;
 566 }
 567
 568 void AliL3Transform::Raw2Local(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time)
 569 {
 570   //Transformation from rawdata to local coordinate system
 571
 572   Int_t slice,slicerow;
 573   Sector2Slice(slice, slicerow, sector, row);
 574
 575   xyz[0]=Row2X(slicerow);
 576   Int_t npads= fNPads[slicerow];
 577   if(sector<fNSectorLow)
 578     xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthLow;
 579   else
 580     xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthUp;
 581   xyz[2]=fZWidth*time-3.*fZSigma;
 582   Int_t sign=1;
 583
 584
 585   if(slice < 18)
 586     sign = 1;
 587   else
 588     sign = -1;
 589
 590   xyz[2]=sign*(250.-xyz[2]);
 591
 592   //Int_t nis=fNSectorLow;
 593   //Int_t nos=fNSectorUp;
 594   //if((sector<nis)/2 || ((sector-nis)<nos/2)) sign=1;
 595
 596 }
 597
 598 void AliL3Transform::Local2Global(Float_t *xyz,Int_t sector,Int_t row)
 599 {
 600   //Transformation to global coordinate system
 601   Int_t slice,slicerow;
 602   Sector2Slice(slice, slicerow, sector, row);
 603   Float_t r=Row2X(slicerow);
 604   Float_t cs = fCos[slice];
 605   Float_t sn = fSin[slice];
 606
 607   xyz[0]=r*cs-xyz[1]*sn;
 608   xyz[1]=r*sn+xyz[1]*cs;
 609   xyz[2]=xyz[2];//global z=local z
 610 }
 611
 612 Double_t AliL3Transform::GetMaxY(Int_t slicerow)
 613 {
 614
 615  if(slicerow < fNRowLow)
 616      return fPadPitchWidthLow*fNPads[slicerow]/2;
 617
 618  else
 619      return fPadPitchWidthUp*fNPads[slicerow]/2;
 620 }
 621
 622 void AliL3Transform::Global2Local(Float_t *xyz,Int_t sector,Bool_t isSlice)
 623 {
 624
 625   Int_t slice;
 626   if(!isSlice)
 627     Sector2Slice(slice, sector);
 628   else
 629     slice = sector;
 630   Float_t cs = fCos[slice];
 631   Float_t sn = fSin[slice];
 632   Float_t x1 = xyz[0]*cs + xyz[1]*sn;
 633   Float_t y1 = -xyz[0]*sn + xyz[1]*cs;
 634   xyz[0] = x1;
 635   xyz[1] = y1;
 636 }
 637
 638 void AliL3Transform::Raw2Global(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time)
 639 {
 640   //Transformation from raw to global coordinates
 641
 642   Raw2Local(xyz,sector,row,pad,time);
 643   Local2Global(xyz,sector,row);
 644 }
 645
 646 void AliL3Transform::Local2Raw(Float_t *xyz,Int_t sector,Int_t row)
 647 {
 648   //Transformation from local coordinates to raw
 649
 650   Int_t slice,slicerow;
 651   Sector2Slice(slice, slicerow, sector, row);
 652
 653   if(sector<fNSectorLow)
 654     xyz[1]=xyz[1]/fPadPitchWidthLow+0.5*(fNPads[slicerow]-1);
 655   else
 656     xyz[1]=xyz[1]/fPadPitchWidthUp+0.5*(fNPads[slicerow]-1);
 657   Int_t sign=-1;
 658   Int_t nis=fNSectorLow;
 659   Int_t nos=fNSectorUp;
 660
 661   if ((sector<nis/2) || ((sector-nis)<nos/2)) sign=1;
 662   xyz[2]=250-sign*xyz[2];
 663   xyz[2]=(xyz[2]+3.*fZSigma)/fZWidth;
 664 }
 665
 666 void AliL3Transform::Global2Raw(Float_t *xyz,Int_t sector,Int_t row)
 667 {
 668   //Transformation from global coordinates to raw.
 669
 670   Global2Local(xyz,sector);
 671   Local2Raw(xyz,sector,row);
 672
 673 }