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 }