[u/mrichter/AliRoot.git] / HLT / src / AliL3Transform.cxx

/* $Id$
// Author: Anders Vestbo <mailto:vestbo@fi.uib.no>, Uli Frankenfeld <mailto:franken@fi.uib.no>
// -- Copyright &copy AV
// changes done by Constantin Loizides <mailto:loizides@ikf.physik.uni-frankfurt.de>
*/

#include "AliL3Logging.h"
#include "AliL3Transform.h"
#include <math.h>

//_____________________________________________________________
// AliL3Transform
//
// Transformation class for ALICE TPC.

ClassImp(AliL3Transform);


AliL3Transform::AliL3Transform()
{
  //constructor
  
  LOG(AliL3Log::kWarning,"AliL3Transform::AliL3Transform()","Init")
    <<"You are using the _old_ version of parameter inits. V3.06 only! Use the constructor AliL3Transform(char*) instead."<<ENDLOG;
  Init();
}

AliL3Transform::AliL3Transform(const Char_t* pathname){
  //constructor
  Init(pathname);
}

AliL3Transform::~AliL3Transform(){
}

/// old init function used by Anders for AliRoot version 3.06
void AliL3Transform::Init()
{
  fVersion = 0;
  
  //sector:
  fNTimeBins = 446;
  fNRowLow = 64;
  fNRowUp = 112;
  fNSectorLow = 36;
  fNSectorUp = 36;
  fNSector = 72;
  fPadPitchWidthLow = 0.400000;
  fPadPitchWidthUp = 0.600000;
  fZWidth = 0.56599998474121093750;
  fZSigma = 0.22880849748219134199;

  //slices:
  fNSlice = 36;
  fNRow = 176;
  fPi = 3.141592653589793;

  fNRotShift=0.5;
  for(Int_t i=0;i<36;i++){
    fCos[i] = cos( (2*fPi/18) * (i+fNRotShift) );
    fSin[i] = sin( (2*fPi/18) * (i+fNRotShift) );
  }

  fX[0] = 84.570007324218750;
  fX[1] = 85.320007324218750;
  fX[2] = 86.070007324218750;
  fX[3] = 86.820007324218750;
  fX[4] = 87.570007324218750;
  fX[5] = 88.320007324218750;
  fX[6] = 89.070007324218750;
  fX[7] = 89.820007324218750;
  fX[8] = 90.570007324218750;
  fX[9] = 91.320007324218750;
  fX[10] = 92.070007324218750;
  fX[11] = 92.820007324218750;
  fX[12] = 93.570007324218750;
  fX[13] = 94.320007324218750;
  fX[14] = 95.070007324218750;
  fX[15] = 95.820007324218750;
  fX[16] = 96.570007324218750;
  fX[17] = 97.320007324218750;
  fX[18] = 98.070007324218750;
  fX[19] = 98.820007324218750;
  fX[20] = 99.570007324218750;
  fX[21] = 100.320007324218750;
  fX[22] = 101.070007324218750;
  fX[23] = 101.820007324218750;
  fX[24] = 102.570007324218750;
  fX[25] = 103.320007324218750;
  fX[26] = 104.070007324218750;
  fX[27] = 104.820007324218750;
  fX[28] = 105.570007324218750;
  fX[29] = 106.320007324218750;
  fX[30] = 107.070007324218750;
  fX[31] = 107.820007324218750;
  fX[32] = 108.570007324218750;
  fX[33] = 109.320007324218750;
  fX[34] = 110.070007324218750;
  fX[35] = 110.820007324218750;
  fX[36] = 111.570007324218750;
  fX[37] = 112.320007324218750;
  fX[38] = 113.070007324218750;
  fX[39] = 113.820007324218750;
  fX[40] = 114.570007324218750;
  fX[41] = 115.320007324218750;
  fX[42] = 116.070007324218750;
  fX[43] = 116.820007324218750;
  fX[44] = 117.570007324218750;
  fX[45] = 118.320007324218750;
  fX[46] = 119.070007324218750;
  fX[47] = 119.820007324218750;
  fX[48] = 120.570007324218750;
  fX[49] = 121.320007324218750;
  fX[50] = 122.070007324218750;
  fX[51] = 122.820007324218750;
  fX[52] = 123.570007324218750;
  fX[53] = 124.320007324218750;
  fX[54] = 125.070007324218750;
  fX[55] = 125.820007324218750;
  fX[56] = 126.570007324218750;
  fX[57] = 127.320007324218750;
  fX[58] = 128.070007324218750;
  fX[59] = 128.820007324218750;
  fX[60] = 129.570007324218750;
  fX[61] = 130.320007324218750;
  fX[62] = 131.070007324218750;
  fX[63] = 131.820007324218750;
  fX[64] = 135.054992675781250;
  fX[65] = 136.054992675781250;
  fX[66] = 137.054992675781250;
  fX[67] = 138.054992675781250;
  fX[68] = 139.054992675781250;
  fX[69] = 140.054992675781250;
  fX[70] = 141.054992675781250;
  fX[71] = 142.054992675781250;
  fX[72] = 143.054992675781250;
  fX[73] = 144.054992675781250;
  fX[74] = 145.054992675781250;
  fX[75] = 146.054992675781250;
  fX[76] = 147.054992675781250;
  fX[77] = 148.054992675781250;
  fX[78] = 149.054992675781250;
  fX[79] = 150.054992675781250;
  fX[80] = 151.054992675781250;
  fX[81] = 152.054992675781250;
  fX[82] = 153.054992675781250;
  fX[83] = 154.054992675781250;
  fX[84] = 155.054992675781250;
  fX[85] = 156.054992675781250;
  fX[86] = 157.054992675781250;
  fX[87] = 158.054992675781250;
  fX[88] = 159.054992675781250;
  fX[89] = 160.054992675781250;
  fX[90] = 161.054992675781250;
  fX[91] = 162.054992675781250;
  fX[92] = 163.054992675781250;
  fX[93] = 164.054992675781250;
  fX[94] = 165.054992675781250;
  fX[95] = 166.054992675781250;
  fX[96] = 167.054992675781250;
  fX[97] = 168.054992675781250;
  fX[98] = 169.054992675781250;
  fX[99] = 170.054992675781250;
  fX[100] = 171.054992675781250;
  fX[101] = 172.054992675781250;
  fX[102] = 173.054992675781250;
  fX[103] = 174.054992675781250;
  fX[104] = 175.054992675781250;
  fX[105] = 176.054992675781250;
  fX[106] = 177.054992675781250;
  fX[107] = 178.054992675781250;
  fX[108] = 179.054992675781250;
  fX[109] = 180.054992675781250;
  fX[110] = 181.054992675781250;
  fX[111] = 182.054992675781250;
  fX[112] = 183.054992675781250;
  fX[113] = 184.054992675781250;
  fX[114] = 185.054992675781250;
  fX[115] = 186.054992675781250;
  fX[116] = 187.054992675781250;
  fX[117] = 188.054992675781250;
  fX[118] = 189.054992675781250;
  fX[119] = 190.054992675781250;
  fX[120] = 191.054992675781250;
  fX[121] = 192.054992675781250;
  fX[122] = 193.054992675781250;
  fX[123] = 194.054992675781250;
  fX[124] = 195.054992675781250;
  fX[125] = 196.054992675781250;
  fX[126] = 197.054992675781250;
  fX[127] = 198.054992675781250;
  fX[128] = 199.054992675781250;
  fX[129] = 200.054992675781250;
  fX[130] = 201.054992675781250;
  fX[131] = 202.054992675781250;
  fX[132] = 203.054992675781250;
  fX[133] = 204.054992675781250;
  fX[134] = 205.054992675781250;
  fX[135] = 206.054992675781250;
  fX[136] = 207.054992675781250;
  fX[137] = 208.054992675781250;
  fX[138] = 209.054992675781250;
  fX[139] = 210.054992675781250;
  fX[140] = 211.054992675781250;
  fX[141] = 212.054992675781250;
  fX[142] = 213.054992675781250;
  fX[143] = 214.054992675781250;
  fX[144] = 215.054992675781250;
  fX[145] = 216.054992675781250;
  fX[146] = 217.054992675781250;
  fX[147] = 218.054992675781250;
  fX[148] = 219.054992675781250;
  fX[149] = 220.054992675781250;
  fX[150] = 221.054992675781250;
  fX[151] = 222.054992675781250;
  fX[152] = 223.054992675781250;
  fX[153] = 224.054992675781250;
  fX[154] = 225.054992675781250;
  fX[155] = 226.054992675781250;
  fX[156] = 227.054992675781250;
  fX[157] = 228.054992675781250;
  fX[158] = 229.054992675781250;
  fX[159] = 230.054992675781250;
  fX[160] = 231.054992675781250;
  fX[161] = 232.054992675781250;
  fX[162] = 233.054992675781250;
  fX[163] = 234.054992675781250;
  fX[164] = 235.054992675781250;
  fX[165] = 236.054992675781250;
  fX[166] = 237.054992675781250;
  fX[167] = 238.054992675781250;
  fX[168] = 239.054992675781250;
  fX[169] = 240.054992675781250;
  fX[170] = 241.054992675781250;
  fX[171] = 242.054992675781250;
  fX[172] = 243.054992675781250;
  fX[173] = 244.054992675781250;
  fX[174] = 245.054992675781250;
  fX[175] = 246.054992675781250;
  fNPads[0] = 67;
  fNPads[1] = 67;
  fNPads[2] = 67;
  fNPads[3] = 69;
  fNPads[4] = 69;
  fNPads[5] = 69;
  fNPads[6] = 71;
  fNPads[7] = 71;
  fNPads[8] = 71;
  fNPads[9] = 73;
  fNPads[10] = 73;
  fNPads[11] = 73;
  fNPads[12] = 75;
  fNPads[13] = 75;
  fNPads[14] = 75;
  fNPads[15] = 77;
  fNPads[16] = 77;
  fNPads[17] = 77;
  fNPads[18] = 79;
  fNPads[19] = 79;
  fNPads[20] = 79;
  fNPads[21] = 81;
  fNPads[22] = 81;
  fNPads[23] = 81;
  fNPads[24] = 83;
  fNPads[25] = 83;
  fNPads[26] = 83;
  fNPads[27] = 85;
  fNPads[28] = 85;
  fNPads[29] = 85;
  fNPads[30] = 87;
  fNPads[31] = 87;
  fNPads[32] = 87;
  fNPads[33] = 89;
  fNPads[34] = 89;
  fNPads[35] = 89;
  fNPads[36] = 91;
  fNPads[37] = 91;
  fNPads[38] = 91;
  fNPads[39] = 93;
  fNPads[40] = 93;
  fNPads[41] = 93;
  fNPads[42] = 95;
  fNPads[43] = 95;
  fNPads[44] = 95;
  fNPads[45] = 97;
  fNPads[46] = 97;
  fNPads[47] = 97;
  fNPads[48] = 99;
  fNPads[49] = 99;
  fNPads[50] = 99;
  fNPads[51] = 101;
  fNPads[52] = 101;
  fNPads[53] = 101;
  fNPads[54] = 103;
  fNPads[55] = 103;
  fNPads[56] = 103;
  fNPads[57] = 105;
  fNPads[58] = 105;
  fNPads[59] = 105;
  fNPads[60] = 107;
  fNPads[61] = 107;
  fNPads[62] = 107;
  fNPads[63] = 109;
  fNPads[64] = 73;
  fNPads[65] = 75;
  fNPads[66] = 75;
  fNPads[67] = 75;
  fNPads[68] = 77;
  fNPads[69] = 77;
  fNPads[70] = 77;
  fNPads[71] = 77;
  fNPads[72] = 79;
  fNPads[73] = 79;
  fNPads[74] = 79;
  fNPads[75] = 81;
  fNPads[76] = 81;
  fNPads[77] = 81;
  fNPads[78] = 83;
  fNPads[79] = 83;
  fNPads[80] = 83;
  fNPads[81] = 83;
  fNPads[82] = 85;
  fNPads[83] = 85;
  fNPads[84] = 85;
  fNPads[85] = 87;
  fNPads[86] = 87;
  fNPads[87] = 87;
  fNPads[88] = 87;
  fNPads[89] = 89;
  fNPads[90] = 89;
  fNPads[91] = 89;
  fNPads[92] = 91;
  fNPads[93] = 91;
  fNPads[94] = 91;
  fNPads[95] = 93;
  fNPads[96] = 93;
  fNPads[97] = 93;
  fNPads[98] = 93;
  fNPads[99] = 95;
  fNPads[100] = 95;
  fNPads[101] = 95;
  fNPads[102] = 97;
  fNPads[103] = 97;
  fNPads[104] = 97;
  fNPads[105] = 97;
  fNPads[106] = 99;
  fNPads[107] = 99;
  fNPads[108] = 99;
  fNPads[109] = 101;
  fNPads[110] = 101;
  fNPads[111] = 101;
  fNPads[112] = 103;
  fNPads[113] = 103;
  fNPads[114] = 103;
  fNPads[115] = 103;
  fNPads[116] = 105;
  fNPads[117] = 105;
  fNPads[118] = 105;
  fNPads[119] = 107;
  fNPads[120] = 107;
  fNPads[121] = 107;
  fNPads[122] = 107;
  fNPads[123] = 109;
  fNPads[124] = 109;
  fNPads[125] = 109;
  fNPads[126] = 111;
  fNPads[127] = 111;
  fNPads[128] = 111;
  fNPads[129] = 113;
  fNPads[130] = 113;
  fNPads[131] = 113;
  fNPads[132] = 113;
  fNPads[133] = 115;
  fNPads[134] = 115;
  fNPads[135] = 115;
  fNPads[136] = 117;
  fNPads[137] = 117;
  fNPads[138] = 117;
  fNPads[139] = 117;
  fNPads[140] = 119;
  fNPads[141] = 119;
  fNPads[142] = 119;
  fNPads[143] = 121;
  fNPads[144] = 121;
  fNPads[145] = 121;
  fNPads[146] = 123;
  fNPads[147] = 123;
  fNPads[148] = 123;
  fNPads[149] = 123;
  fNPads[150] = 125;
  fNPads[151] = 125;
  fNPads[152] = 125;
  fNPads[153] = 127;
  fNPads[154] = 127;
  fNPads[155] = 127;
  fNPads[156] = 127;
  fNPads[157] = 129;
  fNPads[158] = 129;
  fNPads[159] = 129;
  fNPads[160] = 131;
  fNPads[161] = 131;
  fNPads[162] = 131;
  fNPads[163] = 133;
  fNPads[164] = 133;
  fNPads[165] = 133;
  fNPads[166] = 133;
  fNPads[167] = 135;
  fNPads[168] = 135;
  fNPads[169] = 135;
  fNPads[170] = 137;
  fNPads[171] = 137;
  fNPads[172] = 137;
  fNPads[173] = 137;
  fNPads[174] = 139;
  fNPads[175] = 139;
}

void AliL3Transform::Init(const Char_t* path){
  Init(); //call the old init function to have some initial values

  Char_t *pathname=new Char_t[1024];
  strcpy(pathname,path);
  strcat(pathname,"/l3transform.config");

  FILE *fptr=fopen(pathname,"r");
  if(!fptr){
    LOG(AliL3Log::kWarning,"AliL3Transform::Init","File Open")
      <<"Pointer to Config File \""<<pathname<<"\" 0x0!"<<ENDLOG;
    return;
  }

  Char_t d1[250], d2[100], d3[100];
  Int_t dummy=0;
  Double_t ddummy=0.0;

  while(!feof(fptr)) {
    fscanf(fptr,"%s",d1);

    if(strcmp(d1,"fNTimeBins")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNTimeBins=(Int_t)dummy;}
    else if(strcmp(d1,"fNRowLow")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRowLow=(Int_t)dummy;}    
    else if(strcmp(d1,"fNRowUp")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRowUp=(Int_t)dummy;}
    else if(strcmp(d1,"fNSectorLow")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSectorLow=(Int_t)dummy;}
    else if(strcmp(d1,"fNSectorUp")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSectorUp=(Int_t)dummy;}
    else if(strcmp(d1,"fNSector")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSector=(Int_t)dummy;}
    else if(strcmp(d1,"fPadPitchWidthLow")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPadPitchWidthLow=(Double_t)ddummy;}
    else if(strcmp(d1,"fPadPitchWidthUp")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPadPitchWidthUp=(Double_t)ddummy;}
    else if(strcmp(d1,"fZWidth")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fZWidth=(Double_t)ddummy;}
    else if(strcmp(d1,"fZSigma")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fZSigma=(Double_t)ddummy;}
    else if(strcmp(d1,"fNSlice")==0){fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNSlice=(Int_t)dummy;}
    else if(strcmp(d1,"fNRow")==0){
      fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNRow=(Int_t)dummy;
      if(fNRow!=176){
	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;
      }
    }
    else if(strcmp(d1,"fNRotShift")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fNRotShift=(Double_t)ddummy;}
    else if(strcmp(d1,"fPi")==0){fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fPi=(Double_t)ddummy;}
    else if(strcmp(d1,"fX[0]")==0){
      fscanf(fptr,"%s %lf %s",d2,&ddummy,d3);fX[0]=(Double_t)ddummy;
      for(Int_t i=1;i<fNRow;i++){fscanf(fptr,"%s %s %lf %s",d1,d2,&ddummy,d3);fX[i]=(Double_t)dummy;}
    }
    else if(strcmp(d1,"fNPads[0]")==0){
      fscanf(fptr,"%s %d %s",d2,&dummy,d3);fNPads[0]=(Int_t)dummy;
      for(Int_t i=1;i<fNRow;i++){fscanf(fptr,"%s %s %d %s",d1,d2,&dummy,d3);fNPads[i]=(Int_t)dummy;}
    }
  }

  for(Int_t i=0;i<36;i++){
    fCos[i] = cos( (2*fPi/18) * (i+fNRotShift) );
    fSin[i] = sin( (2*fPi/18) * (i+fNRotShift) );
  }

  fclose(fptr);
  delete pathname;
  fVersion=1; //new version
}


Double_t AliL3Transform::GetEta(Float_t *xyz)
{
  Double_t r3 = sqrt(xyz[0]*xyz[0]+xyz[1]*xyz[1]+xyz[2]*xyz[2]);
  Double_t eta = 0.5 * log((r3+xyz[2])/(r3-xyz[2]));
  return eta;
}

Double_t AliL3Transform::GetEta(Int_t padrow,Int_t pad,Int_t time)
{
  Float_t xyz[3];
  Int_t sector,row;
  Slice2Sector(0,padrow,sector,row);
  Raw2Local(xyz,sector,row,pad,time);
  
  return GetEta(xyz);
}

Double_t AliL3Transform::GetPhi(Float_t *xyz)
{
  
  Double_t phi = atan2(xyz[1],xyz[0]);
  //if(phi<0) phi=phi+2*TMath::Pi();
  return phi;
}

Bool_t AliL3Transform::Slice2Sector(Int_t slice, Int_t slicerow, Int_t & sector, Int_t &row) const{
  if(slicerow<0&&slicerow>=fNRow) return kFALSE;
  if(slice<0||slice>=fNSlice) return kFALSE;

  if(slicerow<fNRowLow){
    sector = slice;
    row    = slicerow;
  }
  else {
    sector = slice+fNSlice;
    row    = slicerow-fNRowLow;
  }
  return kTRUE;
}

Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t  sector) const{
  if(sector<0||sector>=fNSector) return kFALSE;
  if(sector<fNSectorLow) slice = sector;
  else          slice = sector - fNSectorLow;
  return kTRUE;
}

Bool_t AliL3Transform::Sector2Slice(Int_t & slice, Int_t & slicerow,Int_t  sector, Int_t row) const{
  if(sector<0||sector>=fNSector||row<0) return kFALSE;
  if(sector<fNSectorLow){
    if(row>=fNRowLow) return kFALSE;
    slice = sector;
    slicerow = row;
  }
  else{
    if(row>=fNRowUp) return kFALSE;
    slice = sector - fNSectorLow;
    slicerow = row + fNRowLow;
  }
  return kTRUE;
}

Double_t AliL3Transform::Row2X(Int_t slicerow){
  if(slicerow<0||slicerow>=fNRow) return 0;
  return fX[slicerow];
}

void AliL3Transform::Local2Global(Float_t *xyz,Int_t slice)
{
  //Transformation to global coordinate system
  Float_t x0 = xyz[0];
  Float_t y0 = xyz[1];
  Float_t cs,sn;
  cs = fCos[slice];
  sn = fSin[slice];
  xyz[0]=x0*cs-y0*sn;
  xyz[1]=x0*sn+y0*cs;
  xyz[2]=xyz[2];//global z=local z
}

void AliL3Transform::Local2GlobalAngle(Float_t *angle,Int_t slice){
  angle[0] = fmod(angle[0]+(slice+0.5)*(2*fPi/18),2*fPi);
}

void AliL3Transform::Global2LocalAngle(Float_t *angle,Int_t slice){
  angle[0] = angle[0]-(slice+0.5)*(2*fPi/18);
  if(angle[0]<0) angle[0]+=2*fPi;
}

void AliL3Transform::Raw2Local(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time)
{
  //Transformation from rawdata to local coordinate system
  
  Int_t slice,slicerow;
  Sector2Slice(slice, slicerow, sector, row);  

  xyz[0]=Row2X(slicerow); 
  Int_t npads= fNPads[slicerow];
  if(sector<fNSectorLow)
    xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthLow;
  else
    xyz[1]=(pad-0.5*(npads-1))*fPadPitchWidthUp;
  xyz[2]=fZWidth*time-3.*fZSigma;
  Int_t sign=1;
  
  
  if(slice < 18)
    sign = 1;
  else
    sign = -1;
  
  xyz[2]=sign*(250.-xyz[2]);
  
  //Int_t nis=fNSectorLow;
  //Int_t nos=fNSectorUp;
  //if((sector<nis)/2 || ((sector-nis)<nos/2)) sign=1;
  
}

void AliL3Transform::Local2Global(Float_t *xyz,Int_t sector,Int_t row)
{
  //Transformation to global coordinate system
  Int_t slice,slicerow;
  Sector2Slice(slice, slicerow, sector, row);  
  Float_t r=Row2X(slicerow);
  Float_t cs = fCos[slice];
  Float_t sn = fSin[slice];

  xyz[0]=r*cs-xyz[1]*sn;
  xyz[1]=r*sn+xyz[1]*cs;
  xyz[2]=xyz[2];//global z=local z
}

Double_t AliL3Transform::GetMaxY(Int_t slicerow)
{
 
 if(slicerow < fNRowLow)
     return fPadPitchWidthLow*fNPads[slicerow]/2; 
 
 else
     return fPadPitchWidthUp*fNPads[slicerow]/2;
}

void AliL3Transform::Global2Local(Float_t *xyz,Int_t sector,Bool_t isSlice)
{
  
  Int_t slice;
  if(!isSlice)
    Sector2Slice(slice, sector);  
  else
    slice = sector;
  Float_t cs = fCos[slice];
  Float_t sn = fSin[slice];
  Float_t x1 = xyz[0]*cs + xyz[1]*sn;
  Float_t y1 = -xyz[0]*sn + xyz[1]*cs;
  xyz[0] = x1;
  xyz[1] = y1;
}

void AliL3Transform::Raw2Global(Float_t *xyz,Int_t sector,Int_t row,Float_t pad,Float_t time)
{
  //Transformation from raw to global coordinates
  
  Raw2Local(xyz,sector,row,pad,time);
  Local2Global(xyz,sector,row);
}

void AliL3Transform::Local2Raw(Float_t *xyz,Int_t sector,Int_t row)
{
  //Transformation from local coordinates to raw
  
  Int_t slice,slicerow;
  Sector2Slice(slice, slicerow, sector, row);  
   
  if(sector<fNSectorLow)
    xyz[1]=xyz[1]/fPadPitchWidthLow+0.5*(fNPads[slicerow]-1);
  else
    xyz[1]=xyz[1]/fPadPitchWidthUp+0.5*(fNPads[slicerow]-1);
  Int_t sign=-1;
  Int_t nis=fNSectorLow;
  Int_t nos=fNSectorUp;
 
  if ((sector<nis/2) || ((sector-nis)<nos/2)) sign=1; 
  xyz[2]=250-sign*xyz[2];
  xyz[2]=(xyz[2]+3.*fZSigma)/fZWidth;
}

void AliL3Transform::Global2Raw(Float_t *xyz,Int_t sector,Int_t row)
{
  //Transformation from global coordinates to raw. 

  Global2Local(xyz,sector);
  Local2Raw(xyz,sector,row);

}
Commit	Line	Data
a4e75a99	1	/* $Id$
b661165c	2	// Author: Anders Vestbo <mailto:vestbo@fi.uib.no>, Uli Frankenfeld <mailto:franken@fi.uib.no>
a4e75a99	3	// -- Copyright &copy AV
	4	// changes done by Constantin Loizides <mailto:loizides@ikf.physik.uni-frankfurt.de>
	5	*/
108615fc	6
	7	#include "AliL3Logging.h"
	8	#include "AliL3Transform.h"
108615fc	9	#include <math.h>
b661165c	10
b661165c	11	//_____________________________________________________________
108615fc	12	// AliL3Transform
	13	//
	14	// Transformation class for ALICE TPC.
	15
	16	ClassImp(AliL3Transform);
	17
	18
f59eed2d	19	AliL3Transform::AliL3Transform()
f59eed2d	20	{
108615fc	21	//constructor
f59eed2d	22
f59eed2d	23	LOG(AliL3Log::kWarning,"AliL3Transform::AliL3Transform()","Init")
89693986	24	<<"You are using the _old_ version of parameter inits. V3.06 only! Use the constructor AliL3Transform(char*) instead."<<ENDLOG;
108615fc	25	Init();
	26	}
	27
f59eed2d	28	AliL3Transform::AliL3Transform(const Char_t* pathname){
a4e75a99	29	//constructor
	30	Init(pathname);
	31	}
108615fc	32
	33	AliL3Transform::~AliL3Transform(){
	34	}
	35
a4e75a99	36	/// old init function used by Anders for AliRoot version 3.06
f59eed2d	37	void AliL3Transform::Init()
f59eed2d	38	{
89693986	39	fVersion = 0;
b25c64e5	40
a4982538	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;
a4e75a99	57
a4e75a99	58	fNRotShift=0.5;
40f2887b	59	for(Int_t i=0;i<36;i++){
a4e75a99	60	fCos[i] = cos( (2fPi/18) (i+fNRotShift) );
a4e75a99	61	fSin[i] = sin( (2fPi/18) (i+fNRotShift) );
40f2887b	62	}
a4982538	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	}
108615fc	417
f59eed2d	418	void AliL3Transform::Init(const Char_t* path){
a4e75a99	419	Init(); //call the old init function to have some initial values
a4e75a99	420
f59eed2d	421	Char_t *pathname=new Char_t[1024];
a4e75a99	422	strcpy(pathname,path);
	423	strcat(pathname,"/l3transform.config");
	424
	425	FILE *fptr=fopen(pathname,"r");
	426	if(!fptr){
b25c64e5	427	LOG(AliL3Log::kWarning,"AliL3Transform::Init","File Open")
b25c64e5	428	<<"Pointer to Config File \""<<pathname<<"\" 0x0!"<<ENDLOG;
a4e75a99	429	return;
	430	}
	431
f59eed2d	432	Char_t d1[250], d2[100], d3[100];
	433	Int_t dummy=0;
	434	Double_t ddummy=0.0;
a4e75a99	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( (2fPi/18) (i+fNRotShift) );
	470	fSin[i] = sin( (2fPi/18) (i+fNRotShift) );
	471	}
	472
	473	fclose(fptr);
	474	delete pathname;
89693986	475	fVersion=1; //new version
a4e75a99	476	}
a4e75a99	477
89693986	478
108615fc	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
3a1c3849	486	Double_t AliL3Transform::GetEta(Int_t padrow,Int_t pad,Int_t time)
bd7e82e0	487	{
bd7e82e0	488	Float_t xyz[3];
3a1c3849	489	Int_t sector,row;
	490	Slice2Sector(0,padrow,sector,row);
	491	Raw2Local(xyz,sector,row,pad,time);
bd7e82e0	492
	493	return GetEta(xyz);
	494	}
	495
108615fc	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
108615fc	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]=x0cs-y0sn;
	555	xyz[1]=x0sn+y0cs;
	556	xyz[2]=xyz[2];//global z=local z
	557	}
	558
	559	void AliL3Transform::Local2GlobalAngle(Float_t *angle,Int_t slice){
40f2887b	560	angle[0] = fmod(angle[0]+(slice+0.5)(2fPi/18),2*fPi);
108615fc	561	}
	562
	563	void AliL3Transform::Global2LocalAngle(Float_t *angle,Int_t slice){
40f2887b	564	angle[0] = angle[0]-(slice+0.5)(2fPi/18);
108615fc	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]=fZWidthtime-3.fZSigma;
1c225117	582	Int_t sign=1;
	583
	584
0e0b1e59	585	if(slice < 18)
	586	sign = 1;
	587	else
	588	sign = -1;
0e0b1e59	589
1c225117	590	xyz[2]=sign*(250.-xyz[2]);
0e0b1e59	591
	592	//Int_t nis=fNSectorLow;
	593	//Int_t nos=fNSectorUp;
	594	//if((sector<nis)/2 \|\| ((sector-nis)<nos/2)) sign=1;
	595
108615fc	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]=rcs-xyz[1]sn;
	608	xyz[1]=rsn+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
2342214c	622	void AliL3Transform::Global2Local(Float_t *xyz,Int_t sector,Bool_t isSlice)
108615fc	623	{
2342214c	624
108615fc	625	Int_t slice;
2342214c	626	if(!isSlice)
	627	Sector2Slice(slice, sector);
	628	else
	629	slice = sector;
108615fc	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	}