| 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 // G4MSSteppingAction implementation
27 //
28 // Author: M.Asai, 5 May 2006
29 // -------------------------------------------
30
31 #include "G4MSSteppingAction.hh"
32
33 #include "G4LogicalVolume.hh"
34 #include "G4Material.hh"
35 #include "G4Region.hh"
36 #include "G4Step.hh"
37 #include "G4VPhysicalVolume.hh"
38
39 #include <map>
40
41 // -------------------------------------------
42 void G4MSSteppingAction::Initialize(G4bool rSe
43 {
44 regionSensitive = rSens;
45 theRegion = reg;
46 length = 0.;
47 x0 = 0.;
48 lambda = 0.;
49 shape_mat_info_v.clear();
50 }
51
52 // -------------------------------------------
53 void G4MSSteppingAction::UserSteppingAction(co
54 {
55 G4StepPoint* preStepPoint = aStep->GetPreSte
56 G4Region* region = preStepPoint->GetPhysical
57
58 if (regionSensitive && (region != theRegion)
59
60 G4double stlen = aStep->GetStepLength();
61 const G4Material* material = preStepPoint->G
62 length += stlen;
63 x0 += stlen / (material->GetRadlen());
64 lambda += stlen / (material->GetNuclearInter
65
66 // store information per step (1 geantino st
67 {
68 shape_mat_info_v.push_back({});
69
70 shape_mat_info_t & thisMaterialInfo = shap
71
72 // calculate average mass number and atomi
73 {
74 const std::vector<const G4Element*> * El
75 for( auto & element : *ElementVector_ptr
76 {
77 thisMaterialInfo.aveA += element->GetA
78 thisMaterialInfo.aveZ += element->GetZ
79 }
80 thisMaterialInfo.aveA /= ElementVector_p
81 thisMaterialInfo.aveZ /= ElementVector_p
82 }
83 thisMaterialInfo.density = mate
84 thisMaterialInfo.radiation_length = mate
85 thisMaterialInfo.interaction_length = mate
86 thisMaterialInfo.thickness = aS
87 thisMaterialInfo.integrated_thickness = le
88 thisMaterialInfo.lambda = stlen
89 thisMaterialInfo.x0 = stlen
90 thisMaterialInfo.integrated_lambda = lambd
91 thisMaterialInfo.integrated_x0 = x0;
92 thisMaterialInfo.entry_point = preStepPoin
93 thisMaterialInfo.exit_point = aStep->GetP
94 thisMaterialInfo.material_name = material-
95 }
96
97 }
98
99 void G4MSSteppingAction::PrintEachMaterialVerb
100 {
101
102 G4int colwidth = 11;
103 G4int matname_colwidth = 15;
104
105 oss << G4endl<< G4endl;
106 oss << " Material Atomic propertie
107 oss << G4endl;
108
109
110 std::ios::fmtflags os_flags (oss.flags());
111 for( auto & matInfo : shape_mat_info_v)
112 {
113 oss << std::setw(matname_colwidth) << std
114
115 oss << std::setw(colwidth) << std::fixed
116 oss << std::setw(colwidth) << std::fixed
117 oss << std::setw(colwidth) << std::scient
118 oss << std::setw(colwidth) << std::scient
119 oss << std::setw(colwidth) << std::scient
120 oss << std::setw(colwidth) << std::scient
121 oss << std::setw(colwidth) << std::scient
122 oss << std::setw(colwidth) << std::scient
123 oss << std::setw(colwidth) << std::scient
124 oss << std::setw(colwidth) << " ";
125 oss << std::scientific << std::right <<
126 oss << std::scientific << std::left <<
127 oss << ", " << matInfo.entry_point.y()/
128 oss << ", " << matInfo.entry_point.z()/
129 oss << std::setw(colwidth) << " ";
130 oss << std::scientific << std::right <<
131 oss << std::scientific << std::left <<
132 oss << ", " << matInfo.exit_point.y()/CL
133 oss << ", " << matInfo.exit_point.z()/CL
134 oss << G4endl;
135 oss << G4endl;
136 }
137 oss.flags(os_flags); // Restore original st
138 }
139
140 void G4MSSteppingAction::PrintIntegratedMateri
141 {
142 // create database (db) of material name (ke
143 std::map<G4String, shape_mat_info_t> mat_db;
144 // accumulate information for each material
145 for(auto & matInfo : shape_mat_info_v)
146 {
147 G4String key = matInfo.material_name;
148 if( 0 == mat_db.count( key ) )
149 {
150 mat_db[key] = shape_mat_info_t{};
151 }
152
153 mat_db[key].x0 += matInfo.x0;
154 mat_db[key].thickness += matInfo.thickness
155 mat_db[key].lambda += matInfo.lambda;
156 }
157
158 std::ios::fmtflags os_flags (oss.flags());
159 oss << std::scientific << std::setprecision(
160 for(auto & [key,mat] : mat_db)
161 oss << '\t' << key
162 << '\t'<< mat.thickness/CLHEP::mm
163 << '\t'<< mat.x0
164 << '\t'<< mat.lambda;
165 oss.flags(os_flags); // Restore original st
166 return;
167 }
168