#ifndef ALIALIGNOBJ_H
#define ALIALIGNOBJ_H
-/*************************************************************************
- * AliAlignObj: alignment base class for the storage of alignment *
- * information for a single volume, that is a translation, a rotation *
- * and a the identity of the volume itself in form of a TGeo path and *
- * as a unique integer identifier *
- *************************************************************************/
-#include "TObject.h"
-#include "TString.h"
-#include "TGeoMatrix.h"
+//************************************************************************
+// AliAlignObj: alignment base class for the storage of the alignment *
+// constants for a single volume: *
+// - a displacement (a shift and a rotation) either as *
+// - the 6 doubles which identify it or as *
+// - the matrix which identifies it *
+// - the identity of the volume itself in form of a symbolic volume *
+// name for alignable volumes, in form of a TGeo path otherwise, *
+// and as a unique integer identifier *
+//************************************************************************
+#include <TObject.h>
+#include <TString.h>
+#include <Rtypes.h>
+#include "AliGeomManager.h"
+#include <TMatrixD.h>
+#include <TMatrixDSym.h>
+
+
+class AliTrackPoint;
+class AliTrackPointArray;
+class TGeoMatrix;
class AliAlignObj : public TObject {
public:
AliAlignObj();
+ AliAlignObj(const char* symname, UShort_t voluid);
+ AliAlignObj(const char* symname, UShort_t voluid, Double_t* cmat);
+ AliAlignObj(const char* symname, AliGeomManager::ELayerID detId, Int_t modId);
AliAlignObj(const AliAlignObj& theAlignObj);
AliAlignObj& operator= (const AliAlignObj& theAlignObj);
+ AliAlignObj& operator*=(const AliAlignObj& theAlignObj);
virtual ~AliAlignObj();
//Setters
virtual void SetTranslation(Double_t x, Double_t y, Double_t z) = 0;
virtual void SetTranslation(const TGeoMatrix& m) = 0;
+ virtual Bool_t SetLocalTranslation(Double_t x, Double_t y, Double_t z);
+ virtual Bool_t SetLocalTranslation(const TGeoMatrix& m);
virtual void SetRotation(Double_t psi, Double_t theta, Double_t phi) = 0;
virtual Bool_t SetRotation(const TGeoMatrix& m) = 0;
+ virtual Bool_t SetLocalRotation(Double_t psi, Double_t theta, Double_t phi);
+ virtual Bool_t SetLocalRotation(const TGeoMatrix& m);
virtual void SetPars(Double_t x, Double_t y, Double_t z, Double_t psi,
- Double_t theta, Double_t phi) = 0;
- virtual void SetMatrix(const TGeoMatrix& m) = 0;
- void SetVolPath(const TString& volpath) {fVolPath=volpath;}
- void SetVolUID(const UShort_t voluid) {fVolUID=voluid;}
+ Double_t theta, Double_t phi);
+ virtual Bool_t SetLocalPars(Double_t x, Double_t y, Double_t z,
+ Double_t psi, Double_t theta, Double_t phi);
+ virtual Bool_t SetMatrix(const TGeoMatrix& m);
+ virtual Bool_t SetLocalMatrix(const TGeoMatrix& m);
+ void SetSymName(const TString& symname) {fVolPath=symname;}
+ void SetVolUID(UShort_t voluid) {fVolUID=voluid;}
+ void SetVolUID(AliGeomManager::ELayerID layerId, Int_t modId);
+ void SetCorrMatrix(Double_t *cov);
+ void SetCorrMatrix(TMatrixDSym& mcov);
+ Bool_t SetFromLocalCov(Double_t *lCov);
+ Bool_t SetFromLocalCov(TMatrixDSym& lCov);
//Getters
- const char *GetVolPath() const {return fVolPath.Data();}
+ const char *GetSymName() const {return fVolPath.Data();}
UShort_t GetVolUID() const {return fVolUID;}
+ void GetVolUID(AliGeomManager::ELayerID &layerId, Int_t &modId) const;
virtual void GetTranslation(Double_t* tr) const=0;
+ virtual Bool_t GetLocalTranslation(Double_t* tr) const;
virtual Bool_t GetAngles(Double_t* angles) const=0;
- virtual void GetPars(Double_t transl[], Double_t rot[]) const=0;
+ virtual Bool_t GetLocalAngles(Double_t* angles) const;
+ virtual Bool_t GetPars(Double_t transl[], Double_t angles[]) const;
+ virtual Bool_t GetLocalPars(Double_t transl[], Double_t angles[]) const;
virtual void GetMatrix(TGeoHMatrix& m) const=0;
+ virtual Bool_t GetLocalMatrix(TGeoHMatrix& m) const;
+ void GetCovMatrix(Double_t *cov) const;
+ void GetCovMatrix(TMatrixDSym& mcov) const;
+ Bool_t GetJacobian(TMatrixD& mJ) const;
+ Bool_t GetLocalCovMatrix(Double_t *cov) const;
+ Bool_t GetLocalCovMatrix(TMatrixDSym& lCov) const;
+
+ Bool_t IsSortable() const {return kTRUE;}
+ Int_t GetLevel() const;
+ virtual Int_t Compare(const TObject* obj) const;
+
+ virtual AliAlignObj& Inverse() const=0;
+
+ void Transform(AliTrackPoint &p, Bool_t copycov=kFALSE) const;
+ void Transform(AliTrackPointArray &array) const;
void Print(Option_t *) const;
+ static Int_t LayerSize(Int_t layerId) {return AliGeomManager::LayerSize(layerId);}
+ static const char* LayerName(Int_t layerId) {return AliGeomManager::LayerName(layerId);}
+
+ static UShort_t LayerToVolUID(AliGeomManager::ELayerID layerId, Int_t modId) {return AliGeomManager::LayerToVolUID(layerId, modId);}
+ static UShort_t LayerToVolUID(Int_t layerId, Int_t modId) {return AliGeomManager::LayerToVolUID(layerId, modId);}
+ static AliGeomManager::ELayerID VolUIDToLayer(UShort_t voluid, Int_t &modId) {return AliGeomManager::VolUIDToLayer(voluid, modId);}
+ static AliGeomManager::ELayerID VolUIDToLayer(UShort_t voluid) {return AliGeomManager::VolUIDToLayer(voluid);}
+ static UShort_t LayerToVolUIDSafe(AliGeomManager::ELayerID layerId, Int_t modId) {return AliGeomManager::LayerToVolUIDSafe(layerId, modId);}
+ static UShort_t LayerToVolUIDSafe(Int_t layerId, Int_t modId) {return AliGeomManager::LayerToVolUIDSafe(layerId, modId);}
+ static AliGeomManager::ELayerID VolUIDToLayerSafe(UShort_t voluid, Int_t &modId) {return AliGeomManager::VolUIDToLayerSafe(voluid, modId);}
+ static AliGeomManager::ELayerID VolUIDToLayerSafe(UShort_t voluid) {return AliGeomManager::VolUIDToLayerSafe(voluid);}
+
+ static const char* SymName(UShort_t voluid) {return AliGeomManager::SymName(voluid);}
+ static const char* SymName(AliGeomManager::ELayerID layerId, Int_t modId) {return AliGeomManager::SymName(layerId, modId);}
+
+ Bool_t ApplyToGeometry(Bool_t ovlpcheck=kFALSE);
+ static Bool_t GetFromGeometry(const char *symname, AliAlignObj &alobj) {return AliGeomManager::GetFromGeometry(symname, alobj);}
+
+ static AliAlignObj* GetAlignObj(UShort_t voluid) {return AliGeomManager::GetAlignObj(voluid);}
+ static AliAlignObj* GetAlignObj(AliGeomManager::ELayerID layerId, Int_t modId) {return AliGeomManager::GetAlignObj(layerId, modId);}
+
+ static Bool_t GetOrigGlobalMatrix(const char *symname, TGeoHMatrix &m) {return AliGeomManager::GetOrigGlobalMatrix(symname, m);}
+
protected:
void AnglesToMatrix(const Double_t *angles, Double_t *rot) const;
Bool_t MatrixToAngles(const Double_t *rot, Double_t *angles) const;
//Volume identifiers
- TString fVolPath; // Volume path inside TGeo geometry
- UShort_t fVolUID; // Unique volume ID
-
- ClassDef(AliAlignObj, 1)
-};
+ TString fVolPath; // Symbolic volume name; in case could coincide with
+ // the volume path inside TGeo geometry (for non-alignable volumes)
+ UShort_t fVolUID; // Unique volume ID
-/*****************************************************************************
- * AliAlignObjAngles: derived alignment class storing alignment information *
- * for a single volume in form of three doubles for the translation *
- * and three doubles for the rotation expressed with the euler angles *
- * in the xyz-convention (http://mathworld.wolfram.com/EulerAngles.html), *
- * also known as roll, pitch, yaw. PLEASE NOTE THE ANGLES SIGNS ARE *
- * INVERSE WITH RESPECT TO THIS REFERENCE!!! In this way the representation*
- * is fully consistent with the TGeo Rotation methods. *
- *****************************************************************************/
-class AliAlignObjAngles : public AliAlignObj{
- public:
- AliAlignObjAngles();
- AliAlignObjAngles(const AliAlignObjAngles& theAlignObj);
- AliAlignObjAngles& operator= (const AliAlignObjAngles& theAlignObj);
- virtual ~AliAlignObjAngles();
-
- //Setters
- virtual void SetTranslation(Double_t x, Double_t y, Double_t z){
- fTranslation[0]=x; fTranslation[1]=y; fTranslation[2]=z;}
- virtual void SetTranslation(const TGeoMatrix& m);
- virtual void SetRotation(Double_t psi, Double_t theta, Double_t phi){
- fRotation[0]=psi; fRotation[1]=theta; fRotation[2]=phi;}
- virtual Bool_t SetRotation(const TGeoMatrix& m);
- virtual void SetMatrix(const TGeoMatrix& m);
- virtual void SetPars(Double_t x, Double_t y, Double_t z, Double_t psi,
- Double_t theta, Double_t phi){
- fTranslation[0]=x; fTranslation[1]=y; fTranslation[2]=z;
- fRotation[0]=psi; fRotation[1]=theta; fRotation[2]=phi;}
-
- //Getters
- virtual void GetTranslation(Double_t *tr) const {
- tr[0] = fTranslation[0]; tr[1] = fTranslation[1]; tr[2] = fTranslation[2];}
- virtual Bool_t GetAngles(Double_t* angles) const {
- angles[0] = fRotation[0]; angles[1] = fRotation[1];
- angles[2] = fRotation[2]; return kTRUE;}
- virtual void GetPars(Double_t tr[], Double_t angles[]) const;
- virtual void GetMatrix(TGeoHMatrix& m) const;
-
- protected:
- Double_t fTranslation[3]; // Translation vector
- Double_t fRotation[3]; // Roll-pitch-yaw angles
-
- ClassDef(AliAlignObjAngles, 1)
-};
-
-/**************************************************************************
- * AliAlignObjMatrix: derived alignment class storing alignment *
- * information for a single volume in form of TGeoHMatrix, which *
- * includes the information for a translation, a rotation and a scale *
- *************************************************************************/
-class AliAlignObjMatrix : public AliAlignObj {
- public:
- AliAlignObjMatrix();
- AliAlignObjMatrix(const AliAlignObjMatrix& theAlignObj);
- AliAlignObjMatrix& operator= (const AliAlignObjMatrix& theAlignObj);
- virtual ~AliAlignObjMatrix();
-
- //Setters
- virtual void SetTranslation(Double_t x, Double_t y, Double_t z);
- virtual void SetTranslation(const TGeoMatrix& m);
- virtual void SetRotation(Double_t psi, Double_t theta, Double_t phi);
- virtual Bool_t SetRotation(const TGeoMatrix& m);
- virtual void SetMatrix(const TGeoMatrix& m);
- virtual void SetPars(Double_t x, Double_t y, Double_t z, Double_t psi,
- Double_t theta, Double_t phi);
- //Getters
- virtual void GetTranslation(Double_t* tr) const;
- virtual Bool_t GetAngles(Double_t* angles) const;
- virtual void GetPars(Double_t tr[], Double_t rot[]) const;
- virtual void GetMatrix(TGeoHMatrix& m) const;
-
- protected:
- TGeoHMatrix fMatrix; // Transformation matrix
-
- ClassDef(AliAlignObjMatrix, 1)
+ Double32_t fDiag[6]; // diagonal elements of the correlation matrix for shifts and rotations (dx,dy,dz,dpsi,dtheta,dphi)
+ Double32_t fODia[15]; // [-1, 1,8] off-diagonal elements (in 8 bit precision) of the correlation matrix
+
+ ClassDef(AliAlignObj, 3)
};
#endif