1 #include "AliTPCTransform.h"
3 #include "AliTPCCalROC.h"
4 #include "AliTPCcalibDB.h"
5 #include "AliTPCParam.h"
9 cdb=AliCDBManager::Instance()
10 cdb->SetDefaultStorage("local:///u/mmager/mycalib1")
11 c=AliTPCcalibDB::Instance()
13 Double_t x[]={1.0,2.0,3.0}
16 trafo.Transform(x,i,0,1)
19 AliTPCTransform::AliTPCTransform() {
21 for (Int_t i=0;i<18;++i) {
22 Double_t alpha=TMath::DegToRad()*(10.+20.*(i%18));
23 fSins[i]=TMath::Sin(alpha);
24 fCoss[i]=TMath::Cos(alpha);
28 AliTPCTransform::~AliTPCTransform() {
31 void AliTPCTransform::Transform(Double_t *x,Int_t *i,UInt_t time,
32 Int_t coordinateType) {
37 // output: x[0] - x (all in the rotated global coordinate frame)
40 Int_t row=TMath::Nint(x[1]);
41 Int_t pad=TMath::Nint(x[0]);
43 AliTPCcalibDB* const calib=AliTPCcalibDB::Instance();
46 //ugly: calib->SetRun(time);
49 //TODO: x[2]-=calib->GetPadTime0()->GetCalROC(sector)->GetValue(row,pad);
53 // TODO: use a map or parametrisation!
56 Pad2RotatedGlobal(pad,row,x);
59 //TODO: calib->GetParameters()->GetClusterMatrix(sector)->LocalToMaster(x,xx);
61 RotatedGlobal2Global(sector,x);
64 calib->GetExB()->Correct(x,xx);
66 Global2RotatedGlobal(sector,xx);
68 x[0]=xx[0];x[1]=xx[1];x[2]=xx[2];
71 inline void AliTPCTransform::Pad2RotatedGlobal(Int_t pad,Int_t row,Double_t *x)
74 AliTPCROC::Instance()->GetPositionLocal(0,row,pad,tmp);
79 //TODO rotation in the right direction?
80 inline void AliTPCTransform::RotatedGlobal2Global(Int_t sector,Double_t *x)
83 GetCosAndSin(sector,cos,sin);
85 x[0]= cos*tmp+sin*x[1];
86 x[1]=-sin*tmp+cos*x[1];
89 inline void AliTPCTransform::Global2RotatedGlobal(Int_t sector,Double_t *x)
92 GetCosAndSin(sector,cos,sin);
94 x[0]= cos*tmp-sin*x[1];
95 x[1]= sin*tmp+cos*x[1];
98 inline void AliTPCTransform::GetCosAndSin(Int_t sector,Double_t &cos,
99 Double_t &sin) const {
100 cos=fCoss[sector%18];
101 sin=fSins[sector%18];
104 ClassImp(AliTPCTransform)