1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
7 /// \class AliMUONGeometryTransformer
8 /// \brief Top container class for geometry transformations
10 /// Geometry transformations can be filled in these ways:
11 /// - by geometry builder when geometry is built via builders
12 /// - from Root geometry file (*.root) or Root geometry manager
13 /// - from ASCII file (*.dat)
14 /// If geometry is loaded from a file, the list of aligned volume paths
15 /// has to be read first from volpaths.dat file.
17 /// \author Ivana Hrivnacova, IPN Orsay
19 #ifndef ALI_MUON_GEOMETRY_TRANSFORMER_H
20 #define ALI_MUON_GEOMETRY_TRANSFORMER_H
23 #include <TObjArray.h>
24 #include <TGeoMatrix.h>
26 class AliMUONGeometryModuleTransformer;
27 class AliMUONGeometryDetElement;
32 class AliMUONGeometryTransformer : public TObject
35 AliMUONGeometryTransformer(Bool_t isOwner, const TString& detName = "MUON");
36 AliMUONGeometryTransformer();
37 virtual ~AliMUONGeometryTransformer();
40 void AddModuleTransformer(AliMUONGeometryModuleTransformer* transformer);
41 void AddMisAlignModule(Int_t moduleId, const TGeoHMatrix& matrix);
42 void AddMisAlignDetElement(Int_t detElemId, const TGeoHMatrix& matrix);
44 void AddAlignableVolumes() const;
45 TClonesArray* CreateZeroAlignmentData() const;
49 Bool_t ReadGeometryData(const TString& volPathFileName,
50 const TString& transformFileName);
51 Bool_t ReadGeometryData(const TString& volPathFileName,
52 TGeoManager* geoManager);
54 Bool_t WriteGeometryData(const TString& volPathFileName,
55 const TString& transformFileName,
56 const TString& misalignFileName = "") const;
58 Bool_t WriteVolumePaths(const TString& fileName) const;
59 Bool_t WriteTransformations(const TString& fileName) const;
60 Bool_t WriteMisAlignmentData(const TString& fileName) const;
62 // Transformation methods
64 void Global2Local(Int_t detElemId,
65 Float_t xg, Float_t yg, Float_t zg,
66 Float_t& xl, Float_t& yl, Float_t& zl) const;
67 void Global2Local(Int_t detElemId,
68 Double_t xg, Double_t yg, Double_t zg,
69 Double_t& xl, Double_t& yl, Double_t& zl) const;
71 void Local2Global(Int_t detElemId,
72 Float_t xl, Float_t yl, Float_t zl,
73 Float_t& xg, Float_t& yg, Float_t& zg) const;
74 void Local2Global(Int_t detElemId,
75 Double_t xl, Double_t yl, Double_t zl,
76 Double_t& xg, Double_t& yg, Double_t& zg) const;
80 Int_t GetNofModuleTransformers() const;
81 const AliMUONGeometryModuleTransformer* GetModuleTransformer(
82 Int_t index, Bool_t warn = true) const;
84 const AliMUONGeometryModuleTransformer* GetModuleTransformerByDEId(
85 Int_t detElemId, Bool_t warn = true) const;
87 const AliMUONGeometryDetElement* GetDetElement(
88 Int_t detElemId, Bool_t warn = true) const;
90 const TClonesArray* GetMisAlignmentData() const;
92 Bool_t HasDE(Int_t detElemId) const;
95 AliMUONGeometryTransformer(const AliMUONGeometryTransformer& right);
96 AliMUONGeometryTransformer& operator = (const AliMUONGeometryTransformer& right);
100 AliMUONGeometryModuleTransformer* GetModuleTransformerNonConst(
101 Int_t index, Bool_t warn = true) const;
103 TGeoHMatrix GetTransform(
104 Double_t x, Double_t y, Double_t z,
105 Double_t a1, Double_t a2, Double_t a3,
106 Double_t a4, Double_t a5, Double_t a6) const;
108 void FillModuleVolPath(Int_t moduleId, const TString& volPath);
109 void FillDetElemVolPath(Int_t detElemId, const TString& volPath);
111 void FillModuleTransform(Int_t moduleId,
112 Double_t x, Double_t y, Double_t z,
113 Double_t a1, Double_t a2, Double_t a3,
114 Double_t a4, Double_t a5, Double_t a6);
115 void FillDetElemTransform(Int_t id,
116 Double_t x, Double_t y, Double_t z,
117 Double_t a1, Double_t a2, Double_t a3,
118 Double_t a4, Double_t a5, Double_t a6);
120 Bool_t ReadVolPaths(ifstream& in);
121 TString ReadModuleTransforms(ifstream& in);
122 TString ReadDetElemTransforms(ifstream& in);
123 Bool_t LoadTransforms(TGeoManager* tgeoManager);
125 Bool_t ReadVolPaths(const TString& fileName);
126 Bool_t ReadTransformations(const TString& fileName);
127 Bool_t ReadTransformations2(const TString& fileName);
129 void WriteTransform(ofstream& out, const TGeoMatrix* transform) const;
130 void WriteModuleVolPaths(ofstream& out) const;
131 void WriteDetElemVolPaths(ofstream& out) const;
132 void WriteModuleTransforms(ofstream& out) const;
133 void WriteDetElemTransforms(ofstream& out) const;
135 TString GetModuleSymName(Int_t moduleId) const;
136 TString GetDESymName(Int_t detElemId) const;
139 TString fDetectorName; /// < Detector name
140 TObjArray* fModuleTransformers; ///< array of module transformers
141 TClonesArray* fMisAlignArray; ///< array of misalignment data
143 ClassDef(AliMUONGeometryTransformer,2) // Geometry parametrisation
148 /// Return the number of contained module transformers
149 inline Int_t AliMUONGeometryTransformer::GetNofModuleTransformers() const
150 { return fModuleTransformers->GetEntriesFast(); }
152 /// Return the array of misalignment data
153 inline const TClonesArray* AliMUONGeometryTransformer::GetMisAlignmentData() const
154 { return fMisAlignArray; }
156 #endif //ALI_MUON_GEOMETRY_TRANSFORMER_H