| 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 = 25;
59 const G4String id[] = {"0", "1", "2", "3"
60 "9", "10", "11", "12
61 "18", "19", "20", "21
62 const G4String title[] = {
63 "dummy", // 0
64 "total energy deposit", // 1
65 "Edep (MeV/mm) along absorber", // 2
66 "total kinetic energy flow", // 3
67 "energy spectrum of gamma at creation", /
68 "energy spectrum of e+- at creation", //
69 "energy spectrum of neutrons at creation",
70 "energy spectrum of protons at creation",
71 "energy spectrum of deuterons at creation"
72 "energy spectrum of alphas at creation",
73 "energy spectrum of all others ions at cre
74 "energy spectrum of all others baryons at
75 "energy spectrum of all others mesons at c
76 "energy spectrum of all others leptons (ne
77 "energy spectrum of emerging gamma", // 1
78 "energy spectrum of emerging e+-", // 15
79 "energy spectrum of emerging neutrons", /
80 "energy spectrum of emerging protons", //
81 "energy spectrum of emerging deuterons",
82 "energy spectrum of emerging alphas", //
83 "energy spectrum of all others emerging io
84 "energy spectrum of all others emerging ba
85 "energy spectrum of all others emerging me
86 "energy spectrum of all others emerging le
87 "total energy released : edep + eflow" //
88 };
89
90 // Default values (to be reset via /analysis
91 G4int nbins = 100;
92 G4double vmin = 0.;
93 G4double vmax = 100.;
94
95 // Create all histograms as inactivated
96 // as we have not yet set nbins, vmin, vmax
97 for (G4int k = 0; k < kMaxHisto; k++) {
98 G4int ih = analysisManager->CreateH1(id[k]
99 analysisManager->SetH1Activation(ih, false
100 }
101 }
102
103 //....oooOO0OOooo........oooOO0OOooo........oo
104