| 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 /// \file HistoManager.cc
27 /// \brief Implementation of the HistoManager
28 //
29 //
30 //....oooOO0OOooo........oooOO0OOooo........oo
31 //....oooOO0OOooo........oooOO0OOooo........oo
32
33 #include "HistoManager.hh"
34
35 #include "G4UnitsTable.hh"
36
37 //....oooOO0OOooo........oooOO0OOooo........oo
38
39 HistoManager::HistoManager()
40 {
41 Book();
42 }
43
44 //....oooOO0OOooo........oooOO0OOooo........oo
45
46 void HistoManager::Book()
47 {
48 // Create or get analysis manager
49 // The choice of analysis technology is done
50 // in HistoManager.hh
51 G4AnalysisManager* analysisManager = G4Analy
52 analysisManager->SetDefaultFileType("root");
53 analysisManager->SetFileName(fFileName);
54 analysisManager->SetVerboseLevel(1);
55 analysisManager->SetActivation(true); // en
56
57 // Define histograms start values
58 const G4int kMaxHisto = 14;
59 const G4String id[] = {"0", "1", "2", "3", "
60 const G4String title[] = {
61 "dummy", // 0
62 "total energy deposit", // 1
63 "dummy", // 2
64 "total kinetic energy flow", // 3
65 "energy spectrum of emerging gamma", // 4
66 "energy spectrum of emerging e+-", // 5
67 "energy spectrum of emerging neutrons", /
68 "energy spectrum of emerging protons", //
69 "energy spectrum of emerging deuterons",
70 "energy spectrum of emerging alphas", //
71 "energy spectrum of all others emerging io
72 "energy spectrum of all others emerging ba
73 "energy spectrum of all others emerging me
74 "energy spectrum of all others emerging le
75 };
76
77 // Default values (to be reset via /analysis
78 G4int nbins = 100;
79 G4double vmin = 0.;
80 G4double vmax = 100.;
81
82 // Create all histograms as inactivated
83 // as we have not yet set nbins, vmin, vmax
84 for (G4int k = 0; k < kMaxHisto; k++) {
85 G4int ih = analysisManager->CreateH1(id[k]
86 analysisManager->SetH1Activation(ih, false
87 }
88 }
89
90 //....oooOO0OOooo........oooOO0OOooo........oo
91