This program computes the dose deposited in an ionization chamber by a monoenergetic photon beam. The geometry of the chamber satisfies the conditions of charged particle equilibrium. Hence, under idealized conditions, the ratio of the dose deposited over the beam energy fluence must be equal to the mass_energy_transfer coefficient of the wall material.
E.Poon and al, Phys. Med. Biol. 50 (2005) 681
I.Kawrakow, Med. Phys. 27-3 (2000) 499
The chamber is modelized as a cylinder with a cavity in it.
6 parameters define the geometry :
Wall and cavity must be made of the same material, but with different density
All above parameters can be redifined via the UI commands built in DetectorMessenger class
----------------- | | | wall | | ----- | | | | | | | <-+-----+--- cavity ------> | | | | ------> | | | | beam -------------------------------- cylinder axis ------> | | | | ------> | | | | | | | | | | | | | ----- | | | | | -----------------
Monoenergetic incident photon beam is uniformly distribued, perpendicular to the flat end of the chamber. The beam radius can be controled with an UI command built in PrimaryGeneratorMessenger; the default is full wall chamber radius.
Beam regeneration : after each Compton interaction, the scattered photon is reset to its initial state, energy and direction. Consequently, interaction sites are uniformly distribued within the wall material.
This modification must be done in the ParticleChange of the final state of the Compton scattering interaction. Therefore, a specific model (MyKleinNishinaCompton) is assigned to the ComptonScattering process in PhysicsList. MyKleinNishinaCompton inherites from G4KleinNishinaCompton; only the function SampleSecondaries() is overwritten.
The program computes the dose deposited in the cavity and the ratio Dose/Beam_energy_fluence. This ratio is compared to the mass_energy_transfer coefficient of the wall material.
The mass_energy_transfer coefficient needs :
The program needs high statistic to reach precision on the computed dose. The UI command /testem/event/printModulo allows to survey the convergence of the kineticEnergy and dose calculations.
In addition, to increase the program efficiency, the secondary particles which have no chance to reach the cavity are immediately killed (see StackinAction). This feature can be switched off by an UI command (see StackingMessenger).
The simplest way to study the effect of e- tracking parameters on dose deposition is to use the command /testem/stepMax.
The physics lists contains the standard electromagnetic processes, with few modifications listed here.
fanoCavity has several predefined 1D histograms :
The histograms are managed by G4AnalysisManager class and its messenger. The histos can be individually activated with the command :
/analysis/h1/set id nbBins valMin valMax unit
where unit is the desired unit for the histo (MeV or keV, deg or mrad, etc..)
One can control the name of the histograms file with the command:
/analysis/setFileName name (default fanoCavity)
It is possible to choose the format of the histogram file : root (default), hdf5, xml, csv, by changing the default file type in HistoManager.cc
It is also possible to print selected histograms on an ascii file:
/analysis/h1/setAscii id
All selected histos will be written on a file name.ascii (default fanocavity)
% fanoCavity run01.mac
% fanoCavity .... Idle> type your commands .... Idle> exit