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