| 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 // Author: Ivana Hrivnacova, 10/08/2022 (ivan
28
29 #include "G4THnToolsManager.hh"
30
31 #include "tools/histo/h3d"
32
33 using namespace G4Analysis;
34
35 #include <fstream>
36
37 // Specialization for H3 type
38
39
40 //____________________________________________
41 template <>
42 tools::histo::h3d* G4THnToolsManager<kDim3, to
43 const G4String& title,
44 const std::array<G4HnDimension, kDim3>& bins
45 const std::array<G4HnDimensionInformation, k
46 {
47 // Apply hn information to bins
48 auto newXBins(bins[kX]);
49 Update(newXBins, hnInfo[kX]);
50 auto newYBins(bins[kY]);
51 Update(newYBins, hnInfo[kY]);
52 auto newZBins(bins[kZ]);
53 Update(newZBins, hnInfo[kZ]);
54
55 if ((hnInfo[kX].fBinScheme == G4BinScheme::k
56 (hnInfo[kY].fBinScheme == G4BinScheme::k
57 (hnInfo[kZ].fBinScheme == G4BinScheme::k
58 return new tools::histo::h3d(title,
59 newXBins.fNBins, newXBins.fMi
60 newYBins.fNBins, newYBins.fMi
61 newZBins.fNBins, newZBins.fMi
62 }
63
64 return new tools::histo::h3d(title,
65 newXBins.fEdges, newYBins.fEdge
66 }
67
68 //____________________________________________
69 template <>
70 void G4THnToolsManager<kDim3, tools::histo::h3
71 tools::histo::h3d* ht,
72 const std::array<G4HnDimension, kDim3>& bins
73 const std::array<G4HnDimensionInformation, k
74 {
75 // Apply hn information to bins
76 auto newXBins(bins[kX]);
77 Update(newXBins, hnInfo[kX]);
78 auto newYBins(bins[kY]);
79 Update(newYBins, hnInfo[kY]);
80 auto newZBins(bins[kZ]);
81 Update(newZBins, hnInfo[kZ]);
82
83 if ((hnInfo[kX].fBinScheme == G4BinScheme::k
84 (hnInfo[kY].fBinScheme == G4BinScheme::k
85 (hnInfo[kZ].fBinScheme == G4BinScheme::k
86 ht->configure(
87 newXBins.fNBins, newXBins.fMinValue,
88 newYBins.fNBins, newYBins.fMinValue,
89 newZBins.fNBins, newZBins.fMinValue,
90 return;
91 }
92
93 ht->configure(newXBins.fEdges, newYBins.fEdg
94 }
95
96 //____________________________________________
97 template <>
98 G4bool G4THnToolsManager<kDim3, tools::histo::
99 tools::histo::h3d* ht, const G4HnInformation
100 std::array<G4double, kDim3>& value, G4double
101 {
102 const auto& xInfo = hnInformation.GetHnDimen
103 const auto& yInfo = hnInformation.GetHnDimen
104 const auto& zInfo = hnInformation.GetHnDimen
105
106 // Apply hn information to value
107 Update(value[kX], xInfo);
108 Update(value[kY], yInfo);
109 Update(value[kZ], zInfo);
110
111 // Fill updated value
112 ht->fill(value[kX], value[kY], value[kZ], we
113
114 return true;
115 }
116
117 //____________________________________________
118 template <>
119 G4bool G4THnToolsManager<kDim3, tools::histo::
120 std::ofstream& output)
121 {
122 // Write selected objects on ASCII file
123 // According to the implementation by Michel M
124 // extended examples.
125
126 // Do nothing if no histograms are selected
127 if ( ! GetHnManager()->IsAscii() ) return tr
128
129 // Write h3 histograms
130 auto id = GetHnManager()->GetFirstId();
131 for (const auto& [h3, info] : *GetTHnVector(
132
133 if ( (h3 == nullptr) || (! info->GetAscii(
134 // skip writing
135 // if h3 was deleted or writing ascii is
136 id++;
137 continue;
138 }
139
140 Message(kVL3, "write on ascii", "h3d", inf
141
142 output << "\n 3D histogram " << id++ << "
143 << "\n \n \t \t \t X \t\t Y
144
145 for (G4int j=0; j< G4int(h3->axis_x().bins
146 for (G4int k=0; k< G4int(h3->axis_y().bi
147 for (G4int l=0; l< G4int(h3->axis_y().
148 output << " " << j << "\t" << k <<
149 << h3->axis_x().bin_center(j)
150 << h3->axis_y().bin_center(k)
151 << h3->axis_y().bin_center(l)
152 << h3->bin_height(j, k, l) <<
153 }
154 }
155 }
156 }
157
158 return output.good();
159 }
160