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 // Authors: Susanna Guatelli and Francesco Rom 26 // Authors: Susanna Guatelli and Francesco Romano 27 // susanna@uow.edu.au, francesco.romano@ct.inf 27 // susanna@uow.edu.au, francesco.romano@ct.infn.it 28 // 28 // 29 29 30 #include "RunAction.hh" 30 #include "RunAction.hh" 31 #include "AnalysisManager.hh" 31 #include "AnalysisManager.hh" 32 #include "G4SystemOfUnits.hh" 32 #include "G4SystemOfUnits.hh" 33 #include "G4UnitsTable.hh" 33 #include "G4UnitsTable.hh" 34 #include "G4ios.hh" 34 #include "G4ios.hh" 35 #include "G4Run.hh" 35 #include "G4Run.hh" 36 36 37 #ifdef ANALYSIS_USE 37 #ifdef ANALYSIS_USE 38 RunAction::RunAction(AnalysisManager* analysis 38 RunAction::RunAction(AnalysisManager* analysis, DetectorMessenger* detector) 39 { 39 { 40 40 41 analysisMan = analysis; 41 analysisMan = analysis; 42 detectorMess = detector; 42 detectorMess = detector; 43 43 44 } 44 } 45 #else 45 #else 46 RunAction::RunAction(DetectorMessenger* detect 46 RunAction::RunAction(DetectorMessenger* detector) 47 { 47 { 48 detectorMess = detector; 48 detectorMess = detector; 49 } 49 } 50 #endif 50 #endif 51 51 52 RunAction::~RunAction() 52 RunAction::~RunAction() 53 { } 53 { } 54 54 55 void RunAction::BeginOfRunAction(const G4Run* 55 void RunAction::BeginOfRunAction(const G4Run* aRun) 56 { 56 { 57 G4int run_number = aRun->GetRunID(); 57 G4int run_number = aRun->GetRunID(); 58 G4cout << "### Run " << run_number << " star 58 G4cout << "### Run " << run_number << " start." << G4endl; 59 59 60 #ifdef ANALYSIS_USE 60 #ifdef ANALYSIS_USE 61 G4bool additionalOutput; 61 G4bool additionalOutput; 62 G4String detectorType = detectorMess -> GetD 62 G4String detectorType = detectorMess -> GetDetectorType(); 63 if( detectorType == "DiamondTelescope" ) add 63 if( detectorType == "DiamondTelescope" ) additionalOutput = true; 64 else additionalOutput = false; 64 else additionalOutput = false; 65 65 66 // Create ROOT file, histograms and ntuple 66 // Create ROOT file, histograms and ntuple 67 analysisMan -> book(additionalOutput); 67 analysisMan -> book(additionalOutput); 68 #else 68 #else 69 G4cout << "Output to file disabled: no file 69 G4cout << "Output to file disabled: no file will be generated. " 70 << "To enable it, re-run cmake with the -DWI 70 << "To enable it, re-run cmake with the -DWITH_ANALYSIS_USE=ON flag." << G4endl; 71 #endif 71 #endif 72 72 73 // check if there are pending changes in the 73 // check if there are pending changes in the geometry 74 if( detectorMess -> AreTherePendingChanges() 74 if( detectorMess -> AreTherePendingChanges() ) 75 { 75 { 76 G4cout << "WARNING: pending changes to the 76 G4cout << "WARNING: pending changes to the geometry found during BeginOfRunAction. " 77 << "These will be ignored during the curre 77 << "These will be ignored during the current run. " 78 << "Use /geometrySetup/applyChanges if you 78 << "Use /geometrySetup/applyChanges if you want to apply them in your next run." << G4endl; 79 } 79 } 80 80 81 // print the current geometry setup 81 // print the current geometry setup 82 G4cout << "Simulating the " << detectorMess- 82 G4cout << "Simulating the " << detectorMess->GetDetectorType() << " detector " 83 << "with dimensions: " << detectorMess->GetD 83 << "with dimensions: " << detectorMess->GetDetectorSizeWidth()/um << " um width, " 84 << detectorMess->GetDetectorSizeThickness()/ 84 << detectorMess->GetDetectorSizeThickness()/um << " um thickness." << G4endl; 85 } 85 } 86 86 87 void RunAction::EndOfRunAction(const G4Run* aR 87 void RunAction::EndOfRunAction(const G4Run* aRun) 88 { 88 { 89 G4cout << "Number of events = " << aRun->Get 89 G4cout << "Number of events = " << aRun->GetNumberOfEvent() << G4endl; 90 90 91 #ifdef ANALYSIS_USE 91 #ifdef ANALYSIS_USE 92 // Close the output ROOT file with the results 92 // Close the output ROOT file with the results 93 analysisMan -> finish(); 93 analysisMan -> finish(); 94 #endif 94 #endif 95 } 95 } 96 96 97 97