The NEURON example allows for the modelling of neuron cell irradiation, including physical, physico-chemical and chemical processes (eg. production of oxidative radical species in the vicinity of neurons). It uses realistic geometrical models of neurons generated from a standardized SWC file representing neuron morphology.
A typical neuron cell is composed of a cell body (soma), a single axon, a dendritic tree, and thousands of dendritic spines. In the example, individual compartments of a neuron cell are simulated by spherical and cylindrical voxels. The soma is represented by combination of several spheres, while the dendritic tree is described with combinations of cylinders. Each voxel is represented as interconnection of two tracing points of the neuron model.
A standardized neuromorphometric format (SWC) is an output file representing individual neuron morphology generated by digitally tracing tools based on 3D confocal microscopy images. In the SWC file, different numerical markers (e.g. from 1 to 6) describe different types of tracing points:
In order to simulate a neural network, user can create his own file containing a combination of several individual neurons (see NeuralNETWORK.dat sample file describing a network of 10 pyramidal neurons).
Geant4-DNA models are activated in the neuron model, which is declared as a G4Region. Geant4 condensed EM models are used outside the neuron structure.
The example package contains:
To run the example: see HOW TO RUN THE EXAMPLE
of this README To simulation output: see SIMULATION OUTPUT AND RESULT ANALYSIS
of this README
The code can be compiled with cmake. It works in MT mode.
The geometry is cube (World) made of galactic material. Before computation, user loads a standardized SWC file of a neuron and generates a bounding volume and a homogeneous spherical medium of liquid water. Dimensions of the target volume are automatically extrapolated using SWC file describing 3D coordinates of a neuron. The homogeneous medium contains volumes of neuronal cell and a bounding slice. The side cube (World) is again represented as overall dimensions of neuronal cell that is equal to the diameter of the homogeneous medium.
The construction of whole geometry of neuron morphology is set in the DetectorConstruction class.
User can choose between single-neuron simulation and modelling a neural network. Single-neuron simulation is set by default. To switch simulation to neural network, the following command should be used:
./neuron -network FileName.dat
The primary kinematic consists of a single particle starting at the random positions on the sphere surface. Then, the particle beam is directed towards the bounding slice volume, and traverses the individual neurons (default option). The type of the particle and its energy are set in the PrimaryGeneratorAction class, and can be changed via the G4 build-in commands of G4ParticleGun class.
We included the following options for particle directions:
a) Particles are directed to "square" on the XY plane of bounding slice (or YZ, XZ)
./neuron -mac myMacro.mac -sXY
b) Particles are directed to "disk" on the XY plane of bounding slice (or YZ, XZ)
./neuron -mac myMacro.mac -dXY
c) Particles are directed towards the bounding slice (default option)
./neuron -mac myMacro.mac
The following options of physical and chemical processes are included: Default Livermore physics
./neuron -mac myMacro.mac
a) Livermore + DNAphysics with extended Rudd ionisation model.
./neuron -mac myMacro.mac -dnaliv
b) Livermore + DNAPhysics + DNAChemistry
./neuron -mac myMacro.mac -dnachemON
c) Combination of DNA- and Livermore- physics with hadronic physics.
./neuron -mac myMacro.mac -dnahad
NOTE, that it requires more memory or computing resources when chemistry is ON (b) and also long computational time when dnaphysics activated. Conversely, it can works faster when default.
Look at the src/PhyscisList.cc files.
To get help, run:
./neuron -h
In visualization and interactive mode, run:
./neuron -gui
( OGL used by default)
or you may use your own visualization driver, for instance:
./neuron -vis "DAWNFILE"
"GranuleCell-Nr2.CNG.swc" is the default file and it should be placed into same directory as the executable. You can download it here: http://neuromorpho.org/neuron_info.jsp?neuron_name=GranuleCell-Nr2 You can change neuron`s file name using the following command:
./neuron -gui -swc newFileName.swc
In batch mode , run:
./neuron(.exe) [-mac neuron.in] [-mt numberofThreads] ./neuron -mac ../neuron.in -mt 3 > neuron.out
To get visualization, make sure to uncomment the #/control/execute vis.mac line in the macro. User can start a visualization of the chemical track evolution in time and space using SetEndTime (default-10 ps) and SetVerbose setting in src/ActionInitialization.cc file.
The simulation outputs appears in terminal display.
The main output results are stored in OutputPerEvent.out file, containing for each event. Dend3DEdep.out, Axon3DEdep.out and Soma3DEdep.out files for given dose:
This file can be easily analysed using for example the provided ROOT macro file plotDend.C; to do so:
If you have any questions or wish to notify of updates and/or modification please contact: batmu.nosp@m.nkh@.nosp@m.jinr..nosp@m.ru, dem@j.nosp@m.inr..nosp@m.ru