Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
I am running dnadamage1 example (extented-medical-dna-dnadamage1)
Voxel model irradiated with alpha particles with energy 1 MeV. About 4000 particles shots.
Obtained result reveals a significant uderestimation of Indirect Strand Breaks.
Total direct strand breaks SSBd=299843
total indirect strad breaks SSBi=15792
So, the number of indirect SB 20-times smaller, than number of direct SB. But as we know from literature, for alpha particles with 1 MeV enegry, it should be approximatly equal to the amound of direct SB.
This is interesting! This example is done for a short chromatical fiber. It could be this reason. Maybe you want to try moleculedna or dsbandrepair examples ( advanced/dna)?
I have also tried dsbandrepair-example and find out the same underestimation for indirect SB.
Here is the output of dsbandrepair geant 11.2.2 example. Only 1 option changed from original - I have choose G4emDNA_option2 instead of G4emDNA.
Alpha particles 10 MeV, 5000 ions shot
#=========================== Strand Breaks ============================#
Name of Cell Nucleus: Undefined
Volume of Cell Nucleus: 7.32358e-16 (m3)
Mass Density of Cell Nucleus: 1000 (kg/m3)
Mass of Cell Nucleus: 7.32358e-13 (kg)
Energy deposited in Cell Nucleus: 1.67964e+08 (eV)
Dose delivered in Cell Nucleus: 36.6955 (Gy)
Minimum Distance between two clusters: 10 (bp)
Number of basepairs in Cell Nucleus: 6.23368e+09 (bp)
Threshold Energy for direct damage selection: 17.5 (eV)
Propability for indirect damage selection: 40 (%)
#======================================================================#
#=========================== Strand Breaks ============================#
Name of Cell Nucleus: Undefined
Volume of Cell Nucleus: 7.32358e-16 (m3)
Mass Density of Cell Nucleus: 1000 (kg/m3)
Mass of Cell Nucleus: 7.32358e-13 (kg)
Energy deposited in Cell Nucleus: 4.62892e+08 (eV)
Dose delivered in Cell Nucleus: 101.129 (Gy)
Minimum Distance between two clusters: 10 (bp)
Number of basepairs in Cell Nucleus: 6.23368e+09 (bp)
Threshold Energy for direct damage selection: 17.5 (eV)
Propability for indirect damage selection: 40 (%)
#======================================================================#
This is interesting! This example is done for a short chromatical fiber. It could be this reason. Maybe you want to try moleculedna or dsbandrepair examples ( advanced/dna)?
Thank you for your reply. dnadamage1 is the simplest example developed for counting direct-indirect strand breaks with Voxel_straight model, so first of all i have run the simplest example
Thanks for your result! In this example, if some DNA molecules are involved in direct damage, they are no longer engaged in indirect damage. At this LET, this may have happened a lot. Can you cite some papers for equal damage ?