// name for alignable volumes, in form of a TGeo path otherwise, *
// and as a unique integer identifier *
//************************************************************************
-#include "TObject.h"
-#include "TString.h"
-#include "TGeoMatrix.h"
+#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:
- enum ELayerID{kInvalidLayer=0,
- kFirstLayer=1,
- kSPD1=1, kSPD2=2,
- kSDD1=3, kSDD2=4,
- kSSD1=5, kSSD2=6,
- kTPC1=7, kTPC2=8,
- kTRD1=9, kTRD2=10, kTRD3=11, kTRD4=12, kTRD5=13, kTRD6=14,
- kTOF=15,
- kPHOS1=16, kPHOS2=17,
- kHMPID=18,
- kMUON=19,
- kLastLayer=20};
AliAlignObj();
AliAlignObj(const char* symname, UShort_t voluid);
- AliAlignObj(const char* symname, ELayerID detId, Int_t modId);
+ 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 Bool_t SetLocalMatrix(const TGeoMatrix& m);
void SetSymName(const TString& symname) {fVolPath=symname;}
void SetVolUID(UShort_t voluid) {fVolUID=voluid;}
- void SetVolUID(ELayerID layerId, Int_t modId);
+ 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 *GetSymName() const {return fVolPath.Data();}
UShort_t GetVolUID() const {return fVolUID;}
- void GetVolUID(ELayerID &layerId, Int_t &modId) const;
+ 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 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 AliAlignObj& Inverse() const=0;
- void Transform(AliTrackPoint &p) const;
+ 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);
- static const char* LayerName(Int_t layerId);
+ 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(ELayerID layerId, Int_t modId);
- static UShort_t LayerToVolUID(Int_t layerId, Int_t modId);
- static ELayerID VolUIDToLayer(UShort_t voluid, Int_t &modId);
- static ELayerID VolUIDToLayer(UShort_t voluid);
+ 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);
- static const char* SymName(ELayerID layerId, Int_t modId);
+ 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();
- static Bool_t GetFromGeometry(const char *symname, AliAlignObj &alobj);
+ 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);
- static AliAlignObj* GetAlignObj(ELayerID layerId, Int_t modId);
+ 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);
+ 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;
- static void InitSymNames();
- static void InitAlignObjFromGeometry();
-
//Volume identifiers
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
- static Int_t fgLayerSize[kLastLayer - kFirstLayer]; // Size of layers
- static const char* fgLayerName[kLastLayer - kFirstLayer]; // Name of layers
-
- static TString* fgVolPath[kLastLayer - kFirstLayer]; // Symbolic volume names
- static AliAlignObj** fgAlignObjs[kLastLayer - kFirstLayer]; // Alignment objects
-
- ClassDef(AliAlignObj, 2)
+ 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