| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
26 // 23 //
>> 24 // $Id: G4Scene.cc,v 1.13 2003/06/16 17:14:15 gunter Exp $
>> 25 // GEANT4 tag $Name: geant4-05-02-patch-01 $
27 // 26 //
28 // 27 //
29 // Scene data John Allison 19th July 1996. 28 // Scene data John Allison 19th July 1996.
30 29
31 #include "G4Scene.hh" 30 #include "G4Scene.hh"
32 31
33 #include "G4Vector3D.hh" 32 #include "G4Vector3D.hh"
34 #include "G4BoundingExtentScene.hh" << 33 #include "G4BoundingSphereScene.hh"
35 #include "G4VisAttributes.hh" 34 #include "G4VisAttributes.hh"
36 #include "G4PhysicalVolumeModel.hh" 35 #include "G4PhysicalVolumeModel.hh"
37 #include "G4TransportationManager.hh" 36 #include "G4TransportationManager.hh"
38 37
39 #include <set> <<
40 <<
41 #define G4warn G4cout <<
42 <<
43 G4Scene::G4Scene (const G4String& name): 38 G4Scene::G4Scene (const G4String& name):
44 fName (name), 39 fName (name),
45 fRefreshAtEndOfEvent(true), << 40 fRefreshAtEndOfEvent(true)
46 fRefreshAtEndOfRun(true), <<
47 fMaxNumberOfKeptEvents(100) <<
48 {} // Note all other data members have default 41 {} // Note all other data members have default initial values.
49 42
50 G4Scene::~G4Scene () {} 43 G4Scene::~G4Scene () {}
51 44
52 namespace { << 45 G4bool G4Scene::AddRunDurationModel (G4VModel* pModel, G4bool warn) {
53 void PrintInvalidModel(const G4VModel* model << 46 G4int i, nModels = fRunDurationModelList.size ();
54 { << 47 for (i = 0; i < nModels; i++) {
55 G4ExceptionDescription ed; << 48 if (pModel -> GetGlobalDescription () ==
56 ed << "Invalid model \"" << model->GetGlob << 49 fRunDurationModelList [i] -> GetGlobalDescription ()) break;
57 << "\".\n Not included in extent calculat <<
58 G4Exception <<
59 ("G4Scene::CalculateExtent", <<
60 "visman0201", JustWarning, ed); <<
61 } <<
62 } <<
63 <<
64 void G4Scene::CalculateExtent () <<
65 { <<
66 G4BoundingExtentScene boundingExtentScene; <<
67 <<
68 for (size_t i = 0; i < fRunDurationModelList <<
69 if (fRunDurationModelList[i].fActive) { <<
70 G4VModel* model = fRunDurationModelList[ <<
71 if (model -> Validate()) { <<
72 const G4VisExtent& thisExtent = model <<
73 if (thisExtent != G4VisExtent::GetNull <<
74 boundingExtentScene.AccrueBoundingEx <<
75 } <<
76 } else { <<
77 PrintInvalidModel(model); <<
78 } <<
79 } <<
80 } 50 }
81 << 51 if (i < nModels) {
82 for (size_t i = 0; i < fEndOfEventModelList. << 52 if (warn) {
83 if (fEndOfEventModelList[i].fActive) { << 53 G4cout << "G4Scene::AddRunDurationModel: model \""
84 G4VModel* model = fEndOfEventModelList[i << 54 << pModel -> GetGlobalDescription ()
85 if (model -> Validate()) { << 55 << "\"\n is already in the run-duration list of scene \""
86 const G4VisExtent& thisExtent = model << 56 << fName
87 if (thisExtent != G4VisExtent::GetNull << 57 << "\"."
88 boundingExtentScene.AccrueBoundingEx << 58 << G4endl;
89 } <<
90 } else { <<
91 PrintInvalidModel(model); <<
92 } <<
93 } 59 }
>> 60 return false;
94 } 61 }
>> 62 fRunDurationModelList.push_back (pModel);
>> 63 CalculateExtent ();
>> 64 return true;
>> 65 }
95 66
96 for (size_t i = 0; i < fEndOfRunModelList.si << 67 void G4Scene::CalculateExtent () {
97 if (fEndOfRunModelList[i].fActive) { << 68 G4int nModels = fRunDurationModelList.size ();
98 G4VModel* model = fEndOfRunModelList[i]. << 69 G4BoundingSphereScene boundingSphereScene;
99 if (model -> Validate()) { << 70 for (G4int i = 0; i < nModels; i++) {
100 const G4VisExtent& thisExtent = model << 71 const G4VisExtent& thisExtent =
101 if (thisExtent != G4VisExtent::GetNull << 72 fRunDurationModelList[i] -> GetExtent ();
102 boundingExtentScene.AccrueBoundingEx << 73 G4Point3D thisCentre = thisExtent.GetExtentCentre ();
103 } << 74 G4double thisRadius = thisExtent.GetExtentRadius ();
104 } else { << 75 thisCentre.transform (fRunDurationModelList[i] -> GetTransformation ());
105 PrintInvalidModel(model); << 76 boundingSphereScene.AccrueBoundingSphere (thisCentre, thisRadius);
106 } <<
107 } <<
108 } 77 }
109 << 78 fExtent = boundingSphereScene.GetBoundingSphereExtent ();
110 fExtent = boundingExtentScene.GetBoundingExt <<
111 fStandardTargetPoint = fExtent.GetExtentCent 79 fStandardTargetPoint = fExtent.GetExtentCentre ();
112 if (fExtent.GetExtentRadius() <= 0.) { <<
113 G4Exception <<
114 ("G4Scene::CalculateExtent", <<
115 "visman0202", JustWarning, <<
116 "Scene has no extent. Please activate or <<
117 "\nThe camera needs to have something to <<
118 "\nAdd a volume. (You may need \"/r <<
119 "\nOr use \"/vis/scene/add/extent\" <<
120 "\n\"/vis/scene/list\" to see list of mod <<
121 } <<
122 } 80 }
123 81
124 G4bool G4Scene::AddWorldIfEmpty (G4bool warn) 82 G4bool G4Scene::AddWorldIfEmpty (G4bool warn) {
125 G4bool successful = true; 83 G4bool successful = true;
126 if (IsEmpty ()) { 84 if (IsEmpty ()) {
127 successful = false; 85 successful = false;
128 G4VPhysicalVolume* pWorld = 86 G4VPhysicalVolume* pWorld =
129 G4TransportationManager::GetTransportati 87 G4TransportationManager::GetTransportationManager ()
130 -> GetNavigatorForTracking () -> GetWorl 88 -> GetNavigatorForTracking () -> GetWorldVolume ();
131 if (pWorld) { 89 if (pWorld) {
132 const G4VisAttributes* pVisAttribs = 90 const G4VisAttributes* pVisAttribs =
133 pWorld -> GetLogicalVolume () -> GetVisAttri 91 pWorld -> GetLogicalVolume () -> GetVisAttributes ();
134 if (!pVisAttribs || pVisAttribs -> IsVis 92 if (!pVisAttribs || pVisAttribs -> IsVisible ()) {
135 if (warn) { 93 if (warn) {
136 G4warn << << 94 G4cout <<
137 "Your \"world\" has no vis attributes or 95 "Your \"world\" has no vis attributes or is marked as visible."
138 "\n For a better view of the contents, 96 "\n For a better view of the contents, mark the world as"
139 " invisible, e.g.," 97 " invisible, e.g.,"
140 "\n myWorldLogicalVol ->" 98 "\n myWorldLogicalVol ->"
141 " SetVisAttributes (G4VisAttributes::GetIn << 99 " SetVisAttributes (G4VisAttributes::Invisible);"
142 << G4endl; 100 << G4endl;
143 } 101 }
144 } 102 }
145 successful = AddRunDurationModel (new G4 103 successful = AddRunDurationModel (new G4PhysicalVolumeModel (pWorld));
146 // Note: default depth and no modeling p 104 // Note: default depth and no modeling parameters.
147 if (successful) { 105 if (successful) {
148 if (warn) { 106 if (warn) {
149 G4warn << << 107 G4cout <<
150 "G4Scene::AddWorldIfEmpty: The scene had n << 108 "G4Scene::AddWorldIfEmpty: The scene was empty,"
151 "\n \"world\" has been added."; << 109 "\n \"world\" has been added.";
152 G4warn << G4endl; << 110 G4cout << G4endl;
153 } 111 }
154 } 112 }
155 } 113 }
156 } 114 }
157 return successful; 115 return successful;
158 } 116 }
159 117
160 G4bool G4Scene::AddRunDurationModel (G4VModel* <<
161 { <<
162 std::vector<Model>::const_iterator i; <<
163 for (i = fRunDurationModelList.begin (); <<
164 i != fRunDurationModelList.end (); ++i) <<
165 if (pModel -> GetGlobalDescription () == <<
166 i->fpModel->GetGlobalDescription ()) break; <<
167 } <<
168 if (i != fRunDurationModelList.end ()) { <<
169 if (warn) { <<
170 G4warn << "G4Scene::AddRunDurationModel: <<
171 << pModel -> GetGlobalDescription () <<
172 << "\"\n is already in the run-duration <<
173 << fName <<
174 << "\"." <<
175 << G4endl; <<
176 } <<
177 return false; <<
178 } <<
179 fRunDurationModelList.push_back (Model(pMode <<
180 CalculateExtent (); <<
181 return true; <<
182 } <<
183 <<
184 G4bool G4Scene::AddEndOfEventModel (G4VModel* 118 G4bool G4Scene::AddEndOfEventModel (G4VModel* pModel, G4bool warn) {
185 std::size_t i, nModels = fEndOfEventModelLis << 119 G4int i, nModels = fEndOfEventModelList.size ();
186 for (i = 0; i < nModels; ++i) { << 120 for (i = 0; i < nModels; i++) {
187 if (pModel -> GetGlobalDescription () == 121 if (pModel -> GetGlobalDescription () ==
188 fEndOfEventModelList[i].fpModel -> GetGlobal << 122 fEndOfEventModelList [i] -> GetGlobalDescription ()) break;
189 } 123 }
190 if (i < nModels) { 124 if (i < nModels) {
>> 125 delete fEndOfEventModelList[i];
>> 126 fEndOfEventModelList[i] = pModel;
191 if (warn) { 127 if (warn) {
192 G4warn << "G4Scene::AddEndOfEventModel: << 128 G4cout << "G4Scene::AddEndOfEventModel: a model \""
193 << pModel -> GetGlobalDescription () 129 << pModel -> GetGlobalDescription ()
194 << "\"\n is already in the end-of-even << 130 << "\"\n is already in the run-duration list of scene \""
195 << fName << "\"." << 131 << fName <<
>> 132 "\".\n The old model has been deleted; this new model replaces it."
196 << G4endl; 133 << G4endl;
197 } 134 }
198 return false; 135 return false;
199 } 136 }
200 fEndOfEventModelList.push_back (Model(pModel << 137 fEndOfEventModelList.push_back (pModel);
201 CalculateExtent (); <<
202 return true; 138 return true;
203 } 139 }
204 140
205 G4bool G4Scene::AddEndOfRunModel (G4VModel* pM << 141 void G4Scene::Clear () {
206 std::size_t i, nModels = fEndOfRunModelList. << 142 size_t i;
207 for (i = 0; i < nModels; ++i) { << 143 for (i = 0; i < fRunDurationModelList.size(); ++i) {
208 if (pModel -> GetGlobalDescription () == << 144 delete fRunDurationModelList[i];
209 fEndOfRunModelList[i].fpModel -> GetGlobalDe <<
210 } 145 }
211 if (i < nModels) { << 146 for (i = 0; i < fEndOfEventModelList.size(); ++i) {
212 if (warn) { << 147 delete fEndOfEventModelList[i];
213 G4warn << "G4Scene::AddEndOfRunModel: a <<
214 << pModel -> GetGlobalDescription () <<
215 << "\"\n is already in the end-of-run <<
216 << fName << "\"." <<
217 << G4endl; <<
218 } <<
219 return false; <<
220 } 148 }
221 fEndOfRunModelList.push_back (pModel); <<
222 CalculateExtent (); <<
223 return true; <<
224 } 149 }
225 150
226 std::ostream& operator << (std::ostream& os, c << 151 std::ostream& operator << (std::ostream& os, const G4Scene& s) {
227 152
228 size_t i; 153 size_t i;
229 154
230 os << "Scene data:"; 155 os << "Scene data:";
231 156
232 os << "\n Run-duration model list:"; 157 os << "\n Run-duration model list:";
233 if (scene.fRunDurationModelList.size () == 0 << 158 for (i = 0; i < s.fRunDurationModelList.size (); i++) {
234 os << " none"; << 159 os << "\n " << *(s.fRunDurationModelList[i]);
235 } <<
236 for (i = 0; i < scene.fRunDurationModelList. <<
237 if (scene.fRunDurationModelList[i].fActive <<
238 else os << "\n Inactive: "; <<
239 os << *(scene.fRunDurationModelList[i].fpM <<
240 } 160 }
241 161
242 os << "\n End-of-event model list:"; 162 os << "\n End-of-event model list:";
243 if (scene.fEndOfEventModelList.size () == 0) << 163 for (i = 0; i < s.fEndOfEventModelList.size (); i++) {
244 os << " none"; << 164 os << "\n " << *(s.fEndOfEventModelList[i]);
245 } <<
246 for (i = 0; i < scene.fEndOfEventModelList.s <<
247 if (scene.fEndOfEventModelList[i].fActive) <<
248 else os << "\n Inactive: "; <<
249 os << *(scene.fEndOfEventModelList[i].fpMo <<
250 } 165 }
251 166
252 os << "\n End-of-run model list:"; << 167 os << "\n Extent or bounding box: " << s.fExtent;
253 if (scene.fEndOfRunModelList.size () == 0) { <<
254 os << " none"; <<
255 } <<
256 for (i = 0; i < scene.fEndOfRunModelList.siz <<
257 if (scene.fEndOfRunModelList[i].fActive) o <<
258 else os << "\n Inactive: "; <<
259 os << *(scene.fEndOfRunModelList[i].fpMode <<
260 } <<
261 168
262 os << "\n Overall extent or bounding box: " << 169 os << "\n Standard target point: " << s.fStandardTargetPoint;
263 <<
264 os << "\n Standard target point: " << scen <<
265 170
266 os << "\n End of event action set to \""; 171 os << "\n End of event action set to \"";
267 if (scene.fRefreshAtEndOfEvent) os << "refre << 172 if (s.fRefreshAtEndOfEvent) os << "refresh";
268 else { <<
269 os << "accumulate (maximum number of kept <<
270 if (scene.fMaxNumberOfKeptEvents >= 0) os <<
271 else os << "unlimited"; <<
272 os << ")"; <<
273 } <<
274 <<
275 os << "\n End of run action set to \""; <<
276 if (scene.fRefreshAtEndOfRun) os << "refresh <<
277 else os << "accumulate"; 173 else os << "accumulate";
278 os << "\""; 174 os << "\"";
279 175
280 return os; 176 return os;
281 } 177 }
282 178
283 G4bool G4Scene::operator != (const G4Scene& sc << 179 G4bool G4Scene::operator != (const G4Scene& s) const {
284 if ( 180 if (
285 (fRunDurationModelList.size () != 181 (fRunDurationModelList.size () !=
286 scene.fRunDurationModelList.size ()) << 182 s.fRunDurationModelList.size ()) ||
287 (fEndOfEventModelList.size () != << 183 (fExtent != s.fExtent) ||
288 scene.fEndOfEventModelList.size ()) << 184 !(fStandardTargetPoint == s.fStandardTargetPoint) ||
289 (fEndOfRunModelList.size () != << 185 fRefreshAtEndOfEvent != s.fRefreshAtEndOfEvent
290 scene.fEndOfRunModelList.size ()) <<
291 (fExtent != scene.fExtent) <<
292 !(fStandardTargetPoint == scene.fStandar <<
293 fRefreshAtEndOfEvent != scene.fRefresh <<
294 fRefreshAtEndOfRun != scene.fRefresh <<
295 fMaxNumberOfKeptEvents != scene.fMaxNumb <<
296 ) return true; 186 ) return true;
297 187
298 /* A complete comparison should, perhaps, in <<
299 individual models, but it is not easy to <<
300 all models. Also, it would be unfeasible <<
301 implement opeerator!= if we ever get roun <<
302 user-defined models. Moreover, there is <<
303 objects, apart from changing fRefreshAtEn <<
304 models are concerned, all you can ever do <<
305 on size (above) is enough. <<
306 <<
307 for (size_t i = 0; i < fRunDurationModelList 188 for (size_t i = 0; i < fRunDurationModelList.size (); i++) {
308 if (fRunDurationModelList[i] != scene.fRun << 189 if (fRunDurationModelList[i] != s.fRunDurationModelList[i])
309 return true; <<
310 } <<
311 <<
312 for (size_t i = 0; i < fEndOfEventModelList. <<
313 if (fEndOfEventModelList[i] != scene.fEndO <<
314 return true; <<
315 } <<
316 <<
317 for (size_t i = 0; i < fEndOfRunModelList.si <<
318 if (fEndOfRunModelList[i] != scene.fEndOfR <<
319 return true; 190 return true;
320 } 191 }
321 */ <<
322 192
323 return false; 193 return false;
324 } 194 }
325 195