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