89 if(!IsMaster())
return;
102 (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
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 <<
"============================================================="
129 G4cout <<
" Number of event processed : "<< aRun->GetNumberOfEvent() << G4endl;
130 G4cout <<
"============================================================="
132 G4cout << std::setw( 8) <<
"#Z Cell#";
133 G4cout << std::setw(16) << totEdep->GetName();
134 G4cout << std::setw(16) << proEdep->GetName();
135 G4cout << std::setw(12) << proNstep->GetName();
136 G4cout << std::setw(21) << passCFx->GetName();
137 G4cout << std::setw(20) << cFx->GetName();
138 G4cout << std::setw(20) << surfFx->GetName();
139 G4cout << std::setw(20) << gCurr00->GetName();
140 G4cout << std::setw(20) << gCurr01->GetName();
141 G4cout << std::setw(20) << gCurr02->GetName();
142 G4cout << std::setw(20) << gCurr03->GetName()
148 for ( iz = 0; iz <
fNz; iz++){
149 G4double* totED = (*totEdep)[
CopyNo(ix,iy,iz)];
150 G4double* proED = (*proEdep)[
CopyNo(ix,iy,iz)];
151 G4double* proNS = (*proNstep)[
CopyNo(ix,iy,iz)];
152 G4double* pasCF = (*passCFx)[
CopyNo(ix,iy,iz)];
153 G4double* cF = (*cFx)[
CopyNo(ix,iy,iz)];
154 G4double* sfx = (*surfFx)[
CopyNo(ix,iy,iz)];
155 G4double* gcur0 = (*gCurr00)[
CopyNo(ix,iy,iz)];
156 G4double* gcur1 = (*gCurr01)[
CopyNo(ix,iy,iz)];
157 G4double* gcur2 = (*gCurr02)[
CopyNo(ix,iy,iz)];
158 G4double* gcur3 = (*gCurr03)[
CopyNo(ix,iy,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);
163 if ( !cF ) cF =
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++){
188 G4double* totED = (*totEdep)[
CopyNo(ix,iy,iz)];
189 if ( !totED ) totED =
new G4double(0.0);
190 file << ix <<
" "<<iy<<
" "<<iz<<
" "<< *totED/MeV << G4endl;