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/G01/load_gdml.cc 26 /// \file persistency/gdml/G01/load_gdml.cc 27 /// \brief Main program of the persistency/gdm 27 /// \brief Main program of the persistency/gdml/G01 example 28 // 28 // 29 // 29 // 30 // 30 // 31 // 31 // 32 // ------------------------------------------- 32 // -------------------------------------------------------------- 33 // GEANT 4 - load_gdml 33 // GEANT 4 - load_gdml 34 // 34 // 35 // ------------------------------------------- 35 // -------------------------------------------------------------- 36 36 37 #include "FTFP_BERT.hh" << 37 #include <vector> 38 #include "G01ActionInitialization.hh" << 39 #include "G01DetectorConstruction.hh" << 40 #include "G01PrimaryGeneratorAction.hh" << 41 38 42 #include "G4GDMLParser.hh" << 43 #include "G4LogicalVolumeStore.hh" << 44 #include "G4RunManagerFactory.hh" << 45 #include "G4TransportationManager.hh" << 46 #include "G4Types.hh" 39 #include "G4Types.hh" 47 #include "G4UIExecutive.hh" << 40 >> 41 #ifdef G4MULTITHREADED >> 42 #include "G4MTRunManager.hh" >> 43 #else >> 44 #include "G4RunManager.hh" >> 45 #endif >> 46 48 #include "G4UImanager.hh" 47 #include "G4UImanager.hh" >> 48 >> 49 #include "G4LogicalVolumeStore.hh" >> 50 #include "G4TransportationManager.hh" >> 51 >> 52 #include "G01PrimaryGeneratorAction.hh" >> 53 #include "G01DetectorConstruction.hh" >> 54 #include "G01ActionInitialization.hh" >> 55 >> 56 #include "FTFP_BERT.hh" >> 57 49 #include "G4VisExecutive.hh" 58 #include "G4VisExecutive.hh" >> 59 #include "G4UIExecutive.hh" 50 60 51 #include <vector> << 61 #include "G4GDMLParser.hh" 52 62 53 void print_aux(const G4GDMLAuxListType* auxInf << 63 void print_aux(const G4GDMLAuxListType* auxInfoList, G4String prepend="|") 54 { 64 { 55 for (std::vector<G4GDMLAuxStructType>::const << 65 for(std::vector<G4GDMLAuxStructType>::const_iterator 56 iaux != auxInfoList->end(); iaux++) << 66 iaux = auxInfoList->begin(); iaux != auxInfoList->end(); iaux++ ) 57 { << 67 { 58 G4String str = iaux->type; << 68 G4String str=iaux->type; 59 G4String val = iaux->value; << 69 G4String val=iaux->value; 60 G4String unit = iaux->unit; << 70 G4String unit=iaux->unit; 61 71 62 G4cout << prepend << str << " : " << val < << 72 G4cout << prepend << str << " : " << val << " " << unit << G4endl; 63 73 64 if (iaux->auxList) print_aux(iaux->auxList << 74 if (iaux->auxList) print_aux(iaux->auxList, prepend + "|"); 65 } << 75 } 66 return; 76 return; 67 } 77 } 68 78 69 // ------------------------------------------- 79 // -------------------------------------------------------------- 70 80 71 int main(int argc, char** argv) << 81 int main(int argc,char **argv) 72 { 82 { 73 G4cout << G4endl; << 83 G4cout << G4endl; 74 G4cout << "Usage: load_gdml <intput_gdml_fil << 84 G4cout << "Usage: load_gdml <intput_gdml_file:mandatory>" 75 << " <output_gdml_file:optional>" << << 85 << " <output_gdml_file:optional>" << G4endl; 76 G4cout << G4endl; << 86 G4cout << G4endl; 77 << 87 78 if (argc < 2) { << 88 if (argc<2) 79 G4cout << "Error! Mandatory input file is << 89 { 80 G4cout << G4endl; << 90 G4cout << "Error! Mandatory input file is not specified!" << G4endl; 81 return -1; << 91 G4cout << G4endl; 82 } << 92 return -1; 83 << 93 } 84 G4GDMLParser parser; << 94 85 << 95 G4GDMLParser parser; 86 // Uncomment the following if wish to avoid << 96 87 // parser.SetStripFlag(false); << 97 // Uncomment the following if wish to avoid names stripping 88 << 98 // parser.SetStripFlag(false); 89 // Uncomment the following and set a string << 99 90 // schema filename if wishing to use alterna << 100 parser.SetOverlapCheck(true); 91 // parser.SetImportSchema(""); << 101 parser.Read(argv[1]); 92 << 102 93 parser.SetOverlapCheck(true); << 103 if (argc>4) 94 parser.Read(argv[1]); << 104 { 95 << 105 G4cout << "Error! Too many arguments!" << G4endl; 96 if (argc > 4) { << 106 G4cout << G4endl; 97 G4cout << "Error! Too many arguments!" << << 107 return -1; 98 G4cout << G4endl; << 108 } 99 return -1; << 109 100 } << 110 #ifdef G4MULTITHREADED 101 << 111 G4MTRunManager* runManager = new G4MTRunManager; 102 auto* runManager = G4RunManagerFactory::Crea << 112 #else 103 << 113 G4RunManager* runManager = new G4RunManager; 104 runManager->SetUserInitialization(new G01Det << 114 #endif 105 runManager->SetUserInitialization(new FTFP_B << 115 106 runManager->SetUserInitialization(new G01Act << 116 runManager->SetUserInitialization(new G01DetectorConstruction( 107 << 117 parser.GetWorldVolume())); 108 runManager->Initialize(); << 118 runManager->SetUserInitialization(new FTFP_BERT); 109 << 119 runManager->SetUserInitialization(new G01ActionInitialization()); 110 // Initialize visualization << 120 111 G4VisManager* visManager = new G4VisExecutiv << 121 runManager->Initialize(); 112 visManager->Initialize(); << 122 113 << 123 // Initialize visualization 114 // Get the pointer to the User Interface man << 124 G4VisManager* visManager = new G4VisExecutive; 115 G4UImanager* UImanager = G4UImanager::GetUIp << 125 visManager->Initialize(); 116 << 126 117 //////////////////////////////////////////// << 127 // Get the pointer to the User Interface manager 118 // << 128 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 119 // Example how to retrieve Auxiliary Informa << 129 120 // << 130 /////////////////////////////////////////////////////////////////////// 121 << 131 // 122 G4cout << std::endl; << 132 // Example how to retrieve Auxiliary Information 123 << 133 // 124 const G4LogicalVolumeStore* lvs = G4LogicalV << 134 125 std::vector<G4LogicalVolume*>::const_iterato << 135 G4cout << std::endl; 126 for (lvciter = lvs->begin(); lvciter != lvs- << 136 127 G4GDMLAuxListType auxInfo = parser.GetVolu << 137 const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance(); 128 << 138 std::vector<G4LogicalVolume*>::const_iterator lvciter; 129 if (auxInfo.size() > 0) << 139 for( lvciter = lvs->begin(); lvciter != lvs->end(); lvciter++ ) 130 G4cout << "Auxiliary Information is foun << 140 { 131 << G4endl; << 141 G4GDMLAuxListType auxInfo = parser.GetVolumeAuxiliaryInformation(*lvciter); 132 << 142 133 print_aux(&auxInfo); << 143 if (auxInfo.size()>0) 134 } << 144 G4cout << "Auxiliary Information is found for Logical Volume : " 135 << 145 << (*lvciter)->GetName() << G4endl; 136 // now the 'global' auxiliary info << 146 137 G4cout << std::endl; << 147 print_aux(&auxInfo); 138 G4cout << "Global auxiliary info:" << std::e << 148 } 139 G4cout << std::endl; << 149 140 << 150 // now the 'global' auxiliary info 141 print_aux(parser.GetAuxList()); << 151 G4cout << std::endl; 142 << 152 G4cout << "Global auxiliary info:" << std::endl; 143 G4cout << std::endl; << 153 G4cout << std::endl; 144 << 154 145 // << 155 print_aux(parser.GetAuxList()); 146 // End of Auxiliary Information block << 156 147 // << 157 G4cout << std::endl; 148 //////////////////////////////////////////// << 158 149 << 159 // 150 runManager->BeamOn(0); << 160 // End of Auxiliary Information block 151 << 161 // 152 // example of writing out << 162 //////////////////////////////////////////////////////////////////////// 153 << 163 154 if (argc >= 3) { << 164 155 /* << 165 runManager->BeamOn(0); 156 G4GDMLAuxStructType mysubaux = {"mysu << 166 157 G4GDMLAuxListType* myauxlist = new G4 << 167 // example of writing out 158 myauxlist->push_back(mysubaux); << 168 159 << 169 if (argc>=3) 160 G4GDMLAuxStructType myaux = {"mytype" << 170 { 161 parser.AddAuxiliary(myaux); << 171 /* 162 << 172 G4GDMLAuxStructType mysubaux = {"mysubtype", "mysubvalue", "mysubunit", 0}; 163 << 173 G4GDMLAuxListType* myauxlist = new G4GDMLAuxListType(); 164 // example of setting auxiliary info << 174 myauxlist->push_back(mysubaux); 165 // (can be set for any volume) << 175 166 << 176 G4GDMLAuxStructType myaux = {"mytype", "myvalue", "myunit", myauxlist}; 167 G4GDMLAuxStructType mylocalaux = {"so << 177 parser.AddAuxiliary(myaux); 168 << 178 169 parser.AddVolumeAuxiliary(mylocalaux, << 179 170 G4TransportationManager::GetTranspo << 180 // example of setting auxiliary info for world volume 171 ->GetNavigatorForTracking()->GetWor << 181 // (can be set for any volume) 172 */ << 182 173 << 183 G4GDMLAuxStructType mylocalaux = {"sometype", "somevalue", "someunit", 0}; 174 parser.SetRegionExport(true); << 184 175 // parser.SetEnergyCutsExport(true); << 185 parser.AddVolumeAuxiliary(mylocalaux, 176 // parser.SetOutputFileOverwrite(true) << 186 G4TransportationManager::GetTransportationManager() 177 parser.Write(argv[2], G4TransportationMana << 187 ->GetNavigatorForTracking()->GetWorldVolume()->GetLogicalVolume()); 178 ->GetNavigatorForT << 188 */ 179 ->GetWorldVolume() << 189 180 ->GetLogicalVolume << 190 parser.SetRegionExport(true); 181 } << 191 // parser.SetEnergyCutsExport(true); 182 << 192 parser.Write(argv[2], G4TransportationManager::GetTransportationManager() 183 if (argc == 4) // batch mode << 193 ->GetNavigatorForTracking()->GetWorldVolume()->GetLogicalVolume()); 184 { << 194 } 185 G4String command = "/control/execute "; << 195 186 G4String fileName = argv[3]; << 196 187 UImanager->ApplyCommand(command + fileName << 197 if (argc==4) // batch mode 188 } << 198 { 189 else // interactive mode << 199 G4String command = "/control/execute "; 190 { << 200 G4String fileName = argv[3]; 191 G4UIExecutive* ui = new G4UIExecutive(argc << 201 UImanager->ApplyCommand(command+fileName); 192 UImanager->ApplyCommand("/control/execute << 202 } 193 ui->SessionStart(); << 203 else // interactive mode 194 delete ui; << 204 { 195 } << 205 G4UIExecutive* ui = new G4UIExecutive(argc, argv); >> 206 UImanager->ApplyCommand("/control/execute vis.mac"); >> 207 ui->SessionStart(); >> 208 delete ui; >> 209 } 196 210 197 delete visManager; << 211 delete visManager; 198 delete runManager; << 212 delete runManager; 199 213 200 return 0; << 214 return 0; 201 } 215 } 202 216