54 fNx(0), fNy(0), fNz(0)
113 re02Run->
GetHitsMap(
"PhantomSD/chargedPassCellFlux");
115 re02Run->
GetHitsMap(
"PhantomSD/chargedCellFlux");
117 re02Run->
GetHitsMap(
"PhantomSD/chargedSurfFlux");
119 re02Run->
GetHitsMap(
"PhantomSD/gammaSurfCurr000");
121 re02Run->
GetHitsMap(
"PhantomSD/gammaSurfCurr001");
123 re02Run->
GetHitsMap(
"PhantomSD/gammaSurfCurr002");
125 re02Run->
GetHitsMap(
"PhantomSD/gammaSurfCurr003");
127 G4cout <<
"============================================================="
130 G4cout <<
"============================================================="
132 G4cout << std::setw( 8) <<
"#Z Cell#";
148 for ( iz = 0; iz <
fNz; iz++){
159 if ( !totED ) totED =
new G4double(0.0);
160 if ( !proED ) proED =
new G4double(0.0);
161 if ( !proNS ) proNS =
new G4double(0.0);
162 if ( !pasCF ) pasCF =
new G4double(0.0);
164 if ( !sfx ) sfx =
new G4double(0.0);
165 if ( !gcur0 ) gcur0 =
new G4double(0.0);
166 if ( !gcur1 ) gcur1 =
new G4double(0.0);
167 if ( !gcur2 ) gcur2 =
new G4double(0.0);
168 if ( !gcur3 ) gcur3 =
new G4double(0.0);
169 G4cout << std::setw( 6) << iz <<
" "
170 << std::setw(12) <<
G4BestUnit(*totED,
"Energy")
171 << std::setw(12) <<
G4BestUnit(*proED,
"Energy")
172 << std::setw(12) << (*proNS) <<
" "
173 << std::setw(13) << (*pasCF)*
cm*
cm <<
" /cm2"
174 << std::setw(15) << (*cF)*
cm*
cm <<
" /cm2"
175 << std::setw(15) << (*sfx)*
cm*
cm <<
" /cm2"
176 << std::setw(15) << (*gcur0)*
cm*
cm <<
" /cm2"
177 << std::setw(15) << (*gcur1)*
cm*
cm <<
" /cm2"
178 << std::setw(15) << (*gcur2)*
cm*
cm <<
" /cm2"
179 << std::setw(15) << (*gcur3)*
cm*
cm <<
" /cm2"
182 G4cout <<
"============================================="<<
G4endl;
184 std::ofstream
file(
"totED.txt");
185 for ( iz = 0; iz <
fNz; iz++){
186 for ( iy = 0; iy <
fNy; iy++){
187 for ( ix = 0; ix <
fNx; ix++){
189 if ( !totED ) totED =
new G4double(0.0);
190 file << ix <<
" "<<iy<<
" "<<iz<<
" "<< *totED/
MeV <<
G4endl;