8e73b1c9cd6401965bcec4462028136e8a513d73
[u/mrichter/AliRoot.git] / STEER / AliAlignObjAngles.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 //-----------------------------------------------------------------
17 //   Implementation of the alignment object class through
18 //   the concrete representation of alignment object class
19 //   AliAlignObjAngles derived from the base class AliAlignObj
20 //-----------------------------------------------------------------
21
22 #include "AliAlignObj.h"
23 #include "AliAlignObjAngles.h"
24 //#include "AliLog.h"
25
26 ClassImp(AliAlignObjAngles)
27
28 //_____________________________________________________________________________
29 AliAlignObjAngles::AliAlignObjAngles() : AliAlignObj()
30 {
31   // default constructor
32   //
33   fTranslation[0]=fTranslation[1]=fTranslation[2]=0.;
34   fRotation[0]=fRotation[1]=fRotation[2]=0.;
35 }
36
37 //_____________________________________________________________________________
38 AliAlignObjAngles::AliAlignObjAngles(const char* volpath, UShort_t volUId, Double_t x, Double_t y, Double_t z, Double_t psi, Double_t theta, Double_t phi, Bool_t global) throw (const Char_t *) : AliAlignObj(volpath,volUId)
39 {
40   // standard constructor with 3 translation + 3 rotation parameters
41   // If the user explicitly sets the global variable to kFALSE then the
42   // parameters are interpreted as giving the local transformation.
43   // This requires to have a gGeoMenager active instance, otherwise the
44   // constructor will fail (no object created)
45   // 
46   if(global){
47     SetPars(x, y, z, psi, theta, phi);
48   }else{
49     if(!SetLocalPars(x,y,z,psi,theta,phi)) throw "Alignment object creation failed (TGeo instance needed)!\n";
50   }
51 }
52
53 //_____________________________________________________________________________
54 AliAlignObjAngles::AliAlignObjAngles(const char* volpath, UShort_t volUId, TGeoMatrix& m, Bool_t global) throw (const Char_t *) : AliAlignObj(volpath,volUId)
55 {
56   // standard constructor with TGeoMatrix
57   // If the user explicitly sets the global variable to kFALSE then the
58   // parameters are interpreted as giving the local transformation.
59   // This requires to have a gGeoMenager active instance, otherwise the
60   // constructor will fail (no object created)
61   //
62
63   if (!SetMatrix(m)) throw "Alignment object creation failed (can't extract roll-pitch-yall angles from the matrix)!\n";
64
65   if (!global) {
66     if (!SetLocalPars(fTranslation[0],fTranslation[1],fTranslation[2],fRotation[0],fRotation[1],fRotation[2])) throw "Alignment object creation failed (TGeo instance needed)!\n";
67   }
68 }
69
70 //_____________________________________________________________________________
71 AliAlignObjAngles::AliAlignObjAngles(const AliAlignObj& theAlignObj) :
72   AliAlignObj(theAlignObj)
73 {
74   // copy constructor
75   //
76   Double_t tr[3];
77   theAlignObj.GetTranslation(tr);
78   SetTranslation(tr[0],tr[1],tr[2]);
79   Double_t rot[3];
80   if (theAlignObj.GetAngles(rot))
81     SetRotation(rot[0],rot[1],rot[2]);
82 }
83
84 //_____________________________________________________________________________
85 AliAlignObjAngles &AliAlignObjAngles::operator =(const AliAlignObj& theAlignObj)
86 {
87   // assignment operator
88   //
89   if(this==&theAlignObj) return *this;
90   ((AliAlignObj *)this)->operator=(theAlignObj);
91
92   Double_t tr[3];
93   theAlignObj.GetTranslation(tr);
94   SetTranslation(tr[0],tr[1],tr[2]);
95   Double_t rot[3];
96   if (theAlignObj.GetAngles(rot))
97     SetRotation(rot[0],rot[1],rot[2]);
98
99   return *this;
100 }
101
102 //_____________________________________________________________________________
103 AliAlignObjAngles::~AliAlignObjAngles()
104 {
105   // default destructor
106   //
107 }
108
109 //_____________________________________________________________________________
110 void AliAlignObjAngles::SetTranslation(const TGeoMatrix& m)
111 {
112   if(m.IsTranslation()){
113     const Double_t* tr = m.GetTranslation();
114     fTranslation[0]=tr[0];  fTranslation[1]=tr[1]; fTranslation[2]=tr[2];
115   }else{
116     fTranslation[0] = fTranslation[1] = fTranslation[2] = 0.;
117   }
118 }
119
120 //_____________________________________________________________________________
121 Bool_t AliAlignObjAngles::SetRotation(const TGeoMatrix& m)
122 {
123   if(m.IsRotation()){
124     const Double_t* rot = m.GetRotationMatrix();
125     return MatrixToAngles(rot,fRotation);
126   }else{
127     fRotation[0] = fRotation[1] = fRotation[2] = 0.;
128     return kTRUE;
129   }
130 }
131
132 //_____________________________________________________________________________
133 void AliAlignObjAngles::GetMatrix(TGeoHMatrix& m) const
134 {
135   m.SetTranslation(&fTranslation[0]);
136   Double_t rot[9];
137   AnglesToMatrix(fRotation,rot);
138   m.SetRotation(rot);
139 }
140
141 //_____________________________________________________________________________
142 AliAlignObj& AliAlignObjAngles::Inverse() const
143 {
144   // Return a temporary inverse of the alignment
145   // object. This means 'mis
146    static AliAlignObjAngles a;
147    a = *this;
148
149    TGeoHMatrix m;
150    GetMatrix(m);
151    a.SetMatrix(m.Inverse());
152
153    return a;
154 }