Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 /// \file B4/B4b/src/RunAction.cc 28 /// \brief Implementation of the B4b::RunActio 29 30 #include "RunAction.hh" 31 32 #include "RunData.hh" 33 34 #include "G4AnalysisManager.hh" 35 #include "G4Run.hh" 36 #include "G4RunManager.hh" 37 #include "G4SystemOfUnits.hh" 38 #include "G4UnitsTable.hh" 39 #include "globals.hh" 40 41 namespace B4b 42 { 43 44 //....oooOO0OOooo........oooOO0OOooo........oo 45 46 RunAction::RunAction() 47 { 48 // set printing event number per each event 49 G4RunManager::GetRunManager()->SetPrintProgr 50 51 // Create analysis manager 52 // The choice of analysis technology is done 53 // in Analysis.hh 54 auto analysisManager = G4AnalysisManager::In 55 56 // Create directories 57 // analysisManager->SetHistoDirectoryName("h 58 // analysisManager->SetNtupleDirectoryName(" 59 analysisManager->SetVerboseLevel(1); 60 analysisManager->SetNtupleMerging(true); 61 // Note: merging ntuples is available only w 62 63 // Book histograms, ntuple 64 // 65 66 // Creating histograms 67 analysisManager->CreateH1("Eabs", "Edep in a 68 analysisManager->CreateH1("Egap", "Edep in g 69 analysisManager->CreateH1("Labs", "trackL in 70 analysisManager->CreateH1("Lgap", "trackL in 71 72 // Creating ntuple 73 // 74 analysisManager->CreateNtuple("B4", "Edep an 75 analysisManager->CreateNtupleDColumn("Eabs") 76 analysisManager->CreateNtupleDColumn("Egap") 77 analysisManager->CreateNtupleDColumn("Labs") 78 analysisManager->CreateNtupleDColumn("Lgap") 79 analysisManager->FinishNtuple(); 80 } 81 82 //....oooOO0OOooo........oooOO0OOooo........oo 83 84 G4Run* RunAction::GenerateRun() 85 { 86 return (new RunData); 87 } 88 89 //....oooOO0OOooo........oooOO0OOooo........oo 90 91 void RunAction::BeginOfRunAction(const G4Run* 92 { 93 G4cout << "### Run " << run->GetRunID() << " 94 95 // inform the runManager to save random numb 96 // G4RunManager::GetRunManager()->SetRandomN 97 98 // Get analysis manager 99 auto analysisManager = G4AnalysisManager::In 100 101 // Open an output file 102 // 103 G4String fileName = "B4.root"; 104 // Other supported output types: 105 // G4String fileName = "B4.csv"; 106 // G4String fileName = "B4.hdf5"; 107 // G4String fileName = "B4.xml"; 108 analysisManager->OpenFile(fileName); 109 G4cout << "Using " << analysisManager->GetTy 110 } 111 112 //....oooOO0OOooo........oooOO0OOooo........oo 113 114 void RunAction::EndOfRunAction(const G4Run* /* 115 { 116 // print histogram statistics 117 // 118 auto analysisManager = G4AnalysisManager::In 119 if (analysisManager->GetH1(1)) { 120 G4cout << G4endl << " ----> print histogra 121 if (isMaster) { 122 G4cout << "for the entire run " << G4end 123 } 124 else { 125 G4cout << "for the local thread " << G4e 126 } 127 128 G4cout << " EAbs : mean = " << G4BestUnit( 129 << " rms = " << G4BestUnit(analysis 130 131 G4cout << " EGap : mean = " << G4BestUnit( 132 << " rms = " << G4BestUnit(analysis 133 134 G4cout << " LAbs : mean = " << G4BestUnit( 135 << " rms = " << G4BestUnit(analysis 136 137 G4cout << " LGap : mean = " << G4BestUnit( 138 << " rms = " << G4BestUnit(analysis 139 } 140 141 // save histograms & ntuple 142 // 143 analysisManager->Write(); 144 analysisManager->CloseFile(); 145 } 146 147 //....oooOO0OOooo........oooOO0OOooo........oo 148 149 } // namespace B4b 150