34 #include "RunAction.hh"
35 #include "DetectorConstruction.hh"
36 #include "PrimaryGeneratorAction.hh"
37 #include "HistoManager.hh"
53 :fDetector(det),fKinematic(kin),fProcCounter(0),fHistoManager(0),
54 fMateWall(0),fMateCavity(0)
85 fSurfaceCavity =
pi*fCavityRadius*fCavityRadius;
86 fVolumeCavity = fSurfaceCavity*fCavityThickness;
89 fMassCavity = fVolumeCavity*fDensityCavity;
97 fEsecondary = fEsecondary2 = 0.;
102 fPartFlowCavity[0] = fPartFlowCavity[1] = 0;
103 fEnerFlowCavity[0] = fEnerFlowCavity[1] = 0.;
107 fEdepCavity = fEdepCavity2 = fTrkSegmCavity = 0.;
112 fStepWall = fStepWall2 = fStepCavity = fStepCavity2 =0.;
113 fNbStepWall = fNbStepCavity = 0;
118 if ( analysisManager->IsActive() ) {
119 analysisManager->OpenFile();
128 size_t nbProc = fProcCounter->size();
130 while ((i<nbProc)&&((*fProcCounter)[i]->GetName()!=procName)) i++;
131 if (i == nbProc) fProcCounter->push_back(
new OneProcessCount(procName));
133 (*fProcCounter)[i]->Count();
140 if (NbofEvents == 0)
return;
145 if (fNbSec > 0) meanEsecond = fEsecondary/fNbSec;
148 if (fOldEmean > 0.) rateEmean = 100*(meanEsecond/fOldEmean - 1.);
149 fOldEmean = meanEsecond;
159 G4double doseCavity = fEdepCavity/fMassCavity;
160 G4double doseOverBeam = doseCavity*surfaceBeam/(NbofEvents*beamEnergy);
163 if (fOldDose > 0.) rateDose = 100*(doseOverBeam/fOldDose - 1.);
164 fOldDose = doseOverBeam;
166 std::ios::fmtflags mode =
G4cout.flags();
167 G4cout.setf(std::ios::fixed,std::ios::floatfield);
170 G4cout <<
"\n ---> NbofEvents= " << NbofEvents
171 <<
" NbOfelectr= " << fNbSec
172 <<
" Tkin= " <<
G4BestUnit(meanEsecond,
"Energy")
173 <<
" (" << rateEmean <<
" %)"
174 <<
" NbOfelec in cav= " << fPartFlowCavity[0]
175 <<
" Dose/EnFluence= " <<
G4BestUnit(doseOverBeam,
"Surface/Mass")
176 <<
" (" << rateDose <<
" %)"
180 G4cout.setf(mode,std::ios::floatfield);
188 std::ios::fmtflags mode =
G4cout.flags();
189 G4cout.setf(std::ios::fixed,std::ios::floatfield);
192 if (NbofEvents == 0)
return;
201 G4cout <<
"\n ======================== run summary ======================\n";
205 G4cout <<
"\n The run consists of " << NbofEvents <<
" "<< partName <<
" of "
206 <<
G4BestUnit(energy,
"Energy") <<
" through 2*"
207 <<
G4BestUnit(fWallThickness,
"Length") <<
" of "
208 << fMateWall->
GetName() <<
" (density: "
211 G4cout <<
"\n the cavity is "
212 <<
G4BestUnit(fCavityThickness,
"Length") <<
" of "
213 << fMateCavity->
GetName() <<
" (density: "
214 <<
G4BestUnit(fDensityCavity,
"Volumic Mass") <<
"); Mass = "
217 G4cout <<
"\n ============================================================\n";
221 G4cout <<
"\n Process calls frequency --->";
222 for (
size_t i=0; i< fProcCounter->size();i++) {
223 G4String procName = (*fProcCounter)[i]->GetName();
224 G4int count = (*fProcCounter)[i]->GetCounter();
225 G4cout <<
" " << procName <<
"= " << count;
232 G4cout <<
"\n Gamma crossSections in wall material :";
234 for (
size_t i=0; i< fProcCounter->size();i++) {
235 G4String procName = (*fProcCounter)[i]->GetName();
238 procName,fMateWall)/fDensityWall;
239 if (massSigma > 0.) {
241 G4cout <<
" " << procName <<
"= "
250 if (fNbSec == 0)
return;
251 G4double meanEsecond = fEsecondary/fNbSec, meanEsecond2 = fEsecondary2/fNbSec;
252 G4double varianceEsec = meanEsecond2 - meanEsecond*meanEsecond;
254 if (varianceEsec>0.) dToverT = std::sqrt(varianceEsec/fNbSec)/meanEsecond;
260 <<
"\n Mean energy of secondary e- = " <<
G4BestUnit(meanEsecond,
"Energy")
261 <<
" +- " << 100*dToverT <<
" %"
262 <<
" (--> range in wall material = " <<
G4BestUnit(csdaRange,
"Length")
270 G4cout <<
" Mass_energy_transfer coef: "
278 dedxWall /= fDensityWall;
281 dedxCavity /= fDensityCavity;
284 <<
"\n StoppingPower in wall = "
287 <<
G4BestUnit(dedxCavity,
"Energy*Surface/Mass")
293 <<
"\n Charged particle flow in cavity :"
294 <<
"\n Enter --> nbParticles = " << fPartFlowCavity[0]
295 <<
"\t Energy = " <<
G4BestUnit (fEnerFlowCavity[0],
"Energy")
296 <<
"\n Exit --> nbParticles = " << fPartFlowCavity[1]
297 <<
"\t Energy = " <<
G4BestUnit (fEnerFlowCavity[1],
"Energy")
300 if (fPartFlowCavity[0] == 0)
return;
309 if (fNbEventCavity == 0)
return;
310 G4double meanEdep = fEdepCavity/fNbEventCavity;
311 G4double meanEdep2 = fEdepCavity2/fNbEventCavity;
312 G4double varianceEdep = meanEdep2 - meanEdep*meanEdep;
314 if(varianceEdep>0.) dEoverE = std::sqrt(varianceEdep/fNbEventCavity)/meanEdep;
318 G4double doseCavity = fEdepCavity/fMassCavity;
319 G4double doseOverBeam = doseCavity*surfaceBeam/(NbofEvents*
energy);
322 G4double meantrack = fTrkSegmCavity/fPartFlowCavity[0];
326 <<
"\n Total edep in cavity = " <<
G4BestUnit(fEdepCavity,
"Energy")
327 <<
" +- " << 100*dEoverE <<
" %"
328 <<
"\t Total charged trackLength = " <<
G4BestUnit(fTrkSegmCavity,
"Length")
329 <<
" (mean value = " <<
G4BestUnit(meantrack,
"Length") <<
")"
330 <<
"\n Total dose in cavity = " << doseCavity/(
MeV/
mg) <<
" MeV/mg"
331 <<
"\n Dose/EnergyFluence = " <<
G4BestUnit(doseOverBeam,
"Surface/Mass")
336 G4double ratio = doseOverBeam/massTransfCoef;
337 G4double error = ratio*std::sqrt(dEoverE*dEoverE + dToverT*dToverT);
341 <<
"\n (Dose/EnergyFluence)/Mass_energy_transfer = " << ratio
342 <<
" +- " << error <<
G4endl;
346 fStepWall /= fNbStepWall; fStepWall2 /= fNbStepWall;
347 G4double rms = fStepWall2 - fStepWall*fStepWall;
348 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
352 <<
"\n StepSize of ch. tracks in wall = "
354 <<
"\t (nbSteps/track = " << double(fNbStepWall)/fNbSec <<
")";
356 fStepCavity /= fNbStepCavity; fStepCavity2 /= fNbStepCavity;
357 rms = fStepCavity2 - fStepCavity*fStepCavity;
358 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
361 <<
"\n StepSize of ch. tracks in cavity = "
363 <<
"\t (nbSteps/track = " << double(fNbStepCavity)/fPartFlowCavity[0] <<
")";
368 G4cout.setf(mode,std::ios::floatfield);
372 while (fProcCounter->size()>0){
374 fProcCounter->pop_back();
381 if ( analysisManager->IsActive() ) {
382 analysisManager->Write();
383 analysisManager->CloseFile();
std::vector< OneProcessCount * > ProcessesCount
G4Material * GetCavityMaterial()
void BeginOfRunAction(const G4Run *)
G4Material * GetWallMaterial()
const G4String & GetName() const
G4double GetDensity() const
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
void SetRandomNumberStore(G4bool flag)
void CountProcesses(G4String procName)
const G4String & GetParticleName() const
G4double GetDEDX(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
G4double GetWallThickness()
G4double GetCSDARange(G4double kinEnergy, const G4ParticleDefinition *, const G4Material *, const G4Region *r=0)
double precision function energy(A, Z)
G4GLOB_DLL std::ostream G4cout
G4int GetNumberOfEvent() const
void EndOfRunAction(const G4Run *)
ExG4HbookAnalysisManager G4AnalysisManager
G4double GetCavityThickness()
static void showEngineStatus()
static G4RunManager * GetRunManager()
G4ParticleGun * GetParticleGun()
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
void SurveyConvergence(G4int)
G4ParticleDefinition * GetParticleDefinition() const
G4double GetCavityRadius()
static G4Electron * Electron()
G4double GetParticleEnergy() const