52 const G4double milligray = 1.e-3*gray;
53 const G4double microgray = 1.e-6*gray;
54 const G4double nanogray = 1.e-9*gray;
55 const G4double picogray = 1.e-12*gray;
64 accumulableManager->RegisterAccumulable(
fEdep);
65 accumulableManager->RegisterAccumulable(
fEdep2);
85 G4int nofEvents = run->GetNumberOfEvent();
86 if (nofEvents == 0)
return;
90 accumulableManager->Merge();
94 G4double edep =
fEdep.GetValue();
95 G4double edep2 =
fEdep2.GetValue();
97 G4double rms = edep2 - edep*edep/nofEvents;
98 if (rms > 0.) rms = std::sqrt(rms);
else rms = 0.;
101 G4RunManager::GetRunManager()->GetUserDetectorConstruction());
103 G4double dose = edep/mass;
104 G4double rmsDose = rms/mass;
110 G4RunManager::GetRunManager()->GetUserPrimaryGeneratorAction());
114 const G4ParticleGun* particleGun = generatorAction->GetParticleGun();
115 runCondition += particleGun->GetParticleDefinition()->GetParticleName();
116 runCondition +=
" of ";
117 G4double particleEnergy = particleGun->GetParticleEnergy();
118 runCondition +=
G4BestUnit(particleEnergy,
"Energy");
126 <<
"--------------------End of Global Run-----------------------";
131 <<
"--------------------End of Local Run------------------------";
136 <<
" The run consists of " << nofEvents <<
" "<< runCondition
138 <<
" Cumulated dose per run, in scoring volume : "
141 <<
"------------------------------------------------------------"