- Double_t rotmatr[9];
- AliITSgeomTGeo::GetRotation(layer,ladder,detector,rotmatr);
- Double_t sAlpha=TMath::ATan2(rotmatr[1],rotmatr[0])+TMath::Pi();
- sAlpha+=TMath::Pi()/2.;
- if(layer==1) sAlpha+=TMath::Pi();
-
-
- // get the radius of this detector. Procedure taken from the
- // AliITStrackerV2 constructor
- Float_t x=0,y=0,z=0;
- Double_t xyz[3];
- AliITSgeomTGeo::GetTranslation(layer,ladder,detector,xyz);
- x=xyz[0];
- y=xyz[1];
- z=xyz[2];
-
- Double_t fi=TMath::ATan2(rotmatr[1],rotmatr[0])+TMath::Pi();
- fi+=TMath::Pi()/2;
- if (layer==1) fi+=TMath::Pi();
- Double_t cp=TMath::Cos(fi), sp=TMath::Sin(fi);
- Double_t sX=x*cp+y*sp;