75 fA1(0), fA2(0), fA3(0), fAsOK(
kFALSE),
89 fA1(t.fA1), fA2(t.fA2), fA3(t.fA3), fAsOK(t.fAsOK),
90 fUseTrans (t.fUseTrans),
91 fEditTrans(t.fEditTrans),
103 fA1(0), fA2(0), fA3(0), fAsOK(
kFALSE),
106 fEditRotation(
kTRUE),
117 fA1(0), fA2(0), fA3(0), fAsOK(
kFALSE),
120 fEditRotation(
kTRUE),
165 memcpy(
fM, t.
fM,
sizeof(
fM));
166 if (copyAngles && t.
fAsOK) {
179 for(
Int_t i=0; i<16; ++i)
fM[i] = arr[i];
188 for(
Int_t i=0; i<16; ++i)
fM[i] = arr[i];
204 t(i,j) = -s; t(j,i) = s;
220 static const float kFromToEpsilon = 0.000001f;
226 f = (e < 0.0f) ? -e : e;
228 if (f > 1.0f - kFromToEpsilon)
234 x.
fX = (from.
fX > 0.0f) ? from.
fX : -from.
fX;
235 x.
fY = (from.
fY > 0.0f) ? from.
fY : -from.
fY;
236 x.
fZ = (from.
fZ > 0.0f) ? from.
fZ : -from.
fZ;
241 x.
fX = 1.0f; x.
fY = x.
fZ = 0.0f;
243 x.
fZ = 1.0f; x.
fX = x.
fY = 0.0f;
249 x.
fY = 1.0f; x.
fX = x.
fZ = 0.0f;
251 x.
fZ = 1.0f; x.
fX = x.
fY = 0.0f;
258 c1 = 2.0f / u.
Mag2();
259 c2 = 2.0f / v.
Mag2();
260 c3 = c1 * c2 * u.
Dot(v);
262 for (
int i = 0; i < 3; i++) {
263 for (
int j = 0; j < 3; j++) {
264 CM(i, j) = - c1 * u[i] * u[j]
275 Float_t h, hvx, hvz, hvxy, hvxz, hvyz;
283 CM(0, 0) = e + hvx * v.
fX;
284 CM(0, 1) = hvxy - v.
fZ;
285 CM(0, 2) = hvxz + v.
fY;
287 CM(1, 0) = hvxy + v.
fZ;
288 CM(1, 1) = e + h * v.
fY * v.
fY;
289 CM(1, 2) = hvyz - v.
fX;
291 CM(2, 0) = hvxz - v.
fY;
292 CM(2, 1) = hvyz + v.
fX;
293 CM(2, 2) = e + hvz * v.
fZ;
304 for(
int c=0;
c<4; ++
c, col+=4) {
306 for(
int r=0;
r<4; ++
r, ++row)
307 buf[
r] = row[0]*col[0] + row[4]*col[1] + row[8]*col[2] + row[12]*col[3];
308 col[0] = buf[0]; col[1] = buf[1]; col[2] = buf[2]; col[3] = buf[3];
320 for(
int r=0;
r<4; ++
r, ++row) {
322 for(
int c=0;
c<4; ++
c, col+=4)
323 buf[
c] = row[0]*col[0] + row[4]*col[1] + row[8]*col[2] + row[12]*col[3];
324 row[0] = buf[0]; row[4] = buf[1]; row[8] = buf[2]; row[12] = buf[3];
358 fM[
F03] += amount*col[0];
fM[
F13] += amount*col[1];
fM[
F23] += amount*col[2];
367 fM[
F13] += x*fM[1] + y*fM[5] + z*fM[9];
368 fM[
F23] += x*fM[2] + y*fM[6] + z*fM[10];
382 --i1 <<= 2; --i2 <<= 2;
383 for (
int r=0;
r<4; ++
r, ++row) {
384 b1 = cos*row[i1] +
sin*row[i2];
385 b2 = cos*row[i2] -
sin*row[i1];
386 row[i1] = b1; row[i2] = b2;
396 fM[
F03 + --ai] += amount;
422 for(
int c=0;
c<4; ++
c, col+=4) {
423 b1 = cos*col[i1] -
sin*col[i2];
424 b2 = cos*col[i2] +
sin*col[i1];
425 col[i1] = b1; col[i2] = b2;
436 fM[
F03] += amount*vec[0];
437 fM[
F13] += amount*vec[1];
438 fM[
F23] += amount*vec[2];
473 col[0] =
x; col[1] =
y; col[2] = z;
500 v.
SetXYZ(col[0], col[1], col[2]);
564 inline void clamp_angle(
Float_t& a)
576 clamp_angle(a1); clamp_angle(a2); clamp_angle(a3);
604 int n = strspn(pat,
"XxYyZz");
if(n > 3) n = 3;
608 for(
int i=0; i<
n; i++) {
609 if(isupper(pat[i])) a[i] = -a[i];
611 case 'x':
case 'X':
RotateLF(2, 3, a[i]);
break;
612 case 'y':
case 'Y':
RotateLF(3, 1, a[i]);
break;
613 case 'z':
case 'Z':
RotateLF(1, 2, a[i]);
break;
628 if(d>1) d=1;
else if(d<-1) d=-1;
660 return (sx + sy + sz)/3;
792 Double_t t[3] = { v[0], v[1], v[2] };
804 Double_t t[3] = { v[0], v[1], v[2] };
828 c[0] /=
l; c[1] /=
l; c[2] /=
l;
839 const Double_t dp = c[0]*rc[0] + c[1]*rc[1] + c[2]*rc[2];
840 c[0] -= rc[0]*dp; c[1] -= rc[1]*dp; c[2] -= rc[2]*dp;
887 const Double_t det3_012_012 = fM[
F00]*det2_12_12 - fM[
F01]*det2_12_02 + fM[
F02]*det2_12_01;
888 const Double_t det3_012_013 = fM[
F00]*det2_12_13 - fM[
F01]*det2_12_03 + fM[
F03]*det2_12_01;
889 const Double_t det3_012_023 = fM[
F00]*det2_12_23 - fM[
F02]*det2_12_03 + fM[
F03]*det2_12_02;
890 const Double_t det3_012_123 = fM[
F01]*det2_12_23 - fM[
F02]*det2_12_13 + fM[
F03]*det2_12_12;
891 const Double_t det3_013_012 = fM[
F00]*det2_13_12 - fM[
F01]*det2_13_02 + fM[
F02]*det2_13_01;
892 const Double_t det3_013_013 = fM[
F00]*det2_13_13 - fM[
F01]*det2_13_03 + fM[
F03]*det2_13_01;
893 const Double_t det3_013_023 = fM[
F00]*det2_13_23 - fM[
F02]*det2_13_03 + fM[
F03]*det2_13_02;
894 const Double_t det3_013_123 = fM[
F01]*det2_13_23 - fM[
F02]*det2_13_13 + fM[
F03]*det2_13_12;
895 const Double_t det3_023_012 = fM[
F00]*det2_23_12 - fM[
F01]*det2_23_02 + fM[
F02]*det2_23_01;
896 const Double_t det3_023_013 = fM[
F00]*det2_23_13 - fM[
F01]*det2_23_03 + fM[
F03]*det2_23_01;
897 const Double_t det3_023_023 = fM[
F00]*det2_23_23 - fM[
F02]*det2_23_03 + fM[
F03]*det2_23_02;
898 const Double_t det3_023_123 = fM[
F01]*det2_23_23 - fM[
F02]*det2_23_13 + fM[
F03]*det2_23_12;
899 const Double_t det3_123_012 = fM[
F10]*det2_23_12 - fM[
F11]*det2_23_02 + fM[
F12]*det2_23_01;
900 const Double_t det3_123_013 = fM[
F10]*det2_23_13 - fM[
F11]*det2_23_03 + fM[
F13]*det2_23_01;
901 const Double_t det3_123_023 = fM[
F10]*det2_23_23 - fM[
F12]*det2_23_03 + fM[
F13]*det2_23_02;
902 const Double_t det3_123_123 = fM[
F11]*det2_23_23 - fM[
F12]*det2_23_13 + fM[
F13]*det2_23_12;
905 const Double_t det = fM[
F00]*det3_123_123 - fM[
F01]*det3_123_023 +
906 fM[
F02]*det3_123_013 - fM[
F03]*det3_123_012;
909 throw(eh +
"matrix is singular.");
912 const Double_t oneOverDet = 1.0/det;
913 const Double_t mn1OverDet = - oneOverDet;
915 fM[
F00] = det3_123_123 * oneOverDet;
916 fM[
F01] = det3_023_123 * mn1OverDet;
917 fM[
F02] = det3_013_123 * oneOverDet;
918 fM[
F03] = det3_012_123 * mn1OverDet;
920 fM[
F10] = det3_123_023 * mn1OverDet;
921 fM[
F11] = det3_023_023 * oneOverDet;
922 fM[
F12] = det3_013_023 * mn1OverDet;
923 fM[
F13] = det3_012_023 * oneOverDet;
925 fM[
F20] = det3_123_013 * oneOverDet;
926 fM[
F21] = det3_023_013 * mn1OverDet;
927 fM[
F22] = det3_013_013 * oneOverDet;
928 fM[
F23] = det3_012_013 * mn1OverDet;
930 fM[
F30] = det3_123_012 * mn1OverDet;
931 fM[
F31] = det3_023_012 * oneOverDet;
932 fM[
F32] = det3_013_012 * mn1OverDet;
933 fM[
F33] = det3_012_012 * oneOverDet;
942 void TEveTrans::Streamer(
TBuffer &R__b)
958 for(
Int_t i=0; i<4; ++i, ++row)
959 printf(
"%8.3f %8.3f %8.3f | %8.3f\n", row[0], row[4], row[8], row[12]);
969 s.setf(std::ios::fixed, std::ios::floatfield);
971 for(
Int_t i=1; i<=4; i++)
972 for(
Int_t j=1; j<=4; j++)
973 s << t(i,j) << ((j==4) ?
"\n" :
"\t");
1001 m[0] = r[0]*s[0]; m[1] = r[3]*s[0]; m[2] = r[6]*s[0]; m[3] = 0;
1002 m[4] = r[1]*s[1]; m[5] = r[4]*s[1]; m[6] = r[7]*s[1]; m[7] = 0;
1003 m[8] = r[2]*s[2]; m[9] = r[5]*s[2]; m[10] = r[8]*s[2]; m[11] = 0;
1004 m[12] = t[0]; m[13] = t[1]; m[14] = t[2]; m[15] = 1;
1008 m[0] = r[0]; m[1] = r[3]; m[2] = r[6]; m[3] = 0;
1009 m[4] = r[1]; m[5] = r[4]; m[6] = r[7]; m[7] = 0;
1010 m[8] = r[2]; m[9] = r[5]; m[10] = r[8]; m[11] = 0;
1011 m[12] = t[0]; m[13] = t[1]; m[14] = t[2]; m[15] = 1;
1029 r[0] = m[0]/s[0]; r[3] = m[1]/s[0]; r[6] = m[2]/s[0]; m += 4;
1030 r[1] = m[0]/s[1]; r[4] = m[1]/s[1]; r[7] = m[2]/s[1]; m += 4;
1031 r[2] = m[0]/s[2]; r[5] = m[1]/s[2]; r[8] = m[2]/s[2]; m += 4;
1032 t[0] = m[0]; t[1] = m[1]; t[2] = m[2];
1037 r[0] = 1; r[3] = 0; r[6] = 0;
1038 r[1] = 0; r[4] = 1; r[7] = 0;
1039 r[2] = 0; r[5] = 0; r[8] = 1;
1040 s[0] = s[1] = s[2] = 1;
1041 t[0] = t[1] = t[2] = 0;
1056 m[0] =
fM[0]; m[1] =
fM[4]; m[2] =
fM[8]; m[3] =
fM[3];
1057 m[4] =
fM[1]; m[5] =
fM[5]; m[6] =
fM[9]; m[7] =
fM[7];
1058 m[8] =
fM[2]; m[9] =
fM[6]; m[10] =
fM[10]; m[11] =
fM[11];
1059 m[12] =
fM[12]; m[13] =
fM[13]; m[14] =
fM[14]; m[15] =
fM[15];
1078 if (s < low || s > high)
return kTRUE;
1080 if (s < low || s > high)
return kTRUE;
1082 if (s < low || s > high)
return kTRUE;
TEveTrans is a 4x4 transformation matrix for homogeneous coordinates stored internally in a column-ma...
void MultLeft(const TEveTrans &t)
Multiply from left: this = t * this.
TEveTrans()
Default constructor.
Double_t CM(Int_t i, Int_t j) const
void OrtoNorm3()
Orto-norm columns 1 to 3.
TVector3 GetBaseVec(Int_t b) const
Get base-vector with index b.
void MultRight(const TEveTrans &t)
Multiply from right: this = this * t.
void SetScale(Double_t sx, Double_t sy, Double_t sz)
Set scaling.
virtual const Double_t * GetRotationMatrix() const =0
void SetScaleX(Double_t sx)
Change x scaling.
Buffer base class used for serializing objects.
Double_t fLocalMaster[16]
void SetScaleZ(Double_t sz)
Change z scaling.
TT Dot(const TEveVectorT &a) const
void UnitRot()
Reset rotation part of the matrix to unity.
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual const Double_t * GetScale() const
Double_t Unscale()
Remove scaling, make all base vectors of unit length.
TEveTrans operator*(const TEveTrans &t)
Copy, multiply from right and return product.
void Move3LF(Double_t x, Double_t y, Double_t z)
General move in local-frame.
void SetRotByAngles(Float_t a1, Float_t a2, Float_t a3)
Double_t Orto3Column(Int_t col, Int_t ref)
Orto-norm 3-vector in column col with respect to column ref.
void Move(const TEveTrans &a, Int_t ai, Double_t amount)
Move in a's coord-system along axis-index ai.
Bool_t IsScale(Double_t low=0.9, Double_t high=1.1) const
Test if the transformation is a scale.
void SetXYZ(Double_t x, Double_t y, Double_t z)
void ZeroTrans(Double_t w=1.0)
Reset matrix to zero, only the perspective scaling is set to w (1 by default).
if(pyself &&pyself!=Py_None)
virtual const Double_t * GetRotationMatrix() const
Double_t ATan2(Double_t, Double_t)
void Scale(Double_t sx, Double_t sy, Double_t sz)
Scale matrix. Translation part untouched.
void SetScaleY(Double_t sy)
Change y scaling.
TEveVectorT Cross(const TEveVectorT &a) const
void SetupRotation(Int_t i, Int_t j, Double_t f)
Setup the matrix as an elementary rotation.
void SetFrom(Double_t *carr)
void RotatePF(Int_t i1, Int_t i2, Double_t amount)
Rotate in parent frame. Does optimised version of MultLeft.
void SetBuffer3D(TBuffer3D &buff)
Fill transformation part TBuffer3D core section.
void TransposeRotationPart()
Transpose 3x3 rotation sub-matrix.
void SetBaseVec(Int_t b, Double_t x, Double_t y, Double_t z)
Set base-vector with index b.
void SetRotByAnyAngles(Float_t a1, Float_t a2, Float_t a3, const char *pat)
Sets Rotation part as given by angles a1, a1, a3 and pattern pat.
Generic 3D primitive description class.
void GetRotAngles(Float_t *x) const
Get Cardan rotation angles (pattern xYz above).
TVector3 Multiply(const TVector3 &v, Double_t w=1) const
Multiply vector and return it.
void SetTrans(const TEveTrans &t, Bool_t copyAngles=kTRUE)
Set matrix from another,.
std::ostream & operator<<(std::ostream &s, const TEveTrans &t)
Print to std::ostream.
void GetScale(Double_t &sx, Double_t &sy, Double_t &sz) const
Deduce scales from sizes of base vectors.
void MovePF(Int_t ai, Double_t amount)
Move in parent-frame along axis index ai.
void SetPos(Double_t x, Double_t y, Double_t z)
Set position (base-vec 4).
void MoveLF(Int_t ai, Double_t amount)
Move in local-frame along axis with index ai.
virtual void Print(Option_t *option="") const
Print in reasonable format.
ClassImp(TMCParticle) void TMCParticle printf(": p=(%7.3f,%7.3f,%9.3f) ;", fPx, fPy, fPz)
Double_t Invert()
Invert matrix.
void RotateLF(Int_t i1, Int_t i2, Double_t amount)
Rotate in local frame. Does optimised version of MultRight.
void GetXYZ(Double_t *carray) const
Mother of all ROOT objects.
void RotateIP(TVector3 &v) const
Rotate vector in-place. Translation is NOT applied.
void SetupFromToVec(const TEveVector &from, const TEveVector &to)
A function for creating a rotation matrix that rotates a vector called "from" into another vector cal...
void SetGeoHMatrix(TGeoHMatrix &mat)
Set TGeoHMatrix mat.
void Move3(const TEveTrans &a, Double_t x, Double_t y, Double_t z)
General move in a's coord-system.
Exception class thrown by TEve classes and macros.
TEveVectorT & Sub(const TEveVectorT &a, const TEveVectorT &b)
void Rotate(const TEveTrans &a, Int_t i1, Int_t i2, Double_t amount)
Rotate in a's coord-system, rotating base vector with index i1 into i2.
void MultiplyIP(TVector3 &v, Double_t w=1) const
Multiply vector in-place.
Double_t Sqrt(Double_t x)
void SetFromArray(const Double_t arr[16])
Set matrix from Double_t array.
Double_t Norm3Column(Int_t col)
Norm 3-vector in column col.
virtual const Double_t * GetTranslation() const
virtual const Double_t * GetTranslation() const =0
virtual const Double_t * GetScale() const =0
void UnitTrans()
Reset matrix to unity.
void Move3PF(Double_t x, Double_t y, Double_t z)
General move in parent-frame.