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 // 26 // 27 // 27 // 28 // GEANT4 Visualization Manager - John Allison 28 // GEANT4 Visualization Manager - John Allison 02/Jan/1996. 29 // Michael Kelsey 31 Jan 2019 -- Add new comm 29 // Michael Kelsey 31 Jan 2019 -- Add new command for electric field 30 30 31 #include "G4VisManager.hh" 31 #include "G4VisManager.hh" 32 32 33 #include "G4VisCommands.hh" 33 #include "G4VisCommands.hh" 34 #include "G4VisCommandsCompound.hh" 34 #include "G4VisCommandsCompound.hh" 35 #include "G4VisCommandsGeometry.hh" 35 #include "G4VisCommandsGeometry.hh" 36 #include "G4VisCommandsGeometrySet.hh" 36 #include "G4VisCommandsGeometrySet.hh" 37 #include "G4VisCommandsMultithreading.hh" 37 #include "G4VisCommandsMultithreading.hh" 38 #include "G4VisCommandsSet.hh" 38 #include "G4VisCommandsSet.hh" 39 #include "G4VisCommandsScene.hh" 39 #include "G4VisCommandsScene.hh" 40 #include "G4VisCommandsSceneAdd.hh" 40 #include "G4VisCommandsSceneAdd.hh" 41 #include "G4VisCommandsPlotter.hh" << 42 #include "G4VisCommandsSceneHandler.hh" 41 #include "G4VisCommandsSceneHandler.hh" 43 #include "G4VisCommandsTouchable.hh" 42 #include "G4VisCommandsTouchable.hh" 44 #include "G4VisCommandsTouchableSet.hh" 43 #include "G4VisCommandsTouchableSet.hh" 45 #include "G4VisCommandsViewer.hh" 44 #include "G4VisCommandsViewer.hh" 46 #include "G4VisCommandsViewerDefault.hh" 45 #include "G4VisCommandsViewerDefault.hh" 47 #include "G4VisCommandsViewerSet.hh" 46 #include "G4VisCommandsViewerSet.hh" 48 #include "G4UImanager.hh" 47 #include "G4UImanager.hh" 49 #include "G4VisStateDependent.hh" 48 #include "G4VisStateDependent.hh" 50 #include "G4UIdirectory.hh" 49 #include "G4UIdirectory.hh" 51 #include "G4VGraphicsSystem.hh" 50 #include "G4VGraphicsSystem.hh" 52 #include "G4VSceneHandler.hh" 51 #include "G4VSceneHandler.hh" 53 #include "G4VViewer.hh" 52 #include "G4VViewer.hh" 54 #include "G4VPhysicalVolume.hh" 53 #include "G4VPhysicalVolume.hh" 55 #include "G4LogicalVolume.hh" 54 #include "G4LogicalVolume.hh" 56 #include "G4VSolid.hh" 55 #include "G4VSolid.hh" 57 #include "G4Vector3D.hh" 56 #include "G4Vector3D.hh" 58 #include "G4Point3D.hh" 57 #include "G4Point3D.hh" 59 #include "G4RotationMatrix.hh" 58 #include "G4RotationMatrix.hh" 60 #include "G4Polyline.hh" 59 #include "G4Polyline.hh" 61 #include "G4Polyhedron.hh" 60 #include "G4Polyhedron.hh" 62 #include "G4NullModel.hh" 61 #include "G4NullModel.hh" 63 #include "G4ModelingParameters.hh" 62 #include "G4ModelingParameters.hh" 64 #include "G4TransportationManager.hh" 63 #include "G4TransportationManager.hh" 65 #include "G4VisCommandModelCreate.hh" 64 #include "G4VisCommandModelCreate.hh" 66 #include "G4VisCommandsListManager.hh" 65 #include "G4VisCommandsListManager.hh" 67 #include "G4VisModelManager.hh" 66 #include "G4VisModelManager.hh" 68 #include "G4VModelFactory.hh" 67 #include "G4VModelFactory.hh" 69 #include "G4VisFilterManager.hh" 68 #include "G4VisFilterManager.hh" 70 #include "G4VTrajectoryModel.hh" 69 #include "G4VTrajectoryModel.hh" 71 #include "G4TrajectoryDrawByCharge.hh" 70 #include "G4TrajectoryDrawByCharge.hh" 72 #include "Randomize.hh" 71 #include "Randomize.hh" 73 #include "G4RunManager.hh" 72 #include "G4RunManager.hh" 74 #include "G4RunManagerFactory.hh" 73 #include "G4RunManagerFactory.hh" 75 #include "G4EventManager.hh" 74 #include "G4EventManager.hh" 76 #include "G4Run.hh" 75 #include "G4Run.hh" 77 #include "G4Event.hh" 76 #include "G4Event.hh" 78 #include <map> 77 #include <map> 79 #include <set> 78 #include <set> 80 #include <vector> 79 #include <vector> 81 #include <sstream> 80 #include <sstream> 82 81 83 #include "G4Threading.hh" << 82 #ifdef G4MULTITHREADED 84 #include "G4AutoLock.hh" << 83 # include "G4Threading.hh" 85 #include "G4GeometryWorkspace.hh" // no_geant4 << 84 # include "G4AutoLock.hh" 86 #include "G4SolidsWorkspace.hh" << 85 # include "G4GeometryWorkspace.hh" 87 #include <deque> << 86 # include "G4SolidsWorkspace.hh" 88 #include <typeinfo> << 87 # include <deque> 89 #include <chrono> << 88 # include <typeinfo> 90 #include <thread> << 89 # include <chrono> 91 #include <utility> << 90 # include <thread> 92 << 91 #endif 93 #define G4warn G4cout << 94 << 95 // Local threading related variables << 96 namespace { << 97 G4bool mtRunInProgress = false; << 98 std::deque<const G4Event*> mtVisEventQueue; << 99 G4Thread* mtVisSubThread = 0; << 100 [[maybe_unused]] G4Mutex mtVisSubThreadMutex << 101 // G4Mutex visBeginOfRunMutex = G4MUTEX_INI << 102 // G4Mutex visBeginOfEventMutex = G4MUTEX_I << 103 G4Mutex visEndOfEventMutex = G4MUTEX_INITIAL << 104 // G4Mutex visEndOfRunMutex = G4MUTEX_INITI << 105 G4bool isSubEventRunManagerType = false; << 106 G4bool isValidViewForRun = false; << 107 G4bool isFakeRun = false; << 108 } << 109 92 110 // Class statics << 111 G4VisManager* G4VisManager::fpInstance = 0; 93 G4VisManager* G4VisManager::fpInstance = 0; 112 94 113 G4VisManager::Verbosity G4VisManager::fVerbosi 95 G4VisManager::Verbosity G4VisManager::fVerbosity = G4VisManager::warnings; 114 96 115 G4VisManager::G4VisManager (const G4String& ve << 97 G4VisManager::G4VisManager (const G4String& verbosityString): 116 : fVerbose (1) << 98 fVerbose (1), 117 , fDefaultGraphicsSystemName("OGL") / << 99 fInitialised (false), 118 , fDefaultXGeometryString ("600x600-0+0") / << 100 fpGraphicsSystem (0), 119 , fDefaultGraphicsSystemBasis ("G4VisManager i << 101 fpScene (0), 120 , fDefaultXGeometryStringBasis("G4VisManager i << 102 fpSceneHandler (0), 121 , fInitialised (false) << 103 fpViewer (0), 122 , fpGraphicsSystem (0) << 104 fpStateDependent (0), 123 , fpScene (0) << 105 fEventRefreshing (false), 124 , fpSceneHandler (0) << 106 fTransientsDrawnThisRun (false), 125 , fpViewer (0) << 107 fTransientsDrawnThisEvent (false), 126 , fpStateDependent (0) << 108 fNoOfEventsDrawnThisRun (0), 127 , fEventRefreshing (false) << 109 fNKeepRequests (0), 128 , fTransientsDrawnThisRun (false) << 110 fEventKeepingSuspended (false), 129 , fTransientsDrawnThisEvent (false) << 111 fDrawEventOnlyIfToBeKept (false), 130 , fNoOfEventsDrawnThisRun (0) << 112 fpRequestedEvent (0), 131 , fNKeepForPostProcessingRequests (0) << 113 fReviewingKeptEvents (false), 132 , fNKeepTheEventRequests (0) << 114 fAbortReviewKeptEvents (false), 133 , fEventKeepingSuspended (false) << 115 fIsDrawGroup (false), 134 , fDrawEventOnlyIfToBeKept (false) << 116 fDrawGroupNestingDepth (0), 135 , fpRequestedEvent (0) << 117 fIgnoreStateChanges (false) 136 , fReviewingKeptEvents (false) << 118 #ifdef G4MULTITHREADED 137 , fAbortReviewKeptEvents (false) << 138 , fReviewingPlots (false) << 139 , fAbortReviewPlots (false) << 140 , fIsDrawGroup (false) << 141 , fDrawGroupNestingDepth (0) << 142 , fIgnoreStateChanges (false) << 143 , fMaxEventQueueSize (100) 119 , fMaxEventQueueSize (100) 144 , fWaitOnEventQueueFull (true) 120 , fWaitOnEventQueueFull (true) 145 // All other objects use default constructors. << 121 #endif >> 122 // All other objects use default constructors. 146 { 123 { 147 fpTrajDrawModelMgr = new G4VisModelManager<G 124 fpTrajDrawModelMgr = new G4VisModelManager<G4VTrajectoryModel>("/vis/modeling/trajectories"); 148 fpTrajFilterMgr = new G4VisFilterManager<G4V 125 fpTrajFilterMgr = new G4VisFilterManager<G4VTrajectory>("/vis/filtering/trajectories"); 149 fpHitFilterMgr = new G4VisFilterManager<G4VH 126 fpHitFilterMgr = new G4VisFilterManager<G4VHit>("/vis/filtering/hits"); 150 fpDigiFilterMgr = new G4VisFilterManager<G4V 127 fpDigiFilterMgr = new G4VisFilterManager<G4VDigi>("/vis/filtering/digi"); 151 128 152 VerbosityGuidanceStrings.push_back 129 VerbosityGuidanceStrings.push_back 153 ("Simple graded message scheme - digit or 130 ("Simple graded message scheme - digit or string (1st character defines):"); 154 VerbosityGuidanceStrings.push_back 131 VerbosityGuidanceStrings.push_back 155 (" 0) quiet, // Nothing is printe 132 (" 0) quiet, // Nothing is printed."); 156 VerbosityGuidanceStrings.push_back 133 VerbosityGuidanceStrings.push_back 157 (" 1) startup, // Startup and endup 134 (" 1) startup, // Startup and endup messages are printed..."); 158 VerbosityGuidanceStrings.push_back 135 VerbosityGuidanceStrings.push_back 159 (" 2) errors, // ...and errors..." 136 (" 2) errors, // ...and errors..."); 160 VerbosityGuidanceStrings.push_back 137 VerbosityGuidanceStrings.push_back 161 (" 3) warnings, // ...and warnings.. 138 (" 3) warnings, // ...and warnings..."); 162 VerbosityGuidanceStrings.push_back 139 VerbosityGuidanceStrings.push_back 163 (" 4) confirmations, // ...and confirming 140 (" 4) confirmations, // ...and confirming messages..."); 164 VerbosityGuidanceStrings.push_back 141 VerbosityGuidanceStrings.push_back 165 (" 5) parameters, // ...and parameters 142 (" 5) parameters, // ...and parameters of scenes and views..."); 166 VerbosityGuidanceStrings.push_back 143 VerbosityGuidanceStrings.push_back 167 (" 6) all // ...and everything 144 (" 6) all // ...and everything available."); 168 145 169 if (fpInstance) { 146 if (fpInstance) { 170 G4Exception 147 G4Exception 171 ("G4VisManager::G4VisManager", 148 ("G4VisManager::G4VisManager", 172 "visman0001", FatalException, 149 "visman0001", FatalException, 173 "Attempt to Construct more than one Vis 150 "Attempt to Construct more than one VisManager"); 174 } 151 } 175 152 176 fpInstance = this; 153 fpInstance = this; 177 SetConcreteInstance(this); 154 SetConcreteInstance(this); 178 155 179 fpStateDependent = new G4VisStateDependent ( 156 fpStateDependent = new G4VisStateDependent (this); 180 // No need to delete this; G4StateManager do 157 // No need to delete this; G4StateManager does this. 181 158 182 fVerbosity = GetVerbosityValue(verbosityStri 159 fVerbosity = GetVerbosityValue(verbosityString); 183 if (fVerbosity >= startup) { 160 if (fVerbosity >= startup) { 184 G4cout 161 G4cout 185 << "Visualization Manager instantiating with 162 << "Visualization Manager instantiating with verbosity \"" 186 << VerbosityString(fVerbosity) 163 << VerbosityString(fVerbosity) 187 << "\"..." << G4endl; 164 << "\"..." << G4endl; 188 } 165 } 189 166 190 // Note: The specific graphics systems must 167 // Note: The specific graphics systems must be instantiated in a 191 // higher level library to avoid circular de 168 // higher level library to avoid circular dependencies. Also, 192 // some specifically need additional externa 169 // some specifically need additional external libararies that the 193 // user must supply. Therefore we ask the u 170 // user must supply. Therefore we ask the user to implement 194 // RegisterGraphicsSystems() and RegisterMod 171 // RegisterGraphicsSystems() and RegisterModelFactories() 195 // in a subclass. We have to wait for the s 172 // in a subclass. We have to wait for the subclass to instantiate 196 // so RegisterGraphicsSystems() cannot be ca 173 // so RegisterGraphicsSystems() cannot be called from this 197 // constructor; it is called from Initialise 174 // constructor; it is called from Initialise(). So we ask the 198 // user: 175 // user: 199 // (a) to write a subclass and implement 176 // (a) to write a subclass and implement RegisterGraphicsSystems() 200 // and RegisterModelFactories(). See 177 // and RegisterModelFactories(). See 201 // visualization/include/G4VisExecutiv 178 // visualization/include/G4VisExecutive.hh/icc as an example. 202 // (b) instantiate the subclass. 179 // (b) instantiate the subclass. 203 // (c) invoke the Initialise() method of t 180 // (c) invoke the Initialise() method of the subclass. 204 // For example: 181 // For example: 205 // ... 182 // ... 206 // // Instantiate and initialise Visualiza 183 // // Instantiate and initialise Visualization Manager. 207 // G4VisManager* visManager = new G4VisExe 184 // G4VisManager* visManager = new G4VisExecutive; 208 // visManager -> SetVerboseLevel (Verbose) 185 // visManager -> SetVerboseLevel (Verbose); 209 // visManager -> Initialise (); 186 // visManager -> Initialise (); 210 // // (Don't forget to delete visManager;) 187 // // (Don't forget to delete visManager;) 211 // ... 188 // ... 212 189 213 // Make top level command directory... 190 // Make top level command directory... 214 // vis commands should *not* be broadcast to << 191 // Vis commands should *not* be broadcast to threads (2nd argument). 215 G4bool propagateToWorkers; << 192 G4UIcommand* directory = new G4UIdirectory ("/vis/",false); 216 auto directory = new G4UIdirectory ("/vis/", << 217 directory -> SetGuidance ("Visualization com 193 directory -> SetGuidance ("Visualization commands."); 218 // Request commands in name order << 219 directory -> Sort(); // Ordering propagates << 220 fDirectoryList.push_back (directory); 194 fDirectoryList.push_back (directory); 221 195 222 // Instantiate *basic* top level commands so 196 // Instantiate *basic* top level commands so that they can be used 223 // immediately after instantiation of the vi 197 // immediately after instantiation of the vis manager. Other top 224 // level and lower level commands are instan 198 // level and lower level commands are instantiated later in 225 // RegisterMessengers. 199 // RegisterMessengers. 226 G4VVisCommand::SetVisManager (this); // Set 200 G4VVisCommand::SetVisManager (this); // Sets shared pointer 227 RegisterMessenger(new G4VisCommandVerbose); 201 RegisterMessenger(new G4VisCommandVerbose); 228 RegisterMessenger(new G4VisCommandInitialize 202 RegisterMessenger(new G4VisCommandInitialize); 229 } 203 } 230 204 231 G4VisManager::~G4VisManager() << 205 G4VisManager::~G4VisManager () { 232 { << 206 fpInstance = 0; 233 G4UImanager* UImanager = G4UImanager::GetUIp << 207 size_t i; 234 UImanager->SetCoutDestination(nullptr); << 235 std::size_t i; << 236 for (i = 0; i < fSceneList.size (); ++i) { 208 for (i = 0; i < fSceneList.size (); ++i) { 237 delete fSceneList[i]; 209 delete fSceneList[i]; 238 } 210 } 239 for (i = 0; i < fAvailableSceneHandlers.size 211 for (i = 0; i < fAvailableSceneHandlers.size (); ++i) { 240 if (fAvailableSceneHandlers[i] != NULL) { 212 if (fAvailableSceneHandlers[i] != NULL) { 241 delete fAvailableSceneHandlers[i]; 213 delete fAvailableSceneHandlers[i]; 242 } 214 } 243 } 215 } 244 for (i = 0; i < fAvailableGraphicsSystems.si 216 for (i = 0; i < fAvailableGraphicsSystems.size (); ++i) { 245 if (fAvailableGraphicsSystems[i]) { 217 if (fAvailableGraphicsSystems[i]) { 246 delete fAvailableGraphicsSystems[i]; 218 delete fAvailableGraphicsSystems[i]; 247 } 219 } 248 } 220 } 249 if (fVerbosity >= startup) { 221 if (fVerbosity >= startup) { 250 G4cout << "Graphics systems deleted." << G 222 G4cout << "Graphics systems deleted." << G4endl; 251 G4cout << "Visualization Manager deleting. 223 G4cout << "Visualization Manager deleting..." << G4endl; 252 } 224 } 253 for (i = 0; i < fMessengerList.size (); ++i) 225 for (i = 0; i < fMessengerList.size (); ++i) { 254 delete fMessengerList[i]; 226 delete fMessengerList[i]; 255 } 227 } 256 for (i = 0; i < fDirectoryList.size (); ++i) 228 for (i = 0; i < fDirectoryList.size (); ++i) { 257 delete fDirectoryList[i]; 229 delete fDirectoryList[i]; 258 } 230 } 259 231 260 delete fpDigiFilterMgr; 232 delete fpDigiFilterMgr; 261 delete fpHitFilterMgr; 233 delete fpHitFilterMgr; 262 delete fpTrajFilterMgr; 234 delete fpTrajFilterMgr; 263 delete fpTrajDrawModelMgr; 235 delete fpTrajDrawModelMgr; 264 fpInstance = 0; << 265 } 236 } 266 237 267 G4VisManager* G4VisManager::GetInstance () { 238 G4VisManager* G4VisManager::GetInstance () { 268 if (!fpInstance) { 239 if (!fpInstance) { 269 G4Exception 240 G4Exception 270 ("G4VisManager::GetInstance", 241 ("G4VisManager::GetInstance", 271 "visman0002", FatalException, "VisManag 242 "visman0002", FatalException, "VisManager not yet instantiated"); 272 } 243 } 273 return fpInstance; 244 return fpInstance; 274 } 245 } 275 246 276 void G4VisManager::Initialise () { 247 void G4VisManager::Initialise () { 277 248 278 if (fInitialised && fVerbosity >= warnings) 249 if (fInitialised && fVerbosity >= warnings) { 279 G4warn << "WARNING: G4VisManager::Initiali << 250 G4cout << "WARNING: G4VisManager::Initialise: already initialised." 280 << G4endl; 251 << G4endl; 281 return; 252 return; 282 } 253 } 283 254 284 if (fVerbosity >= startup) { 255 if (fVerbosity >= startup) { 285 G4cout << "Visualization Manager initialis 256 G4cout << "Visualization Manager initialising..." << G4endl; 286 } 257 } 287 258 288 if (fVerbosity >= parameters) { 259 if (fVerbosity >= parameters) { 289 G4cout << 260 G4cout << 290 "\nYou have instantiated your own Visual 261 "\nYou have instantiated your own Visualization Manager, inheriting" 291 "\n G4VisManager and implementing Regis 262 "\n G4VisManager and implementing RegisterGraphicsSystems(), in which" 292 "\n you should, normally, instantiate d 263 "\n you should, normally, instantiate drivers which do not need" 293 "\n external packages or libraries, and 264 "\n external packages or libraries, and, optionally, drivers under" 294 "\n control of environment variables." 265 "\n control of environment variables." 295 "\n Also you should implement RegisterM 266 "\n Also you should implement RegisterModelFactories()." 296 "\n See visualization/management/includ 267 "\n See visualization/management/include/G4VisExecutive.hh/icc, for example." 297 "\n In your main() you will have someth 268 "\n In your main() you will have something like:" 298 "\n G4VisManager* visManager = new G4 269 "\n G4VisManager* visManager = new G4VisExecutive;" 299 "\n visManager -> SetVerboseLevel (Ve 270 "\n visManager -> SetVerboseLevel (Verbose);" 300 "\n visManager -> Initialize ();" 271 "\n visManager -> Initialize ();" 301 "\n (Don't forget to delete visManager; 272 "\n (Don't forget to delete visManager;)" 302 "\n" 273 "\n" 303 << G4endl; 274 << G4endl; 304 } 275 } 305 276 306 if (fVerbosity >= startup) { 277 if (fVerbosity >= startup) { 307 G4cout << "Registering graphics systems... 278 G4cout << "Registering graphics systems..." << G4endl; 308 } 279 } 309 280 310 RegisterGraphicsSystems (); 281 RegisterGraphicsSystems (); 311 282 312 if (fVerbosity >= startup) { 283 if (fVerbosity >= startup) { 313 G4cout << 284 G4cout << 314 "\nYou have successfully registered the 285 "\nYou have successfully registered the following graphics systems." 315 << G4endl; 286 << G4endl; 316 PrintAvailableGraphicsSystems (fVerbosity) 287 PrintAvailableGraphicsSystems (fVerbosity); 317 G4cout << G4endl; 288 G4cout << G4endl; 318 } 289 } 319 290 320 // Make command directories for commands ins 291 // Make command directories for commands instantiated in the 321 // modeling subcategory... 292 // modeling subcategory... 322 G4UIcommand* directory; 293 G4UIcommand* directory; 323 directory = new G4UIdirectory ("/vis/modelin 294 directory = new G4UIdirectory ("/vis/modeling/"); 324 directory -> SetGuidance ("Modeling commands 295 directory -> SetGuidance ("Modeling commands."); 325 fDirectoryList.push_back (directory); 296 fDirectoryList.push_back (directory); 326 directory = new G4UIdirectory ("/vis/modelin 297 directory = new G4UIdirectory ("/vis/modeling/trajectories/"); 327 directory -> SetGuidance ("Trajectory model 298 directory -> SetGuidance ("Trajectory model commands."); 328 fDirectoryList.push_back (directory); 299 fDirectoryList.push_back (directory); 329 directory = new G4UIdirectory ("/vis/modelin 300 directory = new G4UIdirectory ("/vis/modeling/trajectories/create/"); 330 directory -> SetGuidance ("Create trajectory 301 directory -> SetGuidance ("Create trajectory models and messengers."); 331 fDirectoryList.push_back (directory); 302 fDirectoryList.push_back (directory); 332 303 333 // Filtering command directory 304 // Filtering command directory 334 directory = new G4UIdirectory ("/vis/filteri 305 directory = new G4UIdirectory ("/vis/filtering/"); 335 directory -> SetGuidance ("Filtering command 306 directory -> SetGuidance ("Filtering commands."); 336 fDirectoryList.push_back (directory); 307 fDirectoryList.push_back (directory); 337 directory = new G4UIdirectory ("/vis/filteri 308 directory = new G4UIdirectory ("/vis/filtering/trajectories/"); 338 directory -> SetGuidance ("Trajectory filter 309 directory -> SetGuidance ("Trajectory filtering commands."); 339 fDirectoryList.push_back (directory); 310 fDirectoryList.push_back (directory); 340 directory = new G4UIdirectory ("/vis/filteri 311 directory = new G4UIdirectory ("/vis/filtering/trajectories/create/"); 341 directory -> SetGuidance ("Create trajectory 312 directory -> SetGuidance ("Create trajectory filters and messengers."); 342 fDirectoryList.push_back (directory); 313 fDirectoryList.push_back (directory); 343 directory = new G4UIdirectory ("/vis/filteri 314 directory = new G4UIdirectory ("/vis/filtering/hits/"); 344 directory -> SetGuidance ("Hit filtering com 315 directory -> SetGuidance ("Hit filtering commands."); 345 fDirectoryList.push_back (directory); 316 fDirectoryList.push_back (directory); 346 directory = new G4UIdirectory ("/vis/filteri 317 directory = new G4UIdirectory ("/vis/filtering/hits/create/"); 347 directory -> SetGuidance ("Create hit filter 318 directory -> SetGuidance ("Create hit filters and messengers."); 348 fDirectoryList.push_back (directory); 319 fDirectoryList.push_back (directory); 349 directory = new G4UIdirectory ("/vis/filteri 320 directory = new G4UIdirectory ("/vis/filtering/digi/"); 350 directory -> SetGuidance ("Digi filtering co 321 directory -> SetGuidance ("Digi filtering commands."); 351 fDirectoryList.push_back (directory); 322 fDirectoryList.push_back (directory); 352 directory = new G4UIdirectory ("/vis/filteri 323 directory = new G4UIdirectory ("/vis/filtering/digi/create/"); 353 directory -> SetGuidance ("Create digi filte 324 directory -> SetGuidance ("Create digi filters and messengers."); 354 fDirectoryList.push_back (directory); 325 fDirectoryList.push_back (directory); 355 326 356 RegisterMessengers (); 327 RegisterMessengers (); 357 328 358 if (fVerbosity >= startup) { 329 if (fVerbosity >= startup) { 359 G4cout << "Registering model factories..." 330 G4cout << "Registering model factories..." << G4endl; 360 } 331 } 361 332 362 RegisterModelFactories(); 333 RegisterModelFactories(); 363 334 364 if (fVerbosity >= startup) { 335 if (fVerbosity >= startup) { 365 G4cout << 336 G4cout << 366 "\nYou have successfully registered the 337 "\nYou have successfully registered the following model factories." 367 << G4endl; 338 << G4endl; 368 PrintAvailableModels (fVerbosity); 339 PrintAvailableModels (fVerbosity); 369 G4cout << G4endl; 340 G4cout << G4endl; 370 } 341 } 371 342 372 if (fVerbosity >= startup) { 343 if (fVerbosity >= startup) { 373 PrintAvailableUserVisActions (fVerbosity); 344 PrintAvailableUserVisActions (fVerbosity); 374 G4cout << G4endl; 345 G4cout << G4endl; 375 } 346 } 376 347 377 InitialiseG4ColourMap(); 348 InitialiseG4ColourMap(); 378 349 379 if (fVerbosity >= startup) { 350 if (fVerbosity >= startup) { 380 G4cout << 351 G4cout << 381 "Some /vis commands (optionally) take a st 352 "Some /vis commands (optionally) take a string to specify colour." 382 "\n\"/vis/list\" to see available colours. 353 "\n\"/vis/list\" to see available colours." 383 << G4endl; 354 << G4endl; 384 } 355 } 385 356 386 fInitialised = true; 357 fInitialised = true; 387 } 358 } 388 359 389 void G4VisManager::InitialiseG4ColourMap() con 360 void G4VisManager::InitialiseG4ColourMap() const 390 { 361 { 391 G4Colour::InitialiseColourMap(); // Initial 362 G4Colour::InitialiseColourMap(); // Initialises (if not already initialised) 392 363 393 // our forever 65 named colors taken long ti 364 // our forever 65 named colors taken long time ago from X11. 394 // Extracted from g4tools/include/tools/colo 365 // Extracted from g4tools/include/tools/colors 395 // Copyright (C) 2010, Guy Barrand. All righ 366 // Copyright (C) 2010, Guy Barrand. All rights reserved. 396 // See the file tools.license for terms. 367 // See the file tools.license for terms. 397 368 398 #define TOOLS_COLORS_STAT(name,r,g,b) \ 369 #define TOOLS_COLORS_STAT(name,r,g,b) \ 399 G4Colour::AddToMap(#name, G4Colour(r,g,b)); 370 G4Colour::AddToMap(#name, G4Colour(r,g,b)); 400 371 401 //0-9 372 //0-9 402 TOOLS_COLORS_STAT(aquamarine,0.496101F,0.996 373 TOOLS_COLORS_STAT(aquamarine,0.496101F,0.996109F,0.828138F) 403 TOOLS_COLORS_STAT(mediumaquamarine,0.398444F 374 TOOLS_COLORS_STAT(mediumaquamarine,0.398444F,0.800793F,0.664073F) 404 // TOOLS_COLORS_STAT(black,0,0,0) (already 375 // TOOLS_COLORS_STAT(black,0,0,0) (already defined) 405 // TOOLS_COLORS_STAT(blue,0,0,1) (already d 376 // TOOLS_COLORS_STAT(blue,0,0,1) (already defined) 406 TOOLS_COLORS_STAT(cadetblue,0.371099F,0.6171 377 TOOLS_COLORS_STAT(cadetblue,0.371099F,0.617197F,0.62501F) 407 TOOLS_COLORS_STAT(cornflowerblue,0.390631F,0 378 TOOLS_COLORS_STAT(cornflowerblue,0.390631F,0.58204F,0.925795F) 408 TOOLS_COLORS_STAT(darkslateblue,0.281254F,0. 379 TOOLS_COLORS_STAT(darkslateblue,0.281254F,0.238285F,0.542977F) 409 TOOLS_COLORS_STAT(lightblue,0.675792F,0.8437 380 TOOLS_COLORS_STAT(lightblue,0.675792F,0.843763F,0.898451F) 410 TOOLS_COLORS_STAT(lightsteelblue,0.68751F,0. 381 TOOLS_COLORS_STAT(lightsteelblue,0.68751F,0.765637F,0.867201F) 411 TOOLS_COLORS_STAT(mediumblue,0,0,0.800793F) 382 TOOLS_COLORS_STAT(mediumblue,0,0,0.800793F) 412 383 413 //10-19 384 //10-19 414 TOOLS_COLORS_STAT(mediumslateblue,0.480476F, 385 TOOLS_COLORS_STAT(mediumslateblue,0.480476F,0.406256F,0.929702F) 415 TOOLS_COLORS_STAT(midnightblue,0.0976577F,0. 386 TOOLS_COLORS_STAT(midnightblue,0.0976577F,0.0976577F,0.437507F) 416 TOOLS_COLORS_STAT(navyblue,0,0,0.500008F) 387 TOOLS_COLORS_STAT(navyblue,0,0,0.500008F) 417 TOOLS_COLORS_STAT(navy,0,0,0.500008F) 388 TOOLS_COLORS_STAT(navy,0,0,0.500008F) 418 TOOLS_COLORS_STAT(skyblue,0.527352F,0.8047F, 389 TOOLS_COLORS_STAT(skyblue,0.527352F,0.8047F,0.917983F) 419 TOOLS_COLORS_STAT(slateblue,0.414069F,0.3515 390 TOOLS_COLORS_STAT(slateblue,0.414069F,0.351568F,0.800793F) 420 TOOLS_COLORS_STAT(steelblue,0.273442F,0.5078 391 TOOLS_COLORS_STAT(steelblue,0.273442F,0.50782F,0.703136F) 421 TOOLS_COLORS_STAT(coral,0.996109F,0.496101F, 392 TOOLS_COLORS_STAT(coral,0.996109F,0.496101F,0.312505F) 422 // TOOLS_COLORS_STAT(cyan,0,1,1) (already d 393 // TOOLS_COLORS_STAT(cyan,0,1,1) (already defined) 423 TOOLS_COLORS_STAT(firebrick,0.695323F,0.1328 394 TOOLS_COLORS_STAT(firebrick,0.695323F,0.132815F,0.132815F) 424 395 425 //20-29 396 //20-29 426 // TOOLS_COLORS_STAT(brown,0.644541F,0.1640 397 // TOOLS_COLORS_STAT(brown,0.644541F,0.164065F,0.164065F) (already defined) 427 TOOLS_COLORS_STAT(gold,0.996109F,0.839857F,0 398 TOOLS_COLORS_STAT(gold,0.996109F,0.839857F,0) 428 TOOLS_COLORS_STAT(goldenrod,0.851575F,0.6445 399 TOOLS_COLORS_STAT(goldenrod,0.851575F,0.644541F,0.125002F) 429 // TOOLS_COLORS_STAT(green,0,1,0) (already 400 // TOOLS_COLORS_STAT(green,0,1,0) (already defined) 430 TOOLS_COLORS_STAT(darkgreen,0,0.390631F,0) 401 TOOLS_COLORS_STAT(darkgreen,0,0.390631F,0) 431 TOOLS_COLORS_STAT(darkolivegreen,0.332036F,0 402 TOOLS_COLORS_STAT(darkolivegreen,0.332036F,0.417975F,0.183597F) 432 TOOLS_COLORS_STAT(forestgreen,0.132815F,0.54 403 TOOLS_COLORS_STAT(forestgreen,0.132815F,0.542977F,0.132815F) 433 TOOLS_COLORS_STAT(limegreen,0.195315F,0.8007 404 TOOLS_COLORS_STAT(limegreen,0.195315F,0.800793F,0.195315F) 434 TOOLS_COLORS_STAT(mediumseagreen,0.234379F,0 405 TOOLS_COLORS_STAT(mediumseagreen,0.234379F,0.699229F,0.441413F) 435 TOOLS_COLORS_STAT(mediumspringgreen,0,0.9765 406 TOOLS_COLORS_STAT(mediumspringgreen,0,0.976577F,0.601572F) 436 407 437 //30-39 408 //30-39 438 TOOLS_COLORS_STAT(palegreen,0.593759F,0.9804 409 TOOLS_COLORS_STAT(palegreen,0.593759F,0.980484F,0.593759F) 439 TOOLS_COLORS_STAT(seagreen,0.17969F,0.542977 410 TOOLS_COLORS_STAT(seagreen,0.17969F,0.542977F,0.339849F) 440 TOOLS_COLORS_STAT(springgreen,0,0.996109F,0. 411 TOOLS_COLORS_STAT(springgreen,0,0.996109F,0.496101F) 441 TOOLS_COLORS_STAT(yellowgreen,0.601572F,0.80 412 TOOLS_COLORS_STAT(yellowgreen,0.601572F,0.800793F,0.195315F) 442 TOOLS_COLORS_STAT(darkslategrey,0.183597F,0. 413 TOOLS_COLORS_STAT(darkslategrey,0.183597F,0.308598F,0.308598F) 443 TOOLS_COLORS_STAT(dimgrey,0.410163F,0.410163 414 TOOLS_COLORS_STAT(dimgrey,0.410163F,0.410163F,0.410163F) 444 TOOLS_COLORS_STAT(lightgrey,0.824231F,0.8242 415 TOOLS_COLORS_STAT(lightgrey,0.824231F,0.824231F,0.824231F) 445 // TOOLS_COLORS_STAT(grey,0.750011F,0.75001 416 // TOOLS_COLORS_STAT(grey,0.750011F,0.750011F,0.750011F) (already defined) 446 TOOLS_COLORS_STAT(khaki,0.937514F,0.898451F, 417 TOOLS_COLORS_STAT(khaki,0.937514F,0.898451F,0.546883F) 447 // TOOLS_COLORS_STAT(magenta,1,0,1) (alread 418 // TOOLS_COLORS_STAT(magenta,1,0,1) (already defined) 448 419 449 //40-49 420 //40-49 450 TOOLS_COLORS_STAT(maroon,0.68751F,0.187503F, 421 TOOLS_COLORS_STAT(maroon,0.68751F,0.187503F,0.375006F) 451 TOOLS_COLORS_STAT(orange,0.996109F,0.644541F 422 TOOLS_COLORS_STAT(orange,0.996109F,0.644541F,0) 452 TOOLS_COLORS_STAT(orchid,0.851575F,0.437507F 423 TOOLS_COLORS_STAT(orchid,0.851575F,0.437507F,0.83595F) 453 TOOLS_COLORS_STAT(darkorchid,0.597665F,0.195 424 TOOLS_COLORS_STAT(darkorchid,0.597665F,0.195315F,0.796887F) 454 TOOLS_COLORS_STAT(mediumorchid,0.726574F,0.3 425 TOOLS_COLORS_STAT(mediumorchid,0.726574F,0.332036F,0.824231F) 455 TOOLS_COLORS_STAT(pink,0.996109F,0.750011F,0 426 TOOLS_COLORS_STAT(pink,0.996109F,0.750011F,0.792981F) 456 TOOLS_COLORS_STAT(plum,0.863294F,0.62501F,0. 427 TOOLS_COLORS_STAT(plum,0.863294F,0.62501F,0.863294F) 457 // TOOLS_COLORS_STAT(red,1,0,0) (already de 428 // TOOLS_COLORS_STAT(red,1,0,0) (already defined) 458 TOOLS_COLORS_STAT(indianred,0.800793F,0.3593 429 TOOLS_COLORS_STAT(indianred,0.800793F,0.35938F,0.35938F) 459 TOOLS_COLORS_STAT(mediumvioletred,0.777356F, 430 TOOLS_COLORS_STAT(mediumvioletred,0.777356F,0.0820325F,0.519539F) 460 431 461 //50-59 432 //50-59 462 TOOLS_COLORS_STAT(orangered,0.996109F,0.2695 433 TOOLS_COLORS_STAT(orangered,0.996109F,0.269535F,0) 463 TOOLS_COLORS_STAT(violetred,0.812512F,0.1250 434 TOOLS_COLORS_STAT(violetred,0.812512F,0.125002F,0.562509F) 464 TOOLS_COLORS_STAT(salmon,0.976577F,0.500008F 435 TOOLS_COLORS_STAT(salmon,0.976577F,0.500008F,0.445319F) 465 TOOLS_COLORS_STAT(sienna,0.62501F,0.320317F, 436 TOOLS_COLORS_STAT(sienna,0.62501F,0.320317F,0.175784F) 466 TOOLS_COLORS_STAT(tan,0.820325F,0.703136F,0. 437 TOOLS_COLORS_STAT(tan,0.820325F,0.703136F,0.546883F) 467 TOOLS_COLORS_STAT(thistle,0.843763F,0.746105 438 TOOLS_COLORS_STAT(thistle,0.843763F,0.746105F,0.843763F) 468 TOOLS_COLORS_STAT(turquoise,0.250004F,0.8750 439 TOOLS_COLORS_STAT(turquoise,0.250004F,0.875013F,0.812512F) 469 TOOLS_COLORS_STAT(darkturquoise,0,0.8047F,0. 440 TOOLS_COLORS_STAT(darkturquoise,0,0.8047F,0.816419F) 470 TOOLS_COLORS_STAT(mediumturquoise,0.281254F, 441 TOOLS_COLORS_STAT(mediumturquoise,0.281254F,0.816419F,0.796887F) 471 TOOLS_COLORS_STAT(violet,0.929702F,0.50782F, 442 TOOLS_COLORS_STAT(violet,0.929702F,0.50782F,0.929702F) 472 443 473 //60-64 444 //60-64 474 TOOLS_COLORS_STAT(blueviolet,0.539071F,0.167 445 TOOLS_COLORS_STAT(blueviolet,0.539071F,0.167971F,0.882826F) 475 TOOLS_COLORS_STAT(wheat,0.957046F,0.867201F, 446 TOOLS_COLORS_STAT(wheat,0.957046F,0.867201F,0.699229F) 476 // TOOLS_COLORS_STAT(white,1,1,1) (already 447 // TOOLS_COLORS_STAT(white,1,1,1) (already defined) 477 // TOOLS_COLORS_STAT(yellow,1,1,0) (already 448 // TOOLS_COLORS_STAT(yellow,1,1,0) (already defined) 478 TOOLS_COLORS_STAT(greenyellow,0.675792F,0.99 449 TOOLS_COLORS_STAT(greenyellow,0.675792F,0.996109F,0.18359F) 479 450 480 #undef TOOLS_COLORS_STAT 451 #undef TOOLS_COLORS_STAT 481 } 452 } 482 453 483 void G4VisManager::RegisterMessengers () { 454 void G4VisManager::RegisterMessengers () { 484 455 485 // Instantiate individual messengers/command 456 // Instantiate individual messengers/commands (often - but not 486 // always - one command per messenger). 457 // always - one command per messenger). 487 458 488 G4UIcommand* directory; 459 G4UIcommand* directory; 489 460 490 directory = new G4UIdirectory ("/vis/geometr 461 directory = new G4UIdirectory ("/vis/geometry/"); 491 directory -> SetGuidance("Operations on vis 462 directory -> SetGuidance("Operations on vis attributes of Geant4 geometry."); 492 fDirectoryList.push_back (directory); 463 fDirectoryList.push_back (directory); 493 RegisterMessenger(new G4VisCommandGeometryLi 464 RegisterMessenger(new G4VisCommandGeometryList); 494 RegisterMessenger(new G4VisCommandGeometryRe 465 RegisterMessenger(new G4VisCommandGeometryRestore); 495 466 496 directory = new G4UIdirectory ("/vis/geometr 467 directory = new G4UIdirectory ("/vis/geometry/set/"); 497 directory -> SetGuidance("Set vis attributes 468 directory -> SetGuidance("Set vis attributes of Geant4 geometry."); 498 fDirectoryList.push_back (directory); 469 fDirectoryList.push_back (directory); 499 RegisterMessenger(new G4VisCommandGeometrySe 470 RegisterMessenger(new G4VisCommandGeometrySetColour); 500 RegisterMessenger(new G4VisCommandGeometrySe 471 RegisterMessenger(new G4VisCommandGeometrySetDaughtersInvisible); 501 RegisterMessenger(new G4VisCommandGeometrySe 472 RegisterMessenger(new G4VisCommandGeometrySetLineStyle); 502 RegisterMessenger(new G4VisCommandGeometrySe 473 RegisterMessenger(new G4VisCommandGeometrySetLineWidth); 503 RegisterMessenger(new G4VisCommandGeometrySe 474 RegisterMessenger(new G4VisCommandGeometrySetForceAuxEdgeVisible); 504 RegisterMessenger(new G4VisCommandGeometrySe 475 RegisterMessenger(new G4VisCommandGeometrySetForceCloud); 505 RegisterMessenger(new G4VisCommandGeometrySe 476 RegisterMessenger(new G4VisCommandGeometrySetForceLineSegmentsPerCircle); 506 RegisterMessenger(new G4VisCommandGeometrySe 477 RegisterMessenger(new G4VisCommandGeometrySetForceSolid); 507 RegisterMessenger(new G4VisCommandGeometrySe 478 RegisterMessenger(new G4VisCommandGeometrySetForceWireframe); 508 RegisterMessenger(new G4VisCommandGeometrySe 479 RegisterMessenger(new G4VisCommandGeometrySetVisibility); 509 480 >> 481 #ifdef G4MULTITHREADED 510 directory = new G4UIdirectory ("/vis/multith 482 directory = new G4UIdirectory ("/vis/multithreading/"); 511 directory -> SetGuidance("Commands unique to 483 directory -> SetGuidance("Commands unique to multithreading mode."); 512 fDirectoryList.push_back (directory); 484 fDirectoryList.push_back (directory); 513 RegisterMessenger(new G4VisCommandMultithrea 485 RegisterMessenger(new G4VisCommandMultithreadingActionOnEventQueueFull); 514 RegisterMessenger(new G4VisCommandMultithrea 486 RegisterMessenger(new G4VisCommandMultithreadingMaxEventQueueSize); >> 487 #endif 515 488 516 directory = new G4UIdirectory ("/vis/set/"); 489 directory = new G4UIdirectory ("/vis/set/"); 517 directory -> SetGuidance 490 directory -> SetGuidance 518 ("Set quantities for use in future command 491 ("Set quantities for use in future commands where appropriate."); 519 fDirectoryList.push_back (directory); 492 fDirectoryList.push_back (directory); 520 RegisterMessenger(new G4VisCommandSetArrow3D 493 RegisterMessenger(new G4VisCommandSetArrow3DLineSegmentsPerCircle); 521 RegisterMessenger(new G4VisCommandSetColour) 494 RegisterMessenger(new G4VisCommandSetColour); 522 RegisterMessenger(new G4VisCommandSetExtentF 495 RegisterMessenger(new G4VisCommandSetExtentForField); 523 RegisterMessenger(new G4VisCommandSetLineWid 496 RegisterMessenger(new G4VisCommandSetLineWidth); 524 RegisterMessenger(new G4VisCommandSetTextCol 497 RegisterMessenger(new G4VisCommandSetTextColour); 525 RegisterMessenger(new G4VisCommandSetTextLay 498 RegisterMessenger(new G4VisCommandSetTextLayout); 526 RegisterMessenger(new G4VisCommandSetTextSiz 499 RegisterMessenger(new G4VisCommandSetTextSize); 527 RegisterMessenger(new G4VisCommandSetTouchab 500 RegisterMessenger(new G4VisCommandSetTouchable); 528 RegisterMessenger(new G4VisCommandSetVolumeF 501 RegisterMessenger(new G4VisCommandSetVolumeForField); 529 502 530 directory = new G4UIdirectory ("/vis/scene/" 503 directory = new G4UIdirectory ("/vis/scene/"); 531 directory -> SetGuidance ("Operations on Gea 504 directory -> SetGuidance ("Operations on Geant4 scenes."); 532 fDirectoryList.push_back (directory); 505 fDirectoryList.push_back (directory); 533 RegisterMessenger(new G4VisCommandSceneActiv 506 RegisterMessenger(new G4VisCommandSceneActivateModel); 534 RegisterMessenger(new G4VisCommandSceneCreat 507 RegisterMessenger(new G4VisCommandSceneCreate); 535 RegisterMessenger(new G4VisCommandSceneEndOf 508 RegisterMessenger(new G4VisCommandSceneEndOfEventAction); 536 RegisterMessenger(new G4VisCommandSceneEndOf 509 RegisterMessenger(new G4VisCommandSceneEndOfRunAction); 537 RegisterMessenger(new G4VisCommandSceneList) 510 RegisterMessenger(new G4VisCommandSceneList); 538 RegisterMessenger(new G4VisCommandSceneNotif 511 RegisterMessenger(new G4VisCommandSceneNotifyHandlers); 539 RegisterMessenger(new G4VisCommandSceneRemov << 540 RegisterMessenger(new G4VisCommandSceneSelec 512 RegisterMessenger(new G4VisCommandSceneSelect); 541 RegisterMessenger(new G4VisCommandSceneShowE 513 RegisterMessenger(new G4VisCommandSceneShowExtents); 542 514 543 directory = new G4UIdirectory ("/vis/scene/a 515 directory = new G4UIdirectory ("/vis/scene/add/"); 544 directory -> SetGuidance ("Add model to curr 516 directory -> SetGuidance ("Add model to current scene."); 545 fDirectoryList.push_back (directory); 517 fDirectoryList.push_back (directory); 546 RegisterMessenger(new G4VisCommandSceneAddAr 518 RegisterMessenger(new G4VisCommandSceneAddArrow); 547 RegisterMessenger(new G4VisCommandSceneAddAr 519 RegisterMessenger(new G4VisCommandSceneAddArrow2D); 548 RegisterMessenger(new G4VisCommandSceneAddAx 520 RegisterMessenger(new G4VisCommandSceneAddAxes); 549 RegisterMessenger(new G4VisCommandSceneAddDa 521 RegisterMessenger(new G4VisCommandSceneAddDate); 550 RegisterMessenger(new G4VisCommandSceneAddDi 522 RegisterMessenger(new G4VisCommandSceneAddDigis); 551 RegisterMessenger(new G4VisCommandSceneAddEv 523 RegisterMessenger(new G4VisCommandSceneAddEventID); 552 RegisterMessenger(new G4VisCommandSceneAddEx 524 RegisterMessenger(new G4VisCommandSceneAddExtent); 553 RegisterMessenger(new G4VisCommandSceneAddEl 525 RegisterMessenger(new G4VisCommandSceneAddElectricField); 554 RegisterMessenger(new G4VisCommandSceneAddFr 526 RegisterMessenger(new G4VisCommandSceneAddFrame); 555 RegisterMessenger(new G4VisCommandSceneAddGP 527 RegisterMessenger(new G4VisCommandSceneAddGPS); 556 RegisterMessenger(new G4VisCommandSceneAddHi 528 RegisterMessenger(new G4VisCommandSceneAddHits); 557 RegisterMessenger(new G4VisCommandSceneAddLi 529 RegisterMessenger(new G4VisCommandSceneAddLine); 558 RegisterMessenger(new G4VisCommandSceneAddLi 530 RegisterMessenger(new G4VisCommandSceneAddLine2D); 559 RegisterMessenger(new G4VisCommandSceneAddLo << 560 RegisterMessenger(new G4VisCommandSceneAddLo 531 RegisterMessenger(new G4VisCommandSceneAddLogicalVolume); 561 RegisterMessenger(new G4VisCommandSceneAddLo 532 RegisterMessenger(new G4VisCommandSceneAddLogo); 562 RegisterMessenger(new G4VisCommandSceneAddLo 533 RegisterMessenger(new G4VisCommandSceneAddLogo2D); 563 RegisterMessenger(new G4VisCommandSceneAddMa 534 RegisterMessenger(new G4VisCommandSceneAddMagneticField); 564 RegisterMessenger(new G4VisCommandSceneAddPl << 565 RegisterMessenger(new G4VisCommandSceneAddPS 535 RegisterMessenger(new G4VisCommandSceneAddPSHits); 566 RegisterMessenger(new G4VisCommandSceneAddSc 536 RegisterMessenger(new G4VisCommandSceneAddScale); 567 RegisterMessenger(new G4VisCommandSceneAddTe 537 RegisterMessenger(new G4VisCommandSceneAddText); 568 RegisterMessenger(new G4VisCommandSceneAddTe 538 RegisterMessenger(new G4VisCommandSceneAddText2D); 569 RegisterMessenger(new G4VisCommandSceneAddTr 539 RegisterMessenger(new G4VisCommandSceneAddTrajectories); 570 RegisterMessenger(new G4VisCommandSceneAddUs 540 RegisterMessenger(new G4VisCommandSceneAddUserAction); 571 RegisterMessenger(new G4VisCommandSceneAddVo 541 RegisterMessenger(new G4VisCommandSceneAddVolume); 572 542 573 RegisterMessenger(new G4VisCommandPlotterCre << 574 RegisterMessenger(new G4VisCommandPlotterSet << 575 RegisterMessenger(new G4VisCommandPlotterAdd << 576 RegisterMessenger(new G4VisCommandPlotterAdd << 577 RegisterMessenger(new G4VisCommandPlotterAdd << 578 RegisterMessenger(new G4VisCommandPlotterCle << 579 RegisterMessenger(new G4VisCommandPlotterCle << 580 RegisterMessenger(new G4VisCommandPlotterLis << 581 RegisterMessenger(new G4VisCommandPlotterAdd << 582 RegisterMessenger(new G4VisCommandPlotterAdd << 583 << 584 directory = new G4UIdirectory ("/vis/sceneHa 543 directory = new G4UIdirectory ("/vis/sceneHandler/"); 585 directory -> SetGuidance ("Operations on Gea 544 directory -> SetGuidance ("Operations on Geant4 scene handlers."); 586 fDirectoryList.push_back (directory); 545 fDirectoryList.push_back (directory); 587 RegisterMessenger(new G4VisCommandSceneHandl 546 RegisterMessenger(new G4VisCommandSceneHandlerAttach); 588 RegisterMessenger(new G4VisCommandSceneHandl 547 RegisterMessenger(new G4VisCommandSceneHandlerCreate); 589 RegisterMessenger(new G4VisCommandSceneHandl 548 RegisterMessenger(new G4VisCommandSceneHandlerList); 590 RegisterMessenger(new G4VisCommandSceneHandl 549 RegisterMessenger(new G4VisCommandSceneHandlerSelect); 591 550 592 directory = new G4UIdirectory ("/vis/touchab 551 directory = new G4UIdirectory ("/vis/touchable/"); 593 directory -> SetGuidance ("Operations on tou 552 directory -> SetGuidance ("Operations on touchables."); 594 fDirectoryList.push_back (directory); 553 fDirectoryList.push_back (directory); 595 RegisterMessenger(new G4VisCommandsTouchable 554 RegisterMessenger(new G4VisCommandsTouchable); 596 555 597 directory = new G4UIdirectory ("/vis/touchab 556 directory = new G4UIdirectory ("/vis/touchable/set/"); 598 directory -> SetGuidance ("Set vis attribute 557 directory -> SetGuidance ("Set vis attributes of current touchable."); 599 fDirectoryList.push_back (directory); 558 fDirectoryList.push_back (directory); 600 RegisterMessenger(new G4VisCommandsTouchable 559 RegisterMessenger(new G4VisCommandsTouchableSet); 601 560 602 directory = new G4UIdirectory ("/vis/viewer/ 561 directory = new G4UIdirectory ("/vis/viewer/"); 603 directory -> SetGuidance ("Operations on Gea 562 directory -> SetGuidance ("Operations on Geant4 viewers."); 604 fDirectoryList.push_back (directory); 563 fDirectoryList.push_back (directory); 605 RegisterMessenger(new G4VisCommandViewerAddC 564 RegisterMessenger(new G4VisCommandViewerAddCutawayPlane); 606 RegisterMessenger(new G4VisCommandViewerCent 565 RegisterMessenger(new G4VisCommandViewerCentreOn); 607 RegisterMessenger(new G4VisCommandViewerChan 566 RegisterMessenger(new G4VisCommandViewerChangeCutawayPlane); 608 RegisterMessenger(new G4VisCommandViewerClea 567 RegisterMessenger(new G4VisCommandViewerClear); 609 RegisterMessenger(new G4VisCommandViewerClea 568 RegisterMessenger(new G4VisCommandViewerClearCutawayPlanes); 610 RegisterMessenger(new G4VisCommandViewerClea 569 RegisterMessenger(new G4VisCommandViewerClearTransients); 611 RegisterMessenger(new G4VisCommandViewerClea 570 RegisterMessenger(new G4VisCommandViewerClearVisAttributesModifiers); 612 RegisterMessenger(new G4VisCommandViewerClon 571 RegisterMessenger(new G4VisCommandViewerClone); 613 RegisterMessenger(new G4VisCommandViewerColo 572 RegisterMessenger(new G4VisCommandViewerColourByDensity); 614 RegisterMessenger(new G4VisCommandViewerCopy 573 RegisterMessenger(new G4VisCommandViewerCopyViewFrom); 615 RegisterMessenger(new G4VisCommandViewerCrea 574 RegisterMessenger(new G4VisCommandViewerCreate); 616 RegisterMessenger(new G4VisCommandViewerDoll 575 RegisterMessenger(new G4VisCommandViewerDolly); 617 RegisterMessenger(new G4VisCommandViewerFlus 576 RegisterMessenger(new G4VisCommandViewerFlush); 618 RegisterMessenger(new G4VisCommandViewerInte 577 RegisterMessenger(new G4VisCommandViewerInterpolate); 619 RegisterMessenger(new G4VisCommandViewerList 578 RegisterMessenger(new G4VisCommandViewerList); 620 RegisterMessenger(new G4VisCommandViewerPan) 579 RegisterMessenger(new G4VisCommandViewerPan); 621 RegisterMessenger(new G4VisCommandViewerRebu 580 RegisterMessenger(new G4VisCommandViewerRebuild); 622 RegisterMessenger(new G4VisCommandViewerRefr 581 RegisterMessenger(new G4VisCommandViewerRefresh); 623 RegisterMessenger(new G4VisCommandViewerRese 582 RegisterMessenger(new G4VisCommandViewerReset); 624 RegisterMessenger(new G4VisCommandViewerRese << 625 RegisterMessenger(new G4VisCommandViewerSave 583 RegisterMessenger(new G4VisCommandViewerSave); 626 RegisterMessenger(new G4VisCommandViewerScal 584 RegisterMessenger(new G4VisCommandViewerScale); 627 RegisterMessenger(new G4VisCommandViewerSele 585 RegisterMessenger(new G4VisCommandViewerSelect); 628 RegisterMessenger(new G4VisCommandViewerUpda 586 RegisterMessenger(new G4VisCommandViewerUpdate); 629 RegisterMessenger(new G4VisCommandViewerZoom 587 RegisterMessenger(new G4VisCommandViewerZoom); 630 588 631 directory = new G4UIdirectory ("/vis/viewer/ 589 directory = new G4UIdirectory ("/vis/viewer/default/"); 632 directory -> SetGuidance("Set default values 590 directory -> SetGuidance("Set default values for future viewers."); 633 fDirectoryList.push_back (directory); 591 fDirectoryList.push_back (directory); 634 RegisterMessenger(new G4VisCommandViewerDefa 592 RegisterMessenger(new G4VisCommandViewerDefaultHiddenEdge); 635 RegisterMessenger(new G4VisCommandViewerDefa 593 RegisterMessenger(new G4VisCommandViewerDefaultStyle); 636 594 637 directory = new G4UIdirectory ("/vis/viewer/ 595 directory = new G4UIdirectory ("/vis/viewer/set/"); 638 directory -> SetGuidance ("Set view paramete 596 directory -> SetGuidance ("Set view parameters of current viewer."); 639 fDirectoryList.push_back (directory); 597 fDirectoryList.push_back (directory); 640 RegisterMessenger(new G4VisCommandsViewerSet 598 RegisterMessenger(new G4VisCommandsViewerSet); 641 599 642 // *Basic* top level commands were instantia 600 // *Basic* top level commands were instantiated in the constructor 643 // so that they can be used immediately afte 601 // so that they can be used immediately after instantiation of the 644 // vis manager. Other top level commands, i 602 // vis manager. Other top level commands, including "compound commands" 645 // (i.e., commands that invoke other command 603 // (i.e., commands that invoke other commands) are instantiated here. 646 604 647 RegisterMessenger(new G4VisCommandAbortRevie 605 RegisterMessenger(new G4VisCommandAbortReviewKeptEvents); 648 RegisterMessenger(new G4VisCommandAbortRevie << 649 RegisterMessenger(new G4VisCommandDrawOnlyTo 606 RegisterMessenger(new G4VisCommandDrawOnlyToBeKeptEvents); 650 RegisterMessenger(new G4VisCommandDrawTree); 607 RegisterMessenger(new G4VisCommandDrawTree); 651 RegisterMessenger(new G4VisCommandDrawView); 608 RegisterMessenger(new G4VisCommandDrawView); 652 RegisterMessenger(new G4VisCommandDrawLogica 609 RegisterMessenger(new G4VisCommandDrawLogicalVolume); 653 RegisterMessenger(new G4VisCommandDrawVolume 610 RegisterMessenger(new G4VisCommandDrawVolume); 654 RegisterMessenger(new G4VisCommandEnable); 611 RegisterMessenger(new G4VisCommandEnable); 655 RegisterMessenger(new G4VisCommandList); 612 RegisterMessenger(new G4VisCommandList); 656 RegisterMessenger(new G4VisCommandOpen); 613 RegisterMessenger(new G4VisCommandOpen); 657 RegisterMessenger(new G4VisCommandPlot); << 658 RegisterMessenger(new G4VisCommandReviewKept 614 RegisterMessenger(new G4VisCommandReviewKeptEvents); 659 RegisterMessenger(new G4VisCommandReviewPlot << 660 RegisterMessenger(new G4VisCommandSpecify); 615 RegisterMessenger(new G4VisCommandSpecify); 661 616 662 // List manager commands 617 // List manager commands 663 RegisterMessenger(new G4VisCommandListManage 618 RegisterMessenger(new G4VisCommandListManagerList< G4VisModelManager<G4VTrajectoryModel> > 664 (fpTrajDrawModelMgr, fpTrajDrawModelMg 619 (fpTrajDrawModelMgr, fpTrajDrawModelMgr->Placement())); 665 RegisterMessenger(new G4VisCommandListManage 620 RegisterMessenger(new G4VisCommandListManagerSelect< G4VisModelManager<G4VTrajectoryModel> > 666 (fpTrajDrawModelMgr, fpTrajDrawModelMg 621 (fpTrajDrawModelMgr, fpTrajDrawModelMgr->Placement())); 667 622 668 // Trajectory filter manager commands 623 // Trajectory filter manager commands 669 RegisterMessenger(new G4VisCommandListManage 624 RegisterMessenger(new G4VisCommandListManagerList< G4VisFilterManager<G4VTrajectory> > 670 (fpTrajFilterMgr, fpTrajFi 625 (fpTrajFilterMgr, fpTrajFilterMgr->Placement())); 671 RegisterMessenger(new G4VisCommandManagerMod 626 RegisterMessenger(new G4VisCommandManagerMode< G4VisFilterManager<G4VTrajectory> > 672 (fpTrajFilterMgr, fpTrajFi 627 (fpTrajFilterMgr, fpTrajFilterMgr->Placement())); 673 628 674 // Hit filter manager commands 629 // Hit filter manager commands 675 RegisterMessenger(new G4VisCommandListManage 630 RegisterMessenger(new G4VisCommandListManagerList< G4VisFilterManager<G4VHit> > 676 (fpHitFilterMgr, fpHitFilt 631 (fpHitFilterMgr, fpHitFilterMgr->Placement())); 677 RegisterMessenger(new G4VisCommandManagerMod 632 RegisterMessenger(new G4VisCommandManagerMode< G4VisFilterManager<G4VHit> > 678 (fpHitFilterMgr, fpHitFilt 633 (fpHitFilterMgr, fpHitFilterMgr->Placement())); 679 634 680 // Digi filter manager commands 635 // Digi filter manager commands 681 RegisterMessenger(new G4VisCommandListManage 636 RegisterMessenger(new G4VisCommandListManagerList< G4VisFilterManager<G4VDigi> > 682 (fpDigiFilterMgr, fpDigiFi 637 (fpDigiFilterMgr, fpDigiFilterMgr->Placement())); 683 RegisterMessenger(new G4VisCommandManagerMod 638 RegisterMessenger(new G4VisCommandManagerMode< G4VisFilterManager<G4VDigi> > 684 (fpDigiFilterMgr, fpDigiFi 639 (fpDigiFilterMgr, fpDigiFilterMgr->Placement())); 685 } 640 } 686 641 687 #include <tools/histo/h1d> << 688 #include <tools/histo/h2d> << 689 << 690 namespace { << 691 struct PlotResults { << 692 std::size_t fNumberOfPlots = 0; << 693 std::size_t fTotalEntries = 0; << 694 }; << 695 template <typename HT> // tools::histo::h1d << 696 PlotResults ResultsOfHnPlots(const G4String& << 697 PlotResults plotResults; << 698 auto ui = G4UImanager::GetUIpointer(); << 699 auto keepControlVerbose = ui->GetVerboseLe << 700 ui->SetVerboseLevel(0); << 701 auto status = ui->ApplyCommand("/analysis/ << 702 ui->SetVerboseLevel(keepControlVerbose); << 703 if(status==G4UIcommandStatus::fCommandSucc << 704 G4String hexString = ui->GetCurrentValue << 705 if(hexString.size()) { << 706 void* ptr; std::istringstream is(hexSt << 707 auto vectorOfPlots = (const std::vecto << 708 for (std::size_t i = 0; i < vectorOfPl << 709 auto plot = (*vectorOfPlots)[i]; << 710 if (plot == nullptr) continue; // I << 711 ++plotResults.fNumberOfPlots; << 712 plotResults.fTotalEntries += plot->e << 713 } << 714 } << 715 } << 716 return plotResults; << 717 } << 718 void PrintListOfPlots() { << 719 std::size_t numberOfPlots = 0; << 720 std::size_t numberOfEntries = 0; << 721 PlotResults h1results = ResultsOfHnPlots<t << 722 numberOfPlots += h1results.fNumberOfPlots; << 723 numberOfEntries += h1results.fTotalEntries << 724 PlotResults h2results = ResultsOfHnPlots<t << 725 numberOfPlots += h2results.fNumberOfPlots; << 726 numberOfEntries += h2results.fTotalEntries << 727 if (numberOfPlots > 0) { << 728 G4warn << "There are histograms that can << 729 if (h1results.fNumberOfPlots > 0) { << 730 G4warn << "\n " << h1results.fNumberO << 731 } << 732 if (h2results.fNumberOfPlots > 0) { << 733 G4warn << "\n " << h2results.fNumberO << 734 } << 735 G4warn << 736 << "\n List them with \"/analysis/list\ << 737 << "\n View them immediately with \"/vi << 738 << G4endl; << 739 if (numberOfEntries == 0) { << 740 G4warn << << 741 " But...there are no entries. To make << 742 "\n plotting in this UI session, use << 743 "\n EndOfRunAction and Reset() in you << 744 << G4endl; << 745 } << 746 } << 747 } << 748 } << 749 << 750 void G4VisManager::Enable() { 642 void G4VisManager::Enable() { 751 if (IsValidView ()) { 643 if (IsValidView ()) { 752 SetConcreteInstance(this); 644 SetConcreteInstance(this); 753 if (fVerbosity >= confirmations) { 645 if (fVerbosity >= confirmations) { 754 G4cout << "G4VisManager::Enable: visuali 646 G4cout << "G4VisManager::Enable: visualization enabled." << G4endl; 755 } 647 } 756 if (fVerbosity >= warnings) { 648 if (fVerbosity >= warnings) { 757 std::size_t nKeptEvents = 0; << 649 G4int nKeptEvents = 0; 758 const G4Run* run = G4RunManager::GetRunM 650 const G4Run* run = G4RunManager::GetRunManager()->GetCurrentRun(); 759 if (run) nKeptEvents = run->GetNumberOfK << 651 if (run) nKeptEvents = run->GetEventVector()->size(); 760 G4String isare("are"),plural("s"); << 761 if (nKeptEvents == 1) {isare = "is"; plu << 762 G4cout << 652 G4cout << 763 "There " << isare << ' ' << nKeptEvents << 653 "There are " << nKeptEvents << " kept events." >> 654 "\n \"/vis/reviewKeptEvents\" to review them one by one." >> 655 "\n \"/vis/viewer/flush\" or \"/vis/viewer/rebuild\" to see them accumulated." 764 << G4endl; 656 << G4endl; 765 if (nKeptEvents > 0) { << 766 G4cout << << 767 " \"/vis/reviewKeptEvents\" to review one b << 768 "\n To see accumulated, \"/vis/enable\", th << 769 << G4endl; << 770 } << 771 PrintListOfPlots(); << 772 } 657 } 773 } 658 } 774 else { 659 else { 775 if (fVerbosity >= warnings) { 660 if (fVerbosity >= warnings) { 776 G4warn << << 661 G4cout << 777 "G4VisManager::Enable: WARNING: visualizatio 662 "G4VisManager::Enable: WARNING: visualization remains disabled for" 778 "\n above reasons. Rectifying with valid v 663 "\n above reasons. Rectifying with valid vis commands will" 779 "\n automatically enable." 664 "\n automatically enable." 780 << G4endl; 665 << G4endl; 781 } 666 } 782 } 667 } 783 } 668 } 784 669 785 void G4VisManager::Disable() { 670 void G4VisManager::Disable() { 786 SetConcreteInstance(0); 671 SetConcreteInstance(0); 787 if (fVerbosity >= confirmations) { 672 if (fVerbosity >= confirmations) { 788 G4cout << 673 G4cout << 789 "G4VisManager::Disable: visualization disa 674 "G4VisManager::Disable: visualization disabled." 790 "\n The pointer returned by GetConcreteIn 675 "\n The pointer returned by GetConcreteInstance will be zero." 791 "\n Note that it will become enabled afte 676 "\n Note that it will become enabled after some valid vis commands." 792 << G4endl; 677 << G4endl; 793 } 678 } 794 if (fVerbosity >= warnings) { 679 if (fVerbosity >= warnings) { 795 G4int currentTrajectoryType = 680 G4int currentTrajectoryType = 796 G4RunManagerKernel::GetRunManagerKernel()- 681 G4RunManagerKernel::GetRunManagerKernel()->GetTrackingManager()->GetStoreTrajectory(); 797 if (currentTrajectoryType > 0) { 682 if (currentTrajectoryType > 0) { 798 G4warn << << 683 G4cout << 799 "You may wish to disable trajectory prod 684 "You may wish to disable trajectory production too:" 800 "\n \"/tracking/storeTrajectory 0\"" 685 "\n \"/tracking/storeTrajectory 0\"" 801 "\nbut don't forget to re-enable with" 686 "\nbut don't forget to re-enable with" 802 "\n \"/vis/enable\"" 687 "\n \"/vis/enable\"" 803 "\n \"/tracking/storeTrajectory " << cu 688 "\n \"/tracking/storeTrajectory " << currentTrajectoryType 804 << "\"\n and maybe \"/vis/viewer/rebuil 689 << "\"\n and maybe \"/vis/viewer/rebuild\"" 805 << G4endl; 690 << G4endl; 806 } 691 } 807 } 692 } 808 } 693 } 809 694 810 const G4GraphicsSystemList& G4VisManager::GetA 695 const G4GraphicsSystemList& G4VisManager::GetAvailableGraphicsSystems () { 811 std::size_t nSystems = fAvailableGraphicsSys << 696 G4int nSystems = fAvailableGraphicsSystems.size (); 812 if (nSystems == 0) { 697 if (nSystems == 0) { 813 if (fVerbosity >= warnings) { 698 if (fVerbosity >= warnings) { 814 G4warn << "G4VisManager::GetAvailableGra << 699 G4cout << "G4VisManager::GetAvailableGraphicsSystems: WARNING: no" 815 "\n graphics system available!" 700 "\n graphics system available!" 816 "\n 1) Did you have environment variables G 701 "\n 1) Did you have environment variables G4VIS_BUILD_xxxx_DRIVER set" 817 "\n when you compiled/built the visualiz 702 "\n when you compiled/built the visualization code?" 818 "\n 2) Did you instantiate your own Visuali 703 "\n 2) Did you instantiate your own Visualization Manager and forget" 819 "\n to implement RegisterGraphicsSystems 704 "\n to implement RegisterGraphicsSystems correctly?" 820 "\n 3) You can register your own graphics s 705 "\n 3) You can register your own graphics system, e.g.," 821 "\n visManager->RegisterGraphicsSystem(n 706 "\n visManager->RegisterGraphicsSystem(new MyGraphicsSystem);)" 822 "\n after instantiating your vis manager 707 "\n after instantiating your vis manager and before" 823 "\n visManager->Initialize()." 708 "\n visManager->Initialize()." 824 << G4endl; 709 << G4endl; 825 } 710 } 826 } 711 } 827 return fAvailableGraphicsSystems; 712 return fAvailableGraphicsSystems; 828 } 713 } 829 714 830 G4bool G4VisManager::RegisterGraphicsSystem (G 715 G4bool G4VisManager::RegisterGraphicsSystem (G4VGraphicsSystem* pSystem) { 831 G4bool happy = true; 716 G4bool happy = true; 832 if (pSystem) { 717 if (pSystem) { 833 fAvailableGraphicsSystems.push_back (pSyst 718 fAvailableGraphicsSystems.push_back (pSystem); 834 if (fVerbosity >= confirmations) { 719 if (fVerbosity >= confirmations) { 835 G4cout << "G4VisManager::RegisterGraphic 720 G4cout << "G4VisManager::RegisterGraphicsSystem: " 836 << pSystem -> GetName (); 721 << pSystem -> GetName (); 837 if (pSystem -> GetNickname () != "") { 722 if (pSystem -> GetNickname () != "") { 838 G4cout << " (" << pSystem -> GetNickname () 723 G4cout << " (" << pSystem -> GetNickname () << ")"; 839 } 724 } 840 G4cout << " registered." << G4endl; 725 G4cout << " registered." << G4endl; 841 } 726 } 842 } 727 } 843 else { 728 else { 844 if (fVerbosity >= errors) { 729 if (fVerbosity >= errors) { 845 G4warn << "G4VisManager::RegisterGraphic << 730 G4cout << "G4VisManager::RegisterGraphicsSystem: null pointer!" 846 << G4endl; 731 << G4endl; 847 } 732 } 848 happy=false; 733 happy=false; 849 } 734 } 850 return happy; 735 return happy; 851 } 736 } 852 737 853 const G4VTrajectoryModel* 738 const G4VTrajectoryModel* 854 G4VisManager::CurrentTrajDrawModel() const 739 G4VisManager::CurrentTrajDrawModel() const 855 { 740 { 856 assert (0 != fpTrajDrawModelMgr); 741 assert (0 != fpTrajDrawModelMgr); 857 742 858 const G4VTrajectoryModel* model = fpTrajDraw 743 const G4VTrajectoryModel* model = fpTrajDrawModelMgr->Current(); 859 744 860 if (0 == model) { 745 if (0 == model) { 861 // No model was registered with the trajec 746 // No model was registered with the trajectory model manager. 862 // Use G4TrajectoryDrawByCharge as a fallb 747 // Use G4TrajectoryDrawByCharge as a fallback. 863 fpTrajDrawModelMgr->Register(new G4Traject 748 fpTrajDrawModelMgr->Register(new G4TrajectoryDrawByCharge("DefaultModel")); 864 if (fVerbosity >= warnings) { 749 if (fVerbosity >= warnings) { 865 G4warn<<"G4VisManager: Using G4Trajector << 750 G4cout<<"G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model."<<G4endl; 866 G4warn<<"See commands in /vis/modeling/t << 751 G4cout<<"See commands in /vis/modeling/trajectories/ for other options."<<G4endl; 867 } 752 } 868 } 753 } 869 754 870 model = fpTrajDrawModelMgr->Current(); 755 model = fpTrajDrawModelMgr->Current(); 871 assert (0 != model); // Should definitely ex 756 assert (0 != model); // Should definitely exist now 872 757 873 return model; 758 return model; 874 } 759 } 875 760 876 void G4VisManager::RegisterModel(G4VTrajectory 761 void G4VisManager::RegisterModel(G4VTrajectoryModel* model) 877 { 762 { 878 fpTrajDrawModelMgr->Register(model); 763 fpTrajDrawModelMgr->Register(model); 879 } 764 } 880 765 881 void 766 void 882 G4VisManager::RegisterModelFactory(G4TrajDrawM 767 G4VisManager::RegisterModelFactory(G4TrajDrawModelFactory* factory) 883 { 768 { 884 fpTrajDrawModelMgr->Register(factory); 769 fpTrajDrawModelMgr->Register(factory); 885 } 770 } 886 771 887 void G4VisManager::RegisterModel(G4VFilter<G4V 772 void G4VisManager::RegisterModel(G4VFilter<G4VTrajectory>* model) 888 { 773 { 889 fpTrajFilterMgr->Register(model); 774 fpTrajFilterMgr->Register(model); 890 } 775 } 891 776 892 void 777 void 893 G4VisManager::RegisterModelFactory(G4TrajFilte 778 G4VisManager::RegisterModelFactory(G4TrajFilterFactory* factory) 894 { 779 { 895 fpTrajFilterMgr->Register(factory); 780 fpTrajFilterMgr->Register(factory); 896 } 781 } 897 782 898 void G4VisManager::RegisterModel(G4VFilter<G4V 783 void G4VisManager::RegisterModel(G4VFilter<G4VHit>* model) 899 { 784 { 900 fpHitFilterMgr->Register(model); 785 fpHitFilterMgr->Register(model); 901 } 786 } 902 787 903 void 788 void 904 G4VisManager::RegisterModelFactory(G4HitFilter 789 G4VisManager::RegisterModelFactory(G4HitFilterFactory* factory) 905 { 790 { 906 fpHitFilterMgr->Register(factory); 791 fpHitFilterMgr->Register(factory); 907 } 792 } 908 793 909 void G4VisManager::RegisterModel(G4VFilter<G4V 794 void G4VisManager::RegisterModel(G4VFilter<G4VDigi>* model) 910 { 795 { 911 fpDigiFilterMgr->Register(model); 796 fpDigiFilterMgr->Register(model); 912 } 797 } 913 798 914 void 799 void 915 G4VisManager::RegisterModelFactory(G4DigiFilte 800 G4VisManager::RegisterModelFactory(G4DigiFilterFactory* factory) 916 { 801 { 917 fpDigiFilterMgr->Register(factory); 802 fpDigiFilterMgr->Register(factory); 918 } 803 } 919 804 920 void G4VisManager::SelectTrajectoryModel(const 805 void G4VisManager::SelectTrajectoryModel(const G4String& model) 921 { 806 { 922 fpTrajDrawModelMgr->SetCurrent(model); 807 fpTrajDrawModelMgr->SetCurrent(model); 923 } 808 } 924 809 925 void G4VisManager::BeginDraw (const G4Transfor 810 void G4VisManager::BeginDraw (const G4Transform3D& objectTransform) 926 { 811 { >> 812 #ifdef G4MULTITHREADED 927 if (G4Threading::IsWorkerThread()) return; 813 if (G4Threading::IsWorkerThread()) return; 928 << 814 #endif 929 fDrawGroupNestingDepth++; 815 fDrawGroupNestingDepth++; 930 if (fDrawGroupNestingDepth > 1) { 816 if (fDrawGroupNestingDepth > 1) { 931 G4Exception 817 G4Exception 932 ("G4VisManager::BeginDraw", 818 ("G4VisManager::BeginDraw", 933 "visman0008", JustWarning, 819 "visman0008", JustWarning, 934 "Nesting detected. It is illegal to nes 820 "Nesting detected. It is illegal to nest Begin/EndDraw." 935 "\n Ignored"); 821 "\n Ignored"); 936 return; 822 return; 937 } 823 } 938 if (IsValidView ()) { 824 if (IsValidView ()) { 939 ClearTransientStoreIfMarked(); 825 ClearTransientStoreIfMarked(); 940 fpSceneHandler -> BeginPrimitives (objectT 826 fpSceneHandler -> BeginPrimitives (objectTransform); 941 fIsDrawGroup = true; 827 fIsDrawGroup = true; 942 } 828 } 943 } 829 } 944 830 945 void G4VisManager::EndDraw () 831 void G4VisManager::EndDraw () 946 { 832 { >> 833 #ifdef G4MULTITHREADED 947 if (G4Threading::IsWorkerThread()) return; 834 if (G4Threading::IsWorkerThread()) return; 948 << 835 #endif 949 fDrawGroupNestingDepth--; 836 fDrawGroupNestingDepth--; 950 if (fDrawGroupNestingDepth != 0) { 837 if (fDrawGroupNestingDepth != 0) { 951 if (fDrawGroupNestingDepth < 0) fDrawGroup 838 if (fDrawGroupNestingDepth < 0) fDrawGroupNestingDepth = 0; 952 return; 839 return; 953 } 840 } 954 if (IsValidView ()) { 841 if (IsValidView ()) { 955 fpSceneHandler -> EndPrimitives (); 842 fpSceneHandler -> EndPrimitives (); 956 } 843 } 957 fIsDrawGroup = false; 844 fIsDrawGroup = false; 958 } 845 } 959 846 960 void G4VisManager::BeginDraw2D (const G4Transf 847 void G4VisManager::BeginDraw2D (const G4Transform3D& objectTransform) 961 { 848 { >> 849 #ifdef G4MULTITHREADED 962 if (G4Threading::IsWorkerThread()) return; 850 if (G4Threading::IsWorkerThread()) return; 963 << 851 #endif 964 fDrawGroupNestingDepth++; 852 fDrawGroupNestingDepth++; 965 if (fDrawGroupNestingDepth > 1) { 853 if (fDrawGroupNestingDepth > 1) { 966 G4Exception 854 G4Exception 967 ("G4VisManager::BeginDraw2D", 855 ("G4VisManager::BeginDraw2D", 968 "visman0009", JustWarning, 856 "visman0009", JustWarning, 969 "Nesting detected. It is illegal to nes 857 "Nesting detected. It is illegal to nest Begin/EndDraw2D." 970 "\n Ignored"); 858 "\n Ignored"); 971 return; 859 return; 972 } 860 } 973 if (IsValidView ()) { 861 if (IsValidView ()) { 974 ClearTransientStoreIfMarked(); 862 ClearTransientStoreIfMarked(); 975 fpSceneHandler -> BeginPrimitives2D (objec 863 fpSceneHandler -> BeginPrimitives2D (objectTransform); 976 fIsDrawGroup = true; 864 fIsDrawGroup = true; 977 } 865 } 978 } 866 } 979 867 980 void G4VisManager::EndDraw2D () 868 void G4VisManager::EndDraw2D () 981 { 869 { >> 870 #ifdef G4MULTITHREADED 982 if (G4Threading::IsWorkerThread()) return; 871 if (G4Threading::IsWorkerThread()) return; 983 << 872 #endif 984 fDrawGroupNestingDepth--; 873 fDrawGroupNestingDepth--; 985 if (fDrawGroupNestingDepth != 0) { 874 if (fDrawGroupNestingDepth != 0) { 986 if (fDrawGroupNestingDepth < 0) fDrawGroup 875 if (fDrawGroupNestingDepth < 0) fDrawGroupNestingDepth = 0; 987 return; 876 return; 988 } 877 } 989 if (IsValidView ()) { 878 if (IsValidView ()) { 990 fpSceneHandler -> EndPrimitives2D (); 879 fpSceneHandler -> EndPrimitives2D (); 991 } 880 } 992 fIsDrawGroup = false; 881 fIsDrawGroup = false; 993 } 882 } 994 883 995 template <class T> void G4VisManager::DrawT 884 template <class T> void G4VisManager::DrawT 996 (const T& graphics_primitive, const G4Transfor 885 (const T& graphics_primitive, const G4Transform3D& objectTransform) { >> 886 #ifdef G4MULTITHREADED 997 if (G4Threading::IsWorkerThread()) return; 887 if (G4Threading::IsWorkerThread()) return; 998 << 888 #endif 999 if (fIsDrawGroup) { 889 if (fIsDrawGroup) { 1000 if (objectTransform != fpSceneHandler->Ge 890 if (objectTransform != fpSceneHandler->GetObjectTransformation()) { 1001 G4Exception 891 G4Exception 1002 ("G4VSceneHandler::DrawT", 892 ("G4VSceneHandler::DrawT", 1003 "visman0010", FatalException, 893 "visman0010", FatalException, 1004 "Different transform detected in Begin/End 894 "Different transform detected in Begin/EndDraw group."); 1005 } 895 } 1006 fpSceneHandler -> AddPrimitive (graphics_ 896 fpSceneHandler -> AddPrimitive (graphics_primitive); 1007 } else { 897 } else { 1008 if (IsValidView ()) { 898 if (IsValidView ()) { 1009 ClearTransientStoreIfMarked(); 899 ClearTransientStoreIfMarked(); 1010 fpSceneHandler -> BeginPrimitives (obje 900 fpSceneHandler -> BeginPrimitives (objectTransform); 1011 fpSceneHandler -> AddPrimitive (graphic 901 fpSceneHandler -> AddPrimitive (graphics_primitive); 1012 fpSceneHandler -> EndPrimitives (); 902 fpSceneHandler -> EndPrimitives (); 1013 } 903 } 1014 } 904 } 1015 } 905 } 1016 906 1017 template <class T> void G4VisManager::DrawT2D 907 template <class T> void G4VisManager::DrawT2D 1018 (const T& graphics_primitive, const G4Transfo 908 (const T& graphics_primitive, const G4Transform3D& objectTransform) { >> 909 #ifdef G4MULTITHREADED 1019 if (G4Threading::IsWorkerThread()) return; 910 if (G4Threading::IsWorkerThread()) return; 1020 << 911 #endif 1021 if (fIsDrawGroup) { 912 if (fIsDrawGroup) { 1022 if (objectTransform != fpSceneHandler->Ge 913 if (objectTransform != fpSceneHandler->GetObjectTransformation()) { 1023 G4Exception 914 G4Exception 1024 ("G4VSceneHandler::DrawT", 915 ("G4VSceneHandler::DrawT", 1025 "visman0011", FatalException, 916 "visman0011", FatalException, 1026 "Different transform detected in Begin/End 917 "Different transform detected in Begin/EndDraw2D group."); 1027 } 918 } 1028 fpSceneHandler -> AddPrimitive (graphics_ 919 fpSceneHandler -> AddPrimitive (graphics_primitive); 1029 } else { 920 } else { 1030 if (IsValidView ()) { 921 if (IsValidView ()) { 1031 ClearTransientStoreIfMarked(); 922 ClearTransientStoreIfMarked(); 1032 fpSceneHandler -> BeginPrimitives2D (ob 923 fpSceneHandler -> BeginPrimitives2D (objectTransform); 1033 fpSceneHandler -> AddPrimitive (graphic 924 fpSceneHandler -> AddPrimitive (graphics_primitive); 1034 fpSceneHandler -> EndPrimitives2D (); 925 fpSceneHandler -> EndPrimitives2D (); 1035 } 926 } 1036 } 927 } 1037 } 928 } 1038 929 1039 void G4VisManager::Draw (const G4Circle& circ 930 void G4VisManager::Draw (const G4Circle& circle, 1040 const G4Transform3D& objectTransform) 931 const G4Transform3D& objectTransform) 1041 { 932 { 1042 DrawT (circle, objectTransform); 933 DrawT (circle, objectTransform); 1043 } 934 } 1044 935 1045 void G4VisManager::Draw (const G4Polyhedron& 936 void G4VisManager::Draw (const G4Polyhedron& polyhedron, 1046 const G4Transform3D& objectTransform) 937 const G4Transform3D& objectTransform) 1047 { 938 { 1048 DrawT (polyhedron, objectTransform); 939 DrawT (polyhedron, objectTransform); 1049 } 940 } 1050 941 1051 void G4VisManager::Draw (const G4Polyline& li 942 void G4VisManager::Draw (const G4Polyline& line, 1052 const G4Transform3D& objectTransform) 943 const G4Transform3D& objectTransform) 1053 { 944 { 1054 DrawT (line, objectTransform); 945 DrawT (line, objectTransform); 1055 } 946 } 1056 947 1057 void G4VisManager::Draw (const G4Polymarker& 948 void G4VisManager::Draw (const G4Polymarker& polymarker, 1058 const G4Transform3D& objectTransform) 949 const G4Transform3D& objectTransform) 1059 { 950 { 1060 DrawT (polymarker, objectTransform); 951 DrawT (polymarker, objectTransform); 1061 } 952 } 1062 953 >> 954 void G4VisManager::Draw (const G4Scale& scale, >> 955 const G4Transform3D& objectTransform) >> 956 { >> 957 DrawT (scale, objectTransform); >> 958 } >> 959 1063 void G4VisManager::Draw (const G4Square& squa 960 void G4VisManager::Draw (const G4Square& square, 1064 const G4Transform3D& objectTransform) 961 const G4Transform3D& objectTransform) 1065 { 962 { 1066 DrawT (square, objectTransform); 963 DrawT (square, objectTransform); 1067 } 964 } 1068 965 1069 void G4VisManager::Draw (const G4Text& text, 966 void G4VisManager::Draw (const G4Text& text, 1070 const G4Transform3D& objectTransform) 967 const G4Transform3D& objectTransform) 1071 { 968 { 1072 DrawT (text, objectTransform); 969 DrawT (text, objectTransform); 1073 } 970 } 1074 971 1075 void G4VisManager::Draw2D (const G4Circle& ci 972 void G4VisManager::Draw2D (const G4Circle& circle, 1076 const G4Transform3D& objectTransform 973 const G4Transform3D& objectTransform) 1077 { 974 { 1078 DrawT2D (circle, objectTransform); 975 DrawT2D (circle, objectTransform); 1079 } 976 } 1080 977 1081 void G4VisManager::Draw2D (const G4Polyhedron 978 void G4VisManager::Draw2D (const G4Polyhedron& polyhedron, 1082 const G4Transform3D& objectTransform 979 const G4Transform3D& objectTransform) 1083 { 980 { 1084 DrawT2D (polyhedron, objectTransform); 981 DrawT2D (polyhedron, objectTransform); 1085 } 982 } 1086 983 1087 void G4VisManager::Draw2D (const G4Polyline& 984 void G4VisManager::Draw2D (const G4Polyline& line, 1088 const G4Transform3D& objectTransform 985 const G4Transform3D& objectTransform) 1089 { 986 { 1090 DrawT2D (line, objectTransform); 987 DrawT2D (line, objectTransform); 1091 } 988 } 1092 989 1093 void G4VisManager::Draw2D (const G4Polymarker 990 void G4VisManager::Draw2D (const G4Polymarker& polymarker, 1094 const G4Transform3D& objectTransform 991 const G4Transform3D& objectTransform) 1095 { 992 { 1096 DrawT2D (polymarker, objectTransform); 993 DrawT2D (polymarker, objectTransform); 1097 } 994 } 1098 995 1099 void G4VisManager::Draw2D (const G4Square& sq 996 void G4VisManager::Draw2D (const G4Square& square, 1100 const G4Transform3D& objectTransform 997 const G4Transform3D& objectTransform) 1101 { 998 { 1102 DrawT2D (square, objectTransform); 999 DrawT2D (square, objectTransform); 1103 } 1000 } 1104 1001 1105 void G4VisManager::Draw2D (const G4Text& text 1002 void G4VisManager::Draw2D (const G4Text& text, 1106 const G4Transform3D& objectTransform 1003 const G4Transform3D& objectTransform) 1107 { 1004 { 1108 DrawT2D (text, objectTransform); 1005 DrawT2D (text, objectTransform); 1109 } 1006 } 1110 1007 1111 void G4VisManager::Draw (const G4VHit& hit) { 1008 void G4VisManager::Draw (const G4VHit& hit) { >> 1009 #ifdef G4MULTITHREADED 1112 if (G4Threading::IsWorkerThread()) return; 1010 if (G4Threading::IsWorkerThread()) return; 1113 << 1011 #endif 1114 if (fIsDrawGroup) { 1012 if (fIsDrawGroup) { 1115 fpSceneHandler -> AddCompound (hit); 1013 fpSceneHandler -> AddCompound (hit); 1116 } else { 1014 } else { 1117 if (IsValidView ()) { 1015 if (IsValidView ()) { 1118 ClearTransientStoreIfMarked(); 1016 ClearTransientStoreIfMarked(); 1119 fpSceneHandler -> AddCompound (hit); 1017 fpSceneHandler -> AddCompound (hit); 1120 } 1018 } 1121 } 1019 } 1122 } 1020 } 1123 1021 1124 void G4VisManager::Draw (const G4VDigi& digi) 1022 void G4VisManager::Draw (const G4VDigi& digi) { >> 1023 #ifdef G4MULTITHREADED 1125 if (G4Threading::IsWorkerThread()) return; 1024 if (G4Threading::IsWorkerThread()) return; 1126 << 1025 #endif 1127 if (fIsDrawGroup) { 1026 if (fIsDrawGroup) { 1128 fpSceneHandler -> AddCompound (digi); 1027 fpSceneHandler -> AddCompound (digi); 1129 } else { 1028 } else { 1130 if (IsValidView ()) { 1029 if (IsValidView ()) { 1131 ClearTransientStoreIfMarked(); 1030 ClearTransientStoreIfMarked(); 1132 fpSceneHandler -> AddCompound (digi); 1031 fpSceneHandler -> AddCompound (digi); 1133 } 1032 } 1134 } 1033 } 1135 } 1034 } 1136 1035 1137 void G4VisManager::Draw (const G4VTrajectory& 1036 void G4VisManager::Draw (const G4VTrajectory& traj) { >> 1037 #ifdef G4MULTITHREADED 1138 if (G4Threading::IsWorkerThread()) return; 1038 if (G4Threading::IsWorkerThread()) return; 1139 << 1039 #endif 1140 // A trajectory needs a trajectories model 1040 // A trajectory needs a trajectories model to provide G4Atts, etc. 1141 static G4TrajectoriesModel trajectoriesMode 1041 static G4TrajectoriesModel trajectoriesModel; 1142 trajectoriesModel.SetCurrentTrajectory(&tra 1042 trajectoriesModel.SetCurrentTrajectory(&traj); 1143 G4RunManager* runManager = G4RunManagerFact 1043 G4RunManager* runManager = G4RunManagerFactory::GetMasterRunManager(); 1144 const G4Run* currentRun = runManager->GetC 1044 const G4Run* currentRun = runManager->GetCurrentRun(); 1145 if (currentRun) { 1045 if (currentRun) { 1146 trajectoriesModel.SetRunID(currentRun->Ge 1046 trajectoriesModel.SetRunID(currentRun->GetRunID()); 1147 } 1047 } 1148 const G4Event* currentEvent = 1048 const G4Event* currentEvent = 1149 G4EventManager::GetEventManager()->GetConst 1049 G4EventManager::GetEventManager()->GetConstCurrentEvent(); 1150 if (currentEvent) { 1050 if (currentEvent) { 1151 trajectoriesModel.SetEventID(currentEvent 1051 trajectoriesModel.SetEventID(currentEvent->GetEventID()); 1152 } 1052 } 1153 if (fIsDrawGroup) { 1053 if (fIsDrawGroup) { 1154 fpSceneHandler -> SetModel (&trajectories 1054 fpSceneHandler -> SetModel (&trajectoriesModel); 1155 fpSceneHandler -> AddCompound (traj); 1055 fpSceneHandler -> AddCompound (traj); 1156 fpSceneHandler -> SetModel (0); 1056 fpSceneHandler -> SetModel (0); 1157 } else { 1057 } else { 1158 if (IsValidView ()) { 1058 if (IsValidView ()) { 1159 ClearTransientStoreIfMarked(); 1059 ClearTransientStoreIfMarked(); 1160 fpSceneHandler -> SetModel (&trajectori 1060 fpSceneHandler -> SetModel (&trajectoriesModel); 1161 fpSceneHandler -> AddCompound (traj); 1061 fpSceneHandler -> AddCompound (traj); 1162 fpSceneHandler -> SetModel (0); 1062 fpSceneHandler -> SetModel (0); 1163 } 1063 } 1164 } 1064 } 1165 } 1065 } 1166 1066 1167 void G4VisManager::Draw (const G4LogicalVolum 1067 void G4VisManager::Draw (const G4LogicalVolume& logicalVol, 1168 const G4VisAttributes& attribs, 1068 const G4VisAttributes& attribs, 1169 const G4Transform3D& objectTransform) 1069 const G4Transform3D& objectTransform) { >> 1070 #ifdef G4MULTITHREADED 1170 if (G4Threading::IsWorkerThread()) return; 1071 if (G4Threading::IsWorkerThread()) return; 1171 << 1072 #endif 1172 // Find corresponding solid. 1073 // Find corresponding solid. 1173 G4VSolid* pSol = logicalVol.GetSolid (); 1074 G4VSolid* pSol = logicalVol.GetSolid (); 1174 Draw (*pSol, attribs, objectTransform); 1075 Draw (*pSol, attribs, objectTransform); 1175 } 1076 } 1176 1077 1177 void G4VisManager::Draw (const G4VSolid& soli 1078 void G4VisManager::Draw (const G4VSolid& solid, 1178 const G4VisAttributes& attribs, 1079 const G4VisAttributes& attribs, 1179 const G4Transform3D& objectTransform) 1080 const G4Transform3D& objectTransform) { >> 1081 #ifdef G4MULTITHREADED 1180 if (G4Threading::IsWorkerThread()) return; 1082 if (G4Threading::IsWorkerThread()) return; 1181 << 1083 #endif 1182 if (fIsDrawGroup) { 1084 if (fIsDrawGroup) { 1183 fpSceneHandler -> PreAddSolid (objectTran 1085 fpSceneHandler -> PreAddSolid (objectTransform, attribs); 1184 solid.DescribeYourselfTo (*fpSceneHandler 1086 solid.DescribeYourselfTo (*fpSceneHandler); 1185 fpSceneHandler -> PostAddSolid (); 1087 fpSceneHandler -> PostAddSolid (); 1186 } else { 1088 } else { 1187 if (IsValidView ()) { 1089 if (IsValidView ()) { 1188 ClearTransientStoreIfMarked(); 1090 ClearTransientStoreIfMarked(); 1189 fpSceneHandler -> PreAddSolid (objectTr 1091 fpSceneHandler -> PreAddSolid (objectTransform, attribs); 1190 solid.DescribeYourselfTo (*fpSceneHandl 1092 solid.DescribeYourselfTo (*fpSceneHandler); 1191 fpSceneHandler -> PostAddSolid (); 1093 fpSceneHandler -> PostAddSolid (); 1192 } 1094 } 1193 } 1095 } 1194 } 1096 } 1195 1097 1196 void G4VisManager::Draw (const G4VPhysicalVol 1098 void G4VisManager::Draw (const G4VPhysicalVolume& physicalVol, 1197 const G4VisAttributes& attribs, 1099 const G4VisAttributes& attribs, 1198 const G4Transform3D& objectTransform) 1100 const G4Transform3D& objectTransform) { >> 1101 #ifdef G4MULTITHREADED 1199 if (G4Threading::IsWorkerThread()) return; 1102 if (G4Threading::IsWorkerThread()) return; 1200 << 1103 #endif 1201 // Note: It is tempting to use a temporary 1104 // Note: It is tempting to use a temporary model here, as for 1202 // trajectories, in order to get at the G4A 1105 // trajectories, in order to get at the G4Atts of the physical 1203 // volume. I tried it (JA). But it's not 1106 // volume. I tried it (JA). But it's not easy to pass the 1204 // vis attributes. Also other aspects of t 1107 // vis attributes. Also other aspects of the model seem not to 1205 // be properly set up. So, the idea has be 1108 // be properly set up. So, the idea has been abandoned for the time 1206 // being. The model pointer will be null. 1109 // being. The model pointer will be null. So when picking there 1207 // will be no G4Atts from this physical vol 1110 // will be no G4Atts from this physical volume. 1208 // 1111 // 1209 // If this is called from DrawHit, for exam 1112 // If this is called from DrawHit, for example, the user may G4Atts to the 1210 // hit and these will be available with "/v 1113 // hit and these will be available with "/vis/scene/add/hits". 1211 // 1114 // 1212 // Find corresponding logical volume and so 1115 // Find corresponding logical volume and solid. 1213 G4LogicalVolume* pLV = physicalVol.GetLogi 1116 G4LogicalVolume* pLV = physicalVol.GetLogicalVolume (); 1214 G4VSolid* pSol = pLV -> GetSolid (); 1117 G4VSolid* pSol = pLV -> GetSolid (); 1215 Draw (*pSol, attribs, objectTransform); 1118 Draw (*pSol, attribs, objectTransform); 1216 } 1119 } 1217 1120 1218 void G4VisManager::DrawGeometry << 1219 (G4VPhysicalVolume* v, const G4Transform3D& t << 1220 // Draws a geometry tree starting at the spec << 1221 { << 1222 auto modelingParameters = fpSceneHandler->C << 1223 auto depth = G4PhysicalVolumeModel::UNLIMIT << 1224 const G4bool useFullExtent = true; << 1225 G4PhysicalVolumeModel aPVModel(v,depth,t,mo << 1226 aPVModel.DescribeYourselfTo(*fpSceneHandler << 1227 delete modelingParameters; << 1228 } << 1229 << 1230 void G4VisManager::CreateSceneHandler (const 1121 void G4VisManager::CreateSceneHandler (const G4String& name) { 1231 if (!fInitialised) Initialise (); 1122 if (!fInitialised) Initialise (); 1232 if (fpGraphicsSystem) { 1123 if (fpGraphicsSystem) { 1233 G4VSceneHandler* pSceneHandler = 1124 G4VSceneHandler* pSceneHandler = 1234 fpGraphicsSystem -> CreateSceneHandler 1125 fpGraphicsSystem -> CreateSceneHandler (name); 1235 if (pSceneHandler) { 1126 if (pSceneHandler) { 1236 fAvailableSceneHandlers.push_back (pSce 1127 fAvailableSceneHandlers.push_back (pSceneHandler); 1237 fpSceneHandler = pSceneHandler; 1128 fpSceneHandler = pSceneHandler; // Make current. 1238 } 1129 } 1239 else { 1130 else { 1240 if (fVerbosity >= errors) { 1131 if (fVerbosity >= errors) { 1241 G4warn << "ERROR in G4VisManager::CreateSce << 1132 G4cout << "ERROR in G4VisManager::CreateSceneHandler during " 1242 << fpGraphicsSystem -> GetName () 1133 << fpGraphicsSystem -> GetName () 1243 << " scene handler creation.\n No a 1134 << " scene handler creation.\n No action taken." 1244 << G4endl; 1135 << G4endl; 1245 } 1136 } 1246 } 1137 } 1247 } 1138 } 1248 else PrintInvalidPointers (); 1139 else PrintInvalidPointers (); 1249 } 1140 } 1250 1141 1251 void G4VisManager::CreateViewer 1142 void G4VisManager::CreateViewer 1252 (const G4String& name, const G4String& XGeome 1143 (const G4String& name, const G4String& XGeometry) 1253 { 1144 { 1254 1145 1255 if (!fInitialised) Initialise (); 1146 if (!fInitialised) Initialise (); 1256 1147 1257 if (!fpSceneHandler) { 1148 if (!fpSceneHandler) { 1258 PrintInvalidPointers (); 1149 PrintInvalidPointers (); 1259 return; 1150 return; 1260 } 1151 } 1261 1152 1262 G4VViewer* p = fpGraphicsSystem -> CreateVi 1153 G4VViewer* p = fpGraphicsSystem -> CreateViewer (*fpSceneHandler, name); 1263 1154 1264 if (!p) { 1155 if (!p) { 1265 if (fVerbosity >= errors) { 1156 if (fVerbosity >= errors) { 1266 G4warn << "ERROR in G4VisManager::Creat << 1157 G4cerr << "ERROR in G4VisManager::CreateViewer: null pointer during " 1267 << fpGraphicsSystem -> GetName () 1158 << fpGraphicsSystem -> GetName () 1268 << " viewer creation.\n No action tak 1159 << " viewer creation.\n No action taken." 1269 << G4endl; 1160 << G4endl; 1270 } 1161 } 1271 return; 1162 return; 1272 } 1163 } 1273 1164 1274 if (p -> GetViewId() < 0) { 1165 if (p -> GetViewId() < 0) { 1275 if (fVerbosity >= errors) { 1166 if (fVerbosity >= errors) { 1276 G4warn << "ERROR in G4VisManager::Creat << 1167 G4cerr << "ERROR in G4VisManager::CreateViewer during " 1277 << fpGraphicsSystem -> GetName () 1168 << fpGraphicsSystem -> GetName () 1278 << " viewer instantiation.\n No actio 1169 << " viewer instantiation.\n No action taken." 1279 << G4endl; 1170 << G4endl; 1280 } 1171 } 1281 return; 1172 return; 1282 } 1173 } 1283 1174 1284 // Viewer is created, now we can set geomet 1175 // Viewer is created, now we can set geometry parameters 1285 // Before 12/2008, it was done in G4VViewer 1176 // Before 12/2008, it was done in G4VViewer.cc but it did not have to be there! 1286 1177 1287 G4ViewParameters initialvp = p -> GetViewPa 1178 G4ViewParameters initialvp = p -> GetViewParameters(); 1288 initialvp.SetXGeometryString(XGeometry); // 1179 initialvp.SetXGeometryString(XGeometry); //parse string and store parameters 1289 p -> SetViewParameters(initialvp); 1180 p -> SetViewParameters(initialvp); 1290 p -> Initialise (); // (Viewer itself may 1181 p -> Initialise (); // (Viewer itself may change view parameters further.) 1291 if (p -> GetViewId() < 0) { 1182 if (p -> GetViewId() < 0) { 1292 if (fVerbosity >= errors) { 1183 if (fVerbosity >= errors) { 1293 G4warn << "ERROR in G4VisManager::Creat << 1184 G4cerr << "ERROR in G4VisManager::CreateViewer during " 1294 << fpGraphicsSystem -> GetName () 1185 << fpGraphicsSystem -> GetName () 1295 << " viewer initialisation.\n No acti 1186 << " viewer initialisation.\n No action taken." 1296 << G4endl; 1187 << G4endl; 1297 } 1188 } 1298 return; 1189 return; 1299 } 1190 } 1300 1191 1301 fpViewer = p; / 1192 fpViewer = p; // Make current. 1302 fpSceneHandler -> AddViewerToList (fpViewer 1193 fpSceneHandler -> AddViewerToList (fpViewer); 1303 fpSceneHandler -> SetCurrentViewer (fpViewe 1194 fpSceneHandler -> SetCurrentViewer (fpViewer); 1304 if (fVerbosity >= confirmations) { 1195 if (fVerbosity >= confirmations) { 1305 G4cout << "G4VisManager::CreateViewer: ne 1196 G4cout << "G4VisManager::CreateViewer: new viewer created." 1306 << G4endl; 1197 << G4endl; 1307 } 1198 } 1308 1199 1309 const G4ViewParameters& vp = fpViewer->GetV 1200 const G4ViewParameters& vp = fpViewer->GetViewParameters(); 1310 if (fVerbosity >= parameters) { 1201 if (fVerbosity >= parameters) { 1311 G4cout << " view parameters are:\n " << 1202 G4cout << " view parameters are:\n " << vp << G4endl; 1312 } 1203 } 1313 1204 1314 if (vp.IsCulling () && vp.IsCullingInvisibl 1205 if (vp.IsCulling () && vp.IsCullingInvisible ()) { 1315 static G4bool warned = false; 1206 static G4bool warned = false; 1316 if (fVerbosity >= confirmations) { 1207 if (fVerbosity >= confirmations) { 1317 if (!warned) { 1208 if (!warned) { 1318 G4cout << 1209 G4cout << 1319 "NOTE: objects with visibility flag set to 1210 "NOTE: objects with visibility flag set to \"false\"" 1320 " will not be drawn!" 1211 " will not be drawn!" 1321 "\n \"/vis/viewer/set/culling global false 1212 "\n \"/vis/viewer/set/culling global false\" to Draw such objects." 1322 "\n Also see other \"/vis/viewer/set\" com 1213 "\n Also see other \"/vis/viewer/set\" commands." 1323 << G4endl; 1214 << G4endl; 1324 warned = true; 1215 warned = true; 1325 } 1216 } 1326 } 1217 } 1327 } 1218 } 1328 if (vp.IsCullingCovered ()) { 1219 if (vp.IsCullingCovered ()) { 1329 static G4bool warned = false; 1220 static G4bool warned = false; 1330 if (fVerbosity >= warnings) { 1221 if (fVerbosity >= warnings) { 1331 if (!warned) { 1222 if (!warned) { 1332 G4warn << << 1223 G4cout << 1333 "WARNING: covered objects in solid mode wil 1224 "WARNING: covered objects in solid mode will not be rendered!" 1334 "\n \"/vis/viewer/set/culling coveredDaugh 1225 "\n \"/vis/viewer/set/culling coveredDaughters false\" to reverse this." 1335 "\n Also see other \"/vis/viewer/set\" com 1226 "\n Also see other \"/vis/viewer/set\" commands." 1336 << G4endl; 1227 << G4endl; 1337 warned = true; 1228 warned = true; 1338 } 1229 } 1339 } 1230 } 1340 } 1231 } 1341 } 1232 } 1342 1233 1343 void G4VisManager::GeometryHasChanged () { 1234 void G4VisManager::GeometryHasChanged () { 1344 if (fVerbosity >= confirmations) { 1235 if (fVerbosity >= confirmations) { 1345 G4cout << "G4VisManager::GeometryHasChang 1236 G4cout << "G4VisManager::GeometryHasChanged() called." << G4endl; 1346 } 1237 } 1347 1238 1348 // Change the world... 1239 // Change the world... 1349 G4VPhysicalVolume* pWorld = 1240 G4VPhysicalVolume* pWorld = 1350 G4TransportationManager::GetTransportatio 1241 G4TransportationManager::GetTransportationManager () 1351 -> GetNavigatorForTracking () -> GetWorld 1242 -> GetNavigatorForTracking () -> GetWorldVolume (); 1352 if (!pWorld) { 1243 if (!pWorld) { 1353 if (fVerbosity >= warnings) { 1244 if (fVerbosity >= warnings) { 1354 G4warn << "WARNING: There is no world v << 1245 G4cout << "WARNING: There is no world volume!" << G4endl; 1355 } 1246 } 1356 } 1247 } 1357 1248 1358 // Check scenes. 1249 // Check scenes. 1359 G4SceneList& sceneList = fSceneList; 1250 G4SceneList& sceneList = fSceneList; 1360 std::size_t iScene, nScenes = sceneList.siz << 1251 G4int iScene, nScenes = sceneList.size (); 1361 for (iScene = 0; iScene < nScenes; ++iScene << 1252 for (iScene = 0; iScene < nScenes; iScene++) { 1362 G4Scene* pScene = sceneList [iScene]; 1253 G4Scene* pScene = sceneList [iScene]; 1363 std::vector<G4Scene::Model>& modelList = 1254 std::vector<G4Scene::Model>& modelList = pScene -> SetRunDurationModelList (); 1364 if (modelList.size ()) { 1255 if (modelList.size ()) { 1365 G4bool modelInvalid; 1256 G4bool modelInvalid; 1366 do { // Remove, if required, one at a 1257 do { // Remove, if required, one at a time. 1367 modelInvalid = false; 1258 modelInvalid = false; 1368 std::vector<G4Scene::Model>::iterator iterM 1259 std::vector<G4Scene::Model>::iterator iterModel; 1369 for (iterModel = modelList.begin(); 1260 for (iterModel = modelList.begin(); 1370 iterModel != modelList.end(); 1261 iterModel != modelList.end(); 1371 ++iterModel) { 1262 ++iterModel) { 1372 modelInvalid = !(iterModel->fpModel->Vali 1263 modelInvalid = !(iterModel->fpModel->Validate(fVerbosity>=warnings)); 1373 if (modelInvalid) { 1264 if (modelInvalid) { 1374 // Model invalid - remove and break. 1265 // Model invalid - remove and break. 1375 if (fVerbosity >= warnings) { 1266 if (fVerbosity >= warnings) { 1376 G4warn << "WARNING: Model \"" << 1267 G4cout << "WARNING: Model \"" 1377 << iterModel->fpModel->GetGlobalDesc 1268 << iterModel->fpModel->GetGlobalDescription () 1378 << 1269 << 1379 "\" is no longer valid - being removed\n 1270 "\" is no longer valid - being removed\n from scene \"" 1380 << pScene -> GetName () << "\"" 1271 << pScene -> GetName () << "\"" 1381 << G4endl; 1272 << G4endl; 1382 } 1273 } 1383 modelList.erase (iterModel); 1274 modelList.erase (iterModel); 1384 break; 1275 break; 1385 } 1276 } 1386 } 1277 } 1387 } while (modelInvalid); 1278 } while (modelInvalid); 1388 1279 1389 if (modelList.size () == 0) { 1280 if (modelList.size () == 0) { 1390 if (fVerbosity >= warnings) { 1281 if (fVerbosity >= warnings) { 1391 G4warn << "WARNING: No run-duration model << 1282 G4cout << "WARNING: No models left in this scene \"" 1392 << pScene -> GetName () 1283 << pScene -> GetName () 1393 << "\"." 1284 << "\"." 1394 << G4endl; 1285 << G4endl; 1395 } 1286 } 1396 if (pWorld) { << 1397 if (fVerbosity >= warnings) { << 1398 G4warn << " Adding current world << 1399 << pScene -> GetName () << 1400 << "\"." << 1401 << G4endl; << 1402 } << 1403 pScene->AddRunDurationModel(new G4P << 1404 // (The above includes a re-calcula << 1405 G4UImanager::GetUIpointer () -> << 1406 ApplyCommand (G4String("/vis/scene/ << 1407 } << 1408 } 1287 } 1409 else { 1288 else { 1410 pScene->CalculateExtent(); // Recalculate << 1289 pScene->CalculateExtent(); 1411 G4UImanager::GetUIpointer () -> 1290 G4UImanager::GetUIpointer () -> 1412 ApplyCommand (G4String("/vis/scene/notify 1291 ApplyCommand (G4String("/vis/scene/notifyHandlers " + pScene->GetName())); 1413 } 1292 } 1414 } 1293 } 1415 } 1294 } 1416 1295 1417 // Check the manager's current scene... 1296 // Check the manager's current scene... 1418 if (fpScene && fpScene -> GetRunDurationMod 1297 if (fpScene && fpScene -> GetRunDurationModelList ().size () == 0) { 1419 if (fVerbosity >= warnings) { 1298 if (fVerbosity >= warnings) { 1420 G4warn << "WARNING: The current scene \ << 1299 G4cout << "WARNING: The current scene \"" 1421 << fpScene -> GetName () 1300 << fpScene -> GetName () 1422 << "\" has no run duration models." 1301 << "\" has no run duration models." 1423 << "\n Use \"/vis/scene/add/vol 1302 << "\n Use \"/vis/scene/add/volume\" or create a new scene." 1424 << G4endl; 1303 << G4endl; 1425 } 1304 } 1426 // Clean up << 1305 fpSceneHandler->ClearTransientStore(); 1427 if (fpSceneHandler) { << 1306 fpSceneHandler->ClearStore(); 1428 fpSceneHandler->ClearTransientStore(); << 1307 fpViewer->NeedKernelVisit(); 1429 fpSceneHandler->ClearStore(); << 1308 fpViewer->SetView(); 1430 if (fpViewer) { << 1309 fpViewer->ClearView(); 1431 fpViewer->NeedKernelVisit(); << 1310 fpViewer->FinishView(); 1432 fpViewer->SetView(); << 1433 fpViewer->ClearView(); << 1434 fpViewer->FinishView(); << 1435 } << 1436 } << 1437 } 1311 } 1438 } 1312 } 1439 1313 1440 void G4VisManager::NotifyHandlers () { 1314 void G4VisManager::NotifyHandlers () { 1441 1315 1442 if (fVerbosity >= confirmations) { 1316 if (fVerbosity >= confirmations) { 1443 G4cout << "G4VisManager::NotifyHandler() 1317 G4cout << "G4VisManager::NotifyHandler() called." << G4endl; 1444 } 1318 } 1445 1319 1446 if (IsValidView()) { 1320 if (IsValidView()) { 1447 1321 1448 // Check scenes. 1322 // Check scenes. 1449 G4SceneList& sceneList = fSceneList; 1323 G4SceneList& sceneList = fSceneList; 1450 std::size_t iScene, nScenes = sceneList.s << 1324 G4int iScene, nScenes = sceneList.size (); 1451 for (iScene = 0; iScene < nScenes; ++iSce << 1325 for (iScene = 0; iScene < nScenes; iScene++) { 1452 G4Scene* pScene = sceneList [iScene]; 1326 G4Scene* pScene = sceneList [iScene]; 1453 std::vector<G4Scene::Model>& modelList 1327 std::vector<G4Scene::Model>& modelList = pScene -> SetRunDurationModelList (); 1454 1328 1455 if (modelList.size ()) { 1329 if (modelList.size ()) { 1456 pScene->CalculateExtent(); 1330 pScene->CalculateExtent(); 1457 G4UImanager::GetUIpointer () -> 1331 G4UImanager::GetUIpointer () -> 1458 ApplyCommand (G4String("/vis/scene/no 1332 ApplyCommand (G4String("/vis/scene/notifyHandlers " + pScene->GetName())); 1459 } 1333 } 1460 } 1334 } 1461 1335 1462 // Check the manager's current scene... 1336 // Check the manager's current scene... 1463 if (fpScene && fpScene -> GetRunDurationM 1337 if (fpScene && fpScene -> GetRunDurationModelList ().size () == 0) { 1464 if (fVerbosity >= warnings) { 1338 if (fVerbosity >= warnings) { 1465 G4warn << "WARNING: The current scene << 1339 G4cout << "WARNING: The current scene \"" 1466 << fpScene -> GetName () 1340 << fpScene -> GetName () 1467 << "\" has no run duration models." 1341 << "\" has no run duration models." 1468 << "\n Use \"/vis/scene/add/volume\" 1342 << "\n Use \"/vis/scene/add/volume\" or create a new scene." 1469 << G4endl; 1343 << G4endl; 1470 } 1344 } 1471 fpSceneHandler->ClearTransientStore(); 1345 fpSceneHandler->ClearTransientStore(); 1472 fpSceneHandler->ClearStore(); 1346 fpSceneHandler->ClearStore(); 1473 fpViewer->NeedKernelVisit(); 1347 fpViewer->NeedKernelVisit(); 1474 fpViewer->SetView(); 1348 fpViewer->SetView(); 1475 fpViewer->ClearView(); 1349 fpViewer->ClearView(); 1476 fpViewer->FinishView(); 1350 fpViewer->FinishView(); 1477 } 1351 } 1478 } 1352 } 1479 } 1353 } 1480 1354 1481 G4bool G4VisManager::FilterTrajectory(const G 1355 G4bool G4VisManager::FilterTrajectory(const G4VTrajectory& trajectory) 1482 { 1356 { 1483 return fpTrajFilterMgr->Accept(trajectory); 1357 return fpTrajFilterMgr->Accept(trajectory); 1484 } 1358 } 1485 1359 1486 G4bool G4VisManager::FilterHit(const G4VHit& 1360 G4bool G4VisManager::FilterHit(const G4VHit& hit) 1487 { 1361 { 1488 return fpHitFilterMgr->Accept(hit); 1362 return fpHitFilterMgr->Accept(hit); 1489 } 1363 } 1490 1364 1491 G4bool G4VisManager::FilterDigi(const G4VDigi 1365 G4bool G4VisManager::FilterDigi(const G4VDigi& digi) 1492 { 1366 { 1493 return fpDigiFilterMgr->Accept(digi); 1367 return fpDigiFilterMgr->Accept(digi); 1494 } 1368 } 1495 1369 1496 void G4VisManager::DispatchToModel(const G4VT 1370 void G4VisManager::DispatchToModel(const G4VTrajectory& trajectory) 1497 { 1371 { 1498 G4bool visible(true); 1372 G4bool visible(true); 1499 1373 1500 // See if trajectory passes filter 1374 // See if trajectory passes filter 1501 G4bool passed = FilterTrajectory(trajectory 1375 G4bool passed = FilterTrajectory(trajectory); 1502 1376 1503 if (!passed) { 1377 if (!passed) { 1504 // Draw invisible trajectory if trajector 1378 // Draw invisible trajectory if trajectory failed filter and 1505 // are filtering in soft mode 1379 // are filtering in soft mode 1506 if (fpTrajFilterMgr->GetMode() == FilterM 1380 if (fpTrajFilterMgr->GetMode() == FilterMode::Soft) visible = false; 1507 else {return;} 1381 else {return;} 1508 } 1382 } 1509 1383 1510 // Go on to draw trajectory 1384 // Go on to draw trajectory 1511 assert (0 != fpTrajDrawModelMgr); 1385 assert (0 != fpTrajDrawModelMgr); 1512 1386 1513 const G4VTrajectoryModel* trajectoryModel = 1387 const G4VTrajectoryModel* trajectoryModel = CurrentTrajDrawModel(); 1514 1388 1515 assert (0 != trajectoryModel); // Should ex 1389 assert (0 != trajectoryModel); // Should exist 1516 1390 1517 if (IsValidView()) { 1391 if (IsValidView()) { 1518 trajectoryModel->Draw(trajectory, visib 1392 trajectoryModel->Draw(trajectory, visible); 1519 } 1393 } 1520 } 1394 } 1521 1395 >> 1396 void G4VisManager::SetUserAction >> 1397 (G4VUserVisAction* pVisAction, >> 1398 const G4VisExtent& extent) { >> 1399 if (fVerbosity >= warnings) { >> 1400 G4cout << >> 1401 "WARNING: SetUserAction is deprecated. Use RegisterRunDurationUserVisAction." >> 1402 << G4endl; >> 1403 } >> 1404 RegisterRunDurationUserVisAction("SetUserAction",pVisAction,extent); >> 1405 } >> 1406 1522 void G4VisManager::RegisterRunDurationUserVis 1407 void G4VisManager::RegisterRunDurationUserVisAction 1523 (const G4String& name, 1408 (const G4String& name, 1524 G4VUserVisAction* pVisAction, 1409 G4VUserVisAction* pVisAction, 1525 const G4VisExtent& extent) { 1410 const G4VisExtent& extent) { 1526 fRunDurationUserVisActions.push_back(UserVi 1411 fRunDurationUserVisActions.push_back(UserVisAction(name,pVisAction)); 1527 if (extent.GetExtentRadius() > 0.) { 1412 if (extent.GetExtentRadius() > 0.) { 1528 fUserVisActionExtents[pVisAction] = exten 1413 fUserVisActionExtents[pVisAction] = extent; 1529 } else { 1414 } else { 1530 if (fVerbosity >= warnings) { 1415 if (fVerbosity >= warnings) { 1531 G4warn << << 1416 G4cout << 1532 "WARNING: No extent set for user vis action 1417 "WARNING: No extent set for user vis action \"" << name << "\"." 1533 << G4endl; 1418 << G4endl; 1534 } 1419 } 1535 } 1420 } 1536 if (fVerbosity >= confirmations) { 1421 if (fVerbosity >= confirmations) { 1537 G4cout 1422 G4cout 1538 << "Run duration user vis action \"" << n 1423 << "Run duration user vis action \"" << name << "\" registered" 1539 << G4endl; 1424 << G4endl; 1540 } 1425 } 1541 } 1426 } 1542 1427 1543 void G4VisManager::RegisterEndOfEventUserVisA 1428 void G4VisManager::RegisterEndOfEventUserVisAction 1544 (const G4String& name, 1429 (const G4String& name, 1545 G4VUserVisAction* pVisAction, 1430 G4VUserVisAction* pVisAction, 1546 const G4VisExtent& extent) { 1431 const G4VisExtent& extent) { 1547 fEndOfEventUserVisActions.push_back(UserVis 1432 fEndOfEventUserVisActions.push_back(UserVisAction(name,pVisAction)); 1548 if (extent.GetExtentRadius() > 0.) { 1433 if (extent.GetExtentRadius() > 0.) { 1549 fUserVisActionExtents[pVisAction] = exten 1434 fUserVisActionExtents[pVisAction] = extent; 1550 } else { 1435 } else { 1551 if (fVerbosity >= warnings) { 1436 if (fVerbosity >= warnings) { 1552 G4warn << << 1437 G4cout << 1553 "WARNING: No extent set for user vis action 1438 "WARNING: No extent set for user vis action \"" << name << "\"." 1554 << G4endl; 1439 << G4endl; 1555 } 1440 } 1556 } 1441 } 1557 if (fVerbosity >= confirmations) { 1442 if (fVerbosity >= confirmations) { 1558 G4cout 1443 G4cout 1559 << "End of event user vis action \"" << n 1444 << "End of event user vis action \"" << name << "\" registered" 1560 << G4endl; 1445 << G4endl; 1561 } 1446 } 1562 } 1447 } 1563 1448 1564 void G4VisManager::RegisterEndOfRunUserVisAct 1449 void G4VisManager::RegisterEndOfRunUserVisAction 1565 (const G4String& name, 1450 (const G4String& name, 1566 G4VUserVisAction* pVisAction, 1451 G4VUserVisAction* pVisAction, 1567 const G4VisExtent& extent) { 1452 const G4VisExtent& extent) { 1568 fEndOfRunUserVisActions.push_back(UserVisAc 1453 fEndOfRunUserVisActions.push_back(UserVisAction(name,pVisAction)); 1569 if (extent.GetExtentRadius() > 0.) { 1454 if (extent.GetExtentRadius() > 0.) { 1570 fUserVisActionExtents[pVisAction] = exten 1455 fUserVisActionExtents[pVisAction] = extent; 1571 } else { 1456 } else { 1572 if (fVerbosity >= warnings) { 1457 if (fVerbosity >= warnings) { 1573 G4warn << << 1458 G4cout << 1574 "WARNING: No extent set for user vis action 1459 "WARNING: No extent set for user vis action \"" << name << "\"." 1575 << G4endl; 1460 << G4endl; 1576 } 1461 } 1577 } 1462 } 1578 if (fVerbosity >= confirmations) { 1463 if (fVerbosity >= confirmations) { 1579 G4cout 1464 G4cout 1580 << "End of run user vis action \"" << nam 1465 << "End of run user vis action \"" << name << "\" registered" 1581 << G4endl; 1466 << G4endl; 1582 } 1467 } 1583 } 1468 } 1584 1469 1585 void G4VisManager::SetCurrentScene (G4Scene* 1470 void G4VisManager::SetCurrentScene (G4Scene* pScene) { 1586 if (pScene != fpScene) { 1471 if (pScene != fpScene) { 1587 // A change of scene. Therefore reset tr 1472 // A change of scene. Therefore reset transients drawn flags. All 1588 // memory of previous transient proceessi 1473 // memory of previous transient proceessing thereby erased... 1589 ResetTransientsDrawnFlags(); 1474 ResetTransientsDrawnFlags(); 1590 } 1475 } 1591 fpScene = pScene; 1476 fpScene = pScene; 1592 } 1477 } 1593 1478 1594 void G4VisManager::SetCurrentGraphicsSystem ( 1479 void G4VisManager::SetCurrentGraphicsSystem (G4VGraphicsSystem* pSystem) { 1595 fpGraphicsSystem = pSystem; 1480 fpGraphicsSystem = pSystem; 1596 if (fVerbosity >= confirmations) { 1481 if (fVerbosity >= confirmations) { 1597 G4cout << "G4VisManager::SetCurrentGraphi 1482 G4cout << "G4VisManager::SetCurrentGraphicsSystem: system now " 1598 << pSystem -> GetName () << G4endl; 1483 << pSystem -> GetName () << G4endl; 1599 } 1484 } 1600 // If current scene handler is of same grap 1485 // If current scene handler is of same graphics system, leave unchanged. 1601 // Else find the most recent scene handler 1486 // Else find the most recent scene handler of same graphics system. 1602 // Or clear pointers. 1487 // Or clear pointers. 1603 if (!(fpSceneHandler && fpSceneHandler -> G 1488 if (!(fpSceneHandler && fpSceneHandler -> GetGraphicsSystem () == pSystem)) { 1604 const G4SceneHandlerList& sceneHandlerLis 1489 const G4SceneHandlerList& sceneHandlerList = fAvailableSceneHandlers; 1605 G4int nSH = (G4int)sceneHandlerList.size << 1490 G4int nSH = sceneHandlerList.size (); // No. of scene handlers. 1606 G4int iSH; 1491 G4int iSH; 1607 for (iSH = nSH - 1; iSH >= 0; iSH--) { 1492 for (iSH = nSH - 1; iSH >= 0; iSH--) { 1608 if (sceneHandlerList [iSH] -> GetGraphi 1493 if (sceneHandlerList [iSH] -> GetGraphicsSystem () == pSystem) break; 1609 } 1494 } 1610 if (iSH >= 0) { 1495 if (iSH >= 0) { 1611 fpSceneHandler = sceneHandlerList [iSH] 1496 fpSceneHandler = sceneHandlerList [iSH]; 1612 if (fVerbosity >= confirmations) { 1497 if (fVerbosity >= confirmations) { 1613 G4cout << " Scene Handler now " 1498 G4cout << " Scene Handler now " 1614 << fpSceneHandler -> GetName () << G 1499 << fpSceneHandler -> GetName () << G4endl; 1615 } 1500 } 1616 if (fpScene != fpSceneHandler -> GetSce 1501 if (fpScene != fpSceneHandler -> GetScene ()) { 1617 fpScene = fpSceneHandler -> GetScene (); 1502 fpScene = fpSceneHandler -> GetScene (); 1618 if (fVerbosity >= confirmations) { 1503 if (fVerbosity >= confirmations) { 1619 G4cout << " Scene now \"" 1504 G4cout << " Scene now \"" 1620 << fpScene -> GetName () << "\"" << G4en 1505 << fpScene -> GetName () << "\"" << G4endl; 1621 } 1506 } 1622 } 1507 } 1623 const G4ViewerList& viewerList = fpScen 1508 const G4ViewerList& viewerList = fpSceneHandler -> GetViewerList (); 1624 if (viewerList.size ()) { 1509 if (viewerList.size ()) { 1625 fpViewer = viewerList [0]; 1510 fpViewer = viewerList [0]; 1626 if (fVerbosity >= confirmations) { 1511 if (fVerbosity >= confirmations) { 1627 G4cout << " Viewer now " << fpViewer -> 1512 G4cout << " Viewer now " << fpViewer -> GetName () << G4endl; 1628 } 1513 } 1629 } 1514 } 1630 else { 1515 else { 1631 fpViewer = 0; 1516 fpViewer = 0; 1632 } 1517 } 1633 } 1518 } 1634 else { 1519 else { 1635 fpSceneHandler = 0; 1520 fpSceneHandler = 0; 1636 fpViewer = 0; 1521 fpViewer = 0; 1637 } 1522 } 1638 } 1523 } 1639 } 1524 } 1640 1525 1641 void G4VisManager::SetCurrentSceneHandler (G4 1526 void G4VisManager::SetCurrentSceneHandler (G4VSceneHandler* pSceneHandler) { 1642 fpSceneHandler = pSceneHandler; 1527 fpSceneHandler = pSceneHandler; 1643 if (fVerbosity >= confirmations) { 1528 if (fVerbosity >= confirmations) { 1644 G4cout << "G4VisManager::SetCurrentSceneH 1529 G4cout << "G4VisManager::SetCurrentSceneHandler: scene handler now \"" 1645 << pSceneHandler -> GetName () << "\"" < 1530 << pSceneHandler -> GetName () << "\"" << G4endl; 1646 } 1531 } 1647 if (fpScene != fpSceneHandler -> GetScene ( 1532 if (fpScene != fpSceneHandler -> GetScene ()) { 1648 fpScene = fpSceneHandler -> GetScene (); 1533 fpScene = fpSceneHandler -> GetScene (); 1649 if (fVerbosity >= confirmations) { 1534 if (fVerbosity >= confirmations) { 1650 G4cout << " Scene now \"" 1535 G4cout << " Scene now \"" 1651 << fpScene -> GetName () << "\"" << G4 1536 << fpScene -> GetName () << "\"" << G4endl; 1652 } 1537 } 1653 } 1538 } 1654 if (fpGraphicsSystem != pSceneHandler -> Ge 1539 if (fpGraphicsSystem != pSceneHandler -> GetGraphicsSystem ()) { 1655 fpGraphicsSystem = pSceneHandler -> GetGr 1540 fpGraphicsSystem = pSceneHandler -> GetGraphicsSystem (); 1656 if (fVerbosity >= confirmations) { 1541 if (fVerbosity >= confirmations) { 1657 G4cout << " Graphics system now \"" 1542 G4cout << " Graphics system now \"" 1658 << fpGraphicsSystem -> GetName () << " 1543 << fpGraphicsSystem -> GetName () << "\"" << G4endl; 1659 } 1544 } 1660 } 1545 } 1661 const G4ViewerList& viewerList = fpSceneHan 1546 const G4ViewerList& viewerList = fpSceneHandler -> GetViewerList (); 1662 std::size_t nViewers = viewerList.size (); << 1547 G4int nViewers = viewerList.size (); 1663 if (nViewers) { 1548 if (nViewers) { 1664 std::size_t iViewer; << 1549 G4int iViewer; 1665 for (iViewer = 0; iViewer < nViewers; ++i << 1550 for (iViewer = 0; iViewer < nViewers; iViewer++) { 1666 if (fpViewer == viewerList [iViewer]) b 1551 if (fpViewer == viewerList [iViewer]) break; 1667 } 1552 } 1668 if (iViewer >= nViewers) { 1553 if (iViewer >= nViewers) { 1669 fpViewer = viewerList [0]; 1554 fpViewer = viewerList [0]; 1670 if (fVerbosity >= confirmations) { 1555 if (fVerbosity >= confirmations) { 1671 G4cout << " Viewer now \"" << fpViewer -> 1556 G4cout << " Viewer now \"" << fpViewer -> GetName () << "\"" 1672 << G4endl; 1557 << G4endl; 1673 } 1558 } 1674 } 1559 } 1675 if (!IsValidView ()) { 1560 if (!IsValidView ()) { 1676 if (fVerbosity >= warnings) { 1561 if (fVerbosity >= warnings) { 1677 G4warn << << 1562 G4cout << 1678 "WARNING: Problem setting scene handler - p 1563 "WARNING: Problem setting scene handler - please report circumstances." 1679 << G4endl; 1564 << G4endl; 1680 } 1565 } 1681 } 1566 } 1682 } 1567 } 1683 else { 1568 else { 1684 fpViewer = 0; 1569 fpViewer = 0; 1685 if (fVerbosity >= warnings) { 1570 if (fVerbosity >= warnings) { 1686 G4warn << << 1571 G4cout << 1687 "WARNING: No viewers for this scene handler 1572 "WARNING: No viewers for this scene handler - please create one." 1688 << G4endl; 1573 << G4endl; 1689 } 1574 } 1690 } 1575 } 1691 } 1576 } 1692 1577 1693 void G4VisManager::SetCurrentViewer (G4VViewe 1578 void G4VisManager::SetCurrentViewer (G4VViewer* pViewer) { 1694 fpViewer = pViewer; 1579 fpViewer = pViewer; 1695 if (fpViewer == nullptr) { << 1696 if (fVerbosity >= confirmations) { << 1697 G4cout << "G4VisManager::SetCurrentView << 1698 << G4endl; << 1699 } << 1700 return; << 1701 } << 1702 if (fVerbosity >= confirmations) { 1580 if (fVerbosity >= confirmations) { 1703 G4cout << "G4VisManager::SetCurrentViewer 1581 G4cout << "G4VisManager::SetCurrentViewer: viewer now " 1704 << pViewer -> GetName () 1582 << pViewer -> GetName () 1705 << G4endl; 1583 << G4endl; 1706 } 1584 } 1707 fpSceneHandler = fpViewer -> GetSceneHandle 1585 fpSceneHandler = fpViewer -> GetSceneHandler (); 1708 if (!fpSceneHandler) { 1586 if (!fpSceneHandler) { 1709 if (fVerbosity >= warnings) { 1587 if (fVerbosity >= warnings) { 1710 G4warn << << 1588 G4cout << 1711 "WARNING: No scene handler for this vie 1589 "WARNING: No scene handler for this viewer - please create one." 1712 << G4endl; 1590 << G4endl; 1713 } 1591 } 1714 return; 1592 return; 1715 } 1593 } 1716 // JA: I don't think we need this. Setview 1594 // JA: I don't think we need this. Setview will be called when needed. 1717 // fpViewer->SetView(); 1595 // fpViewer->SetView(); 1718 fpSceneHandler -> SetCurrentViewer (pViewer 1596 fpSceneHandler -> SetCurrentViewer (pViewer); 1719 fpScene = fpSceneHandler -> GetScene (); 1597 fpScene = fpSceneHandler -> GetScene (); 1720 fpGraphicsSystem = fpSceneHandler -> GetGra 1598 fpGraphicsSystem = fpSceneHandler -> GetGraphicsSystem (); 1721 if (!IsValidView ()) { 1599 if (!IsValidView ()) { 1722 if (fVerbosity >= warnings) { 1600 if (fVerbosity >= warnings) { 1723 G4warn << << 1601 G4cout << 1724 "WARNING: Problem setting viewer - please r 1602 "WARNING: Problem setting viewer - please report circumstances." 1725 << G4endl; 1603 << G4endl; 1726 } 1604 } 1727 } 1605 } 1728 } 1606 } 1729 1607 1730 void G4VisManager::PrintAvailableGraphicsSyst << 1608 void G4VisManager::PrintAvailableGraphicsSystems (Verbosity verbosity) const 1731 (Verbosity verbosity, std::ostream& out) cons << 1732 { 1609 { 1733 out << "Registered graphics systems are:\n" << 1610 G4cout << "Registered graphics systems are:\n"; 1734 if (fAvailableGraphicsSystems.size ()) { 1611 if (fAvailableGraphicsSystems.size ()) { 1735 for (const auto& gs: fAvailableGraphicsSy 1612 for (const auto& gs: fAvailableGraphicsSystems) { 1736 const G4String& name = gs->GetName(); 1613 const G4String& name = gs->GetName(); 1737 const std::vector<G4String>& nicknames 1614 const std::vector<G4String>& nicknames = gs->GetNicknames(); 1738 if (verbosity >= warnings) { << 1615 if (verbosity <= warnings) { 1739 // Brief output 1616 // Brief output 1740 out << " " << name << " ("; << 1617 G4cout << " " << name << " ("; 1741 for (std::size_t i = 0; i < nicknames << 1618 for (size_t i = 0; i < nicknames.size(); ++i) { 1742 if (i != 0) { 1619 if (i != 0) { 1743 out << ", "; << 1620 G4cout << ", "; 1744 } 1621 } 1745 out << nicknames[i]; << 1622 G4cout << nicknames[i]; 1746 } 1623 } 1747 out << ')'; << 1624 G4cout << ')'; 1748 } else { 1625 } else { 1749 // Full output 1626 // Full output 1750 out << *gs; << 1627 G4cout << *gs; 1751 } 1628 } 1752 out << std::endl; << 1629 G4cout << G4endl; 1753 } 1630 } 1754 out << << 1755 "You may choose a graphics system (driver << 1756 "\nthe command \"/vis/open\" or \"/vis/sc << 1757 "\nor you may omit the driver parameter a << 1758 "\n- by argument in the construction of G << 1759 "\n- by environment variable \"G4VIS_DEFA << 1760 "\n- by entry in \"~/.g4session\";" << 1761 "\n- by build flags." << 1762 "\n- Note: This feature is not allowed in << 1763 "\nFor further information see \"examples << 1764 "\nand \"vis.mac\"." << 1765 << std::endl; << 1766 } else { 1631 } else { 1767 out << " NONE!!! None registered - yet! << 1632 G4cout << " NONE!!! None registered - yet! Mmmmm!" << G4endl; 1768 } 1633 } 1769 } 1634 } 1770 1635 1771 void G4VisManager::PrintAvailableModels (Verb 1636 void G4VisManager::PrintAvailableModels (Verbosity verbosity) const 1772 { 1637 { 1773 { 1638 { 1774 //fpTrajDrawModelMgr->Print(G4cout); 1639 //fpTrajDrawModelMgr->Print(G4cout); 1775 G4cout << "Registered model factories:" < 1640 G4cout << "Registered model factories:" << G4endl; 1776 const std::vector<G4VModelFactory<G4VTraj 1641 const std::vector<G4VModelFactory<G4VTrajectoryModel>*>& factoryList = 1777 fpTrajDrawModelMgr->FactoryList(); 1642 fpTrajDrawModelMgr->FactoryList(); 1778 if (factoryList.empty()) G4cout << " Non 1643 if (factoryList.empty()) G4cout << " None" << G4endl; 1779 else { 1644 else { 1780 std::vector<G4VModelFactory<G4VTrajecto 1645 std::vector<G4VModelFactory<G4VTrajectoryModel>*>::const_iterator i; 1781 for (i = factoryList.begin(); i != fact 1646 for (i = factoryList.begin(); i != factoryList.end(); ++i) { 1782 (*i)->Print(G4cout); 1647 (*i)->Print(G4cout); 1783 } 1648 } 1784 } 1649 } 1785 G4cout << "\nRegistered models:" << G4end 1650 G4cout << "\nRegistered models:" << G4endl; 1786 const G4VisListManager<G4VTrajectoryModel 1651 const G4VisListManager<G4VTrajectoryModel>* listManager = 1787 fpTrajDrawModelMgr->ListManager(); 1652 fpTrajDrawModelMgr->ListManager(); 1788 const std::map<G4String, G4VTrajectoryMod 1653 const std::map<G4String, G4VTrajectoryModel*>& modelMap = 1789 listManager->Map(); 1654 listManager->Map(); 1790 if (modelMap.empty()) G4cout << " None" 1655 if (modelMap.empty()) G4cout << " None" << G4endl; 1791 else { 1656 else { 1792 std::map<G4String, G4VTrajectoryModel*> 1657 std::map<G4String, G4VTrajectoryModel*>::const_iterator i; 1793 for (i = modelMap.begin(); i != modelMa 1658 for (i = modelMap.begin(); i != modelMap.end(); ++i) { 1794 G4cout << " " << i->second->Name(); 1659 G4cout << " " << i->second->Name(); 1795 if (i->second == listManager->Current()) G4 1660 if (i->second == listManager->Current()) G4cout << " (Current)"; 1796 G4cout << G4endl; 1661 G4cout << G4endl; 1797 if (verbosity >= parameters) i->second->Pri 1662 if (verbosity >= parameters) i->second->Print(G4cout); 1798 } 1663 } 1799 } 1664 } 1800 } 1665 } 1801 1666 1802 G4cout << G4endl; 1667 G4cout << G4endl; 1803 1668 1804 { 1669 { 1805 //fpTrajFilterMgr->Print(G4cout); 1670 //fpTrajFilterMgr->Print(G4cout); 1806 G4cout << "Registered filter factories:" 1671 G4cout << "Registered filter factories:" << G4endl; 1807 const std::vector<G4VModelFactory<G4VFilt 1672 const std::vector<G4VModelFactory<G4VFilter<G4VTrajectory> >*>& 1808 factoryList = fpTrajFilterMgr->FactoryL 1673 factoryList = fpTrajFilterMgr->FactoryList(); 1809 if (factoryList.empty()) G4cout << " Non 1674 if (factoryList.empty()) G4cout << " None" << G4endl; 1810 else { 1675 else { 1811 std::vector<G4VModelFactory<G4VFilter<G 1676 std::vector<G4VModelFactory<G4VFilter<G4VTrajectory> >*>::const_iterator i; 1812 for (i = factoryList.begin(); i != fact 1677 for (i = factoryList.begin(); i != factoryList.end(); ++i) { 1813 (*i)->Print(G4cout); 1678 (*i)->Print(G4cout); 1814 } 1679 } 1815 } 1680 } 1816 1681 1817 G4cout << "\nRegistered filters:" << G4en 1682 G4cout << "\nRegistered filters:" << G4endl; 1818 const std::vector<G4VFilter<G4VTrajectory 1683 const std::vector<G4VFilter<G4VTrajectory>*>& 1819 filterList = fpTrajFilterMgr->FilterLis 1684 filterList = fpTrajFilterMgr->FilterList(); 1820 if (filterList.empty()) G4cout << " None 1685 if (filterList.empty()) G4cout << " None" << G4endl; 1821 else { 1686 else { 1822 std::vector<G4VFilter<G4VTrajectory>*>: 1687 std::vector<G4VFilter<G4VTrajectory>*>::const_iterator i; 1823 for (i = filterList.begin(); i != filte 1688 for (i = filterList.begin(); i != filterList.end(); ++i) { 1824 G4cout << " " << (*i)->GetName() << G4endl 1689 G4cout << " " << (*i)->GetName() << G4endl; 1825 if (verbosity >= parameters) (*i)->PrintAll 1690 if (verbosity >= parameters) (*i)->PrintAll(G4cout); 1826 } 1691 } 1827 } 1692 } 1828 } 1693 } 1829 } 1694 } 1830 1695 1831 void G4VisManager::PrintAvailableUserVisActio 1696 void G4VisManager::PrintAvailableUserVisActions (Verbosity) const 1832 { 1697 { 1833 G4cout << 1698 G4cout << 1834 "You have successfully registered the fol 1699 "You have successfully registered the following user vis actions." 1835 << G4endl; 1700 << G4endl; 1836 G4cout << "Run Duration User Vis Actions:"; 1701 G4cout << "Run Duration User Vis Actions:"; 1837 if (fRunDurationUserVisActions.empty()) G4c 1702 if (fRunDurationUserVisActions.empty()) G4cout << " none" << G4endl; 1838 else { 1703 else { 1839 G4cout << G4endl; 1704 G4cout << G4endl; 1840 for (std::size_t i = 0; i < fRunDurationU << 1705 for (size_t i = 0; i < fRunDurationUserVisActions.size(); i++) { 1841 const G4String& name = fRunDurationUser 1706 const G4String& name = fRunDurationUserVisActions[i].fName; 1842 G4cout << " " << name << G4endl; 1707 G4cout << " " << name << G4endl; 1843 } 1708 } 1844 } 1709 } 1845 1710 1846 G4cout << "End of Event User Vis Actions:"; 1711 G4cout << "End of Event User Vis Actions:"; 1847 if (fEndOfEventUserVisActions.empty()) G4co 1712 if (fEndOfEventUserVisActions.empty()) G4cout << " none" << G4endl; 1848 else { 1713 else { 1849 G4cout << G4endl; 1714 G4cout << G4endl; 1850 for (std::size_t i = 0; i < fEndOfEventUs << 1715 for (size_t i = 0; i < fEndOfEventUserVisActions.size(); i++) { 1851 const G4String& name = fEndOfEventUserV 1716 const G4String& name = fEndOfEventUserVisActions[i].fName; 1852 G4cout << " " << name << G4endl; 1717 G4cout << " " << name << G4endl; 1853 } 1718 } 1854 } 1719 } 1855 1720 1856 G4cout << "End of Run User Vis Actions:"; 1721 G4cout << "End of Run User Vis Actions:"; 1857 if (fEndOfRunUserVisActions.empty()) G4cout 1722 if (fEndOfRunUserVisActions.empty()) G4cout << " none" << G4endl; 1858 else { 1723 else { 1859 G4cout << G4endl; 1724 G4cout << G4endl; 1860 for (std::size_t i = 0; i < fEndOfRunUser << 1725 for (size_t i = 0; i < fEndOfRunUserVisActions.size(); i++) { 1861 const G4String& name = fEndOfRunUserVis 1726 const G4String& name = fEndOfRunUserVisActions[i].fName; 1862 G4cout << " " << name << G4endl; 1727 G4cout << " " << name << G4endl; 1863 } 1728 } 1864 } 1729 } 1865 } 1730 } 1866 1731 1867 void G4VisManager::PrintAvailableColours (Ver 1732 void G4VisManager::PrintAvailableColours (Verbosity) const { 1868 G4cout << 1733 G4cout << 1869 "Some /vis commands (optionally) take a str << 1734 "Some /vis commands (optionally) take a string to specify colour." 1870 "\nThey are also available in your C++ code << 1735 "\nAvailable colours:\n "; 1871 "\n G4Colour niceColour; // Default - whi << 1736 const std::map<G4String, G4Colour>& map = G4Colour::GetMap(); 1872 "\n G4Colour::GetColour(\"pink\", niceColo << 1737 for (std::map<G4String, G4Colour>::const_iterator i = map.begin(); 1873 "\n logical->SetVisAttributes(niceColour); << 1738 i != map.end();) { 1874 "\nSee G4Colour.hh." << 1739 G4cout << i->first; 1875 "\nAvailable colours"; << 1740 if (++i != map.end()) G4cout << ", "; 1876 for (const auto& i: G4Colour::GetMap()) { << 1877 G4cout << ", " << i.first; << 1878 } 1741 } 1879 G4cout << G4endl; 1742 G4cout << G4endl; 1880 } 1743 } 1881 1744 1882 void G4VisManager::PrintInvalidPointers () co 1745 void G4VisManager::PrintInvalidPointers () const { 1883 if (fVerbosity >= errors) { 1746 if (fVerbosity >= errors) { 1884 G4warn << "ERROR: G4VisManager::PrintInva << 1747 G4cerr << "ERROR: G4VisManager::PrintInvalidPointers:"; 1885 if (!fpGraphicsSystem) { 1748 if (!fpGraphicsSystem) { 1886 G4warn << "\n null graphics system poin << 1749 G4cerr << "\n null graphics system pointer."; 1887 } 1750 } 1888 else { 1751 else { 1889 G4warn << "\n Graphics system is " << << 1752 G4cerr << "\n Graphics system is " << fpGraphicsSystem -> GetName () 1890 << " but:"; 1753 << " but:"; 1891 if (!fpScene) 1754 if (!fpScene) 1892 G4warn << << 1755 G4cerr << 1893 "\n Null scene pointer. Use \"/vis/drawV 1756 "\n Null scene pointer. Use \"/vis/drawVolume\" or" 1894 " \"/vis/scene/create\"."; 1757 " \"/vis/scene/create\"."; 1895 if (!fpSceneHandler) 1758 if (!fpSceneHandler) 1896 G4warn << << 1759 G4cerr << 1897 "\n Null scene handler pointer. Use \"/v 1760 "\n Null scene handler pointer. Use \"/vis/open\" or" 1898 " \"/vis/sceneHandler/create\"."; 1761 " \"/vis/sceneHandler/create\"."; 1899 if (!fpViewer ) 1762 if (!fpViewer ) 1900 G4warn << << 1763 G4cerr << 1901 "\n Null viewer pointer. Use \"/vis/view 1764 "\n Null viewer pointer. Use \"/vis/viewer/create\"."; 1902 } 1765 } 1903 G4warn << G4endl; << 1766 G4cerr << G4endl; 1904 } 1767 } 1905 } 1768 } 1906 1769 >> 1770 #ifdef G4MULTITHREADED >> 1771 >> 1772 namespace { >> 1773 G4bool mtRunInProgress = false; >> 1774 std::deque<const G4Event*> mtVisEventQueue; >> 1775 G4Thread* mtVisSubThread = 0; >> 1776 G4Mutex mtVisSubThreadMutex = G4MUTEX_INITIALIZER; >> 1777 } 1907 1778 1908 G4ThreadFunReturnType G4VisManager::G4VisSubT 1779 G4ThreadFunReturnType G4VisManager::G4VisSubThread(G4ThreadFunArgType p) 1909 { 1780 { 1910 #ifdef G4MULTITHREADED << 1911 G4VisManager* pVisManager = (G4VisManager*) 1781 G4VisManager* pVisManager = (G4VisManager*)p; 1912 G4VSceneHandler* pSceneHandler = pVisManage 1782 G4VSceneHandler* pSceneHandler = pVisManager->GetCurrentSceneHandler(); 1913 if (!pSceneHandler) return 0; 1783 if (!pSceneHandler) return 0; 1914 G4Scene* pScene = pSceneHandler->GetScene() 1784 G4Scene* pScene = pSceneHandler->GetScene(); 1915 if (!pScene) return 0; 1785 if (!pScene) return 0; 1916 G4VViewer* pViewer = pVisManager->GetCurren 1786 G4VViewer* pViewer = pVisManager->GetCurrentViewer(); 1917 if (!pViewer) return 0; 1787 if (!pViewer) return 0; 1918 1788 1919 G4UImanager::GetUIpointer()->SetUpForSpecia 1789 G4UImanager::GetUIpointer()->SetUpForSpecialThread("G4VIS"); 1920 1790 >> 1791 // G4cout << "G4VisManager::G4VisSubThread: thread: " >> 1792 // << G4Threading::G4GetThreadId() << std::endl; >> 1793 1921 // Set up geometry and navigation for a thr 1794 // Set up geometry and navigation for a thread 1922 G4GeometryWorkspace::GetPool()->CreateAndUs 1795 G4GeometryWorkspace::GetPool()->CreateAndUseWorkspace(); 1923 G4SolidsWorkspace::GetPool()->CreateAndUseW 1796 G4SolidsWorkspace::GetPool()->CreateAndUseWorkspace(); 1924 G4Navigator* navigator = G4TransportationMa 1797 G4Navigator* navigator = G4TransportationManager::GetTransportationManager() 1925 ->GetNavigatorFo 1798 ->GetNavigatorForTracking(); 1926 navigator->SetWorldVolume( 1799 navigator->SetWorldVolume( 1927 G4RunManagerFactory::GetMasterRunManagerK 1800 G4RunManagerFactory::GetMasterRunManagerKernel()->GetCurrentWorld()); 1928 1801 1929 pViewer->SwitchToVisSubThread(); 1802 pViewer->SwitchToVisSubThread(); 1930 1803 1931 while (true) { 1804 while (true) { 1932 1805 1933 G4MUTEXLOCK(&mtVisSubThreadMutex); 1806 G4MUTEXLOCK(&mtVisSubThreadMutex); 1934 std::size_t eventQueueSize = mtVisEventQu << 1807 G4int eventQueueSize = mtVisEventQueue.size(); 1935 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1808 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1936 // G4cout << "Event queue size (A): " << 1809 // G4cout << "Event queue size (A): " << eventQueueSize << G4endl; 1937 1810 1938 while (eventQueueSize) { 1811 while (eventQueueSize) { 1939 1812 1940 G4MUTEXLOCK(&mtVisSubThreadMutex); 1813 G4MUTEXLOCK(&mtVisSubThreadMutex); 1941 const G4Event* event = mtVisEventQueue. 1814 const G4Event* event = mtVisEventQueue.front(); 1942 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1815 G4MUTEXUNLOCK(&mtVisSubThreadMutex); >> 1816 // G4int eventID = event->GetEventID(); >> 1817 // G4cout >> 1818 // << "G4VisManager::G4VisSubThread: Vis sub-thread: Dealing with event: >> 1819 // " >> 1820 // << eventID << G4endl; 1943 1821 1944 // Here comes the event drawing 1822 // Here comes the event drawing 1945 pVisManager->SetTransientsDrawnThisEven 1823 pVisManager->SetTransientsDrawnThisEvent(false); 1946 pSceneHandler->SetTransientsDrawnThisEv 1824 pSceneHandler->SetTransientsDrawnThisEvent(false); 1947 1825 1948 // We are about to draw the event (traj 1826 // We are about to draw the event (trajectories, etc.), but first we 1949 // have to clear the previous event(s) 1827 // have to clear the previous event(s) if necessary. If this event 1950 // needs to be drawn afresh, e.g., the 1828 // needs to be drawn afresh, e.g., the first event or any event when 1951 // "accumulate" is not requested, the o 1829 // "accumulate" is not requested, the old event has to be cleared. 1952 // We have postponed this so that, for 1830 // We have postponed this so that, for normal viewers like OGL, the 1953 // previous event(s) stay on screen unt 1831 // previous event(s) stay on screen until this new event comes 1954 // along. For a file-writing viewer th 1832 // along. For a file-writing viewer the geometry has to be drawn. 1955 // See, for example, G4HepRepFileSceneH 1833 // See, for example, G4HepRepFileSceneHandler::ClearTransientStore. 1956 pVisManager->ClearTransientStoreIfMarke 1834 pVisManager->ClearTransientStoreIfMarked(); 1957 1835 1958 // Now draw the event... 1836 // Now draw the event... 1959 pSceneHandler->DrawEvent(event); 1837 pSceneHandler->DrawEvent(event); 1960 ++pVisManager->fNoOfEventsDrawnThisRun; 1838 ++pVisManager->fNoOfEventsDrawnThisRun; 1961 1839 >> 1840 if (pScene->GetRefreshAtEndOfEvent()) { >> 1841 >> 1842 // ShowView guarantees the view is flushed to the screen. It also >> 1843 // triggers other features such picking (if enabled) and allows >> 1844 // file-writing viewers to close the file. >> 1845 pViewer->ShowView(); >> 1846 pSceneHandler->SetMarkForClearingTransientStore(true); >> 1847 >> 1848 } >> 1849 >> 1850 // Testing. >> 1851 // std::this_thread::sleep_for(std::chrono::seconds(5)); >> 1852 1962 // Then pop and release event 1853 // Then pop and release event 1963 G4MUTEXLOCK(&mtVisSubThreadMutex); 1854 G4MUTEXLOCK(&mtVisSubThreadMutex); 1964 mtVisEventQueue.pop_front(); 1855 mtVisEventQueue.pop_front(); 1965 pVisManager->EndOfEventCleanup(event); << 1856 event->PostProcessingFinished(); 1966 eventQueueSize = mtVisEventQueue.size() 1857 eventQueueSize = mtVisEventQueue.size(); 1967 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1858 G4MUTEXUNLOCK(&mtVisSubThreadMutex); >> 1859 // G4cout << "Event queue size (B): " << eventQueueSize << G4endl; 1968 } 1860 } 1969 1861 1970 G4MUTEXLOCK(&mtVisSubThreadMutex); 1862 G4MUTEXLOCK(&mtVisSubThreadMutex); 1971 G4bool runInProgress = mtRunInProgress; << 1863 G4int runInProgress = mtRunInProgress; 1972 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1864 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1973 if (!runInProgress) { 1865 if (!runInProgress) { 1974 // EndOfRun on master thread has signal 1866 // EndOfRun on master thread has signalled end of run. There is 1975 // nothing to draw so... 1867 // nothing to draw so... 1976 break; 1868 break; 1977 } 1869 } 1978 1870 1979 // Run still in progress but nothing to d 1871 // Run still in progress but nothing to draw, so wait a while. 1980 std::this_thread::sleep_for(std::chrono:: 1872 std::this_thread::sleep_for(std::chrono::milliseconds(100)); 1981 } 1873 } 1982 1874 1983 // Inform viewer that we have finished all 1875 // Inform viewer that we have finished all sub-thread drawing 1984 pViewer->DoneWithVisSubThread(); 1876 pViewer->DoneWithVisSubThread(); 1985 pViewer->MovingToMasterThread(); 1877 pViewer->MovingToMasterThread(); 1986 #else << 1878 // G4cout << "G4VisManager::G4VisSubThread: Vis sub-thread: ending" << G4endl; 1987 G4ConsumeParameters(p); << 1988 #endif << 1989 return nullptr; 1879 return nullptr; 1990 } 1880 } 1991 1881 >> 1882 namespace { >> 1883 // G4Mutex visBeginOfRunMutex = G4MUTEX_INITIALIZER; >> 1884 // G4Mutex visBeginOfEventMutex = G4MUTEX_INITIALIZER; >> 1885 G4Mutex visEndOfEventMutex = G4MUTEX_INITIALIZER; >> 1886 // G4Mutex visEndOfRunMutex = G4MUTEX_INITIALIZER; >> 1887 } >> 1888 >> 1889 #endif >> 1890 1992 void G4VisManager::BeginOfRun () 1891 void G4VisManager::BeginOfRun () 1993 { 1892 { 1994 if (fIgnoreStateChanges) return; 1893 if (fIgnoreStateChanges) return; >> 1894 >> 1895 #ifdef G4MULTITHREADED 1995 if (G4Threading::IsWorkerThread()) return; 1896 if (G4Threading::IsWorkerThread()) return; 1996 if (!GetConcreteInstance()) return; << 1897 #endif >> 1898 // G4cout << "G4VisManager::BeginOfRun: thread: " >> 1899 // << G4Threading::G4GetThreadId() << G4endl; 1997 1900 1998 // Check if view is valid << 1901 G4RunManager* runManager = G4RunManagerFactory::GetMasterRunManager(); 1999 // Protect IsValidView() with fpSceneHandle << 2000 isValidViewForRun = fpSceneHandler && IsVal << 2001 if (!isValidViewForRun) return; << 2002 1902 2003 // For a fake run... 1903 // For a fake run... 2004 G4RunManager* runManager = G4RunManagerFact << 2005 G4int nEventsToBeProcessed = runManager->Ge 1904 G4int nEventsToBeProcessed = runManager->GetNumberOfEventsToBeProcessed(); 2006 isFakeRun = (nEventsToBeProcessed == 0); << 1905 if (nEventsToBeProcessed == 0) return; 2007 if (isFakeRun) return; << 2008 << 2009 // Is the run manager a sub-event type? << 2010 G4RunManager::RMType rmType = runManager->G << 2011 isSubEventRunManagerType = << 2012 (rmType == G4RunManager::subEventMasterRM) << 2013 (rmType == G4RunManager::subEventWorkerRM); << 2014 1906 2015 fNKeepForPostProcessingRequests = 0; << 1907 fNKeepRequests = 0; 2016 fNKeepTheEventRequests = 0; << 2017 fEventKeepingSuspended = false; 1908 fEventKeepingSuspended = false; 2018 fTransientsDrawnThisRun = false; 1909 fTransientsDrawnThisRun = false; 2019 if (fpSceneHandler) fpSceneHandler->SetTran 1910 if (fpSceneHandler) fpSceneHandler->SetTransientsDrawnThisRun(false); 2020 fNoOfEventsDrawnThisRun = 0; 1911 fNoOfEventsDrawnThisRun = 0; 2021 1912 2022 // Check to see if the user has created a t 1913 // Check to see if the user has created a trajectory model. If not, create 2023 // a default one. To avoid code duplication 1914 // a default one. To avoid code duplication the following function is used 2024 // and its result (a const G4VTrajectoryMod 1915 // and its result (a const G4VTrajectoryModel*) is thrown away at this point. 2025 // The function is called again later when 1916 // The function is called again later when needed. 2026 CurrentTrajDrawModel(); 1917 CurrentTrajDrawModel(); 2027 1918 2028 // Actions only in MT mode << 1919 #ifdef G4MULTITHREADED >> 1920 // There is a static method G4Threading::IsMultithreadedApplication() >> 1921 // that returns true only if G4MTRunManager is instantiated with MT >> 1922 // installation. Thus method returns false if G4RunManager base class is >> 1923 // instantiated even with the MT installation, or of course with sequential >> 1924 // installation. 2029 if (G4Threading::IsMultithreadedApplication 1925 if (G4Threading::IsMultithreadedApplication()) { 2030 1926 2031 // Inform viewer that we have finished al 1927 // Inform viewer that we have finished all master thread drawing for now... 2032 if (fpViewer) fpViewer->DoneWithMasterThr 1928 if (fpViewer) fpViewer->DoneWithMasterThread(); 2033 1929 2034 // Start vis sub-thread 1930 // Start vis sub-thread >> 1931 // G4cout << "G4VisManager::BeginOfRun: Starting vis sub-thread" << G4endl; 2035 G4MUTEXLOCK(&mtVisSubThreadMutex); 1932 G4MUTEXLOCK(&mtVisSubThreadMutex); 2036 mtRunInProgress = true; 1933 mtRunInProgress = true; 2037 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 1934 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2038 mtVisSubThread = new G4Thread; 1935 mtVisSubThread = new G4Thread; 2039 // Launch vis thread 1936 // Launch vis thread 2040 G4THREADCREATE(mtVisSubThread,G4VisSubThr 1937 G4THREADCREATE(mtVisSubThread,G4VisSubThread,this); 2041 1938 2042 // Tricky things for some viewers (e.g., 1939 // Tricky things for some viewers (e.g., Qt): 2043 // - Launch the vis thread 1940 // - Launch the vis thread 2044 // - Wait for the vis thread to set its Q 1941 // - Wait for the vis thread to set its QThread 2045 // - Then move current QOpenGL context (i 1942 // - Then move current QOpenGL context (if Qt) to this Qthread 2046 // - Go ahead 1943 // - Go ahead 2047 if (fpViewer) fpViewer->MovingToVisSubThr 1944 if (fpViewer) fpViewer->MovingToVisSubThread(); 2048 } 1945 } >> 1946 #endif 2049 } 1947 } 2050 1948 2051 void G4VisManager::BeginOfEvent () 1949 void G4VisManager::BeginOfEvent () 2052 { 1950 { 2053 if (fIgnoreStateChanges) return; 1951 if (fIgnoreStateChanges) return; >> 1952 2054 if (!GetConcreteInstance()) return; 1953 if (!GetConcreteInstance()) return; 2055 if (!isValidViewForRun) return; << 1954 2056 if (isFakeRun) return; << 1955 // G4cout << "G4VisManager::BeginOfEvent: thread: " >> 1956 // << G4Threading::G4GetThreadId() << G4endl; 2057 1957 2058 // Some instructions that should NOT be in 1958 // Some instructions that should NOT be in multithreaded version. 2059 // TODO: protect with G4Threading::IsMultit << 2060 #ifndef G4MULTITHREADED 1959 #ifndef G4MULTITHREADED 2061 // These instructions are in G4VisSubThread 1960 // These instructions are in G4VisSubThread for multithreading. 2062 fTransientsDrawnThisEvent = false; 1961 fTransientsDrawnThisEvent = false; 2063 if (fpSceneHandler) fpSceneHandler->SetTran 1962 if (fpSceneHandler) fpSceneHandler->SetTransientsDrawnThisEvent(false); 2064 #endif 1963 #endif 2065 } 1964 } 2066 1965 2067 // Here begins a sequence of functions that d << 2068 // Sequential/Serial mode: << 2069 // EndOfEvent is invoked by a state change << 2070 // EndOfEvent pulls the event from the Ev << 2071 // EndOfEventKernel draws the event and c << 2072 // (Unless the user him/herself has reque << 2073 // also sets KeepTheEvent for a certain << 2074 // of events. The run manager keeps thes << 2075 // run, so, for example, the vis manager << 2076 // Multithreading/Tasking: << 2077 // EndOfEvent is invoked by a state change << 2078 // EndOfEvent pulls the event from the (w << 2079 // EndOfEventKernel pushes it to the vis << 2080 // are marked KeepForPostProcessing so << 2081 // PostProcessingFinished in EndOfEvent << 2082 // EndOfEventKernel also sets KeepTheEven << 2083 // The vis sub thread pulls events from t << 2084 // calls EndOfEventCleanup, which calls << 2085 // Sub-event parallelism: << 2086 // EndOfEvent is invoked by a state change << 2087 // multithreading, but actually, the even << 2088 // employed on sub-events, so we ignore s << 2089 // G4SubEvtRunManager knows when the even << 2090 // Events are pushed to the Vis Manager by << 2091 // when the event is ready. EventReadyFo << 2092 // processsing continues as for Multithr << 2093 << 2094 void G4VisManager::EndOfEvent () 1966 void G4VisManager::EndOfEvent () 2095 { 1967 { 2096 if (fIgnoreStateChanges) return; 1968 if (fIgnoreStateChanges) return; >> 1969 2097 if (!GetConcreteInstance()) return; 1970 if (!GetConcreteInstance()) return; 2098 if (!isValidViewForRun) return; << 2099 if (isFakeRun) return; << 2100 1971 2101 // For sub-event parallelism, this is not t << 1972 // Don't call IsValidView unless there is a scene handler. This 2102 if (isSubEventRunManagerType) return; << 1973 // avoids WARNING message at end of event and run when the user has >> 1974 // not instantiated a scene handler, e.g., in batch mode. >> 1975 G4bool valid = fpSceneHandler && IsValidView(); >> 1976 if (!valid) return; >> 1977 >> 1978 // G4cout << "G4VisManager::EndOfEvent: thread: " >> 1979 // << G4Threading::G4GetThreadId() << G4endl; 2103 1980 >> 1981 #ifdef G4MULTITHREADED 2104 G4AutoLock al(&visEndOfEventMutex); 1982 G4AutoLock al(&visEndOfEventMutex); >> 1983 // Testing. >> 1984 // std::this_thread::sleep_for(std::chrono::seconds(5)); >> 1985 #endif 2105 1986 2106 G4RunManager* runManager = G4RunManagerFact 1987 G4RunManager* runManager = G4RunManagerFactory::GetMasterRunManager(); 2107 1988 2108 const G4Run* currentRun = runManager->GetCu 1989 const G4Run* currentRun = runManager->GetCurrentRun(); 2109 if (!currentRun) return; 1990 if (!currentRun) return; 2110 1991 2111 // This gets the appropriate event manager << 1992 // This gets the thread-local event manager 2112 G4EventManager* eventManager = G4EventManag 1993 G4EventManager* eventManager = G4EventManager::GetEventManager(); 2113 const G4Event* currentEvent = eventManager- 1994 const G4Event* currentEvent = eventManager->GetConstCurrentEvent(); 2114 if (!currentEvent) return; 1995 if (!currentEvent) return; 2115 1996 2116 // Call kernel << 2117 EndOfEventKernel(currentEvent); << 2118 } << 2119 << 2120 void G4VisManager::EventReadyForVis(const G4E << 2121 // This is invoked by G4SubEvtRunManager. << 2122 // The event is passed to EndOfEventKernel. << 2123 { << 2124 if (fIgnoreStateChanges) return; << 2125 if (!GetConcreteInstance()) return; << 2126 if (!isValidViewForRun) return; << 2127 if (isFakeRun) return; << 2128 << 2129 G4AutoLock al(&visEndOfEventMutex); << 2130 EndOfEventKernel(event); << 2131 } << 2132 << 2133 void G4VisManager::EndOfEventKernel (const G4 << 2134 { << 2135 // Note: we are still subject to the mutex << 2136 << 2137 // Discard event if fDrawEventOnlyIfToBeKep 1997 // Discard event if fDrawEventOnlyIfToBeKept flag is set unless the 2138 // user has requested the event to be kept. 1998 // user has requested the event to be kept. 2139 if (fDrawEventOnlyIfToBeKept) { 1999 if (fDrawEventOnlyIfToBeKept) { 2140 if (!currentEvent->KeepTheEventFlag()) re << 2000 if (!currentEvent->ToBeKept()) return; 2141 } 2001 } 2142 2002 2143 if (G4Threading::IsMultithreadedApplication 2003 if (G4Threading::IsMultithreadedApplication()) { 2144 2004 >> 2005 #ifdef G4MULTITHREADED >> 2006 2145 // Wait if too many events in the queue. 2007 // Wait if too many events in the queue. 2146 G4MUTEXLOCK(&mtVisSubThreadMutex); 2008 G4MUTEXLOCK(&mtVisSubThreadMutex); 2147 std::size_t eventQueueSize = mtVisEventQu << 2009 G4int eventQueueSize = mtVisEventQueue.size(); 2148 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2010 G4MUTEXUNLOCK(&mtVisSubThreadMutex); >> 2011 // G4cout << "Event queue size (1): " << eventQueueSize << G4endl; 2149 2012 2150 G4bool eventQueueFull = false; 2013 G4bool eventQueueFull = false; 2151 while (fMaxEventQueueSize > 0 && (G4int)e << 2014 while (fMaxEventQueueSize > 0 && eventQueueSize >= fMaxEventQueueSize) { 2152 2015 >> 2016 // G4cout << "Event queue size (2): " << eventQueueSize << G4endl; 2153 if (fWaitOnEventQueueFull) { 2017 if (fWaitOnEventQueueFull) { 2154 << 2155 static G4bool warned = false; 2018 static G4bool warned = false; 2156 if (!warned) { 2019 if (!warned) { 2157 G4warn << << 2020 G4cout << 2158 "WARNING: The number of events in t 2021 "WARNING: The number of events in the visualisation queue has exceeded" 2159 "\n the maximum, " 2022 "\n the maximum, " 2160 << fMaxEventQueueSize << 2023 << fMaxEventQueueSize << 2161 ".\n If, during a multithreaded ru 2024 ".\n If, during a multithreaded run, the simulation gets ahead of the" 2162 "\n visualisation by more than thi 2025 "\n visualisation by more than this maximum, the simulation is delayed" 2163 "\n until the vis sub-thread has d 2026 "\n until the vis sub-thread has drawn a few more events and removed them" 2164 "\n from the queue. You may chang 2027 "\n from the queue. You may change this maximum number of events with" 2165 "\n \"/vis/multithreading/maxEvent 2028 "\n \"/vis/multithreading/maxEventQueueSize <N>\", where N is the maximum" 2166 "\n number you wish to allow. N < 2029 "\n number you wish to allow. N <= 0 means \"unlimited\"." 2167 "\n Alternatively you may choose t 2030 "\n Alternatively you may choose to discard events for drawing by setting" 2168 "\n \"/vis/multithreading/actionOn 2031 "\n \"/vis/multithreading/actionOnEventQueueFull discard\"." 2169 "\n To avoid visualisation altoget 2032 "\n To avoid visualisation altogether: \"/vis/disable\"." 2170 "\n And maybe \"/tracking/storeTra 2033 "\n And maybe \"/tracking/storeTrajectories 0\"." 2171 << G4endl; 2034 << G4endl; 2172 warned = true; 2035 warned = true; 2173 } 2036 } 2174 << 2037 // G4cout << "Event queue size (3): " << eventQueueSize << G4endl; 2175 // Wait a while to give event drawing 2038 // Wait a while to give event drawing time to reduce the queue... 2176 std::this_thread::sleep_for(std::chro 2039 std::this_thread::sleep_for(std::chrono::milliseconds(100)); 2177 << 2040 // G4cout << "Event queue size (4): " << eventQueueSize << G4endl; 2178 } else { 2041 } else { 2179 << 2180 static G4bool warned = false; 2042 static G4bool warned = false; 2181 if (!warned) { 2043 if (!warned) { 2182 G4warn << << 2044 G4cout << 2183 "WARNING: The number of events in t 2045 "WARNING: The number of events in the visualisation queue has exceeded" 2184 "\n the maximum, " 2046 "\n the maximum, " 2185 << fMaxEventQueueSize << 2047 << fMaxEventQueueSize << 2186 ".\n Some events have been discard 2048 ".\n Some events have been discarded for drawing. You may change this" 2187 "\n behaviour with \"/vis/multithr 2049 "\n behaviour with \"/vis/multithreading/actionOnEventQueueFull wait\"." 2188 "\n To avoid visualisation altoget 2050 "\n To avoid visualisation altogether: \"/vis/disable\"." 2189 "\n And maybe \"/tracking/storeTra 2051 "\n And maybe \"/tracking/storeTrajectories 0\"." 2190 << G4endl; 2052 << G4endl; 2191 warned = true; 2053 warned = true; 2192 } 2054 } 2193 << 2194 eventQueueFull = true; // Causes eve 2055 eventQueueFull = true; // Causes event to be discarded for drawing. 2195 break; 2056 break; 2196 } 2057 } 2197 2058 2198 G4MUTEXLOCK(&mtVisSubThreadMutex); 2059 G4MUTEXLOCK(&mtVisSubThreadMutex); 2199 eventQueueSize = mtVisEventQueue.size() 2060 eventQueueSize = mtVisEventQueue.size(); 2200 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2061 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2201 } 2062 } 2202 2063 2203 if (!eventQueueFull) { 2064 if (!eventQueueFull) { 2204 << 2205 if (RequiredToBeKeptForVis(currentEvent << 2206 currentEvent->KeepTheEvent(); << 2207 fNKeepTheEventRequests++; // Counts << 2208 } << 2209 << 2210 G4MUTEXLOCK(&mtVisSubThreadMutex); 2065 G4MUTEXLOCK(&mtVisSubThreadMutex); 2211 << 2066 // Keep event for processing and put event on vis event queue 2212 // Keep while post processing (i.e., dr << 2213 // Make sure the event is not deleted b << 2214 currentEvent->KeepForPostProcessing(); 2067 currentEvent->KeepForPostProcessing(); 2215 fNKeepForPostProcessingRequests++; << 2068 if (fpScene->GetRefreshAtEndOfEvent()) { 2216 << 2069 // Keep one event (cannot know which is last so any will do) 2217 // Put event on vis queue << 2070 if (fNKeepRequests == 0) { >> 2071 eventManager->KeepTheCurrentEvent(); >> 2072 fNKeepRequests++; >> 2073 } >> 2074 } 2218 mtVisEventQueue.push_back(currentEvent) 2075 mtVisEventQueue.push_back(currentEvent); 2219 << 2220 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2076 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2221 } 2077 } 2222 2078 >> 2079 // G4MUTEXLOCK(&mtVisSubThreadMutex); >> 2080 // G4int eQS = mtVisEventQueue.size(); >> 2081 // G4MUTEXUNLOCK(&mtVisSubThreadMutex); >> 2082 // G4cout << "Event queue size (5): " << eQS << G4endl; >> 2083 >> 2084 #endif >> 2085 2223 } else { 2086 } else { 2224 2087 2225 // Sequential mode 2088 // Sequential mode 2226 2089 >> 2090 G4int nEventsToBeProcessed = 0; >> 2091 G4int nKeptEvents = 0; >> 2092 G4int eventID = -2; // (If no run manager, triggers ShowView as normal.) >> 2093 if (currentRun) { >> 2094 nEventsToBeProcessed = runManager->GetNumberOfEventsToBeProcessed(); >> 2095 eventID = currentEvent->GetEventID(); >> 2096 const std::vector<const G4Event*>* events = currentRun->GetEventVector(); >> 2097 if (events) nKeptEvents = events->size(); >> 2098 } >> 2099 2227 // We are about to draw the event (trajec 2100 // We are about to draw the event (trajectories, etc.), but first we 2228 // have to clear the previous event(s) if 2101 // have to clear the previous event(s) if necessary. If this event 2229 // needs to be drawn afresh, e.g., the fi 2102 // needs to be drawn afresh, e.g., the first event or any event when 2230 // "accumulate" is not requested, the old 2103 // "accumulate" is not requested, the old event has to be cleared. 2231 // We have postponed this so that, for no 2104 // We have postponed this so that, for normal viewers like OGL, the 2232 // previous event(s) stay on screen until 2105 // previous event(s) stay on screen until this new event comes 2233 // along. For a file-writing viewer the 2106 // along. For a file-writing viewer the geometry has to be drawn. 2234 // See, for example, G4HepRepFileSceneHan 2107 // See, for example, G4HepRepFileSceneHandler::ClearTransientStore. 2235 ClearTransientStoreIfMarked(); 2108 ClearTransientStoreIfMarked(); 2236 2109 2237 // Keep while post processing (i.e., draw << 2238 currentEvent->KeepForPostProcessing(); << 2239 fNKeepForPostProcessingRequests++; << 2240 << 2241 if (RequiredToBeKeptForVis(currentEvent-> << 2242 currentEvent->KeepTheEvent(); << 2243 fNKeepTheEventRequests++; // Counts nu << 2244 } << 2245 << 2246 // Now draw the event... 2110 // Now draw the event... 2247 fpSceneHandler->DrawEvent(currentEvent); 2111 fpSceneHandler->DrawEvent(currentEvent); 2248 ++fNoOfEventsDrawnThisRun; 2112 ++fNoOfEventsDrawnThisRun; 2249 2113 2250 EndOfEventCleanup(currentEvent); << 2114 if (fpScene->GetRefreshAtEndOfEvent()) { 2251 } << 2252 } << 2253 << 2254 void G4VisManager::EndOfEventCleanup(const G4 << 2255 { << 2256 if (fpScene->GetRefreshAtEndOfEvent()) { << 2257 2115 2258 fpSceneHandler->SetMarkForClearingTransie << 2116 // Unless last event (in which case wait end of run)... >> 2117 if (eventID < nEventsToBeProcessed - 1) { >> 2118 // ShowView guarantees the view is flushed to the screen. It also >> 2119 // triggers other features such picking (if enabled) and allows >> 2120 // file-writing viewers to close the file. >> 2121 fpViewer->ShowView(); >> 2122 } else { // Last event... >> 2123 // Keep, but only if user has not kept any... >> 2124 if (nKeptEvents == 0) { >> 2125 eventManager->KeepTheCurrentEvent(); >> 2126 fNKeepRequests++; >> 2127 } >> 2128 } >> 2129 fpSceneHandler->SetMarkForClearingTransientStore(true); 2259 2130 2260 // ShowView guarantees the view is flushe << 2131 } 2261 // triggers other features such as pickin << 2262 // file-writing viewers to close the file << 2263 fpViewer->ShowView(); << 2264 } 2132 } 2265 2133 2266 currentEvent->PostProcessingFinished(); << 2134 // Both modes - sequential and MT 2267 } << 2268 2135 2269 G4bool G4VisManager::RequiredToBeKeptForVis ( << 2136 if (!(fpScene->GetRefreshAtEndOfEvent())) { 2270 // i.e., kept by the run manager at least to << 2271 { << 2272 G4bool requiredToBeKept = false; << 2273 << 2274 G4RunManager* runManager = G4RunManagerFact << 2275 G4int nEventsToBeProcessed = runManager->Ge << 2276 2137 2277 if (fpScene->GetRefreshAtEndOfEvent()) { / << 2138 // Accumulating events... 2278 << 2279 // Keep last event only << 2280 if (eventID == nEventsToBeProcessed - 1) << 2281 requiredToBeKept = true; << 2282 } << 2283 << 2284 } else { // Accumulating events << 2285 2139 2286 G4int maxNumberOfKeptEvents = fpScene->Ge 2140 G4int maxNumberOfKeptEvents = fpScene->GetMaxNumberOfKeptEvents(); 2287 2141 2288 if (maxNumberOfKeptEvents >= 0) { // Eve << 2142 if (maxNumberOfKeptEvents > 0 && 2289 << 2143 fNKeepRequests >= maxNumberOfKeptEvents) { 2290 if (fNKeepTheEventRequests < maxNumberO << 2291 2144 2292 // Not yet reached the limit << 2145 fEventKeepingSuspended = true; 2293 requiredToBeKept = true; << 2146 static G4bool warned = false; 2294 << 2147 if (!warned) { 2295 } else { << 2148 if (fVerbosity >= warnings) { 2296 << 2149 G4cout << 2297 // We have already reached the limit << 2150 "WARNING: G4VisManager::EndOfEvent: Automatic event keeping suspended." 2298 fEventKeepingSuspended = true; << 2151 "\n The number of events exceeds the maximum, " 2299 static G4bool warned = false; << 2152 << maxNumberOfKeptEvents << 2300 if (!warned) { << 2153 ", that may be kept by\n the vis manager." 2301 if (fVerbosity >= warnings) { << 2154 << G4endl; 2302 G4warn << << 2303 "WARNING: G4VisManager::EndOfEven << 2304 << G4endl; << 2305 if (maxNumberOfKeptEvents > 0) { << 2306 G4warn << << 2307 "\n The number of events excee << 2308 << maxNumberOfKeptEvents << << 2309 ", that may be kept by\n the v << 2310 << G4endl; << 2311 } << 2312 } << 2313 warned = true; << 2314 } 2155 } >> 2156 warned = true; 2315 } 2157 } 2316 2158 2317 } else { // Indefinite event keeping (ma << 2159 } else if (maxNumberOfKeptEvents != 0) { 2318 << 2319 requiredToBeKept = true; << 2320 2160 >> 2161 // If not disabled nor suspended. >> 2162 if (GetConcreteInstance() && !fEventKeepingSuspended) { >> 2163 // G4cout << >> 2164 // "Requesting keeping event " << currentEvent->GetEventID() >> 2165 // << G4endl; >> 2166 eventManager->KeepTheCurrentEvent(); >> 2167 fNKeepRequests++; >> 2168 } >> 2169 2321 } 2170 } 2322 } 2171 } 2323 << 2324 return requiredToBeKept; << 2325 } 2172 } 2326 2173 2327 void G4VisManager::EndOfRun () 2174 void G4VisManager::EndOfRun () 2328 { 2175 { 2329 if (fIgnoreStateChanges) return; 2176 if (fIgnoreStateChanges) return; >> 2177 >> 2178 #ifdef G4MULTITHREADED 2330 if (G4Threading::IsWorkerThread()) return; 2179 if (G4Threading::IsWorkerThread()) return; 2331 if (!GetConcreteInstance()) return; << 2180 #endif 2332 if (!isValidViewForRun) return; << 2181 2333 if (isFakeRun) return; << 2182 // G4cout << "G4VisManager::EndOfRun: thread: " >> 2183 // << G4Threading::G4GetThreadId() << G4endl; 2334 2184 2335 G4RunManager* runManager = G4RunManagerFact 2185 G4RunManager* runManager = G4RunManagerFactory::GetMasterRunManager(); 2336 2186 2337 // For a fake run... 2187 // For a fake run... 2338 G4int nEventsToBeProcessed = runManager->Ge 2188 G4int nEventsToBeProcessed = runManager->GetNumberOfEventsToBeProcessed(); 2339 if (nEventsToBeProcessed == 0) return; 2189 if (nEventsToBeProcessed == 0) return; 2340 2190 2341 const G4Run* currentRun = runManager->GetCu 2191 const G4Run* currentRun = runManager->GetCurrentRun(); 2342 if (!currentRun) return; 2192 if (!currentRun) return; 2343 2193 >> 2194 #ifdef G4MULTITHREADED 2344 // G4AutoLock al(&visEndOfRunMutex); ??? 2195 // G4AutoLock al(&visEndOfRunMutex); ??? 2345 if (G4Threading::IsMultithreadedApplication 2196 if (G4Threading::IsMultithreadedApplication()) { 2346 // Reset flag so that sub-thread exits wh 2197 // Reset flag so that sub-thread exits when it has finished processing. 2347 G4MUTEXLOCK(&mtVisSubThreadMutex); 2198 G4MUTEXLOCK(&mtVisSubThreadMutex); 2348 mtRunInProgress = false; 2199 mtRunInProgress = false; 2349 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2200 G4MUTEXUNLOCK(&mtVisSubThreadMutex); 2350 // Wait for sub-thread to finish. 2201 // Wait for sub-thread to finish. 2351 G4THREADJOIN(*mtVisSubThread); 2202 G4THREADJOIN(*mtVisSubThread); 2352 delete mtVisSubThread; 2203 delete mtVisSubThread; 2353 fpViewer->SwitchToMasterThread(); << 2204 if (fpViewer) fpViewer->SwitchToMasterThread(); 2354 } 2205 } >> 2206 #endif 2355 2207 2356 if (G4Threading::IsMultithreadedApplication << 2208 #ifdef G4MULTITHREADED >> 2209 // Print warning about discarded events, if any. >> 2210 // Don't call IsValidView unless there is a scene handler. This >> 2211 // avoids WARNING message from IsValidView() when the user has >> 2212 // not instantiated a scene handler, e.g., in batch mode. >> 2213 if (fpSceneHandler && IsValidView()) { // Events should have been drawn 2357 G4int noOfEventsRequested = runManager->G 2214 G4int noOfEventsRequested = runManager->GetNumberOfEventsToBeProcessed(); 2358 if (fNoOfEventsDrawnThisRun != noOfEvents 2215 if (fNoOfEventsDrawnThisRun != noOfEventsRequested) { 2359 if (!fWaitOnEventQueueFull && fVerbosit 2216 if (!fWaitOnEventQueueFull && fVerbosity >= warnings) { 2360 G4warn << 2217 G4cout 2361 << "WARNING: Number of events drawn t 2218 << "WARNING: Number of events drawn this run, " 2362 << fNoOfEventsDrawnThisRun << ", is d 2219 << fNoOfEventsDrawnThisRun << ", is different to number requested, " 2363 << noOfEventsRequested << 2220 << noOfEventsRequested << 2364 ".\n (This is because you requested 2221 ".\n (This is because you requested \"/vis/multithreading/actionOnEventQueueFull discard\".)" 2365 << G4endl; 2222 << G4endl; 2366 } 2223 } 2367 } 2224 } 2368 } 2225 } >> 2226 #endif 2369 2227 2370 G4int nKeptEvents = currentRun->GetNumberOf << 2228 G4int nKeptEvents = 0; 2371 if (fVerbosity >= warnings && nKeptEvents > << 2229 const std::vector<const G4Event*>* events = currentRun->GetEventVector(); 2372 G4warn << nKeptEvents; << 2230 if (events) nKeptEvents = events->size(); 2373 if (nKeptEvents == 1) G4warn << " event h << 2231 if (fVerbosity >= warnings) { 2374 else G4warn << " events have"; << 2232 G4cout << nKeptEvents; 2375 G4warn << " been kept for refreshing and/ << 2233 if (nKeptEvents == 1) G4cout << " event has"; 2376 if (nKeptEvents != fNKeepTheEventRequests << 2234 else G4cout << " events have"; 2377 if (fNKeepTheEventRequests == 0) { << 2235 G4cout << " been kept for refreshing and/or reviewing." << G4endl; 2378 G4warn << "No keep requests were"; << 2236 if (nKeptEvents != fNKeepRequests) { 2379 } else if (fNKeepTheEventRequests == 1) << 2237 G4cout << " (Note: "; 2380 G4warn << "1 keep request was"; << 2238 if (fNKeepRequests == 0) { >> 2239 G4cout << "No keep requests were"; >> 2240 } else if (fNKeepRequests == 1) { >> 2241 G4cout << "1 keep request was"; 2381 } else { 2242 } else { 2382 G4warn << fNKeepTheEventRequests << " << 2243 G4cout << fNKeepRequests << " keep requests were"; 2383 } 2244 } 2384 G4warn << " made by the vis manager."; << 2245 G4cout << " made by the vis manager."; 2385 if (fNKeepTheEventRequests == 0) { << 2246 if (fNKeepRequests == 0) { 2386 G4warn << << 2247 G4cout << 2387 "\n The kept events are those you ha << 2248 "\n The kept events are those you have asked to be kept in your user action(s).)"; 2388 } else { 2249 } else { 2389 G4warn << << 2250 G4cout << 2390 "\n The same or further events may h << 2251 "\n The same or further events may have been kept by you in your user action(s).)"; 2391 "\n To turn off event keeping by the << 2252 } 2392 "\n or use /vis/scene/endOfEventActi << 2253 G4cout << G4endl; 2393 } << 2254 } 2394 G4warn << G4endl; << 2255 G4cout << 2395 } << 2256 " \"/vis/reviewKeptEvents\" to review them one by one." 2396 G4warn << << 2257 "\n \"/vis/enable\", then \"/vis/viewer/flush\" or \"/vis/viewer/rebuild\" to see them accumulated." 2397 " \"/vis/reviewKeptEvents\" to review one << 2398 "\n To see accumulated, \"/vis/enable\", t << 2399 << G4endl; 2258 << G4endl; 2400 } 2259 } 2401 2260 2402 if (fVerbosity >= warnings) PrintListOfPlot << 2261 // static G4bool warned = false; >> 2262 // if (!valid && fVerbosity >= warnings && !warned) { >> 2263 // G4cout << >> 2264 // " Only useful if before starting the run:" >> 2265 // "\n a) trajectories are stored (\"/vis/scene/add/trajectories [smooth|rich]\"), or" >> 2266 // "\n b) the Draw method of any hits or digis is implemented." >> 2267 // "\n To view trajectories, hits or digis:" >> 2268 // "\n open a viewer, draw a volume, \"/vis/scene/add/trajectories\"" >> 2269 // "\n \"/vis/scene/add/hits\" or \"/vis/scene/add/digitisations\"" >> 2270 // "\n and, possibly, \"/vis/viewer/flush\"." >> 2271 // "\n To see all events: \"/vis/scene/endOfEventAction accumulate\"." >> 2272 // "\n (You may need \"/vis/viewer/flush\" or even \"/vis/viewer/rebuild\".)" >> 2273 // "\n To see events individually: \"/vis/reviewKeptEvents\"." >> 2274 // << G4endl; >> 2275 // warned = true; >> 2276 // } 2403 2277 2404 if (fEventKeepingSuspended && fVerbosity >= 2278 if (fEventKeepingSuspended && fVerbosity >= warnings) { 2405 G4warn << << 2279 G4cout << 2406 "WARNING: G4VisManager::EndOfRun: Automat 2280 "WARNING: G4VisManager::EndOfRun: Automatic event keeping was suspended." >> 2281 "\n The number of events in the run exceeded the maximum, " >> 2282 << fpScene->GetMaxNumberOfKeptEvents() << >> 2283 ", that may be\n kept by the vis manager." << >> 2284 "\n The number of events kept by the vis manager can be changed with" >> 2285 "\n \"/vis/scene/endOfEventAction accumulate <N>\", where N is the" >> 2286 "\n maximum number you wish to allow. N < 0 means \"unlimited\"." 2407 << G4endl; 2287 << G4endl; 2408 if (fpScene->GetMaxNumberOfKeptEvents() > << 2409 G4warn << << 2410 " The number of events in the run exce << 2411 << fpScene->GetMaxNumberOfKeptEvents() << 2412 ", that may be\n kept by the vis manag << 2413 "\n The number of events kept by the v << 2414 "\n \"/vis/scene/endOfEventAction accu << 2415 "\n maximum number you wish to allow. << 2416 << G4endl; << 2417 } << 2418 } 2288 } 2419 2289 >> 2290 // Don't call IsValidView unless there is a scene handler. This >> 2291 // avoids WARNING message at end of event and run when the user has >> 2292 // not instantiated a scene handler, e.g., in batch mode. >> 2293 G4bool valid = fpSceneHandler && IsValidView(); >> 2294 if (GetConcreteInstance() && valid) { >> 2295 // // ???? I can't remember why 2420 // // if (!fpSceneHandler->GetMarkForClear 2296 // // if (!fpSceneHandler->GetMarkForClearingTransientStore()) { 2421 // // is here. It prevents ShowView at en 2297 // // is here. It prevents ShowView at end of run, which seems to be OK 2422 // // for sequential mode, but MT mode see 2298 // // for sequential mode, but MT mode seems to need it (I have not 2423 // // figured out why). ???? JA ???? 2299 // // figured out why). ???? JA ???? 2424 // if (!fpSceneHandler->GetMarkForClearing 2300 // if (!fpSceneHandler->GetMarkForClearingTransientStore()) { 2425 if (fpScene->GetRefreshAtEndOfRun()) { << 2301 if (fpScene->GetRefreshAtEndOfRun()) { 2426 fpSceneHandler->DrawEndOfRunModels(); << 2302 fpSceneHandler->DrawEndOfRunModels(); 2427 // An extra refresh for auto-refresh view << 2303 // An extra refresh for auto-refresh viewers. 2428 // ???? I DON'T WHY THIS IS NECESSARY ??? << 2304 // ???? I DON'T WHY THIS IS NECESSARY ???? JA ???? 2429 if (fpViewer->GetViewParameters().IsAutoR << 2305 if (fpViewer->GetViewParameters().IsAutoRefresh()) { 2430 fpViewer->RefreshView(); << 2306 fpViewer->RefreshView(); 2431 } << 2307 } 2432 // ShowView guarantees the view is flushe << 2308 // ShowView guarantees the view is flushed to the screen. It also 2433 // triggers other features such picking ( << 2309 // triggers other features such picking (if enabled) and allows 2434 // file-writing viewers to close the file << 2310 // file-writing viewers to close the file. 2435 fpViewer->ShowView(); << 2311 fpViewer->ShowView(); 2436 fpSceneHandler->SetMarkForClearingTransie << 2312 fpSceneHandler->SetMarkForClearingTransientStore(true); 2437 } else { << 2313 } else { 2438 if (fpGraphicsSystem->GetFunctionality() << 2314 if (fpGraphicsSystem->GetFunctionality() == 2439 G4VGraphicsSystem::fileWriter) { << 2315 G4VGraphicsSystem::fileWriter) { 2440 if (fVerbosity >= warnings) { << 2316 if (fVerbosity >= warnings) { 2441 G4warn << "\"/vis/viewer/update\" to << 2317 G4cout << "\"/vis/viewer/update\" to close file." << G4endl; >> 2318 } >> 2319 } 2442 } 2320 } 2443 } << 2321 // } 2444 } 2322 } 2445 //} << 2446 fEventRefreshing = false; 2323 fEventRefreshing = false; 2447 } 2324 } 2448 2325 2449 void G4VisManager::ClearTransientStoreIfMarke 2326 void G4VisManager::ClearTransientStoreIfMarked(){ 2450 // Assumes valid view. 2327 // Assumes valid view. 2451 if (fpSceneHandler->GetMarkForClearingTrans 2328 if (fpSceneHandler->GetMarkForClearingTransientStore()) { 2452 fpSceneHandler->SetMarkForClearingTransie 2329 fpSceneHandler->SetMarkForClearingTransientStore(false); 2453 fpSceneHandler->ClearTransientStore(); 2330 fpSceneHandler->ClearTransientStore(); 2454 } 2331 } 2455 // Record if transients drawn. These local 2332 // Record if transients drawn. These local flags are only set 2456 // *after* ClearTransientStore. In the cod 2333 // *after* ClearTransientStore. In the code in G4VSceneHandler 2457 // triggered by ClearTransientStore, use th 2334 // triggered by ClearTransientStore, use these flags so that 2458 // event refreshing is not done too early. 2335 // event refreshing is not done too early. 2459 fTransientsDrawnThisEvent = fpSceneHandler- 2336 fTransientsDrawnThisEvent = fpSceneHandler->GetTransientsDrawnThisEvent(); 2460 fTransientsDrawnThisRun = fpSceneHandler->G 2337 fTransientsDrawnThisRun = fpSceneHandler->GetTransientsDrawnThisRun(); 2461 } 2338 } 2462 2339 2463 void G4VisManager::ResetTransientsDrawnFlags( 2340 void G4VisManager::ResetTransientsDrawnFlags() 2464 { 2341 { 2465 fTransientsDrawnThisRun = false; 2342 fTransientsDrawnThisRun = false; 2466 fTransientsDrawnThisEvent = false; 2343 fTransientsDrawnThisEvent = false; 2467 G4SceneHandlerListConstIterator i; 2344 G4SceneHandlerListConstIterator i; 2468 for (i = fAvailableSceneHandlers.begin(); 2345 for (i = fAvailableSceneHandlers.begin(); 2469 i != fAvailableSceneHandlers.end(); ++ 2346 i != fAvailableSceneHandlers.end(); ++i) { 2470 (*i)->SetTransientsDrawnThisEvent(false); 2347 (*i)->SetTransientsDrawnThisEvent(false); 2471 (*i)->SetTransientsDrawnThisRun(false); 2348 (*i)->SetTransientsDrawnThisRun(false); 2472 } 2349 } 2473 } 2350 } 2474 2351 2475 G4String G4VisManager::ViewerShortName (const 2352 G4String G4VisManager::ViewerShortName (const G4String& viewerName) const { 2476 const G4String& viewerShortName = viewerNam << 2353 G4String viewerShortName (viewerName); 2477 return G4StrUtil::strip_copy(viewerShortNam << 2354 viewerShortName = viewerShortName (0, viewerShortName.find (' ')); >> 2355 return viewerShortName.strip (); 2478 } 2356 } 2479 2357 2480 G4VViewer* G4VisManager::GetViewer (const G4S 2358 G4VViewer* G4VisManager::GetViewer (const G4String& viewerName) const { 2481 G4String viewerShortName = ViewerShortName 2359 G4String viewerShortName = ViewerShortName (viewerName); 2482 std::size_t nHandlers = fAvailableSceneHand << 2360 size_t nHandlers = fAvailableSceneHandlers.size (); 2483 std::size_t iHandler, iViewer; << 2361 size_t iHandler, iViewer; 2484 G4VViewer* viewer = 0; 2362 G4VViewer* viewer = 0; 2485 G4bool found = false; 2363 G4bool found = false; 2486 for (iHandler = 0; iHandler < nHandlers; iH 2364 for (iHandler = 0; iHandler < nHandlers; iHandler++) { 2487 G4VSceneHandler* sceneHandler = fAvailabl 2365 G4VSceneHandler* sceneHandler = fAvailableSceneHandlers [iHandler]; 2488 const G4ViewerList& viewerList = sceneHan 2366 const G4ViewerList& viewerList = sceneHandler -> GetViewerList (); 2489 for (iViewer = 0; iViewer < viewerList.si 2367 for (iViewer = 0; iViewer < viewerList.size (); iViewer++) { 2490 viewer = viewerList [iViewer]; 2368 viewer = viewerList [iViewer]; 2491 if (viewerShortName == viewer -> GetSho 2369 if (viewerShortName == viewer -> GetShortName ()) { 2492 found = true; 2370 found = true; 2493 break; 2371 break; 2494 } 2372 } 2495 } 2373 } 2496 if (found) break; 2374 if (found) break; 2497 } 2375 } 2498 if (found) return viewer; 2376 if (found) return viewer; 2499 else return 0; 2377 else return 0; 2500 } 2378 } 2501 2379 2502 std::vector<G4String> G4VisManager::Verbosity 2380 std::vector<G4String> G4VisManager::VerbosityGuidanceStrings; 2503 2381 2504 G4String G4VisManager::VerbosityString(Verbos 2382 G4String G4VisManager::VerbosityString(Verbosity verbosity) { 2505 G4String rs; 2383 G4String rs; 2506 switch (verbosity) { 2384 switch (verbosity) { 2507 case quiet: rs = "quiet (0)"; break 2385 case quiet: rs = "quiet (0)"; break; 2508 case startup: rs = "startup (1)"; bre 2386 case startup: rs = "startup (1)"; break; 2509 case errors: rs = "errors (2)"; brea 2387 case errors: rs = "errors (2)"; break; 2510 case warnings: rs = "warnings (3)"; br 2388 case warnings: rs = "warnings (3)"; break; 2511 case confirmations: rs = "confirmations (4) 2389 case confirmations: rs = "confirmations (4)"; break; 2512 case parameters: rs = "parameters (5)"; 2390 case parameters: rs = "parameters (5)"; break; 2513 case all: rs = "all (6)"; break; 2391 case all: rs = "all (6)"; break; 2514 } 2392 } 2515 return rs; 2393 return rs; 2516 } 2394 } 2517 2395 2518 void G4VisManager::PrintAvailableVerbosity(st << 2519 { << 2520 os << "Available verbosity options:"; << 2521 for (std::size_t i = 0; i < VerbosityGuidan << 2522 os << '\n' << VerbosityGuidanceStrings[i] << 2523 } << 2524 os << "\nCurrent verbosity: " << G4VisManag << 2525 os << std::endl; << 2526 } << 2527 << 2528 G4VisManager::Verbosity 2396 G4VisManager::Verbosity 2529 G4VisManager::GetVerbosityValue(const G4Strin 2397 G4VisManager::GetVerbosityValue(const G4String& verbosityString) { 2530 G4String ss = G4StrUtil::to_lower_copy(verb << 2398 G4String ss(verbosityString); ss.toLower(); 2531 Verbosity verbosity; 2399 Verbosity verbosity; 2532 if (ss[0] == 'q') verbosity = quiet; << 2400 if (ss(0) == 'q') verbosity = quiet; 2533 else if (ss[0] == 's') verbosity = startup; << 2401 else if (ss(0) == 's') verbosity = startup; 2534 else if (ss[0] == 'e') verbosity = errors; << 2402 else if (ss(0) == 'e') verbosity = errors; 2535 else if (ss[0] == 'w') verbosity = warnings << 2403 else if (ss(0) == 'w') verbosity = warnings; 2536 else if (ss[0] == 'c') verbosity = confirma << 2404 else if (ss(0) == 'c') verbosity = confirmations; 2537 else if (ss[0] == 'p') verbosity = paramete << 2405 else if (ss(0) == 'p') verbosity = parameters; 2538 else if (ss[0] == 'a') verbosity = all; << 2406 else if (ss(0) == 'a') verbosity = all; 2539 else { 2407 else { 2540 // Could be an integer << 2541 G4int intVerbosity; 2408 G4int intVerbosity; 2542 std::istringstream is(ss); 2409 std::istringstream is(ss); 2543 is >> intVerbosity; 2410 is >> intVerbosity; 2544 if (!is) { 2411 if (!is) { 2545 G4warn << "ERROR: G4VisManager::GetVerb << 2412 G4cerr << "ERROR: G4VisManager::GetVerbosityValue: invalid verbosity \"" 2546 << verbosityString << "\"\n"; << 2413 << verbosityString << "\""; 2547 PrintAvailableVerbosity(G4warn); << 2414 for (size_t i = 0; i < VerbosityGuidanceStrings.size(); ++i) { 2548 // Return existing verbosity << 2415 G4cerr << '\n' << VerbosityGuidanceStrings[i]; 2549 return fVerbosity; << 2416 } >> 2417 verbosity = warnings; >> 2418 G4cerr << "\n Returning " << VerbosityString(verbosity) >> 2419 << G4endl; 2550 } 2420 } 2551 else { 2421 else { 2552 verbosity = GetVerbosityValue(intVerbos 2422 verbosity = GetVerbosityValue(intVerbosity); 2553 } 2423 } 2554 } 2424 } 2555 return verbosity; 2425 return verbosity; 2556 } 2426 } 2557 2427 2558 G4VisManager::Verbosity G4VisManager::GetVerb 2428 G4VisManager::Verbosity G4VisManager::GetVerbosityValue(G4int intVerbosity) { 2559 Verbosity verbosity; 2429 Verbosity verbosity; 2560 if (intVerbosity < quiet) verbosity = 2430 if (intVerbosity < quiet) verbosity = quiet; 2561 else if (intVerbosity > all) verbosity = 2431 else if (intVerbosity > all) verbosity = all; 2562 else verbosity = 2432 else verbosity = Verbosity(intVerbosity); 2563 return verbosity; 2433 return verbosity; 2564 } 2434 } 2565 2435 2566 G4VisManager::Verbosity G4VisManager::GetVerb 2436 G4VisManager::Verbosity G4VisManager::GetVerbosity () { 2567 return fVerbosity; 2437 return fVerbosity; 2568 } 2438 } 2569 2439 2570 void G4VisManager::SetVerboseLevel (G4int int 2440 void G4VisManager::SetVerboseLevel (G4int intVerbosity) { 2571 fVerbosity = GetVerbosityValue(intVerbosity 2441 fVerbosity = GetVerbosityValue(intVerbosity); 2572 } 2442 } 2573 2443 2574 void G4VisManager::SetVerboseLevel (const G4S 2444 void G4VisManager::SetVerboseLevel (const G4String& verbosityString) { 2575 fVerbosity = GetVerbosityValue(verbosityStr 2445 fVerbosity = GetVerbosityValue(verbosityString); 2576 } 2446 } 2577 2447 2578 G4bool G4VisManager::IsValidView () { 2448 G4bool G4VisManager::IsValidView () { 2579 2449 2580 if (!fInitialised) Initialise (); 2450 if (!fInitialised) Initialise (); 2581 2451 2582 static G4bool noGSPrinting = true; 2452 static G4bool noGSPrinting = true; 2583 if (!fpGraphicsSystem) { 2453 if (!fpGraphicsSystem) { 2584 // Limit printing - we do not want printi 2454 // Limit printing - we do not want printing if the user simply does 2585 // not want to use graphics, e.g., in bat 2455 // not want to use graphics, e.g., in batch mode. 2586 if (noGSPrinting) { 2456 if (noGSPrinting) { 2587 noGSPrinting = false; 2457 noGSPrinting = false; 2588 if (fVerbosity >= warnings) { 2458 if (fVerbosity >= warnings) { 2589 G4warn << << 2459 G4cout << 2590 "WARNING: G4VisManager::IsValidView(): Atte 2460 "WARNING: G4VisManager::IsValidView(): Attempt to draw when no graphics system" 2591 "\n has been instantiated. Use \"/vis/ope 2461 "\n has been instantiated. Use \"/vis/open\" or \"/vis/sceneHandler/create\"." 2592 "\n Alternatively, to avoid this message, 2462 "\n Alternatively, to avoid this message, suppress instantiation of vis" 2593 "\n manager (G4VisExecutive) and ensure dr 2463 "\n manager (G4VisExecutive) and ensure drawing code is executed only if" 2594 "\n G4VVisManager::GetConcreteInstance() i 2464 "\n G4VVisManager::GetConcreteInstance() is non-zero." 2595 << G4endl; 2465 << G4endl; 2596 } 2466 } 2597 } 2467 } 2598 return false; 2468 return false; 2599 } 2469 } 2600 2470 2601 if ((!fpScene) || (!fpSceneHandler) || (!fp 2471 if ((!fpScene) || (!fpSceneHandler) || (!fpViewer)) { 2602 if (fVerbosity >= errors) { 2472 if (fVerbosity >= errors) { 2603 G4warn << << 2473 G4cerr << 2604 "ERROR: G4VisManager::IsValidView(): Curren 2474 "ERROR: G4VisManager::IsValidView(): Current view is not valid." 2605 << G4endl; 2475 << G4endl; 2606 PrintInvalidPointers (); 2476 PrintInvalidPointers (); 2607 } 2477 } 2608 return false; 2478 return false; 2609 } 2479 } 2610 2480 2611 if (fpScene != fpSceneHandler -> GetScene ( 2481 if (fpScene != fpSceneHandler -> GetScene ()) { 2612 if (fVerbosity >= errors) { 2482 if (fVerbosity >= errors) { 2613 G4warn << "ERROR: G4VisManager::IsValid << 2483 G4cerr << "ERROR: G4VisManager::IsValidView ():"; 2614 if (fpSceneHandler -> GetScene ()) { 2484 if (fpSceneHandler -> GetScene ()) { 2615 G4warn << << 2485 G4cout << 2616 "\n The current scene \"" 2486 "\n The current scene \"" 2617 << fpScene -> GetName () 2487 << fpScene -> GetName () 2618 << "\" is not handled by" 2488 << "\" is not handled by" 2619 "\n the current scene handler \"" 2489 "\n the current scene handler \"" 2620 << fpSceneHandler -> GetName () 2490 << fpSceneHandler -> GetName () 2621 << "\"" 2491 << "\"" 2622 "\n (it currently handles scene \"" 2492 "\n (it currently handles scene \"" 2623 << fpSceneHandler -> GetScene () -> 2493 << fpSceneHandler -> GetScene () -> GetName () 2624 << "\")." 2494 << "\")." 2625 "\n Either:" 2495 "\n Either:" 2626 "\n (a) attach it to the scene handler w 2496 "\n (a) attach it to the scene handler with" 2627 "\n /vis/sceneHandler/attach " 2497 "\n /vis/sceneHandler/attach " 2628 << fpScene -> GetName () 2498 << fpScene -> GetName () 2629 << ", or" 2499 << ", or" 2630 "\n (b) create a new scene handler with 2500 "\n (b) create a new scene handler with " 2631 "\n /vis/sceneHandler/create <graphi 2501 "\n /vis/sceneHandler/create <graphics-system>," 2632 "\n in which case it should pick up 2502 "\n in which case it should pick up the the new scene." 2633 << G4endl; 2503 << G4endl; 2634 } 2504 } 2635 else { 2505 else { 2636 G4warn << "\n Scene handler \"" << 2506 G4cout << "\n Scene handler \"" 2637 << fpSceneHandler -> GetName () 2507 << fpSceneHandler -> GetName () 2638 << "\" has null scene pointer." 2508 << "\" has null scene pointer." 2639 "\n Attach a scene with /vis/sceneHandle 2509 "\n Attach a scene with /vis/sceneHandler/attach [<scene-name>]" 2640 << G4endl; 2510 << G4endl; 2641 } 2511 } 2642 } 2512 } 2643 return false; 2513 return false; 2644 } 2514 } 2645 2515 2646 const G4ViewerList& viewerList = fpSceneHan 2516 const G4ViewerList& viewerList = fpSceneHandler -> GetViewerList (); 2647 if (viewerList.size () == 0) { 2517 if (viewerList.size () == 0) { 2648 if (fVerbosity >= errors) { 2518 if (fVerbosity >= errors) { 2649 G4warn << << 2519 G4cerr << 2650 "ERROR: G4VisManager::IsValidView (): the c 2520 "ERROR: G4VisManager::IsValidView (): the current scene handler\n \"" 2651 << fpSceneHandler -> GetName () 2521 << fpSceneHandler -> GetName () 2652 << "\" has no viewers. Do /vis/viewer 2522 << "\" has no viewers. Do /vis/viewer/create." 2653 << G4endl; 2523 << G4endl; 2654 } 2524 } 2655 return false; 2525 return false; 2656 } 2526 } 2657 2527 2658 G4bool isValid = true; 2528 G4bool isValid = true; 2659 if (fpScene -> IsEmpty ()) { // Add world 2529 if (fpScene -> IsEmpty ()) { // Add world by default if possible... 2660 G4bool warn(fVerbosity >= warnings); 2530 G4bool warn(fVerbosity >= warnings); 2661 G4bool successful = fpScene -> AddWorldIf 2531 G4bool successful = fpScene -> AddWorldIfEmpty (warn); 2662 if (!successful || fpScene -> IsEmpty ()) 2532 if (!successful || fpScene -> IsEmpty ()) { // If still empty... 2663 if (fVerbosity >= errors) { 2533 if (fVerbosity >= errors) { 2664 G4warn << "ERROR: G4VisManager::IsValidView << 2534 G4cerr << "ERROR: G4VisManager::IsValidView ():"; 2665 G4warn << << 2535 G4cerr << 2666 "\n Attempt at some drawing operation wh 2536 "\n Attempt at some drawing operation when scene is empty." 2667 "\n Maybe the geometry has not yet been 2537 "\n Maybe the geometry has not yet been defined." 2668 " Try /run/initialize." 2538 " Try /run/initialize." 2669 "\n Or use \"/vis/scene/add/extent 2539 "\n Or use \"/vis/scene/add/extent\"." 2670 << G4endl; 2540 << G4endl; 2671 } 2541 } 2672 isValid = false; 2542 isValid = false; 2673 } 2543 } 2674 else { 2544 else { 2675 G4UImanager::GetUIpointer()->ApplyComma 2545 G4UImanager::GetUIpointer()->ApplyCommand ("/vis/scene/notifyHandlers"); 2676 if (fVerbosity >= warnings) { 2546 if (fVerbosity >= warnings) { 2677 G4warn << << 2547 G4cout << 2678 "WARNING: G4VisManager: the scene was emp 2548 "WARNING: G4VisManager: the scene was empty, \"world\" has been" 2679 "\n added and the scene handlers notifie 2549 "\n added and the scene handlers notified."; 2680 G4warn << G4endl; << 2550 G4cout << G4endl; 2681 } 2551 } 2682 } 2552 } 2683 } 2553 } 2684 return isValid; 2554 return isValid; 2685 } 2555 } 2686 2556 2687 void 2557 void 2688 G4VisManager::RegisterModelFactories() 2558 G4VisManager::RegisterModelFactories() 2689 { 2559 { 2690 if (fVerbosity >= warnings) { 2560 if (fVerbosity >= warnings) { 2691 G4warn<<"G4VisManager: No model factories << 2561 G4cout<<"G4VisManager: No model factories registered with G4VisManager."<<G4endl; 2692 G4warn<<"G4VisManager::RegisterModelFacto << 2562 G4cout<<"G4VisManager::RegisterModelFactories() should be overridden in derived"<<G4endl; 2693 G4warn<<"class. See G4VisExecutive for an << 2563 G4cout<<"class. See G4VisExecutive for an example."<<G4endl; 2694 } 2564 } 2695 } 2565 } 2696 2566 >> 2567 #ifdef G4MULTITHREADED 2697 void G4VisManager::SetUpForAThread() 2568 void G4VisManager::SetUpForAThread() 2698 { 2569 { 2699 // TODO: protect with G4Threading::IsMultit << 2700 #ifdef G4MULTITHREADED << 2701 new G4VisStateDependent(this); 2570 new G4VisStateDependent(this); 2702 #endif << 2703 } 2571 } >> 2572 #endif 2704 2573 2705 void G4VisManager::IgnoreStateChanges(G4bool 2574 void G4VisManager::IgnoreStateChanges(G4bool val) 2706 { 2575 { 2707 fIgnoreStateChanges = val; 2576 fIgnoreStateChanges = val; 2708 } 2577 } 2709 2578