| 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 G4H1Wrapper.cc
27 /// \brief This allows statistics compatible
28 /// It is fully relying on G4VAnalysis
29 //
30 // Author: G.Hugo, 08 December 2022
31 //
32 // *******************************************
33 //
34 // G4H1Wrapper
35 //
36 /// This allows event-level statistics.
37 /// This is necessary IN THIS SPECIFIC CASE O
38 /// with e.g. FLUKA, MCNP, PENELOPE, SERPENT,
39 ///
40 /// It is fully relying on G4VAnalysisManager
41 /// Collects info within an event, then flush
42 //
43 // *******************************************
44
45 #include "G4H1Wrapper.hh"
46
47 #include "G4RootAnalysisManager.hh"
48
49 //....oooOO0OOooo........oooOO0OOooo........oo
50
51 G4H1Wrapper::G4H1Wrapper(G4VAnalysisManager* c
52 : fAnalysisManager(analysisManager),
53 fHistoIndex(histoIndex),
54 fHisto(dynamic_cast<G4RootAnalysisManager*
55 // fAxis(fHisto->axis()),
56 fNumBins(fHisto->axis().bins()),
57 fIsSet(false),
58 fUnderflow(0.),
59 fOverflow(0.),
60 fSumSquaredEventTotals(0.),
61 fSumSquaredEventInRangeTotals(0.)
62 {
63 fEventData.assign(fNumBins, 0.);
64 }
65
66 // *******************************************
67 // Fill data WITHIN THE EVENT. Can be called a
68 // *******************************************
69 void G4H1Wrapper::Fill(const G4double abscissa
70 {
71 if (!fIsSet) {
72 fIsSet = true;
73 }
74
75 const G4int binIndex = fHisto->axis().coord_
76
77 // underflow
78 // See external/g4tools/include/tools/histo/
79 if (binIndex == tools::histo::axis_UNDERFLOW
80 fUnderflow += weight;
81 }
82 // overflow
83 // See external/g4tools/include/tools/histo/
84 else if (binIndex == tools::histo::axis_OVER
85 fOverflow += weight;
86 }
87 // in range
88 else {
89 fEventData[binIndex] += weight;
90 }
91 }
92
93 // *******************************************
94 // End of event: flush all data collected duri
95 // to G4VAnalysisManager's G4H1, then reset ev
96 // *******************************************
97 void G4H1Wrapper::EndOfEvent()
98 {
99 // Only do the flush if there was at least o
100 if (fIsSet) {
101 G4double eventTotal = 0.;
102 // FLUSH IN RANGE DATA.
103 for (G4int binIndex = 0; binIndex < fNumBi
104 if (fEventData[binIndex] > 0.) {
105 const G4double binKineticEnergy = fHis
106 const G4double binEventWeight = fEvent
107 fAnalysisManager->FillH1(fHistoIndex,
108 eventTotal += binEventWeight;
109
110 fEventData[binIndex] = 0.;
111 }
112 }
113 // Also keep track of the event's total (I
114 const G4double squaredEventTotal = eventTo
115 fSumSquaredEventInRangeTotals += squaredEv
116 fSumSquaredEventTotals += squaredEventTota
117
118 // FLUSH UNDERFLOW DATA.
119 if (fUnderflow > 0.) {
120 fAnalysisManager->FillH1(fHistoIndex, fH
121 fSumSquaredEventTotals += fUnderflow;
122 fUnderflow = 0.;
123 }
124
125 // FLUSH OVERFLOW DATA.
126 if (fOverflow > 0.) {
127 fAnalysisManager->FillH1(fHistoIndex, fH
128 fSumSquaredEventTotals += fOverflow;
129 fOverflow = 0.;
130 }
131
132 fIsSet = false;
133 }
134 }
135
136 //....oooOO0OOooo........oooOO0OOooo........oo
137