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)
74 CLHEP::HepRandom::showEngineStatus();
118 if ( analysisManager->IsActive() ) {
119 analysisManager->OpenFile();
130 while ((i<nbProc)&&((*
fProcCounter)[i]->GetName()!=procName)) i++;
133 (*fProcCounter)[i]->Count();
140 if (NbofEvents == 0)
return;
160 G4double doseOverBeam = doseCavity*surfaceBeam/(NbofEvents*beamEnergy);
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 <<
" %)"
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*"
211 G4cout <<
"\n the cavity is "
212 <<
G4BestUnit(fCavityThickness,
"Length") <<
" of "
214 <<
G4BestUnit(fDensityCavity,
"Volumic Mass") <<
"); Mass = "
217 G4cout <<
"\n ============================================================\n";
221 G4cout <<
"\n Process calls frequency --->";
223 G4String procName = (*fProcCounter)[i]->GetName();
224 G4int count = (*fProcCounter)[i]->GetCounter();
225 G4cout <<
" " << procName <<
"= " << count;
232 G4cout <<
"\n Gamma crossSections in wall material :";
235 G4String procName = (*fProcCounter)[i]->GetName();
239 if (massSigma > 0.) {
241 G4cout <<
" " << procName <<
"= "
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: "
284 <<
"\n StoppingPower in wall = "
287 <<
G4BestUnit(dedxCavity,
"Energy*Surface/Mass")
293 <<
"\n Charged particle flow in cavity :"
312 G4double varianceEdep = meanEdep2 - meanEdep*meanEdep;
314 if(varianceEdep>0.) dEoverE = std::sqrt(varianceEdep/
fNbEventCavity)/meanEdep;
319 G4double doseOverBeam = doseCavity*surfaceBeam/(NbofEvents*
energy);
327 <<
" +- " << 100*dEoverE <<
" %"
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;
348 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
352 <<
"\n StepSize of ch. tracks in wall = "
358 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
361 <<
"\n StepSize of ch. tracks in cavity = "
368 G4cout.setf(mode,std::ios::floatfield);
381 if ( analysisManager->IsActive() ) {
382 analysisManager->Write();
383 analysisManager->CloseFile();
387 CLHEP::HepRandom::showEngineStatus();
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)
DetectorConstruction * fDetector
const G4String & GetParticleName() const
HistoManager * fHistoManager
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)
PrimaryGeneratorAction * fKinematic
G4GLOB_DLL std::ostream G4cout
G4int GetNumberOfEvent() const
void EndOfRunAction(const G4Run *)
ExG4HbookAnalysisManager G4AnalysisManager
G4double GetCavityThickness()
G4long fPartFlowCavity[2]
std::map< G4String, G4int > fProcCounter
static G4RunManager * GetRunManager()
G4ParticleGun * GetParticleGun()
G4double energy(const ThreeVector &p, const G4double m)
G4double ComputeCrossSectionPerVolume(G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, const G4Material *, G4double cut=0.0)
G4double fEnerFlowCavity[2]
void SurveyConvergence(G4int)
G4ParticleDefinition * GetParticleDefinition() const
G4double GetCavityRadius()
static G4Electron * Electron()
Detector construction class to demonstrate various ways of placement.
static PROLOG_HANDLER error
G4double GetParticleEnergy() const
G4double fCavityThickness