From c1e53270f7fdab009da61f6adae844db80489d95 Mon Sep 17 00:00:00 2001 From: sgorbuno Date: Mon, 17 Mar 2014 15:05:05 +0100 Subject: [PATCH] bug fix in inverse alignment calculation --- .../transform/AliHLTTPCFastTransform.cxx | 21 +++++++++++-------- 1 file changed, 12 insertions(+), 9 deletions(-) diff --git a/HLT/TPCLib/transform/AliHLTTPCFastTransform.cxx b/HLT/TPCLib/transform/AliHLTTPCFastTransform.cxx index 6eadb91584f..a39daa83012 100644 --- a/HLT/TPCLib/transform/AliHLTTPCFastTransform.cxx +++ b/HLT/TPCLib/transform/AliHLTTPCFastTransform.cxx @@ -149,18 +149,21 @@ Int_t AliHLTTPCFastTransform::Init( AliTPCTransform *transform, Long_t TimeStam bool AliHLTTPCFastTransform::CalcAdjugateRotation(const Float_t *mA, Float_t *mB, bool bCheck) { // check rotation matrix and adjugate for consistency + // + // ( for a rotation matrix inverse == transpose ) + // - mB[0]= mA[4]*mA[8]-mA[5]*mA[7]; - mB[1]= mA[5]*mA[6]-mA[3]*mA[8]; - mB[2]= mA[3]*mA[7]-mA[4]*mA[6]; + mB[0] = mA[0]; + mB[1] = mA[3]; + mB[2] = mA[6]; - mB[3]= mA[2]*mA[7]-mA[1]*mA[8]; - mB[4]= mA[0]*mA[8]-mA[2]*mA[6]; - mB[5]= mA[2]*mA[6]-mA[0]*mA[7]; + mB[3] = mA[1]; + mB[4] = mA[4]; + mB[5] = mA[7]; - mB[6]= mA[1]*mA[5]-mA[2]*mA[4]; - mB[7]= mA[2]*mA[3]-mA[0]*mA[5]; - mB[8]= mA[0]*mA[4]-mA[1]*mA[3]; + mB[6] = mA[2]; + mB[7] = mA[5]; + mB[8] = mA[8]; if (bCheck) { for (int r=0; r<3; r++) { -- 2.43.0