STEER/AliAlignObjMatrix.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 derived concrete representation of alignment object class:
  19 //   AliAlignObjMatrix derived from the base class AliAlignObj
  20 //-----------------------------------------------------------------
  21
  22 #include "AliAlignObj.h"
  23 #include "AliAlignObjMatrix.h"
  24 //#include "AliLog.h"
  25
  26 ClassImp(AliAlignObjMatrix)
  27
  28 //_____________________________________________________________________________
  29 AliAlignObjMatrix::AliAlignObjMatrix() :
  30   AliAlignObj(),
  31   fMatrix()
  32 {
  33   // Default constructor
  34   //
  35 }
  36
  37 //_____________________________________________________________________________
  38 AliAlignObjMatrix::AliAlignObjMatrix(const char* volpath, UShort_t voluid, Double_t x, Double_t y, Double_t z, Double_t psi, Double_t theta, Double_t phi) :
  39   AliAlignObj(volpath, voluid),
  40   fMatrix()
  41 {
  42   // standard constructor with 3 translation + 3 rotation parameters
  43   //
  44   SetTranslation(x, y, z);
  45   SetRotation(psi, theta, phi);
  46 }
  47
  48 //_____________________________________________________________________________
  49 AliAlignObjMatrix::AliAlignObjMatrix(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 *) :
  50   AliAlignObj(volpath,layerId,volId),
  51   fMatrix()
  52 {
  53   // standard constructor with 3 translation + 3 rotation parameters
  54   // If the user explicitly sets the global variable to kFALSE then the
  55   // parameters are interpreted as giving the local transformation.
  56   // This requires to have a gGeoMenager active instance, otherwise the
  57   // constructor will fail (no object created)
  58   //
  59   if(global){
  60     SetTranslation(x, y, z);
  61     SetRotation(psi, theta, phi);
  62   }else{
  63     if(!SetLocalPars(x,y,z,psi,theta,phi)) throw "Alignment object creation failed (TGeo instance needed)!\n";
  64   }
  65 }
  66
  67
  68 //_____________________________________________________________________________
  69 AliAlignObjMatrix::AliAlignObjMatrix(const char* volpath, UShort_t volUID, TGeoMatrix& m) :
  70   AliAlignObj(volpath,volUID),
  71   fMatrix()
  72 {
  73   // standard constructor with TGeoMatrix
  74   //
  75   SetTranslation(m);
  76   SetRotation(m);
  77 }
  78
  79 //_____________________________________________________________________________
  80 AliAlignObjMatrix::AliAlignObjMatrix(const char* volpath, ELayerID layerId, Int_t volId, TGeoMatrix& m) :
  81   AliAlignObj(volpath,layerId,volId),
  82   fMatrix()
  83 {
  84   // standard constructor with TGeoMatrix
  85   //
  86   SetTranslation(m);
  87   SetRotation(m);
  88 }
  89
  90 //_____________________________________________________________________________
  91 AliAlignObjMatrix::AliAlignObjMatrix(const AliAlignObj& theAlignObj) :
  92   AliAlignObj(theAlignObj),
  93   fMatrix()
  94 {
  95   //copy constructor
  96   //
  97   Double_t tr[3];
  98   theAlignObj.GetTranslation(tr);
  99   SetTranslation(tr[0],tr[1],tr[2]);
 100   Double_t rot[3];
 101   theAlignObj.GetAngles(rot);
 102   SetRotation(rot[0],rot[1],rot[2]);
 103 }
 104
 105 //_____________________________________________________________________________
 106 AliAlignObjMatrix &AliAlignObjMatrix::operator =(const AliAlignObj& theAlignObj)
 107 {
 108   // assignment operator
 109   //
 110   if(this==&theAlignObj) return *this;
 111   ((AliAlignObj *)this)->operator=(theAlignObj);
 112   Double_t tr[3];
 113   theAlignObj.GetTranslation(tr);
 114   SetTranslation(tr[0],tr[1],tr[2]);
 115   Double_t rot[3];
 116   theAlignObj.GetAngles(rot);
 117   SetRotation(rot[0],rot[1],rot[2]);
 118   return *this;
 119 }
 120
 121 //_____________________________________________________________________________
 122 AliAlignObjMatrix::~AliAlignObjMatrix()
 123 {
 124   // Destructor
 125   //
 126 }
 127
 128 //_____________________________________________________________________________
 129 void AliAlignObjMatrix::SetTranslation(Double_t x, Double_t y, Double_t z)
 130 {
 131   Double_t tr[3];
 132   tr[0]=x; tr[1]=y; tr[2]=z;
 133   fMatrix.SetTranslation(tr);
 134 }
 135
 136 //_____________________________________________________________________________
 137 void AliAlignObjMatrix::SetTranslation(const TGeoMatrix& m)
 138 {
 139   const Double_t *tr = m.GetTranslation();
 140   fMatrix.SetTranslation(tr);
 141 }
 142
 143 //_____________________________________________________________________________
 144 void AliAlignObjMatrix::SetRotation(Double_t psi, Double_t theta, Double_t phi)
 145 {
 146   Double_t angles[3] = {psi, theta, phi};
 147   Double_t rot[9];
 148   AnglesToMatrix(angles,rot);
 149   fMatrix.SetRotation(rot);
 150 }
 151
 152 //_____________________________________________________________________________
 153 Bool_t AliAlignObjMatrix::SetRotation(const TGeoMatrix& m)
 154 {
 155   const Double_t* rot = m.GetRotationMatrix();
 156   fMatrix.SetRotation(rot);
 157   return kTRUE;
 158 }
 159
 160 //_____________________________________________________________________________
 161 void AliAlignObjMatrix::SetMatrix(const TGeoMatrix& m)
 162 {
 163   // Set rotation matrix and translation
 164   // using TGeoMatrix
 165   SetTranslation(m);
 166   SetRotation(m);
 167 }
 168
 169 //_____________________________________________________________________________
 170 void AliAlignObjMatrix::SetPars(Double_t x, Double_t y, Double_t z,
 171                        Double_t psi, Double_t theta, Double_t phi)
 172 {
 173   // Set rotation matrix and translation
 174   // using 3 angles and 3 translations
 175   SetTranslation(x,y,z);
 176   SetRotation(psi,theta,phi);
 177 }
 178
 179 //_____________________________________________________________________________
 180 void AliAlignObjMatrix::GetTranslation(Double_t *tr) const
 181 {
 182   // Get Translation from TGeoMatrix
 183   const Double_t* translation = fMatrix.GetTranslation();
 184   tr[0] = translation[0];
 185   tr[1] = translation[1];
 186   tr[2] = translation[2];
 187 }
 188
 189 //_____________________________________________________________________________
 190 Bool_t AliAlignObjMatrix::GetAngles(Double_t *angles) const
 191 {
 192   // Get rotation angles from the TGeoHMatrix
 193   const Double_t* rot = fMatrix.GetRotationMatrix();
 194   return MatrixToAngles(rot,angles);
 195 }
 196
 197 //_____________________________________________________________________________
 198 void AliAlignObjMatrix::GetPars(Double_t tr[], Double_t angles[]) const
 199 {
 200   GetTranslation(tr);
 201   GetAngles(angles);
 202 }
 203
 204 //_____________________________________________________________________________
 205 void AliAlignObjMatrix::GetMatrix(TGeoHMatrix& m) const
 206 {
 207   // Get TGeoHMatrix
 208   //
 209   const Double_t *tr = fMatrix.GetTranslation();
 210   m.SetTranslation(tr);
 211   const Double_t *rot = fMatrix.GetRotationMatrix();
 212   m.SetRotation(rot);
 213 }
 214
 215 //_____________________________________________________________________________
 216 AliAlignObj& AliAlignObjMatrix::Inverse() const
 217 {
 218   // Return a temporary inverse of the alignment
 219   // object. This means 'mis
 220    static AliAlignObjMatrix a;
 221    a = *this;
 222
 223    TGeoHMatrix m;
 224    GetMatrix(m);
 225    a.SetMatrix(m.Inverse());
 226
 227    return a;
 228 }