New constructors added to AliAlignObj and the derived classes (R.Grosso)
[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) : AliAlignObj(volpath,volUID)
39 {
40   // standard constructor with 3 translation + 3 rotation parameters
41   //
42   fTranslation[0]=x; fTranslation[1]=y; fTranslation[2]=z;
43   fRotation[0]=psi; fRotation[1]=theta; fRotation[2]=phi;
44 }
45
46 //_____________________________________________________________________________
47 AliAlignObjAngles::AliAlignObjAngles(const char* volpath, ELayerID layerId, Int_t volId, 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,layerId,volId)
48 {
49   // standard constructor with 3 translation + 3 rotation parameters
50   // If the user explicitly sets the global variable to kFALSE then the
51   // parameters are interpreted as giving the local transformation.
52   // This requires to have a gGeoMenager active instance, otherwise the
53   // constructor will fail (no object created)
54   // 
55   if(global){
56     fTranslation[0]=x; fTranslation[1]=y; fTranslation[2]=z;
57     fRotation[0]=psi; fRotation[1]=theta; fRotation[2]=phi;
58   }else{
59     if(!SetLocalPars(x,y,z,psi,theta,phi)) throw "Alignment object creation failed (TGeo instance needed)!\n";
60   }
61 }
62
63 //_____________________________________________________________________________
64 AliAlignObjAngles::AliAlignObjAngles(const char* volpath, UShort_t voluid, TGeoMatrix& m) : AliAlignObj()
65 {
66   // standard constructor with TGeoMatrix
67   //
68   fVolPath=volpath;
69   fVolUID=voluid;
70   SetTranslation(m);
71   SetRotation(m);
72 }
73
74 //_____________________________________________________________________________
75 AliAlignObjAngles::AliAlignObjAngles(const char* volpath, ELayerID layerId, Int_t volId, TGeoMatrix& m) : AliAlignObj(volpath,layerId,volId)
76 {
77   // standard constructor with TGeoMatrix
78   //
79   SetTranslation(m);
80   SetRotation(m);
81 }
82
83 //_____________________________________________________________________________
84 AliAlignObjAngles::AliAlignObjAngles(const AliAlignObj& theAlignObj) :
85   AliAlignObj(theAlignObj)
86 {
87   // copy constructor
88   //
89   Double_t tr[3];
90   theAlignObj.GetTranslation(tr);
91   SetTranslation(tr[0],tr[1],tr[2]);
92   Double_t rot[3];
93   theAlignObj.GetAngles(rot);
94   SetRotation(rot[0],rot[1],rot[2]);
95 }
96
97 //_____________________________________________________________________________
98 AliAlignObjAngles &AliAlignObjAngles::operator =(const AliAlignObj& theAlignObj)
99 {
100   // assignment operator
101   //
102   if(this==&theAlignObj) return *this;
103   ((AliAlignObj *)this)->operator=(theAlignObj);
104
105   Double_t tr[3];
106   theAlignObj.GetTranslation(tr);
107   SetTranslation(tr[0],tr[1],tr[2]);
108   Double_t rot[3];
109   theAlignObj.GetAngles(rot);
110   SetRotation(rot[0],rot[1],rot[2]);
111   return *this;
112 }
113
114 //_____________________________________________________________________________
115 AliAlignObjAngles::~AliAlignObjAngles()
116 {
117   // default destructor
118   //
119 }
120
121 //_____________________________________________________________________________
122 void AliAlignObjAngles::SetTranslation(const TGeoMatrix& m)
123 {
124   if(m.IsTranslation()){
125     const Double_t* tr = m.GetTranslation();
126     fTranslation[0]=tr[0];  fTranslation[1]=tr[1]; fTranslation[2]=tr[2];
127   }else{
128 //     AliWarning("Argument matrix is not a translation! Setting zero-translation.");
129     fTranslation[0] = fTranslation[1] = fTranslation[2] = 0.;
130   }
131 }
132
133 //_____________________________________________________________________________
134 Bool_t AliAlignObjAngles::SetRotation(const TGeoMatrix& m)
135 {
136   if(m.IsRotation()){
137     const Double_t* rot = m.GetRotationMatrix();
138     return MatrixToAngles(rot,fRotation);
139   }else{
140 //     AliWarning("Argument matrix is not a rotation! Setting yaw-pitch-roll to zero.");
141     fRotation[0] = fRotation[1] = fRotation[2] = 0.;
142     return kTRUE;
143   }
144 }
145
146 //_____________________________________________________________________________
147 void AliAlignObjAngles::SetMatrix(const TGeoMatrix& m)
148 {
149   SetTranslation(m);
150   SetRotation(m);
151 }
152
153 //_____________________________________________________________________________
154 void AliAlignObjAngles::GetPars(Double_t tr[], Double_t angles[]) const
155 {
156   GetTranslation(tr);
157   GetAngles(angles);
158 }
159
160 //_____________________________________________________________________________
161 void AliAlignObjAngles::GetMatrix(TGeoHMatrix& m) const
162 {
163   m.SetTranslation(&fTranslation[0]);
164   Double_t rot[9];
165   AnglesToMatrix(fRotation,rot);
166   m.SetRotation(rot);
167 }
168
169 //_____________________________________________________________________________
170 AliAlignObj& AliAlignObjAngles::Inverse() const
171 {
172   // Return a temporary inverse of the alignment
173   // object. This means 'mis
174    static AliAlignObjAngles a;
175    a = *this;
176
177    TGeoHMatrix m;
178    GetMatrix(m);
179    a.SetMatrix(m.Inverse());
180
181    return a;
182 }