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 /// (this way is used when running simulation and building geometry
14 /// - from Root geometry file (*.root) or ASCII file (*.dat) using
15 /// the method LoadGeometryData(const TString& fileName)
16 /// - from geometry loaded in AliGeomManager using
17 /// the method LoadGeometryData() without arguments
19 /// \author Ivana Hrivnacova, IPN Orsay
21 #ifndef ALI_MUON_GEOMETRY_TRANSFORMER_H
22 #define ALI_MUON_GEOMETRY_TRANSFORMER_H
25 #include <TObjArray.h>
26 #include <TGeoMatrix.h>
28 class AliMUONGeometryModuleTransformer;
29 class AliMUONGeometryDetElement;
36 class AliMUONGeometryTransformer : public TObject
39 AliMUONGeometryTransformer();
40 AliMUONGeometryTransformer(TRootIOCtor* /*ioCtor*/);
41 virtual ~AliMUONGeometryTransformer();
44 void AddModuleTransformer(AliMUONGeometryModuleTransformer* transformer);
45 void AddMisAlignModule(Int_t moduleId, const TGeoHMatrix& matrix, Bool_t bGlobal = kTRUE);
46 void AddMisAlignDetElement(Int_t detElemId, const TGeoHMatrix& matrix, Bool_t bGlobal = kTRUE);
49 void AddAlignableVolumes() const;
50 TClonesArray* CreateZeroAlignmentData() const;
51 void ClearMisAlignmentData();
55 Bool_t LoadTransformations();
56 Bool_t LoadGeometryData(const TString& fileName);
57 Bool_t LoadGeometryData();
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;
79 void SetDetName(const TString& detName);
80 void SetOwner(Bool_t isOwner);
84 Int_t GetNofModuleTransformers() const;
85 const AliMUONGeometryModuleTransformer* GetModuleTransformer(
86 Int_t index, Bool_t warn = true) const;
88 const AliMUONGeometryModuleTransformer* GetModuleTransformerByDEId(
89 Int_t detElemId, Bool_t warn = true) const;
91 const AliMUONGeometryDetElement* GetDetElement(
92 Int_t detElemId, Bool_t warn = true) const;
94 const TClonesArray* GetMisAlignmentData() const;
96 Bool_t HasDE(Int_t detElemId) const;
98 AliMpArea* GetDEArea(Int_t detElemId) const;
102 AliMUONGeometryTransformer(const AliMUONGeometryTransformer& right);
104 AliMUONGeometryTransformer& operator = (const AliMUONGeometryTransformer& right);
108 static const TString& GetDefaultDetectorName();
112 void CreateDEAreas() const;
114 Bool_t LoadMapping() const;
115 AliMUONGeometryModuleTransformer* GetModuleTransformerNonConst(
116 Int_t index, Bool_t warn = true) const;
118 TGeoHMatrix GetTransform(
119 Double_t x, Double_t y, Double_t z,
120 Double_t a1, Double_t a2, Double_t a3,
121 Double_t a4, Double_t a5, Double_t a6) const;
123 void FillModuleTransform(Int_t moduleId,
124 Double_t x, Double_t y, Double_t z,
125 Double_t a1, Double_t a2, Double_t a3,
126 Double_t a4, Double_t a5, Double_t a6);
127 void FillDetElemTransform(Int_t id,
128 Double_t x, Double_t y, Double_t z,
129 Double_t a1, Double_t a2, Double_t a3,
130 Double_t a4, Double_t a5, Double_t a6);
132 TString ReadModuleTransforms(ifstream& in);
133 TString ReadDetElemTransforms(ifstream& in);
134 Bool_t ReadTransformations(const TString& fileName);
136 void WriteTransform(ofstream& out, const TGeoMatrix* transform) const;
137 void WriteModuleTransforms(ofstream& out) const;
138 void WriteDetElemTransforms(ofstream& out) const;
140 TString GetModuleSymName(Int_t moduleId) const;
141 TString GetDESymName(Int_t detElemId) const;
144 TString fDetectorName; ///< Detector name
145 TObjArray* fModuleTransformers; ///< array of module transformers
146 TClonesArray* fMisAlignArray; ///< array of misalignment data
147 mutable AliMpExMap* fDEAreas; ///< areas of detection elements in global coordinates
149 ClassDef(AliMUONGeometryTransformer,4) // Geometry parametrisation
154 /// Return the number of contained module transformers
155 inline Int_t AliMUONGeometryTransformer::GetNofModuleTransformers() const
156 { return fModuleTransformers->GetEntriesFast(); }
158 /// Return the array of misalignment data
159 inline const TClonesArray* AliMUONGeometryTransformer::GetMisAlignmentData() const
160 { return fMisAlignArray; }
162 /// Set detector name
163 inline void AliMUONGeometryTransformer::SetDetName(const TString& detName)
164 { fDetectorName = detName; }
166 /// Set ownership of array module transformers
167 inline void AliMUONGeometryTransformer::SetOwner(Bool_t isOwner)
168 { fModuleTransformers->SetOwner(isOwner); }
170 #endif //ALI_MUON_GEOMETRY_TRANSFORMER_H