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 /// \file persistency/gdml/G04/gdml_det.cc 26 /// \file persistency/gdml/G04/gdml_det.cc 27 /// \brief Main program of the persistency/gdm 27 /// \brief Main program of the persistency/gdml/G04 example 28 // 28 // 29 // 29 // >> 30 // $Id$ 30 // 31 // 31 // 32 // 32 // ------------------------------------------- 33 // -------------------------------------------------------------- 33 // GEANT 4 - gdml_det 34 // GEANT 4 - gdml_det 34 // 35 // 35 // ------------------------------------------- 36 // -------------------------------------------------------------- 36 37 37 #include "FTFP_BERT.hh" << 38 #include <vector> 38 #include "G04ActionInitialization.hh" << 39 >> 40 #include "G4RunManager.hh" >> 41 #include "G4UImanager.hh" >> 42 >> 43 #include "G4LogicalVolumeStore.hh" >> 44 #include "G4TransportationManager.hh" >> 45 #include "G4SDManager.hh" >> 46 >> 47 #include "G04PrimaryGeneratorAction.hh" 39 #include "G04DetectorConstruction.hh" 48 #include "G04DetectorConstruction.hh" >> 49 #include "G04PhysicsList.hh" 40 #include "G04SensitiveDetector.hh" 50 #include "G04SensitiveDetector.hh" 41 51 42 #include "G4GDMLParser.hh" << 52 #ifdef G4VIS_USE 43 #include "G4RunManagerFactory.hh" << 44 #include "G4UIExecutive.hh" << 45 #include "G4UImanager.hh" << 46 #include "G4VisExecutive.hh" 53 #include "G4VisExecutive.hh" >> 54 #endif 47 55 48 #include <vector> << 56 #ifdef G4UI_USE >> 57 #include "G4UIExecutive.hh" >> 58 #endif >> 59 >> 60 #include "G4GDMLParser.hh" 49 61 50 int main(int argc, char** argv) << 62 int main(int argc,char **argv) 51 { 63 { 52 G4cout << G4endl; << 64 G4cout << G4endl; 53 G4cout << "Usage: gdml_det <intput_gdml_file << 65 G4cout << "Usage: gdml_det <intput_gdml_file:mandatory>" 54 G4cout << G4endl; << 66 << G4endl; 55 << 67 G4cout << G4endl; 56 if (argc < 2) { << 68 57 G4cout << "Error! Mandatory input file is << 69 if (argc<2) 58 G4cout << G4endl; << 70 { 59 return -1; << 71 G4cout << "Error! Mandatory input file is not specified!" << G4endl; 60 } << 72 G4cout << G4endl; 61 << 73 return -1; 62 // Detect interactive mode (if only one argu << 74 } 63 // << 75 64 G4UIExecutive* ui = 0; << 76 G4GDMLParser parser; 65 if (argc == 2) { << 77 parser.Read(argv[1]); 66 ui = new G4UIExecutive(argc, argv); << 78 67 } << 79 68 << 80 G4RunManager* runManager = new G4RunManager; 69 G4GDMLParser parser; << 81 70 parser.Read(argv[1]); << 82 runManager->SetUserInitialization(new G04DetectorConstruction( 71 << 83 parser.GetWorldVolume())); 72 auto* runManager = G4RunManagerFactory::Crea << 84 runManager->SetUserInitialization(new G04PhysicsList); 73 << 85 runManager->SetUserAction(new G04PrimaryGeneratorAction); 74 runManager->SetUserInitialization(new G04Det << 86 75 runManager->SetUserInitialization(new FTFP_B << 87 runManager->Initialize(); 76 << 88 77 // User action initialization << 89 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 78 runManager->SetUserInitialization(new G04Act << 90 79 runManager->Initialize(); << 91 //------------------------------------------------ 80 << 92 // Sensitive detectors 81 // Initialize visualization << 93 //------------------------------------------------ 82 G4VisManager* visManager = new G4VisExecutiv << 94 83 visManager->Initialize(); << 95 G4SDManager* SDman = G4SDManager::GetSDMpointer(); 84 << 96 85 // Get the pointer to the User Interface man << 97 G4String trackerChamberSDname = "Tracker"; 86 G4UImanager* UImanager = G4UImanager::GetUIp << 98 G04SensitiveDetector* aTrackerSD = new G04SensitiveDetector(trackerChamberSDname); 87 << 99 SDman->AddNewDetector( aTrackerSD ); 88 // Process macro or start UI session << 100 89 if (!ui) // batch mode << 101 /////////////////////////////////////////////////////////////////////// 90 { << 102 // 91 G4String command = "/control/execute "; << 103 // Example how to retrieve Auxiliary Information for sensitive detector 92 G4String fileName = argv[2]; << 104 // 93 UImanager->ApplyCommand(command + fileName << 105 const G4GDMLAuxMapType* auxmap = parser.GetAuxMap(); 94 } << 106 std::cout << "Found " << auxmap->size() 95 else // interactive mode << 107 << " volume(s) with auxiliary information." 96 { << 108 << G4endl << G4endl; 97 UImanager->ApplyCommand("/control/execute << 109 for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin(); 98 ui->SessionStart(); << 110 iter!=auxmap->end(); iter++) 99 delete ui; << 111 { 100 } << 112 G4cout << "Volume " << ((*iter).first)->GetName() >> 113 << " has the following list of auxiliary information: " >> 114 << G4endl << G4endl; >> 115 for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin(); >> 116 vit!=(*iter).second.end(); vit++) >> 117 { >> 118 std::cout << "--> Type: " << (*vit).type >> 119 << " Value: " << (*vit).value << std::endl; >> 120 } >> 121 } >> 122 G4cout << G4endl; >> 123 >> 124 // The same as above, but now we are looking for >> 125 // sensitive detectors setting them for the volumes >> 126 >> 127 for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin(); >> 128 iter!=auxmap->end(); iter++) >> 129 { >> 130 G4cout << "Volume " << ((*iter).first)->GetName() >> 131 << " has the following list of auxiliary information: " >> 132 << G4endl << G4endl; >> 133 for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin(); >> 134 vit!=(*iter).second.end();vit++) >> 135 { >> 136 if ((*vit).type=="SensDet") >> 137 { >> 138 G4cout << "Attaching sensitive detector " << (*vit).value >> 139 << " to volume " << ((*iter).first)->GetName() >> 140 << G4endl << G4endl; >> 141 >> 142 G4VSensitiveDetector* mydet = SDman->FindSensitiveDetector((*vit).value); >> 143 if(mydet) >> 144 { >> 145 G4LogicalVolume* myvol = (*iter).first; >> 146 myvol->SetSensitiveDetector(mydet); >> 147 } >> 148 else >> 149 { >> 150 G4cout << (*vit).value << " detector not found" << G4endl; >> 151 } >> 152 } >> 153 } >> 154 } >> 155 // >> 156 // End of Auxiliary Information block >> 157 // >> 158 //////////////////////////////////////////////////////////////////////// >> 159 >> 160 #ifdef G4VIS_USE >> 161 G4VisManager* visManager = new G4VisExecutive; >> 162 visManager->Initialize(); >> 163 #endif >> 164 >> 165 if (argc==3) // batch mode >> 166 { >> 167 G4String command = "/control/execute "; >> 168 G4String fileName = argv[2]; >> 169 UImanager->ApplyCommand(command+fileName); >> 170 } >> 171 else // interactive mode >> 172 { >> 173 #ifdef G4UI_USE >> 174 G4UIExecutive* ui = new G4UIExecutive(argc, argv); >> 175 #ifdef G4VIS_USE >> 176 UImanager->ApplyCommand("/control/execute vis.mac"); >> 177 #endif >> 178 ui->SessionStart(); >> 179 delete ui; >> 180 #endif >> 181 } >> 182 >> 183 #ifdef G4VIS_USE >> 184 delete visManager; >> 185 #endif >> 186 delete runManager; 101 187 102 delete visManager; << 188 return 0; 103 delete runManager; << 104 } 189 } 105 190