| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
>> 26 // $Id: G4TrajectoryDrawByAttribute.cc 78838 2014-01-28 08:46:17Z gcosmo $
26 // 27 //
27 // Jane Tinslay August 2006 28 // Jane Tinslay August 2006
28 // 29 //
29 #include "G4TrajectoryDrawByAttribute.hh" 30 #include "G4TrajectoryDrawByAttribute.hh"
30 #include "G4AttDef.hh" 31 #include "G4AttDef.hh"
31 #include "G4AttFilterUtils.hh" 32 #include "G4AttFilterUtils.hh"
32 #include "G4AttUtils.hh" 33 #include "G4AttUtils.hh"
33 #include "G4AttValue.hh" 34 #include "G4AttValue.hh"
34 #include "G4TrajectoryDrawerUtils.hh" 35 #include "G4TrajectoryDrawerUtils.hh"
35 #include "G4VAttValueFilter.hh" 36 #include "G4VAttValueFilter.hh"
36 #include "G4VisTrajContext.hh" 37 #include "G4VisTrajContext.hh"
37 #include "G4VTrajectory.hh" 38 #include "G4VTrajectory.hh"
38 #include <assert.h> 39 #include <assert.h>
39 40
40 G4TrajectoryDrawByAttribute::G4TrajectoryDrawB 41 G4TrajectoryDrawByAttribute::G4TrajectoryDrawByAttribute(const G4String& name, G4VisTrajContext* context)
41 :G4VTrajectoryModel(name, context) 42 :G4VTrajectoryModel(name, context)
42 ,fAttName("") 43 ,fAttName("")
43 ,fFirst(true) 44 ,fFirst(true)
44 ,fWarnedMissingAttribute(false) 45 ,fWarnedMissingAttribute(false)
45 ,filter(0) 46 ,filter(0)
46 {} 47 {}
47 48
48 G4TrajectoryDrawByAttribute::~G4TrajectoryDraw 49 G4TrajectoryDrawByAttribute::~G4TrajectoryDrawByAttribute()
49 { 50 {
50 ContextMap::iterator iter = fContextMap.begi 51 ContextMap::iterator iter = fContextMap.begin();
51 52
52 while (iter != fContextMap.end()) { 53 while (iter != fContextMap.end()) {
53 delete iter->second; 54 delete iter->second;
54 iter++; 55 iter++;
55 } 56 }
56 57
57 delete filter; 58 delete filter;
58 } 59 }
59 60
60 void 61 void
61 G4TrajectoryDrawByAttribute::Draw(const G4VTra 62 G4TrajectoryDrawByAttribute::Draw(const G4VTrajectory& object,
62 const G4bool& /*visible*/) const << 63 const G4bool& visible) const
63 { 64 {
64 // Return if attribute name has not been set 65 // Return if attribute name has not been set. Just print one warning
65 if (fAttName.empty()) { << 66 if (fAttName.isNull()) {
66 67
67 if (!fWarnedMissingAttribute) { 68 if (!fWarnedMissingAttribute) {
68 G4ExceptionDescription ed; 69 G4ExceptionDescription ed;
69 ed<<"Null attribute name"; 70 ed<<"Null attribute name";
70 G4Exception("G4TrajectoryDrawByAttribute 71 G4Exception("G4TrajectoryDrawByAttribute::Draw",
71 "modeling0116", 72 "modeling0116",
72 JustWarning, ed); 73 JustWarning, ed);
73 fWarnedMissingAttribute = true; 74 fWarnedMissingAttribute = true;
74 } 75 }
75 76
76 return; 77 return;
77 } 78 }
78 79
79 // Basically cache data loaded filter for ef 80 // Basically cache data loaded filter for efficiency
80 if (fFirst) { 81 if (fFirst) {
81 82
82 fFirst = false; 83 fFirst = false;
83 84
84 // Get attribute definition 85 // Get attribute definition
85 G4AttDef attDef; 86 G4AttDef attDef;
86 87
87 // Expect definition to exist 88 // Expect definition to exist
88 if (!G4AttUtils::ExtractAttDef(object, fAt 89 if (!G4AttUtils::ExtractAttDef(object, fAttName, attDef)) {
89 static G4bool warnedUnableToExtract = fa 90 static G4bool warnedUnableToExtract = false;
90 if (!warnedUnableToExtract) { 91 if (!warnedUnableToExtract) {
91 G4ExceptionDescription ed; 92 G4ExceptionDescription ed;
92 ed <<"Unable to extract attribute definition << 93 ed <<"Unable to extract attribute definition named "<<fAttName;
93 << "Available attributes:\n" <<
94 << *object.GetAttDefs(); <<
95 G4Exception 94 G4Exception
96 ("G4TrajectoryDrawByAttribute::Draw", 95 ("G4TrajectoryDrawByAttribute::Draw",
97 "modeling0117", JustWarning, ed, ". Inval 96 "modeling0117", JustWarning, ed, ". Invalid attribute name");
>> 97 G4cout << "Available attributes:\n"
>> 98 << *object.GetAttDefs();
98 warnedUnableToExtract = true; 99 warnedUnableToExtract = true;
99 } 100 }
100 return; 101 return;
101 } 102 }
102 103
103 // Get new G4AttValue filter 104 // Get new G4AttValue filter
104 filter = G4AttFilterUtils::GetNewFilter(at 105 filter = G4AttFilterUtils::GetNewFilter(attDef);
105 assert (0 != filter); 106 assert (0 != filter);
106 107
107 // Load both interval and single valued da 108 // Load both interval and single valued data. Single valued data should
108 // override interval data. 109 // override interval data.
109 ContextMap::const_iterator iter = fContext 110 ContextMap::const_iterator iter = fContextMap.begin();
110 111
111 while (iter != fContextMap.end()) { 112 while (iter != fContextMap.end()) {
112 if (iter->first.second == G4TrajectoryDr 113 if (iter->first.second == G4TrajectoryDrawByAttribute::Interval) {
113 filter->LoadIntervalElement(iter->first.firs 114 filter->LoadIntervalElement(iter->first.first);
114 } 115 }
115 else if (iter->first.second == G4Traject 116 else if (iter->first.second == G4TrajectoryDrawByAttribute::SingleValue) {
116 filter->LoadSingleValueElement(iter->first.f 117 filter->LoadSingleValueElement(iter->first.first);
117 } 118 }
118 iter++; 119 iter++;
119 } 120 }
120 } 121 }
121 122
122 // Get attribute value 123 // Get attribute value
123 G4AttValue attVal; 124 G4AttValue attVal;
124 125
125 // Expect value to exist 126 // Expect value to exist
126 if (!G4AttUtils::ExtractAttValue(object, fAt 127 if (!G4AttUtils::ExtractAttValue(object, fAttName, attVal)) {
127 static G4bool warnedUnableToExtract = fals 128 static G4bool warnedUnableToExtract = false;
128 if (!warnedUnableToExtract) { 129 if (!warnedUnableToExtract) {
129 G4ExceptionDescription ed; 130 G4ExceptionDescription ed;
130 ed <<"Unable to extract attribute defini << 131 ed <<"Unable to extract attribute value named "<<fAttName;
131 << "Available attributes:\n" <<
132 << *object.GetAttDefs(); <<
133 G4Exception 132 G4Exception
134 ("G4TrajectoryDrawByAttribute::Draw", 133 ("G4TrajectoryDrawByAttribute::Draw",
135 "modeling0118", JustWarning, ed, ". Invalid 134 "modeling0118", JustWarning, ed, ". Invalid attribute name");
>> 135 G4cout << "Available attributes:\n"
>> 136 << *object.GetAttDefs();
136 warnedUnableToExtract = true; 137 warnedUnableToExtract = true;
137 } 138 }
138 return; 139 return;
139 } 140 }
140 141
141 G4VisTrajContext myContext(GetContext()); 142 G4VisTrajContext myContext(GetContext());
142 G4String key; 143 G4String key;
143 144
144 // If attribute value passes filter, get cor 145 // If attribute value passes filter, get corresponding interval/single value
145 // key loaded into G4AttValue filter. 146 // key loaded into G4AttValue filter.
146 if (filter->GetValidElement(attVal, key)) { 147 if (filter->GetValidElement(attVal, key)) {
147 148
148 // Extract context corresponding to valid 149 // Extract context corresponding to valid key.
149 // Single value match should have override 150 // Single value match should have overriden interval match.
150 ContextMap::const_iterator iter = fContext 151 ContextMap::const_iterator iter = fContextMap.begin();
151 152
152 G4bool gotContext(false); 153 G4bool gotContext(false);
153 154
154 while (!gotContext && (iter != fContextMap 155 while (!gotContext && (iter != fContextMap.end())) {
155 if (iter->first.first == key) { 156 if (iter->first.first == key) {
156 myContext = *(iter->second); 157 myContext = *(iter->second);
157 gotContext = true; 158 gotContext = true;
158 } 159 }
159 iter++; 160 iter++;
160 } 161 }
161 162
162 assert (gotContext); 163 assert (gotContext);
163 } 164 }
>> 165
>> 166 myContext.SetVisible(visible);
164 167
165 if (GetVerbose()) { 168 if (GetVerbose()) {
166 G4cout<<"G4TrajectoryDrawByAttribute drawe 169 G4cout<<"G4TrajectoryDrawByAttribute drawer named "<<Name();
167 G4cout<<", drawing style selected accordin 170 G4cout<<", drawing style selected according to value of attribute "<<fAttName;
168 G4cout<<" : "<<attVal.GetValue()<<". Sele 171 G4cout<<" : "<<attVal.GetValue()<<". Selected context:"<<G4endl;
169 myContext.Print(G4cout); 172 myContext.Print(G4cout);
170 } 173 }
171 174
172 // Draw the trajectory 175 // Draw the trajectory
173 G4TrajectoryDrawerUtils::DrawLineAndPoints(o 176 G4TrajectoryDrawerUtils::DrawLineAndPoints(object, myContext);
174 } 177 }
175 178
176 void 179 void
177 G4TrajectoryDrawByAttribute::Print(std::ostrea 180 G4TrajectoryDrawByAttribute::Print(std::ostream& ostr) const
178 { 181 {
179 ostr<<"G4TrajectoryDrawByAttribute, dumping 182 ostr<<"G4TrajectoryDrawByAttribute, dumping configuration for model named "<< Name() <<":"<<std::endl;;
180 183
181 ostr<<"Default configuration:"<<G4endl; 184 ostr<<"Default configuration:"<<G4endl;
182 GetContext().Print(ostr); 185 GetContext().Print(ostr);
183 186
184 ostr<<"\nAttribute name "<<fAttName<<std::en 187 ostr<<"\nAttribute name "<<fAttName<<std::endl;
185 ostr<<"\nKey<->Context map dump:"<<std::endl 188 ostr<<"\nKey<->Context map dump:"<<std::endl;
186 189
187 ContextMap::const_iterator iter = fContextMa 190 ContextMap::const_iterator iter = fContextMap.begin();
188 191
189 while (iter != fContextMap.end()) { 192 while (iter != fContextMap.end()) {
190 ostr<<"Context for key "<<iter->first.firs 193 ostr<<"Context for key "<<iter->first.first<<":"<<std::endl;
191 iter->second->Print(ostr); 194 iter->second->Print(ostr);
192 195
193 iter++; 196 iter++;
194 } 197 }
195 } 198 }
196 199
197 void 200 void
198 G4TrajectoryDrawByAttribute::Set(const G4Strin 201 G4TrajectoryDrawByAttribute::Set(const G4String& name)
199 { 202 {
200 fAttName = name; 203 fAttName = name;
201 } 204 }
202 205
203 void 206 void
204 G4TrajectoryDrawByAttribute::AddIntervalContex 207 G4TrajectoryDrawByAttribute::AddIntervalContext(const G4String& name, G4VisTrajContext* context)
205 { 208 {
206 // Takes ownership of context 209 // Takes ownership of context
207 std::pair<G4String, Config> myPair(name, G4T 210 std::pair<G4String, Config> myPair(name, G4TrajectoryDrawByAttribute::Interval);
208 211
209 ContextMap::iterator iter = fContextMap.find 212 ContextMap::iterator iter = fContextMap.find(myPair);
210 213
211 if (iter != fContextMap.end()) { 214 if (iter != fContextMap.end()) {
212 G4ExceptionDescription ed; 215 G4ExceptionDescription ed;
213 ed <<"Interval "<< name <<" already exists 216 ed <<"Interval "<< name <<" already exists";
214 G4Exception 217 G4Exception
215 ("G4TrajectoryDrawByAttribute::AddInterv 218 ("G4TrajectoryDrawByAttribute::AddIntervalContext",
216 "modeling0119", FatalErrorInArgument, e 219 "modeling0119", FatalErrorInArgument, ed, ". Invalid interval");
217 } 220 }
218 221
219 fContextMap[myPair] = context; 222 fContextMap[myPair] = context;
220 } 223 }
221 224
222 void 225 void
223 G4TrajectoryDrawByAttribute::AddValueContext(c 226 G4TrajectoryDrawByAttribute::AddValueContext(const G4String& name, G4VisTrajContext* context)
224 { 227 {
225 // Takes ownership of context 228 // Takes ownership of context
226 std::pair<G4String, Config> myPair(name, G4T 229 std::pair<G4String, Config> myPair(name, G4TrajectoryDrawByAttribute::SingleValue);
227 230
228 ContextMap::iterator iter = fContextMap.find 231 ContextMap::iterator iter = fContextMap.find(myPair);
229 232
230 if (iter != fContextMap.end()) { 233 if (iter != fContextMap.end()) {
231 G4ExceptionDescription ed; 234 G4ExceptionDescription ed;
232 ed <<"Single value "<< name <<" already ex 235 ed <<"Single value "<< name <<" already exists";
233 G4Exception 236 G4Exception
234 ("G4TrajectoryDrawByAttribute::AddSingle 237 ("G4TrajectoryDrawByAttribute::AddSingleValueContext",
235 "modeling0120", FatalErrorInArgument, e 238 "modeling0120", FatalErrorInArgument, ed, ". Invalid value");
236 } 239 }
237 240
238 fContextMap[myPair] = context; 241 fContextMap[myPair] = context;
239 } 242 }
240 243