From 93a64d344b0ac48d3333de962f58c7a302adec03 Mon Sep 17 00:00:00 2001 From: Thomas BAUDIER Date: Thu, 27 Jan 2022 15:17:34 +0100 Subject: [PATCH] Add t21_tle benchmark --- .github/workflows/main.yml | 1 + t21_tle/data/GateMaterials.db | 454 +++++++++++++++++++ t21_tle/data/Schneider2000DensitiesTable.txt | 13 + t21_tle/data/Schneider2000MaterialsTable.txt | 36 ++ t21_tle/data/SimpleMaterialsTable.txt | 16 + t21_tle/data/phantom-HU2mat.txt | 42 ++ t21_tle/data/phantom-HUmaterials.db | 425 +++++++++++++++++ t21_tle/data/phantom.mhd | 13 + t21_tle/data/phantom.raw | 1 + t21_tle/mac/beam-photon.mac | 99 ++++ t21_tle/mac/beam-photonCT.mac | 152 +++++++ t21_tle/mac/detectors.mac | 23 + t21_tle/mac/main-setle.mac | 107 +++++ t21_tle/mac/main.mac | 109 +++++ t21_tle/mac/template.mac | 8 + t21_tle/mac/verbose.mac | 13 + t21_tle/mac/visu.mac | 15 + t21_tle/output/.gitkeep | 0 t21_tle/runTest.sh | 2 + 19 files changed, 1529 insertions(+) create mode 100644 t21_tle/data/GateMaterials.db create mode 100644 t21_tle/data/Schneider2000DensitiesTable.txt create mode 100644 t21_tle/data/Schneider2000MaterialsTable.txt create mode 100644 t21_tle/data/SimpleMaterialsTable.txt create mode 100644 t21_tle/data/phantom-HU2mat.txt create mode 100644 t21_tle/data/phantom-HUmaterials.db create mode 100644 t21_tle/data/phantom.mhd create mode 120000 t21_tle/data/phantom.raw create mode 100644 t21_tle/mac/beam-photon.mac create mode 100644 t21_tle/mac/beam-photonCT.mac create mode 100644 t21_tle/mac/detectors.mac create mode 100644 t21_tle/mac/main-setle.mac create mode 100644 t21_tle/mac/main.mac create mode 100644 t21_tle/mac/template.mac create mode 100644 t21_tle/mac/verbose.mac create mode 100644 t21_tle/mac/visu.mac create mode 100644 t21_tle/output/.gitkeep create mode 100755 t21_tle/runTest.sh diff --git a/.github/workflows/main.yml b/.github/workflows/main.yml index f0aabcc..8857620 100644 --- a/.github/workflows/main.yml +++ b/.github/workflows/main.yml @@ -35,6 +35,7 @@ jobs: t17_extended_source, t18_movement_benchmark, t20_digi_adder, + t21_tle, t22_merged_actor] steps: diff --git a/t21_tle/data/GateMaterials.db b/t21_tle/data/GateMaterials.db new file mode 100644 index 0000000..284287f --- /dev/null +++ b/t21_tle/data/GateMaterials.db @@ -0,0 +1,454 @@ +[Elements] +Hydrogen: S= H ; Z= 1. ; A= 1.01 g/mole +Helium: S= He ; Z= 2. ; A= 4.003 g/mole +Lithium: S= Li ; Z= 3. ; A= 6.941 g/mole +Beryllium: S= Be ; Z= 4. ; A= 9.012 g/mole +Boron: S= B ; Z= 5. ; A= 10.811 g/mole +Carbon: S= C ; Z= 6. ; A= 12.01 g/mole +Nitrogen: S= N ; Z= 7. ; A= 14.01 g/mole +Oxygen: S= O ; Z= 8. ; A= 16.00 g/mole +Fluorine: S= F ; Z= 9. ; A= 18.998 g/mole +Neon: S= Ne ; Z= 10. ; A= 20.180 g/mole +Sodium: S= Na ; Z= 11. ; A= 22.99 g/mole +Magnesium: S= Mg ; Z= 12. ; A= 24.305 g/mole +Aluminium: S= Al ; Z= 13. ; A= 26.98 g/mole +Silicon: S= Si ; Z= 14. ; A= 28.09 g/mole +Phosphor: S= P ; Z= 15. ; A= 30.97 g/mole +Sulfur: S= S ; Z= 16. ; A= 32.066 g/mole +Chlorine: S= Cl ; Z= 17. ; A= 35.45 g/mole +Argon: S= Ar ; Z= 18. ; A= 39.95 g/mole +Potassium: S= K ; Z= 19. ; A= 39.098 g/mole +Calcium: S= Ca ; Z= 20. ; A= 40.08 g/mole +Scandium: S= Sc ; Z= 21. ; A= 44.956 g/mole +Titanium: S= Ti ; Z= 22. ; A= 47.867 g/mole +Vandium: S= V ; Z= 23. ; A= 50.942 g/mole +Chromium: S= Cr ; Z= 24. ; A= 51.996 g/mole +Manganese: S= Mn ; Z= 25. ; A= 54.938 g/mole +Iron: S= Fe ; Z= 26. ; A= 55.845 g/mole +Cobalt: S= Co ; Z= 27. ; A= 58.933 g/mole +Nickel: S= Ni ; Z= 28. ; A= 58.693 g/mole +Copper: S= Cu ; Z= 29. ; A= 63.39 g/mole +Zinc: S= Zn ; Z= 30. ; A= 65.39 g/mole +Gallium: S= Ga ; Z= 31. ; A= 69.723 g/mole +Germanium: S= Ge ; Z= 32. ; A= 72.61 g/mole +Yttrium: S= Y ; Z= 39. ; A= 88.91 g/mole +Silver: S= Ag ; Z= 47. ; A= 107.868 g/mole +Cadmium: S= Cd ; Z= 48. ; A= 112.41 g/mole +Tin: S= Sn ; Z= 50. ; A= 118.71 g/mole +Tellurium: S= Te ; Z= 52. ; A= 127.6 g/mole +Iodine: S= I ; Z= 53. ; A= 126.90 g/mole +Cesium: S= Cs ; Z= 55. ; A= 132.905 g/mole +Gadolinium: S= Gd ; Z= 64. ; A= 157.25 g/mole +Lutetium: S= Lu ; Z= 71. ; A= 174.97 g/mole +Tungsten: S= W ; Z= 74. ; A= 183.84 g/mole +Gold: S= Au ; Z= 79. ; A= 196.967 g/mole +Thallium: S= Tl ; Z= 81. ; A= 204.37 g/mole +Lead: S= Pb ; Z= 82. ; A= 207.20 g/mole +Bismuth: S= Bi ; Z= 83. ; A= 208.98 g/mole +Uranium: S= U ; Z= 92. ; A= 238.03 g/mole + +[Materials] +Vacuum: d=0.000001 mg/cm3 ; n=1 + +el: name=Hydrogen ; n=1 + +Aluminium: d=2.7 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +AluminiumEGS: d=2.702 g/cm3 ; n=1 ; state=solid + +el: name=Aluminium ; n=1 + +Uranium: d=18.90 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +Silicon: d=2.33 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +Germanium: d=5.32 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +Yttrium: d=4.47 g/cm3 ; n=1 + +el: name=auto ; n=1 + +Gadolinium: d=7.9 g/cm3 ; n=1 + +el: name=auto ; n=1 + +Lutetium: d=9.84 g/cm3 ; n=1 + +el: name=auto ; n=1 + +Tungsten: d=19.3 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +Lead: d=11.4 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +Bismuth: d=9.75 g/cm3 ; n=1 ; state=solid + +el: name=auto ; n=1 + +NaI: d=3.67 g/cm3; n=2; state=solid + +el: name=Sodium ; n=1 + +el: name=Iodine ; n=1 + +PWO: d=8.28 g/cm3; n=3 ; state=Solid + +el: name=Lead; n=1 + +el: name=Tungsten; n=1 + +el: name=Oxygen; n=4 + +BGO: d=7.13 g/cm3; n= 3 ; state=solid + +el: name=Bismuth; n=4 + +el: name=Germanium; n=3 + +el: name=Oxygen; n=12 + +LSO: d=7.4 g/cm3; n=3 ; state=Solid + +el: name=Lutetium ; n=2 + +el: name=Silicon; n=1 + +el: name=Oxygen; n=5 + +Plexiglass: d=1.19 g/cm3; n=3; state=solid + +el: name=Hydrogen; f=0.080538 + +el: name=Carbon; f=0.599848 + +el: name=Oxygen; f=0.319614 + +GSO: d=6.7 g/cm3; n=3 ; state=Solid + +el: name=Gadolinium ; n=2 + +el: name=Silicon; n=1 + +el: name=Oxygen; n=5 + +LuAP: d=8.34 g/cm3; n=3 ; state=Solid + +el: name=Lutetium ; n=1 + +el: name=Aluminium; n=1 + +el: name=Oxygen; n=3 + +YAP: d=5.55 g/cm3; n=3 ; state=Solid + +el: name=Yttrium ; n=1 + +el: name=Aluminium; n=1 + +el: name=Oxygen; n=3 + +Water: d=1.00 g/cm3; n=2 ; state=liquid + +el: name=Hydrogen ; n=2 + +el: name=Oxygen; n=1 + +Quartz: d=2.2 g/cm3; n=2 ; state=Solid + +el: name=Silicon ; n=1 + +el: name=Oxygen ; n=2 + +Breast: d=1.020 g/cm3 ; n = 8 + +el: name=Oxygen; f=0.5270 + +el: name=Carbon; f=0.3320 + +el: name=Hydrogen ; f=0.1060 + +el: name=Nitrogen; f=0.0300 + +el: name=Sulfur ; f=0.0020 + +el: name=Sodium ; f=0.0010 + +el: name=Phosphor; f=0.0010 + +el: name=Chlorine ; f=0.0010 + +Air: d=1.29 mg/cm3 ; n=4 ; state=gas + +el: name=Nitrogen; f=0.755268 + +el: name=Oxygen; f=0.231781 + +el: name=Argon; f=0.012827 + +el: name=Carbon; f=0.000124 + +Glass: d=2.5 g/cm3; n=4; state=solid + +el: name=Sodium ; f=0.1020 + +el: name=Calcium; f=0.0510 + +el: name=Silicon; f=0.2480 + +el: name=Oxygen; f=0.5990 + +Scinti-C9H10: d=1.032 g/cm3 ; n=2 + +el: name=Carbon; n=9 + +el: name=Hydrogen; n=10 + +LuYAP-70: d=7.1 g/cm3 ; n=4 + +el: name=Lutetium ; n= 7 + +el: name=Yttrium ; n= 3 + +el: name=Aluminium; n=10 + +el: name=Oxygen; n=30 + +LuYAP-80: d=7.5 g/cm3 ; n=4 + +el: name=Lutetium ; n= 8 + +el: name=Yttrium ; n= 2 + +el: name=Aluminium; n=10 + +el: name=Oxygen; n=30 + +Plastic: d=1.18 g/cm3 ; n=3; state=solid + +el: name=Carbon ; n=5 + +el: name=Hydrogen ; n=8 + +el: name=Oxygen ; n=2 + +Biomimic: d=1.05 g/cm3 ; n=3; state=solid + +el: name=Carbon ; n=5 + +el: name=Hydrogen ; n=8 + +el: name=Oxygen ; n=2 + +CZT: d=5.68 g/cm3 ; n=3; state=solid + +el: name=Cadmium ; n=9 + +el: name=Zinc ; n=1 + +el: name=Tellurium ; n=10 + +Lung: d=0.26 g/cm3 ; n=9 + +el: name=Hydrogen ; f=0.103 + +el: name=Carbon ; f=0.105 + +el: name=Nitrogen ; f=0.031 + +el: name=Oxygen ; f=0.749 + +el: name=Sodium ; f=0.002 + +el: name=Phosphor ; f=0.002 + +el: name=Sulfur ; f=0.003 + +el: name=Chlorine ; f=0.003 + +el: name=Potassium ; f=0.002 + +Polyethylene: d=0.96 g/cm3 ; n=2 + +el: name=Hydrogen ; n=2 + +el: name=Carbon ; n=1 + +PVC: d=1.65 g/cm3 ; n=3 ; state=solid + +el: name=Hydrogen ; n=3 + +el: name=Carbon ; n=2 + +el: name=Chlorine ; n=1 + +SS304: d=7.92 g/cm3 ; n=4 ; state=solid + +el: name=Iron ; f=0.695 + +el: name=Chromium ; f=0.190 + +el: name=Nickel ; f=0.095 + +el: name=Manganese ; f=0.020 + +PTFE: d= 2.18 g/cm3 ; n=2 ; state=solid + +el: name=Carbon ; n=1 + +el: name=Fluorine ; n=2 + + +LYSO: d=5.37 g/cm3; n=4 ; state=Solid + +el: name=Lutetium ; f=0.31101534 + +el: name=Yttrium ; f=0.368765605 + +el: name=Silicon; f=0.083209699 + +el: name=Oxygen; f=0.237009356 + +Body: d=1.00 g/cm3 ; n=2 + +el: name=Hydrogen ; f=0.112 + +el: name=Oxygen ; f=0.888 + +Muscle: d=1.05 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.102 + +el: name=Carbon ; f=0.143 + +el: name=Nitrogen ; f=0.034 + +el: name=Oxygen ; f=0.71 + +el: name=Sodium ; f=0.001 + +el: name=Phosphor ; f=0.002 + +el: name=Sulfur ; f=0.003 + +el: name=Chlorine ; f=0.001 + +el: name=Potassium ; f=0.004 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +LungMoby: d=0.30 g/cm3 ; n=6 + +el: name=Hydrogen ; f=0.099 + +el: name=Carbon ; f=0.100 + +el: name=Nitrogen ; f=0.028 + +el: name=Oxygen ; f=0.740 + +el: name=Phosphor ; f=0.001 + +el: name=Calcium ; f=0.032 + +SpineBone: d=1.42 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.063 + +el: name=Carbon ; f=0.261 + +el: name=Nitrogen ; f=0.039 + +el: name=Oxygen ; f=0.436 + +el: name=Sodium ; f=0.001 + +el: name=Magnesium ; f=0.001 + +el: name=Phosphor ; f=0.061 + +el: name=Sulfur ; f=0.003 + +el: name=Chlorine ; f=0.001 + +el: name=Potassium ; f=0.001 + +el: name=Calcium ; f=0.133 + +RibBone: d=1.92 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.034 + +el: name=Carbon ; f=0.155 + +el: name=Nitrogen ; f=0.042 + +el: name=Oxygen ; f=0.435 + +el: name=Sodium ; f=0.001 + +el: name=Magnesium ; f=0.002 + +el: name=Phosphor ; f=0.103 + +el: name=Sulfur ; f=0.003 + +el: name=Calcium ; f=0.225 + +el: name=Scandium ; f=0.0 + +el: name=Titanium ; f=0.0 + +Adipose: d=0.92 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.120 + +el: name=Carbon ; f=0.640 + +el: name=Nitrogen ; f=0.008 + +el: name=Oxygen ; f=0.229 + +el: name=Phosphor ; f=0.002 + +el: name=Calcium ; f=0.001 + +el: name=Scandium ; f=0.0 + +el: name=Titanium ; f=0.0 + +el: name=Vandium , f=0.0 + +el: name=Chromium ; f=0.0 + +el: name=Manganese ; f=0.0 + +Blood: d=1.06 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.102 + +el: name=Carbon ; f=0.11 + +el: name=Nitrogen ; f=0.033 + +el: name=Oxygen ; f=0.745 + +el: name=Sodium ; f=0.001 + +el: name=Phosphor ; f=0.001 + +el: name=Sulfur ; f=0.002 + +el: name=Chlorine ; f=0.003 + +el: name=Potassium ; f=0.002 + +el: name=Iron ; f=0.001 + +el: name=Cobalt ; f=0.0 + +Heart: d=1.05 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.104 + +el: name=Carbon ; f=0.139 + +el: name=Nitrogen ; f=0.029 + +el: name=Oxygen ; f=0.718 + +el: name=Sodium ; f=0.001 + +el: name=Phosphor ; f=0.002 + +el: name=Sulfur ; f=0.002 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.003 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Kidney: d=1.05 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.103 + +el: name=Carbon ; f=0.132 + +el: name=Nitrogen ; f=0.03 + +el: name=Oxygen ; f=0.724 + +el: name=Sodium ; f=0.002 + +el: name=Phosphor ; f=0.002 + +el: name=Sulfur ; f=0.002 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.002 + +el: name=Calcium ; f=0.001 + +el: name=Scandium ; f=0.0 + +Liver: d=1.06 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.102 + +el: name=Carbon ; f=0.139 + +el: name=Nitrogen ; f=0.03 + +el: name=Oxygen ; f=0.716 + +el: name=Sodium ; f=0.002 + +el: name=Phosphor ; f=0.003 + +el: name=Sulfur ; f=0.003 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.003 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Lymph: d=1.03 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.108 + +el: name=Carbon ; f=0.041 + +el: name=Nitrogen ; f=0.011 + +el: name=Oxygen ; f=0.832 + +el: name=Sodium ; f=0.003 + +el: name=Sulfur ; f=0.001 + +el: name=Chlorine ; f=0.004 + +el: name=Argon ; f=0.0 + +el: name=Potassium ; f=0.0 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Pancreas: d=1.04 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.106 + +el: name=Carbon ; f=0.169 + +el: name=Nitrogen ; f=0.022 + +el: name=Oxygen ; f=0.694 + +el: name=Sodium ; f=0.002 + +el: name=Phosphor ; f=0.002 + +el: name=Sulfur ; f=0.001 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.002 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Intestine: d=1.03 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.106 + +el: name=Carbon ; f=0.115 + +el: name=Nitrogen ; f=0.022 + +el: name=Oxygen ; f=0.751 + +el: name=Sodium ; f=0.001 + +el: name=Phosphor ; f=0.001 + +el: name=Sulfur ; f=0.001 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.001 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Skull: d=1.61 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.05 + +el: name=Carbon ; f=0.212 + +el: name=Nitrogen ; f=0.04 + +el: name=Oxygen ; f=0.435 + +el: name=Sodium ; f=0.001 + +el: name=Magnesium ; f=0.002 + +el: name=Phosphor ; f=0.081 + +el: name=Sulfur ; f=0.003 + +el: name=Calcium ; f=0.176 + +el: name=Scandium ; f=0.0 + +el: name=Titanium ; f=0.0 + +Cartilage: d=1.10 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.096 + +el: name=Carbon ; f=0.099 + +el: name=Nitrogen ; f=0.022 + +el: name=Oxygen ; f=0.744 + +el: name=Sodium ; f=0.005 + +el: name=Phosphor ; f=0.022 + +el: name=Sulfur ; f=0.009 + +el: name=Chlorine ; f=0.003 + +el: name=Argon ; f=0.0 + +el: name=Potassium ; f=0.0 + +el: name=Calcium ; f=0.0 + +Brain: d=1.04 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.107 + +el: name=Carbon ; f=0.145 + +el: name=Nitrogen ; f=0.022 + +el: name=Oxygen ; f=0.712 + +el: name=Sodium ; f=0.002 + +el: name=Phosphor ; f=0.004 + +el: name=Sulfur ; f=0.002 + +el: name=Chlorine ; f=0.003 + +el: name=Potassium ; f=0.003 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Spleen: d=1.06 g/cm3 ; n=11 + +el: name=Hydrogen ; f=0.103 + +el: name=Carbon ; f=0.113 + +el: name=Nitrogen ; f=0.032 + +el: name=Oxygen ; f=0.741 + +el: name=Sodium ; f=0.001 + +el: name=Phosphor ; f=0.003 + +el: name=Sulfur ; f=0.002 + +el: name=Chlorine ; f=0.002 + +el: name=Potassium ; f=0.003 + +el: name=Calcium ; f=0.0 + +el: name=Scandium ; f=0.0 + +Testis: d=1.04 g/cm3 ; n=9 + +el: name=Hydrogen ; f=0.106000 + +el: name=Carbon ; f=0.099000 + +el: name=Nitrogen ; f=0.020000 + +el: name=Oxygen ; f=0.766000 + +el: name=Sodium ; f=0.002000 + +el: name=Phosphor ; f=0.001000 + +el: name=Sulfur ; f=0.002000 + +el: name=Chlorine ; f=0.002000 + +el: name=Potassium ; f=0.002000 + +PMMA: d=1.195 g/cm3; n=3 ; state=Solid + +el: name=Hydrogen ; f=0.080541 + +el: name=Carbon ; f=0.599846 + +el: name=Oxygen; f=0.319613 + +Epoxy: d=1.0 g/cm3; n=3; state=Solid + +el: name=Carbon; n=1 + +el: name=Hydrogen; n=1 + +el: name=Oxygen; n=1 + +Carbide: d=15.8 g/cm3; n=2 ; state=Solid + +el: name=Tungsten ; n=1 + +el: name=Carbon ; n=1 diff --git a/t21_tle/data/Schneider2000DensitiesTable.txt b/t21_tle/data/Schneider2000DensitiesTable.txt new file mode 100644 index 0000000..d6c1e19 --- /dev/null +++ b/t21_tle/data/Schneider2000DensitiesTable.txt @@ -0,0 +1,13 @@ +# =================== +# HU density g/cm3 +# =================== +-1000 1.21e-3 +-98 0.93 +-97 0.930486 +14 1.03 +23 1.031 +100 1.119900 +101 1.076200 +1600 1.964200 +3000 2.8 + diff --git a/t21_tle/data/Schneider2000MaterialsTable.txt b/t21_tle/data/Schneider2000MaterialsTable.txt new file mode 100644 index 0000000..993a327 --- /dev/null +++ b/t21_tle/data/Schneider2000MaterialsTable.txt @@ -0,0 +1,36 @@ +[Elements] +Hydrogen Carbon Nitrogen Oxygen Sodium Magnesium Phosphor Sulfur +Chlorine Argon Potassium Calcium +Titanium Copper Zinc Silver Tin +[/Elements] +# =============================================================================== +# HU H C N O Na Mg P S Cl Ar K Ca Ti Cu Zn Ag Sn +# =============================================================================== + -1050 0 0 75.5 23.2 0 0 0 0 0 1.3 0 0 0 0 0 0 0 Air + -950 10.3 10.5 3.1 74.9 0.2 0 0.2 0.3 0.3 0 0.2 0 0 0 0 0 0 Lung + -120 11.6 68.1 0.2 19.8 0.1 0 0 0.1 0.1 0 0 0 0 0 0 0 0 AT_AG_SI1 + -82 11.3 56.7 0.9 30.8 0.1 0 0 0.1 0.1 0 0 0 0 0 0 0 0 AT_AG_SI2 + -52 11.0 45.8 1.5 41.1 0.1 0 0.1 0.2 0.2 0 0 0 0 0 0 0 0 AT_AG_SI3 + -22 10.8 35.6 2.2 50.9 0 0 0.1 0.2 0.2 0 0 0 0 0 0 0 0 AT_AG_SI4 + 8 10.6 28.4 2.6 57.8 0 0 0.1 0.2 0.2 0 0.1 0 0 0 0 0 0 AT_AG_SI5 + 19 10.3 13.4 3.0 72.3 0.2 0 0.2 0.2 0.2 0 0.2 0 0 0 0 0 0 SoftTissus + 80 9.4 20.7 6.2 62.2 0.6 0 0 0.6 0.3 0 0.0 0 0 0 0 0 0 ConnectiveTissue + 120 9.5 45.5 2.5 35.5 0.1 0 2.1 0.1 0.1 0 0.1 4.5 0 0 0 0 0 Marrow_Bone01 + 200 8.9 42.3 2.7 36.3 0.1 0 3.0 0.1 0.1 0 0.1 6.4 0 0 0 0 0 Marrow_Bone02 + 300 8.2 39.1 2.9 37.2 0.1 0 3.9 0.1 0.1 0 0.1 8.3 0 0 0 0 0 Marrow_Bone03 + 400 7.6 36.1 3.0 38.0 0.1 0.1 4.7 0.2 0.1 0 0 10.1 0 0 0 0 0 Marrow_Bone04 + 500 7.1 33.5 3.2 38.7 0.1 0.1 5.4 0.2 0 0 0 11.7 0 0 0 0 0 Marrow_Bone05 + 600 6.6 31.0 3.3 39.4 0.1 0.1 6.1 0.2 0 0 0 13.2 0 0 0 0 0 Marrow_Bone06 + 700 6.1 28.7 3.5 40.0 0.1 0.1 6.7 0.2 0 0 0 14.6 0 0 0 0 0 Marrow_Bone07 + 800 5.6 26.5 3.6 40.5 0.1 0.2 7.3 0.3 0 0 0 15.9 0 0 0 0 0 Marrow_Bone08 + 900 5.2 24.6 3.7 41.1 0.1 0.2 7.8 0.3 0 0 0 17.0 0 0 0 0 0 Marrow_Bone09 + 1000 4.9 22.7 3.8 41.6 0.1 0.2 8.3 0.3 0 0 0 18.1 0 0 0 0 0 Marrow_Bone10 + 1100 4.5 21.0 3.9 42.0 0.1 0.2 8.8 0.3 0 0 0 19.2 0 0 0 0 0 Marrow_Bone11 + 1200 4.2 19.4 4.0 42.5 0.1 0.2 9.2 0.3 0 0 0 20.1 0 0 0 0 0 Marrow_Bone12 + 1300 3.9 17.9 4.1 42.9 0.1 0.2 9.6 0.3 0 0 0 21.0 0 0 0 0 0 Marrow_Bone13 + 1400 3.6 16.5 4.2 43.2 0.1 0.2 10.0 0.3 0 0 0 21.9 0 0 0 0 0 Marrow_Bone14 + 1500 3.4 15.5 4.2 43.5 0.1 0.2 10.3 0.3 0 0 0 22.5 0 0 0 0 0 Marrow_Bone15 + 1640 0 0 0 0 0 0 0 0 0 0 0 0 0 4 2 65 29 AmalgamTooth + 2300 0 0 0 0 0 0 0 0 0 0 0 0 100 0 0 0 0 MetallImplants + 3000 0 0 0 0 0 0 0 0 0 0 0 0 100 0 0 0 0 MetallImplants + diff --git a/t21_tle/data/SimpleMaterialsTable.txt b/t21_tle/data/SimpleMaterialsTable.txt new file mode 100644 index 0000000..f5140a4 --- /dev/null +++ b/t21_tle/data/SimpleMaterialsTable.txt @@ -0,0 +1,16 @@ +# =============================================================================== +[Elements] +Hydrogen Carbon Nitrogen Oxygen Sodium Magnesium Phosphor Sulfur +Chlorine Argon Potassium Calcium +Titanium Copper Zinc Silver Tin +[/Elements] +# =============================================================================== +# HU H C N O Na Mg P S Cl Ar K Ca Ti Cu Zn Ag Sn +# =============================================================================== + -1050 0 0 75.5 23.2 0 0 0 0 0 1.3 0 0 0 0 0 0 0 Air + -950 10.3 10.5 3.1 74.9 0.2 0 0.2 0.3 0.3 0 0.2 0 0 0 0 0 0 Lung + 19 10.3 13.4 3.0 72.3 0.2 0 0.2 0.2 0.2 0 0.2 0 0 0 0 0 0 SoftTissus + 200 5.6 26.5 3.6 40.5 0.1 0.2 7.3 0.3 0 0 0 15.9 0 0 0 0 0 Bone + 1640 0 0 0 0 0 0 0 0 0 0 0 0 0 4 2 65 29 AmalgamTooth + 2300 0 0 0 0 0 0 0 0 0 0 0 0 100 0 0 0 0 MetallImplants + 3000 0 0 0 0 0 0 0 0 0 0 0 0 100 0 0 0 0 MetallImplants diff --git a/t21_tle/data/phantom-HU2mat.txt b/t21_tle/data/phantom-HU2mat.txt new file mode 100644 index 0000000..dfd1df0 --- /dev/null +++ b/t21_tle/data/phantom-HU2mat.txt @@ -0,0 +1,42 @@ +-1050 -950 Air_0 +-950 -852.884 Lung_1 +-852.884 -755.769 Lung_2 +-755.769 -658.653 Lung_3 +-658.653 -561.538 Lung_4 +-561.538 -464.422 Lung_5 +-464.422 -367.306 Lung_6 +-367.306 -270.191 Lung_7 +-270.191 -173.075 Lung_8 +-173.075 -120 Lung_9 +-120 -82 AT_AG_SI1_10 +-82 -52 AT_AG_SI2_11 +-52 -22 AT_AG_SI3_12 +-22 8 AT_AG_SI4_13 +8 19 AT_AG_SI5_14 +19 80 SoftTissus_15 +80 120 ConnectiveTissue_16 +120 200 Marrow_Bone01_17 +200 300 Marrow_Bone02_18 +300 400 Marrow_Bone03_19 +400 500 Marrow_Bone04_20 +500 600 Marrow_Bone05_21 +600 700 Marrow_Bone06_22 +700 800 Marrow_Bone07_23 +800 900 Marrow_Bone08_24 +900 1000 Marrow_Bone09_25 +1000 1100 Marrow_Bone10_26 +1100 1200 Marrow_Bone11_27 +1200 1300 Marrow_Bone12_28 +1300 1400 Marrow_Bone13_29 +1400 1500 Marrow_Bone14_30 +1500 1640 Marrow_Bone15_31 +1640 1807.5 AmalgamTooth_32 +1807.5 1975.01 AmalgamTooth_33 +1975.01 2142.51 AmalgamTooth_34 +2142.51 2300 AmalgamTooth_35 +2300 2467.5 MetallImplants_36 +2467.5 2635.01 MetallImplants_37 +2635.01 2802.51 MetallImplants_38 +2802.51 2970.02 MetallImplants_39 +2970.02 3000 MetallImplants_40 +3000 3001 MetallImplants_41 diff --git a/t21_tle/data/phantom-HUmaterials.db b/t21_tle/data/phantom-HUmaterials.db new file mode 100644 index 0000000..c2421a3 --- /dev/null +++ b/t21_tle/data/phantom-HUmaterials.db @@ -0,0 +1,425 @@ +[Materials] +# Material corresponding to H=[ -1050;-950 ] +Air_0: d=1.21 mg/cm3; n=3; ++el: name=Nitrogen; f=0.755 ++el: name=Oxygen; f=0.232 ++el: name=Argon; f=0.013 + +# Material corresponding to H=[ -950;-852.884 ] +Lung_1: d=102.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -852.884;-755.769 ] +Lung_2: d=202.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -755.769;-658.653 ] +Lung_3: d=302.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -658.653;-561.538 ] +Lung_4: d=402.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -561.538;-464.422 ] +Lung_5: d=502.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -464.422;-367.306 ] +Lung_6: d=602.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -367.306;-270.191 ] +Lung_7: d=702.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -270.191;-173.075 ] +Lung_8: d=802.695 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -173.075;-120 ] +Lung_9: d=880.021 mg/cm3; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.105 ++el: name=Nitrogen; f=0.031 ++el: name=Oxygen; f=0.749 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.003 ++el: name=Chlorine; f=0.003 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ -120;-82 ] +AT_AG_SI1_10: d=926.911 mg/cm3; n=7; ++el: name=Hydrogen; f=0.116 ++el: name=Carbon; f=0.681 ++el: name=Nitrogen; f=0.002 ++el: name=Oxygen; f=0.198 ++el: name=Sodium; f=0.001 ++el: name=Sulfur; f=0.001 ++el: name=Chlorine; f=0.001 + +# Material corresponding to H=[ -82;-52 ] +AT_AG_SI2_11: d=957.382 mg/cm3; n=7; ++el: name=Hydrogen; f=0.113 ++el: name=Carbon; f=0.567 ++el: name=Nitrogen; f=0.009 ++el: name=Oxygen; f=0.308 ++el: name=Sodium; f=0.001 ++el: name=Sulfur; f=0.001 ++el: name=Chlorine; f=0.001 + +# Material corresponding to H=[ -52;-22 ] +AT_AG_SI3_12: d=984.277 mg/cm3; n=8; ++el: name=Hydrogen; f=0.11 ++el: name=Carbon; f=0.458 ++el: name=Nitrogen; f=0.015 ++el: name=Oxygen; f=0.411 ++el: name=Sodium; f=0.001 ++el: name=Phosphor; f=0.001 ++el: name=Sulfur; f=0.002 ++el: name=Chlorine; f=0.002 + +# Material corresponding to H=[ -22;8 ] +AT_AG_SI4_13: d=1.01117 g/cm3 ; n=7; ++el: name=Hydrogen; f=0.108 ++el: name=Carbon; f=0.356 ++el: name=Nitrogen; f=0.022 ++el: name=Oxygen; f=0.509 ++el: name=Phosphor; f=0.001 ++el: name=Sulfur; f=0.002 ++el: name=Chlorine; f=0.002 + +# Material corresponding to H=[ 8;19 ] +AT_AG_SI5_14: d=1.02955 g/cm3 ; n=8; ++el: name=Hydrogen; f=0.106 ++el: name=Carbon; f=0.284 ++el: name=Nitrogen; f=0.026 ++el: name=Oxygen; f=0.578 ++el: name=Phosphor; f=0.001 ++el: name=Sulfur; f=0.002 ++el: name=Chlorine; f=0.002 ++el: name=Potassium; f=0.001 + +# Material corresponding to H=[ 19;80 ] +SoftTissus_15: d=1.0616 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.103 ++el: name=Carbon; f=0.134 ++el: name=Nitrogen; f=0.03 ++el: name=Oxygen; f=0.723 ++el: name=Sodium; f=0.002 ++el: name=Phosphor; f=0.002 ++el: name=Sulfur; f=0.002 ++el: name=Chlorine; f=0.002 ++el: name=Potassium; f=0.002 + +# Material corresponding to H=[ 80;120 ] +ConnectiveTissue_16: d=1.1199 g/cm3 ; n=7; ++el: name=Hydrogen; f=0.094 ++el: name=Carbon; f=0.207 ++el: name=Nitrogen; f=0.062 ++el: name=Oxygen; f=0.622 ++el: name=Sodium; f=0.006 ++el: name=Sulfur; f=0.006 ++el: name=Chlorine; f=0.003 + +# Material corresponding to H=[ 120;200 ] +Marrow_Bone01_17: d=1.11115 g/cm3 ; n=10; ++el: name=Hydrogen; f=0.095 ++el: name=Carbon; f=0.455 ++el: name=Nitrogen; f=0.025 ++el: name=Oxygen; f=0.355 ++el: name=Sodium; f=0.001 ++el: name=Phosphor; f=0.021 ++el: name=Sulfur; f=0.001 ++el: name=Chlorine; f=0.001 ++el: name=Potassium; f=0.001 ++el: name=Calcium; f=0.045 + +# Material corresponding to H=[ 200;300 ] +Marrow_Bone02_18: d=1.16447 g/cm3 ; n=10; ++el: name=Hydrogen; f=0.089 ++el: name=Carbon; f=0.423 ++el: name=Nitrogen; f=0.027 ++el: name=Oxygen; f=0.363 ++el: name=Sodium; f=0.001 ++el: name=Phosphor; f=0.03 ++el: name=Sulfur; f=0.001 ++el: name=Chlorine; f=0.001 ++el: name=Potassium; f=0.001 ++el: name=Calcium; f=0.064 + +# Material corresponding to H=[ 300;400 ] +Marrow_Bone03_19: d=1.22371 g/cm3 ; n=10; ++el: name=Hydrogen; f=0.082 ++el: name=Carbon; f=0.391 ++el: name=Nitrogen; f=0.029 ++el: name=Oxygen; f=0.372 ++el: name=Sodium; f=0.001 ++el: name=Phosphor; f=0.039 ++el: name=Sulfur; f=0.001 ++el: name=Chlorine; f=0.001 ++el: name=Potassium; f=0.001 ++el: name=Calcium; f=0.083 + +# Material corresponding to H=[ 400;500 ] +Marrow_Bone04_20: d=1.28295 g/cm3 ; n=10; ++el: name=Hydrogen; f=0.076 ++el: name=Carbon; f=0.361 ++el: name=Nitrogen; f=0.03 ++el: name=Oxygen; f=0.38 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.001 ++el: name=Phosphor; f=0.047 ++el: name=Sulfur; f=0.002 ++el: name=Chlorine; f=0.001 ++el: name=Calcium; f=0.101 + +# Material corresponding to H=[ 500;600 ] +Marrow_Bone05_21: d=1.34219 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.071 ++el: name=Carbon; f=0.335 ++el: name=Nitrogen; f=0.032 ++el: name=Oxygen; f=0.387 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.001 ++el: name=Phosphor; f=0.054 ++el: name=Sulfur; f=0.002 ++el: name=Calcium; f=0.117 + +# Material corresponding to H=[ 600;700 ] +Marrow_Bone06_22: d=1.40142 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.066 ++el: name=Carbon; f=0.31 ++el: name=Nitrogen; f=0.033 ++el: name=Oxygen; f=0.394 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.001 ++el: name=Phosphor; f=0.061 ++el: name=Sulfur; f=0.002 ++el: name=Calcium; f=0.132 + +# Material corresponding to H=[ 700;800 ] +Marrow_Bone07_23: d=1.46066 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.061 ++el: name=Carbon; f=0.287 ++el: name=Nitrogen; f=0.035 ++el: name=Oxygen; f=0.4 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.001 ++el: name=Phosphor; f=0.067 ++el: name=Sulfur; f=0.002 ++el: name=Calcium; f=0.146 + +# Material corresponding to H=[ 800;900 ] +Marrow_Bone08_24: d=1.5199 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.056 ++el: name=Carbon; f=0.265 ++el: name=Nitrogen; f=0.036 ++el: name=Oxygen; f=0.405 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.073 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.159 + +# Material corresponding to H=[ 900;1000 ] +Marrow_Bone09_25: d=1.57914 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.052 ++el: name=Carbon; f=0.246 ++el: name=Nitrogen; f=0.037 ++el: name=Oxygen; f=0.411 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.078 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.17 + +# Material corresponding to H=[ 1000;1100 ] +Marrow_Bone10_26: d=1.63838 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.049 ++el: name=Carbon; f=0.227 ++el: name=Nitrogen; f=0.038 ++el: name=Oxygen; f=0.416 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.083 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.181 + +# Material corresponding to H=[ 1100;1200 ] +Marrow_Bone11_27: d=1.69762 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.045 ++el: name=Carbon; f=0.21 ++el: name=Nitrogen; f=0.039 ++el: name=Oxygen; f=0.42 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.088 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.192 + +# Material corresponding to H=[ 1200;1300 ] +Marrow_Bone12_28: d=1.75686 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.042 ++el: name=Carbon; f=0.194 ++el: name=Nitrogen; f=0.04 ++el: name=Oxygen; f=0.425 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.092 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.201 + +# Material corresponding to H=[ 1300;1400 ] +Marrow_Bone13_29: d=1.8161 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.039 ++el: name=Carbon; f=0.179 ++el: name=Nitrogen; f=0.041 ++el: name=Oxygen; f=0.429 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.096 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.21 + +# Material corresponding to H=[ 1400;1500 ] +Marrow_Bone14_30: d=1.87534 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.036 ++el: name=Carbon; f=0.165 ++el: name=Nitrogen; f=0.042 ++el: name=Oxygen; f=0.432 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.1 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.219 + +# Material corresponding to H=[ 1500;1640 ] +Marrow_Bone15_31: d=1.94643 g/cm3 ; n=9; ++el: name=Hydrogen; f=0.034 ++el: name=Carbon; f=0.155 ++el: name=Nitrogen; f=0.042 ++el: name=Oxygen; f=0.435 ++el: name=Sodium; f=0.001 ++el: name=Magnesium; f=0.002 ++el: name=Phosphor; f=0.103 ++el: name=Sulfur; f=0.003 ++el: name=Calcium; f=0.225 + +# Material corresponding to H=[ 1640;1807.5 ] +AmalgamTooth_32: d=2.03808 g/cm3 ; n=4; ++el: name=Copper; f=0.04 ++el: name=Zinc; f=0.02 ++el: name=Silver; f=0.65 ++el: name=Tin; f=0.29 + +# Material corresponding to H=[ 1807.5;1975.01 ] +AmalgamTooth_33: d=2.13808 g/cm3 ; n=4; ++el: name=Copper; f=0.04 ++el: name=Zinc; f=0.02 ++el: name=Silver; f=0.65 ++el: name=Tin; f=0.29 + +# Material corresponding to H=[ 1975.01;2142.51 ] +AmalgamTooth_34: d=2.23808 g/cm3 ; n=4; ++el: name=Copper; f=0.04 ++el: name=Zinc; f=0.02 ++el: name=Silver; f=0.65 ++el: name=Tin; f=0.29 + +# Material corresponding to H=[ 2142.51;2300 ] +AmalgamTooth_35: d=2.33509 g/cm3 ; n=4; ++el: name=Copper; f=0.04 ++el: name=Zinc; f=0.02 ++el: name=Silver; f=0.65 ++el: name=Tin; f=0.29 + +# Material corresponding to H=[ 2300;2467.5 ] +MetallImplants_36: d=2.4321 g/cm3 ; n=1; ++el: name=Titanium; f=1 + +# Material corresponding to H=[ 2467.5;2635.01 ] +MetallImplants_37: d=2.5321 g/cm3 ; n=1; ++el: name=Titanium; f=1 + +# Material corresponding to H=[ 2635.01;2802.51 ] +MetallImplants_38: d=2.6321 g/cm3 ; n=1; ++el: name=Titanium; f=1 + +# Material corresponding to H=[ 2802.51;2970.02 ] +MetallImplants_39: d=2.7321 g/cm3 ; n=1; ++el: name=Titanium; f=1 + +# Material corresponding to H=[ 2970.02;3000 ] +MetallImplants_40: d=2.79105 g/cm3 ; n=1; ++el: name=Titanium; f=1 + +# Material corresponding to H=[ 3000;3001 ] +MetallImplants_41: d=2.8 g/cm3 ; n=1; ++el: name=Titanium; f=1 + diff --git a/t21_tle/data/phantom.mhd b/t21_tle/data/phantom.mhd new file mode 100644 index 0000000..df280b1 --- /dev/null +++ b/t21_tle/data/phantom.mhd @@ -0,0 +1,13 @@ +ObjectType = Image +NDims = 3 +BinaryData = True +BinaryDataByteOrderMSB = False +CompressedData = False +TransformMatrix = 1 0 0 0 1 0 0 0 1 +Offset = 0 0 0 +CenterOfRotation = 0 0 0 +AnatomicalOrientation = RAI +ElementSpacing = 1.20312 1.20312 6.11538 +DimSize = 256 256 26 +ElementType = MET_SHORT +ElementDataFile = phantom.raw diff --git a/t21_tle/data/phantom.raw b/t21_tle/data/phantom.raw new file mode 120000 index 0000000..9bb19be --- /dev/null +++ b/t21_tle/data/phantom.raw @@ -0,0 +1 @@ +../../gatebenchmarks_data/t21_tle/data/phantom.raw \ No newline at end of file diff --git a/t21_tle/mac/beam-photon.mac b/t21_tle/mac/beam-photon.mac new file mode 100644 index 0000000..a4fd7ae --- /dev/null +++ b/t21_tle/mac/beam-photon.mac @@ -0,0 +1,99 @@ +#===================================================== +# Type of particle & Energy spectrum +#===================================================== +/gate/source/addSource mybeam gps +#/gate/source/mybeam/setIntensity 10.0 +/gate/source/mybeam/gps/particle gamma +/gate/source/mybeam/gps/ene/type User +/gate/source/mybeam/gps/hist/type energy +/gate/source/mybeam/gps/hist/point 0.008 0 +/gate/source/mybeam/gps/hist/point 0.009 1.22018704408557E-016 +/gate/source/mybeam/gps/hist/point 0.01 0.000000000000505397767065596 +/gate/source/mybeam/gps/hist/point 0.011 0.000000000064392397008331 +/gate/source/mybeam/gps/hist/point 0.012 0.00000000460580972493099 +/gate/source/mybeam/gps/hist/point 0.013 0.000000113273532394405 +/gate/source/mybeam/gps/hist/point 0.014 0.00000182260671953458 +/gate/source/mybeam/gps/hist/point 0.015 0.0000145270770463357 +/gate/source/mybeam/gps/hist/point 0.016 0.0000656837267561139 +/gate/source/mybeam/gps/hist/point 0.017 0.000212992406869583 +/gate/source/mybeam/gps/hist/point 0.018 0.000500524192943408 +/gate/source/mybeam/gps/hist/point 0.019 0.00105103241253155 +/gate/source/mybeam/gps/hist/point 0.02 0.00195666587816552 +/gate/source/mybeam/gps/hist/point 0.021 0.00300350864307093 +/gate/source/mybeam/gps/hist/point 0.022 0.00430927258781235 +/gate/source/mybeam/gps/hist/point 0.023 0.00599095697024167 +/gate/source/mybeam/gps/hist/point 0.024 0.00789829883474465 +/gate/source/mybeam/gps/hist/point 0.025 0.0100771547354759 +/gate/source/mybeam/gps/hist/point 0.026 0.0123781130379142 +/gate/source/mybeam/gps/hist/point 0.027 0.0149153846841967 +/gate/source/mybeam/gps/hist/point 0.028 0.0171821944888603 +/gate/source/mybeam/gps/hist/point 0.029 0.0195586716680301 +/gate/source/mybeam/gps/hist/point 0.03 0.0216075337074845 +/gate/source/mybeam/gps/hist/point 0.031 0.023006569921638 +/gate/source/mybeam/gps/hist/point 0.032 0.0241717797973874 +/gate/source/mybeam/gps/hist/point 0.033 0.0253120717917074 +/gate/source/mybeam/gps/hist/point 0.034 0.0262006023122374 +/gate/source/mybeam/gps/hist/point 0.035 0.0269640173283155 +/gate/source/mybeam/gps/hist/point 0.036 0.0275984428506985 +/gate/source/mybeam/gps/hist/point 0.037 0.028141470923155 +/gate/source/mybeam/gps/hist/point 0.038 0.02849835620382 +/gate/source/mybeam/gps/hist/point 0.039 0.0287555221317396 +/gate/source/mybeam/gps/hist/point 0.04 0.028862319884584 +/gate/source/mybeam/gps/hist/point 0.041 0.0286851685739993 +/gate/source/mybeam/gps/hist/point 0.042 0.0281217183607167 +/gate/source/mybeam/gps/hist/point 0.043 0.0275427721904607 +/gate/source/mybeam/gps/hist/point 0.044 0.0272384966402609 +/gate/source/mybeam/gps/hist/point 0.045 0.0268646327647639 +/gate/source/mybeam/gps/hist/point 0.046 0.0263110827462144 +/gate/source/mybeam/gps/hist/point 0.047 0.0257216146298653 +/gate/source/mybeam/gps/hist/point 0.048 0.0251458728704649 +/gate/source/mybeam/gps/hist/point 0.049 0.0245102994994173 +/gate/source/mybeam/gps/hist/point 0.05 0.0237583629699861 +/gate/source/mybeam/gps/hist/point 0.051 0.0228612714114986 +/gate/source/mybeam/gps/hist/point 0.052 0.0220616035084271 +/gate/source/mybeam/gps/hist/point 0.053 0.0212257783724177 +/gate/source/mybeam/gps/hist/point 0.054 0.020351978576418 +/gate/source/mybeam/gps/hist/point 0.055 0.0194578044666199 +/gate/source/mybeam/gps/hist/point 0.056 0.0200749644320047 +/gate/source/mybeam/gps/hist/point 0.057 0.0207400470791406 +/gate/source/mybeam/gps/hist/point 0.058 0.0212810185894062 +/gate/source/mybeam/gps/hist/point 0.059 0.0217563714176734 +/gate/source/mybeam/gps/hist/point 0.06 0.0184378452739822 +/gate/source/mybeam/gps/hist/point 0.061 0.014898693051531 +/gate/source/mybeam/gps/hist/point 0.062 0.0138241153932648 +/gate/source/mybeam/gps/hist/point 0.063 0.0127098413020113 +/gate/source/mybeam/gps/hist/point 0.064 0.0119687136808393 +/gate/source/mybeam/gps/hist/point 0.065 0.0111278667069219 +/gate/source/mybeam/gps/hist/point 0.066 0.0110022251051654 +/gate/source/mybeam/gps/hist/point 0.067 0.0108997317848136 +/gate/source/mybeam/gps/hist/point 0.068 0.00976967477442051 +/gate/source/mybeam/gps/hist/point 0.069 0.00829148482943069 +/gate/source/mybeam/gps/hist/point 0.07 0.00701182487640557 +/gate/source/mybeam/gps/hist/point 0.071 0.00573847809103627 +/gate/source/mybeam/gps/hist/point 0.072 0.00511380928230323 +/gate/source/mybeam/gps/hist/point 0.073 0.00438812022567142 +/gate/source/mybeam/gps/hist/point 0.074 0.00376247574557339 +/gate/source/mybeam/gps/hist/point 0.075 0.00321745806876471 +/gate/source/mybeam/gps/hist/point 0.076 0.0025820664404223 +/gate/source/mybeam/gps/hist/point 0.077 0.0018611696519656 +/gate/source/mybeam/gps/hist/point 0.078 0.0011988561897333 +/gate/source/mybeam/gps/hist/point 0.079 0.00023865399857969 +/gate/source/mybeam/gps/hist/point 0.08 0.0000104276790214574 +/gate/source/mybeam/gps/hist/point 0.081 0 + + +#===================================================== +# Position and shape +#===================================================== +/gate/source/mybeam/gps/position 0 -1000 0 mm +/gate/source/mybeam/gps/pos/type Point +/gate/source/mybeam/gps/ang/type iso +/gate/source/mybeam/gps/ang/mintheta 0. deg +/gate/source/mybeam/gps/ang/maxtheta 10.3168 deg +/gate/source/mybeam/gps/ang/minphi 0. deg +#/gate/source/mybeam/gps/ang/maxphi 5.3684 deg +/gate/source/mybeam/gps/ang/maxphi 360 deg +/gate/source/mybeam/gps/ang/rot1 1 0 0 +/gate/source/mybeam/gps/ang/rot2 0 0 1 +/gate/source/mybeam/gps/pos/rot1 1 0 0 +/gate/source/mybeam/gps/pos/rot2 0 0 1 diff --git a/t21_tle/mac/beam-photonCT.mac b/t21_tle/mac/beam-photonCT.mac new file mode 100644 index 0000000..e6e5aef --- /dev/null +++ b/t21_tle/mac/beam-photonCT.mac @@ -0,0 +1,152 @@ +#===================================================== +# Type of particle & Energy spectrum +#===================================================== +/gate/source/addSource mybeam gps +#/gate/source/mybeam/setIntensity 10.0 +/gate/source/mybeam/gps/particle gamma +/gate/source/mybeam/gps/ene/type User +/gate/source/mybeam/gps/hist/type energy +/gate/source/mybeam/gps/hist/point 0.008 0 +/gate/source/mybeam/gps/hist/point 0.009 0.00000243111325076991 +/gate/source/mybeam/gps/hist/point 0.01 0.0000168741040906372 +/gate/source/mybeam/gps/hist/point 0.011 0.000062709968414845 +/gate/source/mybeam/gps/hist/point 0.012 0.000123888952832963 +/gate/source/mybeam/gps/hist/point 0.013 0.000216192218515745 +/gate/source/mybeam/gps/hist/point 0.014 0.000472978337810378 +/gate/source/mybeam/gps/hist/point 0.015 0.000813988702075342 +/gate/source/mybeam/gps/hist/point 0.016 0.00158546425123098 +/gate/source/mybeam/gps/hist/point 0.017 0.00250323491753462 +/gate/source/mybeam/gps/hist/point 0.018 0.00329921897292835 +/gate/source/mybeam/gps/hist/point 0.019 0.0042156790501742 +/gate/source/mybeam/gps/hist/point 0.02 0.00510079297864447 +/gate/source/mybeam/gps/hist/point 0.021 0.00607485025440779 +/gate/source/mybeam/gps/hist/point 0.022 0.00736006198849925 +/gate/source/mybeam/gps/hist/point 0.023 0.00873362804109373 +/gate/source/mybeam/gps/hist/point 0.024 0.00989543824367597 +/gate/source/mybeam/gps/hist/point 0.025 0.0110279442639547 +/gate/source/mybeam/gps/hist/point 0.026 0.0119765260445497 +/gate/source/mybeam/gps/hist/point 0.027 0.0129308508558473 +/gate/source/mybeam/gps/hist/point 0.028 0.0137310955528222 +/gate/source/mybeam/gps/hist/point 0.029 0.0145522016861103 +/gate/source/mybeam/gps/hist/point 0.03 0.0151379417320528 +/gate/source/mybeam/gps/hist/point 0.031 0.0157520533313813 +/gate/source/mybeam/gps/hist/point 0.032 0.0161516602796775 +/gate/source/mybeam/gps/hist/point 0.033 0.0165556849438988 +/gate/source/mybeam/gps/hist/point 0.034 0.0168193244131653 +/gate/source/mybeam/gps/hist/point 0.035 0.0170472785549031 +/gate/source/mybeam/gps/hist/point 0.036 0.0171097156066452 +/gate/source/mybeam/gps/hist/point 0.037 0.0171567397383818 +/gate/source/mybeam/gps/hist/point 0.038 0.0171392652176112 +/gate/source/mybeam/gps/hist/point 0.039 0.0171143296655004 +/gate/source/mybeam/gps/hist/point 0.04 0.0169615257802227 +/gate/source/mybeam/gps/hist/point 0.041 0.0168179990983878 +/gate/source/mybeam/gps/hist/point 0.042 0.0166455118343769 +/gate/source/mybeam/gps/hist/point 0.043 0.0164815164020887 +/gate/source/mybeam/gps/hist/point 0.044 0.0163533535545269 +/gate/source/mybeam/gps/hist/point 0.045 0.0162053600710343 +/gate/source/mybeam/gps/hist/point 0.046 0.015961502151966 +/gate/source/mybeam/gps/hist/point 0.047 0.0157163189181201 +/gate/source/mybeam/gps/hist/point 0.048 0.0154464455608258 +/gate/source/mybeam/gps/hist/point 0.049 0.0151701419725739 +/gate/source/mybeam/gps/hist/point 0.05 0.0148686574035486 +/gate/source/mybeam/gps/hist/point 0.051 0.0145637368332481 +/gate/source/mybeam/gps/hist/point 0.052 0.0142883658739137 +/gate/source/mybeam/gps/hist/point 0.053 0.0140086262843867 +/gate/source/mybeam/gps/hist/point 0.054 0.0137566692194555 +/gate/source/mybeam/gps/hist/point 0.055 0.0135001962671343 +/gate/source/mybeam/gps/hist/point 0.056 0.0197551911598232 +/gate/source/mybeam/gps/hist/point 0.057 0.0260310474888254 +/gate/source/mybeam/gps/hist/point 0.058 0.0314083403986518 +/gate/source/mybeam/gps/hist/point 0.059 0.036779350334718 +/gate/source/mybeam/gps/hist/point 0.06 0.0254430985849203 +/gate/source/mybeam/gps/hist/point 0.061 0.0140603135607109 +/gate/source/mybeam/gps/hist/point 0.062 0.012527415220416 +/gate/source/mybeam/gps/hist/point 0.063 0.0109846506478882 +/gate/source/mybeam/gps/hist/point 0.064 0.0106465972081492 +/gate/source/mybeam/gps/hist/point 0.065 0.0102876332463644 +/gate/source/mybeam/gps/hist/point 0.066 0.0141914215522229 +/gate/source/mybeam/gps/hist/point 0.067 0.0181040943756642 +/gate/source/mybeam/gps/hist/point 0.068 0.0148963417566795 +/gate/source/mybeam/gps/hist/point 0.069 0.0116331222599679 +/gate/source/mybeam/gps/hist/point 0.07 0.00936266170310275 +/gate/source/mybeam/gps/hist/point 0.071 0.00708046965616969 +/gate/source/mybeam/gps/hist/point 0.072 0.00680819109798287 +/gate/source/mybeam/gps/hist/point 0.073 0.00651927247647799 +/gate/source/mybeam/gps/hist/point 0.074 0.00647632246053902 +/gate/source/mybeam/gps/hist/point 0.075 0.00645187776575319 +/gate/source/mybeam/gps/hist/point 0.076 0.00630540593996818 +/gate/source/mybeam/gps/hist/point 0.077 0.00614037970729753 +/gate/source/mybeam/gps/hist/point 0.078 0.0060136403459782 +/gate/source/mybeam/gps/hist/point 0.079 0.00582235324641915 +/gate/source/mybeam/gps/hist/point 0.08 0.00579835032322583 +/gate/source/mybeam/gps/hist/point 0.081 0.00569944257223495 +/gate/source/mybeam/gps/hist/point 0.082 0.00540051046129972 +/gate/source/mybeam/gps/hist/point 0.083 0.00513515299139565 +/gate/source/mybeam/gps/hist/point 0.084 0.00501582557568406 +/gate/source/mybeam/gps/hist/point 0.085 0.00489785292618953 +/gate/source/mybeam/gps/hist/point 0.086 0.00479280945863594 +/gate/source/mybeam/gps/hist/point 0.087 0.00468705424796107 +/gate/source/mybeam/gps/hist/point 0.088 0.00454654143010176 +/gate/source/mybeam/gps/hist/point 0.089 0.00440563592638245 +/gate/source/mybeam/gps/hist/point 0.09 0.0043473957047692 +/gate/source/mybeam/gps/hist/point 0.091 0.00429148705544979 +/gate/source/mybeam/gps/hist/point 0.092 0.00418611471348844 +/gate/source/mybeam/gps/hist/point 0.093 0.00408896423172111 +/gate/source/mybeam/gps/hist/point 0.094 0.00396341765370157 +/gate/source/mybeam/gps/hist/point 0.095 0.00383736549263726 +/gate/source/mybeam/gps/hist/point 0.096 0.00368201896641594 +/gate/source/mybeam/gps/hist/point 0.097 0.00352505751959532 +/gate/source/mybeam/gps/hist/point 0.098 0.00348458142457444 +/gate/source/mybeam/gps/hist/point 0.099 0.00344352120933678 +/gate/source/mybeam/gps/hist/point 0.1 0.00327984974288319 +/gate/source/mybeam/gps/hist/point 0.101 0.00311829387150041 +/gate/source/mybeam/gps/hist/point 0.102 0.00310819693632481 +/gate/source/mybeam/gps/hist/point 0.103 0.00310632186134325 +/gate/source/mybeam/gps/hist/point 0.104 0.00293639196399569 +/gate/source/mybeam/gps/hist/point 0.105 0.00276600556933586 +/gate/source/mybeam/gps/hist/point 0.106 0.00264406679265519 +/gate/source/mybeam/gps/hist/point 0.107 0.00252526459428138 +/gate/source/mybeam/gps/hist/point 0.108 0.00243171209670643 +/gate/source/mybeam/gps/hist/point 0.109 0.00233198952255602 +/gate/source/mybeam/gps/hist/point 0.11 0.00223253201136111 +/gate/source/mybeam/gps/hist/point 0.111 0.00214196883489552 +/gate/source/mybeam/gps/hist/point 0.112 0.00211559016224911 +/gate/source/mybeam/gps/hist/point 0.113 0.00209457165159189 +/gate/source/mybeam/gps/hist/point 0.114 0.00197457176134504 +/gate/source/mybeam/gps/hist/point 0.115 0.00184646290808891 +/gate/source/mybeam/gps/hist/point 0.116 0.00178871354380068 +/gate/source/mybeam/gps/hist/point 0.117 0.00172842153856887 +/gate/source/mybeam/gps/hist/point 0.118 0.00166106609643007 +/gate/source/mybeam/gps/hist/point 0.119 0.00158532681117998 +/gate/source/mybeam/gps/hist/point 0.12 0.00145023305818886 +/gate/source/mybeam/gps/hist/point 0.121 0.00129877903055492 +/gate/source/mybeam/gps/hist/point 0.122 0.00125330109939212 +/gate/source/mybeam/gps/hist/point 0.123 0.00118031552373537 +/gate/source/mybeam/gps/hist/point 0.124 0.00115900740725639 +/gate/source/mybeam/gps/hist/point 0.125 0.00113971180580998 +/gate/source/mybeam/gps/hist/point 0.126 0.00101418977065669 +/gate/source/mybeam/gps/hist/point 0.127 0.000888294683936384 +/gate/source/mybeam/gps/hist/point 0.128 0.0000692850022834449 +/gate/source/mybeam/gps/hist/point 0.129 0.000802679349306955 +/gate/source/mybeam/gps/hist/point 0.13 0.0000176175075093751 +/gate/source/mybeam/gps/hist/point 0.131 0.000628989484849764 +/gate/source/mybeam/gps/hist/point 0.132 0.000283663012125928 +/gate/source/mybeam/gps/hist/point 0.133 0.000250387303205714 +/gate/source/mybeam/gps/hist/point 0.134 0.000213110616229289 +/gate/source/mybeam/gps/hist/point 0.135 0.000175829020679614 +/gate/source/mybeam/gps/hist/point 0.136 0.000145202468742517 +/gate/source/mybeam/gps/hist/point 0.137 0.000114571989946819 +/gate/source/mybeam/gps/hist/point 0.138 0.0000759425001826476 +/gate/source/mybeam/gps/hist/point 0.139 0.0000373078073308311 +/gate/source/mybeam/gps/hist/point 0.14 0 + +#===================================================== +# Position and shape +#===================================================== +/gate/source/mybeam/gps/type Surface +/gate/source/mybeam/gps/shape Cylinder +/gate/source/mybeam/gps/radius 1000.0 mm +/gate/source/mybeam/gps/halfz 80.0 mm +/gate/source/mybeam/gps/centre 0.0 0.0 0.0 mm + + diff --git a/t21_tle/mac/detectors.mac b/t21_tle/mac/detectors.mac new file mode 100644 index 0000000..5e78bd0 --- /dev/null +++ b/t21_tle/mac/detectors.mac @@ -0,0 +1,23 @@ +# ---------------------------------------------------- +# the following actor regularly store the current number of +# event/track/step in a file + +/gate/actor/addActor SimulationStatisticActor stat +/gate/actor/stat/saveEveryNSeconds 600 + + +/gate/actor/addActor TLEDoseActor doseDistribution +/gate/actor/doseDistribution/attachTo thorax +/gate/actor/doseDistribution/stepHitType random +/gate/actor/doseDistribution/setPosition 0 0 0 mm +/gate/actor/doseDistribution/setResolution 256 256 26 +/gate/actor/doseDistribution/enableDose true +/gate/actor/doseDistribution/save output/TLE-thorax.mhd + +/gate/actor/addActor DoseActor doseDistribution2 +/gate/actor/doseDistribution2/attachTo thorax +/gate/actor/doseDistribution2/stepHitType random +/gate/actor/doseDistribution2/setPosition 0 0 0 mm +/gate/actor/doseDistribution2/setResolution 256 256 26 +/gate/actor/doseDistribution2/enableDose true +/gate/actor/doseDistribution2/save output/Dose-thorax.mhd \ No newline at end of file diff --git a/t21_tle/mac/main-setle.mac b/t21_tle/mac/main-setle.mac new file mode 100644 index 0000000..29959e5 --- /dev/null +++ b/t21_tle/mac/main-setle.mac @@ -0,0 +1,107 @@ + +#===================================================== +# VERBOSITY +#===================================================== + +/control/execute mac/verbose.mac + +#===================================================== +# VISUALISATION +#===================================================== + +#/control/execute mac/visu.mac + +#===================================================== +# GEOMETRY +#===================================================== + +/gate/geometry/setMaterialDatabase data/GateMaterials.db + +# World +/gate/world/setMaterial Air +/gate/world/geometry/setXLength 1.0 m +/gate/world/geometry/setYLength 2.0 m +/gate/world/geometry/setZLength 1.0 m + +# Generate materials from Hounsfield units +/gate/HounsfieldMaterialGenerator/SetMaterialTable data/SimpleMaterialsTable.txt +/gate/HounsfieldMaterialGenerator/SetMaterialTable data/Schneider2000MaterialsTable.txt +/gate/HounsfieldMaterialGenerator/SetDensityTable data/Schneider2000DensitiesTable.txt +/gate/HounsfieldMaterialGenerator/SetDensityTolerance 0.1 g/cm3 +/gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename data/phantom-HUmaterials.db +/gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename data/phantom-HU2mat.txt +/gate/HounsfieldMaterialGenerator/Generate + +/gate/world/daughters/name phantom +/gate/world/daughters/insert ImageNestedParametrisedVolume + +/gate/geometry/setMaterialDatabase data/phantom-HUmaterials.db +/gate/phantom/geometry/setHUToMaterialFile data/phantom-HU2mat.txt +/gate/phantom/geometry/setImage data/phantom.mhd + +/gate/phantom/geometry/TranslateTheImageAtThisIsoCenter 154 145 80 mm + + +#===================================================== +# PHYSICS +#===================================================== + +/gate/physics/addPhysicsList emlivermore +/gate/physics/addAtomDeexcitation + +/gate/physics/Gamma/SetCutInRegion phantom 0.1 mm +/gate/physics/Electron/SetCutInRegion phantom 0.1 m +/gate/physics/Positron/SetCutInRegion phantom 0.1 m + +/gate/physics/SetMinKineticEnergyInRegion phantom 100 keV +/gate/physics/ActivateSpecialCuts e- + +#===================================================== +# DETECTORS +#===================================================== + +# ---------------------------------------------------- +# the following actor regularly store the current number of +# event/track/step in a file + +/gate/actor/addActor SimulationStatisticActor stat +/gate/actor/stat/saveEveryNSeconds 600 + +/gate/actor/addActor SETLEDoseActor setle +/gate/actor/setle/attachTo phantom +/gate/actor/setle/setVoxelSize 2 2 2 mm +/gate/actor/setle/enableHybridino true +/gate/actor/setle/setPrimaryMultiplicity 200 +/gate/actor/setle/setSecondaryMultiplicity 400 +/gate/actor/setle/enableDose true +/gate/actor/setle/enableUncertaintyDose true +/gate/actor/setle/enablePrimaryDose true +/gate/actor/setle/enablePrimaryUncertaintyDose false +/gate/actor/setle/enableSecondaryDose true +/gate/actor/setle/enableSecondaryUncertaintyDose false +/gate/actor/setle/save output/dose-setle.mhd + +# Set the names of the outputs +/gate/actor/stat/save output/stat.txt + +#===================================================== +# INITIALISATION +#===================================================== + +/gate/run/initialize + +#===================================================== +# BEAMS +#===================================================== + +/control/execute mac/beam-photon.mac + +#===================================================== +# START BEAMS +#===================================================== + +/gate/random/setEngineName MersenneTwister +/gate/random/setEngineSeed auto + +/gate/application/setTotalNumberOfPrimaries 500 +/gate/application/start diff --git a/t21_tle/mac/main.mac b/t21_tle/mac/main.mac new file mode 100644 index 0000000..0ec43ed --- /dev/null +++ b/t21_tle/mac/main.mac @@ -0,0 +1,109 @@ + +#===================================================== +# VERBOSITY +#===================================================== + +/control/execute mac/verbose.mac + +#===================================================== +# VISUALISATION +#===================================================== + +#/control/execute mac/visu.mac + +#===================================================== +# GEOMETRY +#===================================================== + +/gate/geometry/setMaterialDatabase data/GateMaterials.db + +# World +/gate/world/setMaterial Air +/gate/world/geometry/setXLength 1.0 m +/gate/world/geometry/setYLength 2.0 m +/gate/world/geometry/setZLength 1.0 m + +# Generate materials from Hounsfield units +/gate/HounsfieldMaterialGenerator/SetMaterialTable data/SimpleMaterialsTable.txt +/gate/HounsfieldMaterialGenerator/SetMaterialTable data/Schneider2000MaterialsTable.txt +/gate/HounsfieldMaterialGenerator/SetDensityTable data/Schneider2000DensitiesTable.txt +/gate/HounsfieldMaterialGenerator/SetDensityTolerance 0.1 g/cm3 +/gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename data/phantom-HUmaterials.db +/gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename data/phantom-HU2mat.txt +/gate/HounsfieldMaterialGenerator/Generate + +/gate/world/daughters/name phantom +/gate/world/daughters/insert ImageNestedParametrisedVolume + +/gate/geometry/setMaterialDatabase data/phantom-HUmaterials.db +/gate/phantom/geometry/setHUToMaterialFile data/phantom-HU2mat.txt +/gate/phantom/geometry/setImage data/phantom.mhd + +/gate/phantom/geometry/TranslateTheImageAtThisIsoCenter 154 145 80 mm + + +#===================================================== +# PHYSICS +#===================================================== + +/gate/physics/addPhysicsList emlivermore +/gate/physics/addAtomDeexcitation + +/gate/physics/Gamma/SetCutInRegion phantom 0.1 mm +/gate/physics/Electron/SetCutInRegion phantom 0.1 m +/gate/physics/Positron/SetCutInRegion phantom 0.1 m + +/gate/physics/SetMinKineticEnergyInRegion phantom 100 keV +/gate/physics/ActivateSpecialCuts e- + +#===================================================== +# DETECTORS +#===================================================== + +# ---------------------------------------------------- +# the following actor regularly store the current number of +# event/track/step in a file + +/gate/actor/addActor SimulationStatisticActor stat +/gate/actor/stat/saveEveryNSeconds 600 + +/gate/actor/addActor DoseActor dose +/gate/actor/dose/attachTo phantom +/gate/actor/dose/stepHitType random +/gate/actor/dose/setVoxelSize 2 2 2 mm +/gate/actor/dose/enableDose true +/gate/actor/dose/enableUncertaintyDose true +/gate/actor/dose/save output/output.mhd + +/gate/actor/addActor TLEDoseActor tle +/gate/actor/tle/attachTo phantom +/gate/actor/tle/stepHitType random +/gate/actor/tle/setVoxelSize 2 2 2 mm +/gate/actor/tle/enableDose true +/gate/actor/tle/enableUncertaintyDose true +/gate/actor/tle/save output/output-tle.mhd + +# Set the names of the outputs +/gate/actor/stat/save output/stat.txt + +#===================================================== +# INITIALISATION +#===================================================== + +/gate/run/initialize + +#===================================================== +# BEAMS +#===================================================== + +/control/execute mac/beam-photon.mac + +#===================================================== +# START BEAMS +#===================================================== + +/gate/random/setEngineName MersenneTwister +/gate/random/setEngineSeed auto + +/gate/application/setTotalNumberOfPrimaries 5000000 +/gate/application/start diff --git a/t21_tle/mac/template.mac b/t21_tle/mac/template.mac new file mode 100644 index 0000000..4c8138a --- /dev/null +++ b/t21_tle/mac/template.mac @@ -0,0 +1,8 @@ + +/gate/random/setEngineSeed auto +/gate/application/noGlobalOutput + +#/gate/application/enableTrackTimeStudy * +/gate/application/setTotalNumberOfPrimaries +/gate/application/start + diff --git a/t21_tle/mac/verbose.mac b/t21_tle/mac/verbose.mac new file mode 100644 index 0000000..fe48693 --- /dev/null +++ b/t21_tle/mac/verbose.mac @@ -0,0 +1,13 @@ +/gate/verbose Physic 1 +/gate/verbose Cuts 1 +/gate/verbose SD 0 +/gate/verbose Actions 0 +/gate/verbose Actor 0 +/gate/verbose Step 0 +/gate/verbose Error 1 +/gate/verbose Warning 1 +/gate/verbose Output 0 +/gate/verbose Beam 0 +/gate/verbose Volume 0 +/gate/verbose Image 0 +/gate/verbose Geometry 0 diff --git a/t21_tle/mac/visu.mac b/t21_tle/mac/visu.mac new file mode 100644 index 0000000..e2265be --- /dev/null +++ b/t21_tle/mac/visu.mac @@ -0,0 +1,15 @@ + +/vis/open OGLIQt +/vis/drawVolume +/vis/viewer/flush +/tracking/storeTrajectory 1 +/vis/scene/add/trajectories +/vis/scene/endOfEventAction accumulate + +/vis/modeling/trajectories/create/drawByParticleID + +/vis/scene/add/axes 0 0 0 500 mm +/vis/scene/add/text 10 0 0 cm 20 0 0 X +/vis/scene/add/text 0 10 0 cm 20 0 0 Y +/vis/scene/add/text 0 0 10 cm 20 0 0 Z +/vis/viewer/set/auxiliaryEdge true diff --git a/t21_tle/output/.gitkeep b/t21_tle/output/.gitkeep new file mode 100644 index 0000000..e69de29 diff --git a/t21_tle/runTest.sh b/t21_tle/runTest.sh new file mode 100755 index 0000000..e1ec8bd --- /dev/null +++ b/t21_tle/runTest.sh @@ -0,0 +1,2 @@ +Gate mac/main.mac +