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 B5/src/RunAction.cc 28 /// \brief Implementation of the B5::RunAction 29 30 #include "RunAction.hh" 31 32 #include "EventAction.hh" 33 34 #include "G4AnalysisManager.hh" 35 36 namespace B5 37 { 38 39 //....oooOO0OOooo........oooOO0OOooo........oo 40 41 RunAction::RunAction(EventAction* eventAction) 42 { 43 // Create the generic analysis manager 44 auto analysisManager = G4AnalysisManager::In 45 analysisManager->SetDefaultFileType("root"); 46 // If the filename extension is not provided 47 // will be used for all files specified with 48 analysisManager->SetVerboseLevel(1); 49 50 // Default settings 51 analysisManager->SetNtupleMerging(true); 52 // Note: merging ntuples is available only w 53 analysisManager->SetFileName("B5"); 54 55 // Book histograms, ntuple 56 57 // Creating 1D histograms 58 analysisManager->CreateH1("Chamber1", "Drift 59 analysisManager->CreateH1("Chamber2", "Drift 60 61 // Creating 2D histograms 62 analysisManager->CreateH2("Chamber1 XY", "Dr 63 50, -1000., 1000, 64 analysisManager->CreateH2("Chamber2 XY", "Dr 65 50, -1500., 1500, 66 67 // Creating ntuple 68 if (fEventAction) { 69 analysisManager->CreateNtuple("B5", "Hits" 70 analysisManager->CreateNtupleIColumn("Dc1H 71 analysisManager->CreateNtupleIColumn("Dc2H 72 analysisManager->CreateNtupleDColumn("ECEn 73 analysisManager->CreateNtupleDColumn("HCEn 74 analysisManager->CreateNtupleDColumn("Time 75 analysisManager->CreateNtupleDColumn("Time 76 analysisManager // column Id = 6 77 ->CreateNtupleDColumn("ECEnergyVector", 78 analysisManager // column Id = 7 79 ->CreateNtupleDColumn("HCEnergyVector", 80 analysisManager->FinishNtuple(); 81 } 82 83 // Set ntuple output file 84 analysisManager->SetNtupleFileName(0, "B5ntu 85 } 86 87 //....oooOO0OOooo........oooOO0OOooo........oo 88 89 void RunAction::BeginOfRunAction(const G4Run* 90 { 91 // inform the runManager to save random numb 92 // G4RunManager::GetRunManager()->SetRandomN 93 94 // Get analysis manager 95 auto analysisManager = G4AnalysisManager::In 96 97 // Reset histograms from previous run 98 analysisManager->Reset(); 99 100 // Open an output file 101 // The default file name is set in RunAction 102 // it can be overwritten in a macro 103 analysisManager->OpenFile(); 104 } 105 106 //....oooOO0OOooo........oooOO0OOooo........oo 107 108 void RunAction::EndOfRunAction(const G4Run* /* 109 { 110 // save histograms & ntuple 111 // 112 auto analysisManager = G4AnalysisManager::In 113 analysisManager->Write(); 114 analysisManager->CloseFile(false); 115 // Keep content of histos so that they are p 116 // The content will be reset at start of the 117 } 118 119 //....oooOO0OOooo........oooOO0OOooo........oo 120 121 } // namespace B5 122