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 /// \file medical/DICOM/DICOM.cc 27 /// \file medical/DICOM/DICOM.cc 28 /// \brief Main program of the medical/DICOM e 28 /// \brief Main program of the medical/DICOM example 29 // 29 // 30 // The code was written by : 30 // The code was written by : 31 // *Louis Archambault louis.archambault 31 // *Louis Archambault louis.archambault@phy.ulaval.ca, 32 // *Luc Beaulieu beaulieu@phy.ulaval.ca 32 // *Luc Beaulieu beaulieu@phy.ulaval.ca 33 // +Vincent Hubert-Tremblay at tigre.2@sy 33 // +Vincent Hubert-Tremblay at tigre.2@sympatico.ca 34 // 34 // 35 // 35 // 36 // *Centre Hospitalier Universitaire de Quebec 36 // *Centre Hospitalier Universitaire de Quebec (CHUQ), 37 // Hotel-Dieu de Quebec, departement de Radio- 37 // Hotel-Dieu de Quebec, departement de Radio-oncologie 38 // 11 cote du palais. Quebec, QC, Canada, G1R 38 // 11 cote du palais. Quebec, QC, Canada, G1R 2J6 39 // tel (418) 525-4444 #6720 39 // tel (418) 525-4444 #6720 40 // fax (418) 691 5268 40 // fax (418) 691 5268 41 // 41 // 42 // + Université Laval, Québec (QC) Canada << 42 // + Universit� Laval, Qu�bec (QC) Canada 43 // ******************************************* 43 // ******************************************************* 44 #include "DicomActionInitialization.hh" << 45 #include "DicomNestedParamDetectorConstruction << 46 #include "DicomPartialDetectorConstruction.hh" << 47 #include "DicomRegularDetectorConstruction.hh" << 48 << 49 #include "G4GenericPhysicsList.hh" << 50 #include "G4RunManagerFactory.hh" << 51 #include "G4Types.hh" 44 #include "G4Types.hh" >> 45 >> 46 #ifdef G4MULTITHREADED >> 47 #include "G4MTRunManager.hh" >> 48 #else >> 49 #include "G4RunManager.hh" >> 50 #endif >> 51 >> 52 #include "globals.hh" 52 #include "G4UImanager.hh" 53 #include "G4UImanager.hh" 53 #include "Randomize.hh" 54 #include "Randomize.hh" 54 #include "globals.hh" << 55 >> 56 #include "G4GenericPhysicsList.hh" >> 57 >> 58 #include "DicomRegularDetectorConstruction.hh" >> 59 #include "DicomNestedParamDetectorConstruction.hh" >> 60 #include "DicomPartialDetectorConstruction.hh" >> 61 >> 62 #include "DicomActionInitialization.hh" 55 63 56 #ifdef G4_DCMTK 64 #ifdef G4_DCMTK 57 # include "DicomFileMgr.hh" << 65 # include "DicomFileMgr.hh" 58 #else 66 #else 59 # include "DicomHandler.hh" << 67 # include "DicomHandler.hh" 60 #endif 68 #endif 61 69 62 #include "DicomIntersectVolume.hh" 70 #include "DicomIntersectVolume.hh" 63 #include "QGSP_BIC.hh" 71 #include "QGSP_BIC.hh" 64 #include "Shielding.hh" << 72 #include "G4tgrMessenger.hh" 65 73 66 #include "G4UIExecutive.hh" << 67 #include "G4VisExecutive.hh" 74 #include "G4VisExecutive.hh" 68 #include "G4tgrMessenger.hh" << 75 #include "G4UIExecutive.hh" >> 76 >> 77 #include "Shielding.hh" 69 78 70 //....oooOO0OOooo........oooOO0OOooo........oo 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 71 80 72 int main(int argc, char** argv) << 81 int main(int argc,char** argv) 73 { 82 { 74 // Instantiate G4UIExecutive if interactive 83 // Instantiate G4UIExecutive if interactive mode 75 G4UIExecutive* ui = nullptr; 84 G4UIExecutive* ui = nullptr; 76 if (argc == 1) { << 85 if ( argc == 1 ) { 77 ui = new G4UIExecutive(argc, argv); 86 ui = new G4UIExecutive(argc, argv); 78 } 87 } 79 88 80 new G4tgrMessenger; 89 new G4tgrMessenger; 81 char* part = std::getenv("DICOM_PARTIAL_PARA << 90 char* part = std::getenv( "DICOM_PARTIAL_PARAM" ); 82 G4bool bPartial = FALSE; 91 G4bool bPartial = FALSE; 83 if (part && G4String(part) == "1") { << 92 if( part && G4String(part) == "1" ) { 84 bPartial = TRUE; 93 bPartial = TRUE; 85 } 94 } 86 95 87 CLHEP::HepRandom::setTheEngine(new CLHEP::Ra 96 CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine); 88 CLHEP::HepRandom::setTheSeed(G4long(24534575 97 CLHEP::HepRandom::setTheSeed(G4long(24534575684783)); 89 G4long seeds[2]; 98 G4long seeds[2]; 90 seeds[0] = G4long(534524575674523); 99 seeds[0] = G4long(534524575674523); 91 seeds[1] = G4long(526345623452457); 100 seeds[1] = G4long(526345623452457); 92 CLHEP::HepRandom::setTheSeeds(seeds); 101 CLHEP::HepRandom::setTheSeeds(seeds); 93 102 94 // Construct the default run manager 103 // Construct the default run manager >> 104 #ifdef G4MULTITHREADED 95 char* nthread_c = std::getenv("DICOM_NTHREAD 105 char* nthread_c = std::getenv("DICOM_NTHREADS"); 96 106 97 unsigned nthreads = 4; 107 unsigned nthreads = 4; 98 unsigned env_threads = 0; 108 unsigned env_threads = 0; 99 109 100 if (nthread_c) { << 110 if(nthread_c) {env_threads=unsigned(G4UIcommand::ConvertToDouble(nthread_c));} 101 env_threads = unsigned(G4UIcommand::Conver << 111 if(env_threads > 0) {nthreads=env_threads;} 102 } << 103 if (env_threads > 0) { << 104 nthreads = env_threads; << 105 } << 106 112 107 auto* runManager = G4RunManagerFactory::Crea << 113 G4MTRunManager* runManager = new G4MTRunManager; 108 runManager->SetNumberOfThreads(nthreads); 114 runManager->SetNumberOfThreads(nthreads); 109 115 >> 116 G4cout << "\n\n\tDICOM running in multithreaded mode with " << nthreads >> 117 << " threads\n\n" << G4endl; >> 118 >> 119 >> 120 #else >> 121 G4RunManager* runManager = new G4RunManager; >> 122 G4cout << "\n\n\tDICOM running in serial mode\n\n" << G4endl; >> 123 >> 124 #endif >> 125 110 DicomDetectorConstruction* theGeometry = 0; 126 DicomDetectorConstruction* theGeometry = 0; 111 127 112 #ifdef G4_DCMTK 128 #ifdef G4_DCMTK 113 DicomFileMgr* theFileMgr = 0; 129 DicomFileMgr* theFileMgr = 0; 114 #else 130 #else 115 DicomHandler* dcmHandler = 0; 131 DicomHandler* dcmHandler = 0; 116 #endif 132 #endif 117 133 118 if (!bPartial) { << 134 if( !bPartial ){ 119 #ifdef G4_DCMTK 135 #ifdef G4_DCMTK 120 136 121 theFileMgr = DicomFileMgr::GetInstance(); 137 theFileMgr = DicomFileMgr::GetInstance(); 122 theFileMgr->Convert("Data.dat"); 138 theFileMgr->Convert("Data.dat"); 123 139 124 #else 140 #else 125 // Treatment of DICOM images before creati 141 // Treatment of DICOM images before creating the G4runManager 126 dcmHandler = DicomHandler::Instance(); 142 dcmHandler = DicomHandler::Instance(); 127 dcmHandler->CheckFileFormat(); 143 dcmHandler->CheckFileFormat(); 128 #endif 144 #endif 129 145 130 // Initialisation of physics, geometry, pr 146 // Initialisation of physics, geometry, primary particles ... 131 char* nest = std::getenv("DICOM_NESTED_PAR << 147 char* nest = std::getenv( "DICOM_NESTED_PARAM" ); 132 if (nest && G4String(nest) == "1") { << 148 if( nest && G4String(nest) == "1" ) { 133 theGeometry = new DicomNestedParamDetect 149 theGeometry = new DicomNestedParamDetectorConstruction(); 134 } << 150 } else { 135 else { << 136 theGeometry = new DicomRegularDetectorCo 151 theGeometry = new DicomRegularDetectorConstruction(); 137 } 152 } 138 } << 153 } else { 139 else { << 140 theGeometry = new DicomPartialDetectorCons 154 theGeometry = new DicomPartialDetectorConstruction(); 141 } 155 } 142 runManager->SetUserInitialization(theGeometr 156 runManager->SetUserInitialization(theGeometry); 143 157 144 // std::vector<G4String>* MyConstr = new 158 // std::vector<G4String>* MyConstr = new std::vector<G4String>; 145 // MyConstr->push_back("G4EmStandardPhysi 159 // MyConstr->push_back("G4EmStandardPhysics"); 146 // G4VModularPhysicsList* phys = new G4Ge 160 // G4VModularPhysicsList* phys = new G4GenericPhysicsList(MyConstr); 147 G4VModularPhysicsList* phys = new Shielding( 161 G4VModularPhysicsList* phys = new Shielding(); 148 runManager->SetUserInitialization(phys); 162 runManager->SetUserInitialization(phys); 149 163 150 // Set user action classes 164 // Set user action classes 151 runManager->SetUserInitialization(new DicomA 165 runManager->SetUserInitialization(new DicomActionInitialization()); 152 166 153 runManager->Initialize(); 167 runManager->Initialize(); 154 168 155 new DicomIntersectVolume(); 169 new DicomIntersectVolume(); 156 170 157 // visualisation manager 171 // visualisation manager 158 G4VisManager* visManager = new G4VisExecutiv 172 G4VisManager* visManager = new G4VisExecutive; 159 visManager->Initialize(); 173 visManager->Initialize(); 160 174 161 G4UImanager* UImanager = G4UImanager::GetUIp 175 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 162 176 163 if (ui) { << 177 if (ui) 164 UImanager->ApplyCommand("/control/execute << 178 { 165 ui->SessionStart(); << 179 UImanager->ApplyCommand("/control/execute vis.mac"); 166 delete ui; << 180 ui->SessionStart(); 167 } << 181 delete ui; 168 else { << 182 } 169 G4String command = "/control/execute "; << 183 else 170 G4String fileName = argv[1]; << 184 { 171 UImanager->ApplyCommand(command + fileName << 185 G4String command = "/control/execute "; 172 } << 186 G4String fileName = argv[1]; >> 187 UImanager->ApplyCommand(command+fileName); >> 188 } 173 189 174 delete visManager; 190 delete visManager; 175 delete runManager; 191 delete runManager; 176 192 177 if (!bPartial) { << 193 if( !bPartial ) { 178 #ifdef G4_DCMTK 194 #ifdef G4_DCMTK 179 delete theFileMgr; 195 delete theFileMgr; 180 #endif 196 #endif 181 } 197 } 182 198 183 return 0; 199 return 0; 184 } 200 } >> 201 185 202