]>
Commit | Line | Data |
---|---|---|
75678bb3 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * SigmaEffect_thetadegrees * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpeateose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | // $Id$ | |
17 | // | |
18 | // Class AliMUONGeometry | |
19 | // ---------------------------- | |
20 | // Manager class for geometry construction via geometry builders. | |
21 | // | |
22 | // Author: Ivana Hrivnacova, IPN Orsay | |
23 | ||
24 | #include <iostream> | |
25 | ||
26 | #include <TObjArray.h> | |
27 | #include <Riostream.h> | |
28 | #include <TSystem.h> | |
29 | ||
30 | #include "AliMUONGeometry.h" | |
31 | #include "AliMUONGeometryTransformer.h" | |
32 | #include "AliMUONGeometryModule.h" | |
33 | #include "AliMUONGeometrySVMap.h" | |
e5df8ea1 | 34 | #include "AliMUONGeometryStore.h" |
75678bb3 | 35 | #include "AliLog.h" |
36 | ||
37 | ||
38 | ClassImp(AliMUONGeometry) | |
39 | ||
40 | //______________________________________________________________________________ | |
41 | AliMUONGeometry::AliMUONGeometry(Bool_t isOwner) | |
42 | : TObject(), | |
43 | fModules(0), | |
44 | fTransformer(0) | |
45 | ||
46 | { | |
47 | /// Standard constructor | |
48 | ||
49 | // Create array for geometry modules | |
50 | fModules = new TObjArray(); | |
51 | fModules->SetOwner(isOwner); | |
52 | ||
53 | // Geometry parametrisation | |
54 | fTransformer = new AliMUONGeometryTransformer(false); | |
55 | } | |
56 | ||
57 | //______________________________________________________________________________ | |
58 | AliMUONGeometry::AliMUONGeometry() | |
59 | : TObject(), | |
60 | fModules(0), | |
61 | fTransformer(0) | |
62 | { | |
63 | /// Default constructor | |
64 | } | |
65 | ||
66 | //______________________________________________________________________________ | |
67 | AliMUONGeometry::AliMUONGeometry(const AliMUONGeometry& right) | |
68 | : TObject(right) | |
69 | { | |
70 | /// Copy constructor (not implemented) | |
71 | ||
72 | AliFatal("Copy constructor not provided."); | |
73 | } | |
74 | ||
75 | //______________________________________________________________________________ | |
76 | AliMUONGeometry::~AliMUONGeometry() | |
77 | { | |
78 | /// Destructor | |
79 | ||
80 | delete fModules; | |
81 | delete fTransformer; | |
82 | } | |
83 | ||
84 | //______________________________________________________________________________ | |
85 | AliMUONGeometry& | |
86 | AliMUONGeometry::operator=(const AliMUONGeometry& right) | |
87 | { | |
88 | /// Assignement operator (not implemented) | |
89 | ||
90 | // check assignement to self | |
91 | if (this == &right) return *this; | |
92 | ||
93 | AliFatal("Assignement operator not provided."); | |
94 | ||
95 | return *this; | |
96 | } | |
97 | ||
98 | // | |
99 | // private methods | |
100 | // | |
101 | ||
102 | //______________________________________________________________________________ | |
103 | TString AliMUONGeometry::ComposePath(const TString& volName, | |
104 | Int_t copyNo) const | |
105 | { | |
106 | // Compose path from given volName and copyNo | |
107 | // --- | |
108 | ||
109 | TString path(volName); | |
110 | path += "."; | |
111 | path += copyNo; | |
112 | ||
113 | return path; | |
114 | } | |
115 | ||
116 | //______________________________________________________________________________ | |
117 | void AliMUONGeometry::FillData3(const TString& sensVolumePath, | |
118 | Int_t detElemId) | |
119 | { | |
120 | // Fill the mapping of the sensitive volume path to the detection element. | |
121 | // --- | |
122 | ||
123 | // Module Id | |
e5df8ea1 | 124 | Int_t moduleId = AliMUONGeometryStore::GetModuleId(detElemId); |
75678bb3 | 125 | |
126 | // Get module | |
127 | AliMUONGeometryModule* module | |
128 | = (AliMUONGeometryModule*)fModules->At(moduleId); | |
129 | ||
130 | if ( !module ) { | |
131 | AliWarningStream() | |
132 | << "Geometry module for det element " << detElemId << " not defined." | |
133 | << endl; | |
134 | return; | |
135 | } | |
136 | ||
137 | // Get module sensitive volumes map | |
138 | AliMUONGeometrySVMap* svMap = module->GetSVMap(); | |
139 | ||
140 | // Map the sensitive volume to detection element | |
141 | svMap->Add(sensVolumePath, detElemId); | |
142 | } | |
143 | ||
144 | //______________________________________________________________________________ | |
145 | TString AliMUONGeometry::ReadData3(ifstream& in) | |
146 | { | |
147 | // Reads SV maps from a file | |
148 | // Returns true, if reading finished correctly. | |
149 | // --- | |
150 | ||
151 | TString key("SV"); | |
152 | while ( key == TString("SV") ) { | |
153 | ||
154 | // Input data | |
155 | TString volumePath; | |
156 | Int_t detElemId; | |
157 | ||
158 | in >> volumePath; | |
159 | in >> detElemId; | |
160 | ||
161 | //cout << "volumePath=" << volumePath << " " | |
162 | // << "detElemId=" << detElemId | |
163 | // << endl; | |
164 | ||
165 | // Fill data | |
166 | FillData3(volumePath, detElemId); | |
167 | ||
168 | // Go to next line | |
169 | in >> key; | |
170 | } | |
171 | ||
172 | return key; | |
173 | } | |
174 | ||
175 | //______________________________________________________________________________ | |
176 | void AliMUONGeometry::WriteData3(ofstream& out) const | |
177 | { | |
178 | // Writes association of sensitive volumes and detection elements | |
179 | // from the sensitive volume map | |
180 | // --- | |
181 | ||
182 | for (Int_t i=0; i<fModules->GetEntriesFast(); i++) { | |
183 | AliMUONGeometryModule* geometry | |
184 | = (AliMUONGeometryModule*)fModules->At(i); | |
185 | AliMUONGeometrySVMap* svMap | |
186 | = geometry->GetSVMap(); | |
187 | ||
188 | svMap->WriteMap(out); | |
189 | out << endl; | |
190 | } | |
191 | } | |
192 | ||
193 | // | |
194 | // public functions | |
195 | // | |
196 | ||
197 | //_____________________________________________________________________________ | |
198 | void AliMUONGeometry::AddModule(AliMUONGeometryModule* module) | |
199 | { | |
200 | /// Add the geometrymodule to the array | |
201 | ||
202 | fModules->Add(module); | |
203 | ||
204 | if (module) | |
205 | fTransformer->AddModuleTransformer(module->GetTransformer()); | |
206 | } | |
207 | ||
208 | //______________________________________________________________________________ | |
209 | Bool_t | |
210 | AliMUONGeometry::ReadSVMap(const TString& fileName) | |
211 | { | |
212 | // Reads the sensitive volume maps from a file | |
213 | // Returns true, if reading finished correctly. | |
214 | // --- | |
215 | ||
216 | // No reading | |
217 | // if builder is not associated with any geometry module | |
218 | if (fModules->GetEntriesFast() == 0) return false; | |
219 | ||
220 | // File path | |
221 | TString filePath = gSystem->Getenv("ALICE_ROOT"); | |
222 | filePath += "/MUON/data/"; | |
223 | filePath += fileName; | |
224 | ||
225 | // Open input file | |
226 | ifstream in(filePath, ios::in); | |
227 | if (!in) { | |
228 | cerr << filePath << endl; | |
229 | AliFatal("File not found."); | |
230 | return false; | |
231 | } | |
232 | ||
233 | TString key; | |
234 | in >> key; | |
235 | while ( !in.eof() ) { | |
236 | if (key == TString("SV")) | |
237 | key = ReadData3(in); | |
238 | else { | |
239 | AliFatal(Form("%s key not recognized", key.Data())); | |
240 | return false; | |
241 | } | |
242 | } | |
243 | ||
244 | return true; | |
245 | } | |
246 | ||
247 | //______________________________________________________________________________ | |
248 | Bool_t | |
249 | AliMUONGeometry::WriteSVMap(const TString& fileName) const | |
250 | { | |
251 | // Writes sensitive volume map into a file | |
252 | // Returns true, if writing finished correctly. | |
253 | // --- | |
254 | ||
255 | // No writing | |
256 | // if builder is not associated with any geometry module | |
257 | if (fModules->GetEntriesFast() == 0) return false; | |
258 | ||
259 | // File path | |
260 | TString filePath = gSystem->Getenv("ALICE_ROOT"); | |
261 | filePath += "/MUON/data/"; | |
262 | filePath += fileName; | |
263 | ||
264 | // Open input file | |
265 | ofstream out(filePath, ios::out); | |
266 | if (!out) { | |
267 | cerr << filePath << endl; | |
268 | AliError("File not found."); | |
269 | return false; | |
270 | } | |
271 | #if !defined (__DECCXX) | |
272 | out.setf(std::ios::fixed); | |
273 | #endif | |
274 | WriteData3(out); | |
275 | ||
276 | return true; | |
277 | } | |
278 | ||
279 | //_____________________________________________________________________________ | |
280 | const AliMUONGeometryModule* | |
281 | AliMUONGeometry::GetModule(Int_t index, Bool_t warn) const | |
282 | { | |
283 | /// Return the geometry module specified by index | |
284 | ||
285 | if (index < 0 || index >= fModules->GetEntriesFast()) { | |
286 | if (warn) { | |
287 | AliWarningStream() | |
288 | << "Index: " << index << " outside limits" << std::endl; | |
289 | } | |
290 | return 0; | |
291 | } | |
292 | ||
293 | return (const AliMUONGeometryModule*) fModules->At(index); | |
294 | } | |
295 | ||
296 | //_____________________________________________________________________________ | |
297 | const AliMUONGeometryModule* | |
298 | AliMUONGeometry::GetModuleByDEId(Int_t detElemId, Bool_t warn) const | |
299 | { | |
300 | /// Return the geometry module specified by index | |
301 | ||
302 | // Get module index | |
e5df8ea1 | 303 | Int_t index = AliMUONGeometryStore::GetModuleId(detElemId); |
75678bb3 | 304 | |
305 | return GetModule(index, warn); | |
306 | } |