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: gdml_det.cc 82285 2014-06-13 14:56:02Z gcosmo $ 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 "FTFP_BERT.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 FTFP_BERT); 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 = 87 << 99 new G04SensitiveDetector(trackerChamberSDname); 88 // Process macro or start UI session << 100 SDman->AddNewDetector( aTrackerSD ); 89 if (!ui) // batch mode << 101 90 { << 102 /////////////////////////////////////////////////////////////////////// 91 G4String command = "/control/execute "; << 103 // 92 G4String fileName = argv[2]; << 104 // Example how to retrieve Auxiliary Information for sensitive detector 93 UImanager->ApplyCommand(command + fileName << 105 // 94 } << 106 const G4GDMLAuxMapType* auxmap = parser.GetAuxMap(); 95 else // interactive mode << 107 G4cout << "Found " << auxmap->size() 96 { << 108 << " volume(s) with auxiliary information." 97 UImanager->ApplyCommand("/control/execute << 109 << G4endl << G4endl; 98 ui->SessionStart(); << 110 for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin(); 99 delete ui; << 111 iter!=auxmap->end(); iter++) 100 } << 112 { >> 113 G4cout << "Volume " << ((*iter).first)->GetName() >> 114 << " has the following list of auxiliary information: " >> 115 << G4endl << G4endl; >> 116 for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin(); >> 117 vit!=(*iter).second.end(); vit++) >> 118 { >> 119 G4cout << "--> Type: " << (*vit).type >> 120 << " Value: " << (*vit).value << G4endl; >> 121 } >> 122 } >> 123 G4cout << G4endl; >> 124 >> 125 // The same as above, but now we are looking for >> 126 // sensitive detectors setting them for the volumes >> 127 >> 128 for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin(); >> 129 iter!=auxmap->end(); iter++) >> 130 { >> 131 G4cout << "Volume " << ((*iter).first)->GetName() >> 132 << " has the following list of auxiliary information: " >> 133 << G4endl << G4endl; >> 134 for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin(); >> 135 vit!=(*iter).second.end();vit++) >> 136 { >> 137 if ((*vit).type=="SensDet") >> 138 { >> 139 G4cout << "Attaching sensitive detector " << (*vit).value >> 140 << " to volume " << ((*iter).first)->GetName() >> 141 << G4endl << G4endl; >> 142 >> 143 G4VSensitiveDetector* mydet = >> 144 SDman->FindSensitiveDetector((*vit).value); >> 145 if(mydet) >> 146 { >> 147 G4LogicalVolume* myvol = (*iter).first; >> 148 myvol->SetSensitiveDetector(mydet); >> 149 } >> 150 else >> 151 { >> 152 G4cout << (*vit).value << " detector not found" << G4endl; >> 153 } >> 154 } >> 155 } >> 156 } >> 157 // >> 158 // End of Auxiliary Information block >> 159 // >> 160 //////////////////////////////////////////////////////////////////////// >> 161 >> 162 #ifdef G4VIS_USE >> 163 G4VisManager* visManager = new G4VisExecutive; >> 164 visManager->Initialize(); >> 165 #endif >> 166 >> 167 if (argc==3) // batch mode >> 168 { >> 169 G4String command = "/control/execute "; >> 170 G4String fileName = argv[2]; >> 171 UImanager->ApplyCommand(command+fileName); >> 172 } >> 173 else // interactive mode >> 174 { >> 175 #ifdef G4UI_USE >> 176 G4UIExecutive* ui = new G4UIExecutive(argc, argv); >> 177 #ifdef G4VIS_USE >> 178 UImanager->ApplyCommand("/control/execute vis.mac"); >> 179 #endif >> 180 ui->SessionStart(); >> 181 delete ui; >> 182 #endif >> 183 } >> 184 >> 185 #ifdef G4VIS_USE >> 186 delete visManager; >> 187 #endif >> 188 delete runManager; 101 189 102 delete visManager; << 190 return 0; 103 delete runManager; << 104 } 191 } 105 192