Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 /// \file B4/B4c/src/RunAction.cc 28 /// \brief Implementation of the B4::RunAction class 29 30 #include "RunAction.hh" 31 32 #include "G4AnalysisManager.hh" 33 #include "G4RunManager.hh" 34 #include "G4SystemOfUnits.hh" 35 #include "G4UnitsTable.hh" 36 #include "globals.hh" 37 38 namespace B4 39 { 40 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 42 43 RunAction::RunAction() 44 { 45 // set printing event number per each event 46 G4RunManager::GetRunManager()->SetPrintProgress(1); 47 48 // Create analysis manager 49 // The choice of the output format is done via the specified 50 // file extension. 51 auto analysisManager = G4AnalysisManager::Instance(); 52 53 // Create directories 54 // analysisManager->SetHistoDirectoryName("histograms"); 55 // analysisManager->SetNtupleDirectoryName("ntuple"); 56 analysisManager->SetVerboseLevel(1); 57 analysisManager->SetNtupleMerging(true); 58 // Note: merging ntuples is available only with Root output 59 60 // Book histograms, ntuple 61 // 62 63 // Creating histograms 64 analysisManager->CreateH1("Eabs", "Edep in absorber", 110, 0., 330 * MeV); 65 analysisManager->CreateH1("Egap", "Edep in gap", 100, 0., 30 * MeV); 66 analysisManager->CreateH1("Labs", "trackL in absorber", 100, 0., 50 * cm); 67 analysisManager->CreateH1("Lgap", "trackL in gap", 100, 0., 50 * cm); 68 69 // Creating ntuple 70 // 71 analysisManager->CreateNtuple("B4", "Edep and TrackL"); 72 analysisManager->CreateNtupleDColumn("Eabs"); 73 analysisManager->CreateNtupleDColumn("Egap"); 74 analysisManager->CreateNtupleDColumn("Labs"); 75 analysisManager->CreateNtupleDColumn("Lgap"); 76 analysisManager->FinishNtuple(); 77 } 78 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 80 81 void RunAction::BeginOfRunAction(const G4Run* /*run*/) 82 { 83 // inform the runManager to save random number seed 84 // G4RunManager::GetRunManager()->SetRandomNumberStore(true); 85 86 // Get analysis manager 87 auto analysisManager = G4AnalysisManager::Instance(); 88 89 // Open an output file 90 // 91 G4String fileName = "B4.root"; 92 // Other supported output types: 93 // G4String fileName = "B4.csv"; 94 // G4String fileName = "B4.hdf5"; 95 // G4String fileName = "B4.xml"; 96 analysisManager->OpenFile(fileName); 97 G4cout << "Using " << analysisManager->GetType() << G4endl; 98 } 99 100 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 101 102 void RunAction::EndOfRunAction(const G4Run* /*run*/) 103 { 104 // print histogram statistics 105 // 106 auto analysisManager = G4AnalysisManager::Instance(); 107 if (analysisManager->GetH1(1)) { 108 G4cout << G4endl << " ----> print histograms statistic "; 109 if (isMaster) { 110 G4cout << "for the entire run " << G4endl << G4endl; 111 } 112 else { 113 G4cout << "for the local thread " << G4endl << G4endl; 114 } 115 116 G4cout << " EAbs : mean = " << G4BestUnit(analysisManager->GetH1(0)->mean(), "Energy") 117 << " rms = " << G4BestUnit(analysisManager->GetH1(0)->rms(), "Energy") << G4endl; 118 119 G4cout << " EGap : mean = " << G4BestUnit(analysisManager->GetH1(1)->mean(), "Energy") 120 << " rms = " << G4BestUnit(analysisManager->GetH1(1)->rms(), "Energy") << G4endl; 121 122 G4cout << " LAbs : mean = " << G4BestUnit(analysisManager->GetH1(2)->mean(), "Length") 123 << " rms = " << G4BestUnit(analysisManager->GetH1(2)->rms(), "Length") << G4endl; 124 125 G4cout << " LGap : mean = " << G4BestUnit(analysisManager->GetH1(3)->mean(), "Length") 126 << " rms = " << G4BestUnit(analysisManager->GetH1(3)->rms(), "Length") << G4endl; 127 } 128 129 // save histograms & ntuple 130 // 131 analysisManager->Write(); 132 analysisManager->CloseFile(); 133 } 134 135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 136 137 } // namespace B4 138