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 // GEANT 4 - Underground Dark Matter Detecto 28 // GEANT 4 - Underground Dark Matter Detector Advanced Example 29 // 29 // 30 // For information related to this code c 30 // For information related to this code contact: Alex Howard 31 // e-mail: alexander.howard@cern.ch 31 // e-mail: alexander.howard@cern.ch 32 // ------------------------------------------- 32 // -------------------------------------------------------------- 33 // Comments 33 // Comments 34 // 34 // 35 // Underground Advanced example mai 35 // Underground Advanced example main program 36 // by A. Howard and H. Araujo << 36 // by A. Howard and H. Araujo 37 // (27th November 2001) 37 // (27th November 2001) 38 // 38 // 39 // main program 39 // main program 40 // ------------------------------------------- 40 // -------------------------------------------------------------- 41 41 42 #include "G4Types.hh" << 42 #include "G4RunManager.hh" 43 #include "G4RunManagerFactory.hh" << 44 << 45 #include "G4UImanager.hh" 43 #include "G4UImanager.hh" >> 44 46 #include "Randomize.hh" 45 #include "Randomize.hh" 47 46 >> 47 #ifdef G4VIS_USE 48 #include "G4VisExecutive.hh" 48 #include "G4VisExecutive.hh" >> 49 #endif >> 50 >> 51 #ifdef G4UI_USE 49 #include "G4UIExecutive.hh" 52 #include "G4UIExecutive.hh" 50 #include "G4AnalysisManager.hh" << 53 #endif 51 54 >> 55 #include "DMXAnalysisManager.hh" 52 #include "DMXDetectorConstruction.hh" 56 #include "DMXDetectorConstruction.hh" 53 #include "DMXPhysicsList.hh" 57 #include "DMXPhysicsList.hh" 54 #include "DMXActionInitializer.hh" << 58 #include "DMXPrimaryGeneratorAction.hh" >> 59 #include "DMXRunAction.hh" >> 60 #include "DMXEventAction.hh" >> 61 #include "DMXSteppingAction.hh" >> 62 #include "DMXStackingAction.hh" >> 63 >> 64 #include <vector> 55 65 56 int main(int argc,char** argv) { 66 int main(int argc,char** argv) { 57 67 >> 68 // choose the Random engine >> 69 CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine); >> 70 58 // Construct the default run manager 71 // Construct the default run manager 59 auto* runManager = G4RunManagerFactory::Crea << 72 G4RunManager * runManager = new G4RunManager; 60 G4int nThreads = 4; << 61 runManager->SetNumberOfThreads(nThreads); << 62 73 63 // set mandatory initialization classes 74 // set mandatory initialization classes 64 runManager->SetUserInitialization(new DMXDet << 75 DMXDetectorConstruction* detector = new DMXDetectorConstruction; >> 76 runManager->SetUserInitialization(detector); 65 runManager->SetUserInitialization(new DMXPhy 77 runManager->SetUserInitialization(new DMXPhysicsList); 66 runManager->SetUserInitialization(new DMXAct << 78 67 << 79 #ifdef G4VIS_USE 68 // visualization manager 80 // visualization manager 69 G4VisManager* visManager = new G4VisExecutiv 81 G4VisManager* visManager = new G4VisExecutive; 70 visManager->Initialize(); 82 visManager->Initialize(); >> 83 #endif >> 84 >> 85 // output environment variables: >> 86 #ifdef G4ANALYSIS_USE >> 87 G4cout << G4endl << G4endl << G4endl >> 88 << " User Environment " << G4endl >> 89 << " Using AIDA 3.2.1 analysis " << G4endl; >> 90 #else >> 91 G4cout << G4endl << G4endl << G4endl >> 92 << " User Environment " << G4endl >> 93 << " G4ANALYSIS_USE environment variable not set, NO ANALYSIS " >> 94 << G4endl; >> 95 #endif 71 96 72 #ifdef DMXENV_GPS_USE 97 #ifdef DMXENV_GPS_USE 73 G4cout << " Using GPS and not DMX gun " << G 98 G4cout << " Using GPS and not DMX gun " << G4endl; 74 #else 99 #else 75 G4cout << " Using the DMX gun " << G4endl; 100 G4cout << " Using the DMX gun " << G4endl; 76 #endif 101 #endif 77 << 102 >> 103 // set user action classes >> 104 DMXPrimaryGeneratorAction* DMXGenerator = new DMXPrimaryGeneratorAction; >> 105 runManager->SetUserAction(DMXGenerator); >> 106 // runManager->SetUserAction(new DMXPrimaryGeneratorAction); >> 107 // RunAction is inherited by EventAction for output filenames - will all >> 108 // change when implement proper analysis manager? >> 109 DMXRunAction* DMXRun = new DMXRunAction; >> 110 runManager->SetUserAction(DMXRun); >> 111 DMXEventAction* eventAction = new DMXEventAction(DMXRun,DMXGenerator); >> 112 runManager->SetUserAction(eventAction); >> 113 // eventAction is inherited by SteppingAction in order to switch colour >> 114 // flag: >> 115 runManager->SetUserAction(new DMXSteppingAction(eventAction)); >> 116 runManager->SetUserAction(new DMXStackingAction); 78 117 79 //Initialize G4 kernel 118 //Initialize G4 kernel 80 runManager->Initialize(); 119 runManager->Initialize(); 81 << 120 82 // get the pointer to the User Interface man << 121 // get the pointer to the User Interface manager 83 G4UImanager* UImanager = G4UImanager::GetUIp << 122 G4UImanager* UImanager = G4UImanager::GetUIpointer(); 84 123 85 // Define UI session for interactive mode. 124 // Define UI session for interactive mode. 86 if(argc == 1) 125 if(argc == 1) 87 { 126 { >> 127 #ifdef G4UI_USE 88 G4UIExecutive* ui = new G4UIExecutive(ar 128 G4UIExecutive* ui = new G4UIExecutive(argc, argv); 89 UImanager->ApplyCommand("/control/execut << 129 #ifdef G4VIS_USE >> 130 UImanager->ApplyCommand("/control/execute initInter.mac"); >> 131 #endif 90 ui->SessionStart(); 132 ui->SessionStart(); 91 delete ui; 133 delete ui; >> 134 #endif 92 } 135 } 93 // Batch mode 136 // Batch mode 94 else 137 else 95 { << 138 { 96 G4String command = "/control/execute "; 139 G4String command = "/control/execute "; 97 G4String fileName = argv[1]; 140 G4String fileName = argv[1]; 98 UImanager->ApplyCommand(command+fileName 141 UImanager->ApplyCommand(command+fileName); 99 } 142 } 100 143 101 //Close-out analysis: << 144 // job termination 102 // Save histograms << 145 #ifdef G4ANALYSIS_USE 103 G4AnalysisManager* man = G4AnalysisManager:: << 146 DMXAnalysisManager::getInstance()->Finish(); 104 man->Write(); << 147 G4cout << "Analysis files closed" << G4endl; 105 man->CloseFile(); << 148 #endif 106 149 >> 150 #ifdef G4VIS_USE 107 if(visManager) delete visManager; 151 if(visManager) delete visManager; 108 << 152 #endif 109 delete runManager; 153 delete runManager; 110 154 111 return 0; 155 return 0; 112 } 156 } 113 157 114 158