88{
89 if(!IsMaster()) return;
90
91
93
95
96
97
98
99
102 (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
104
105
106
107
108
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");
126
127 G4cout << "============================================================="
128 << G4endl;
129 G4cout << " Number of event processed : "<< aRun->GetNumberOfEvent() << G4endl;
130 G4cout << "============================================================="
131 << G4endl;
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()
143 << G4endl;
146 G4int iz;
147
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"
180 << G4endl;
181 }
182 G4cout << "============================================="<<G4endl;
183
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;
191
192 }
193 }
194 }
195 file.close();
196}
Uer detector construction class.
void GetNumberOfSegmentsInPhantom(G4int &nx, G4int &ny, G4int &nz) const
G4int CopyNo(G4int ix, G4int iy, G4int iz)
G4THitsMap< G4double > * GetHitsMap(G4int i)