Bug fix. Removed delete statement
[u/mrichter/AliRoot.git] / STEER / AliAlignObj.h
index 7a6bd38..adb147b 100644 (file)
 //      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);
@@ -59,12 +53,16 @@ class AliAlignObj : public TObject {
   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;
@@ -73,6 +71,11 @@ class AliAlignObj : public TObject {
   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;
@@ -80,48 +83,48 @@ class AliAlignObj : public TObject {
 
   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(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 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 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);}
 
-  static const char* SymName(UShort_t voluid);
-  static const char* SymName(ELayerID layerId, Int_t modId);
+  Bool_t ApplyToGeometry(Bool_t ovlpcheck=kFALSE);
+  static Bool_t   GetFromGeometry(const char *symname, AliAlignObj &alobj) {return AliGeomManager::GetFromGeometry(symname, alobj);}
 
-  Bool_t ApplyToGeometry();
-  static Bool_t   GetFromGeometry(const char *symname, AliAlignObj &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 AliAlignObj* GetAlignObj(UShort_t voluid);
-  static AliAlignObj* GetAlignObj(ELayerID layerId, Int_t 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;
 
-  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