catsim-test 0.0.5__cp313-cp313-manylinux_2_31_x86_64.whl

gecatsim/examples/cfg/Protocol_Sample_Helical.cfg

@@ -0,0 +1,33 @@
+protocol.scanTypes = [1, 1, 1, 1]           # flags for airscan, offset scan, phantom scan, prep
+
+# Table and gantry
+protocol.scanTrajectory = "Gantry_Helical"  # name of the function that defines the scanning trajectory and model
+protocol.viewsPerRotation = 500             # total numbers of view per rotation
+protocol.viewCount = 2500                   # total number of views in scan
+protocol.startViewId = 0                    # index of the first view in the scan
+protocol.stopViewId = protocol.startViewId+protocol.viewCount-1 # index of the last view in the scan
+protocol.airViewCount = 1                   # number of views averaged for air scan
+protocol.offsetViewCount = 1                # number of views averaged for offset scan
+protocol.rotationTime = 1.0                 # gantry rotation period (in seconds)
+protocol.rotationDirection = 1              # gantry rotation direction (1=CW, -1 CCW, seen from table foot-end)
+protocol.startAngle = 0                     # relative to vertical y-axis (n degrees)
+protocol.tableSpeed = 5.116                     # speed of table translation along positive z-axis (in mm/sec)
+protocol.startZ = -12.79                         # start z-position of table
+protocol.tiltAngle = 0                      # gantry tilt angle towards negative z-axis (in degrees)
+protocol.wobbleDistance = 0.0               # focalspot wobble distance
+protocol.focalspotOffset = [0, 0, 0]        # focalspot position offset
+
+# X-ray tube technique and filtration
+protocol.mA = 200                           # tube current (in mA)
+protocol.spectrumCallback = "Spectrum"      # name of function that reads and models the X-ray spectrum
+protocol.spectrumFilename = "xcist_kVp120_tar7_bin1.dat" # name of the spectrum file
+protocol.spectrumScaling = 1                # scaling factor, works for both mono- and poly-chromatic spectra
+protocol.spectrumUnit_mm = 1;               # Is the spectrum file in units of photons/sec/mm^2/<current>?
+protocol.spectrumUnit_mA = 1;               # Is the spectrum file in units of photons/sec/<area>/mA?
+protocol.bowtie = "large.txt"               # name of the bowtie file (or [])
+protocol.filterCallback = "Xray_Filter"     # name of function to compute additional filtration
+protocol.flatFilter = ['Al',3.0]            # additional filtration - materials and thicknesses (in mm)
+protocol.dutyRatio = 1.0                    # tube ON time fraction (for pulsed tubes)
+
+# Pre- and Post-log processes
+protocol.maxPrep = 9                        # set the upper limit of prep, can be used for simple low signal correction

gecatsim/examples/cfg/Protocol_Sample_axial.cfg

@@ -0,0 +1,33 @@
+protocol.scanTypes = [1, 1, 1, 1]           # flags for airscan, offset scan, phantom scan, prep
+
+# Table and gantry
+protocol.scanTrajectory = "Gantry_Helical"  # name of the function that defines the scanning trajectory and model
+protocol.viewsPerRotation = 1000            # total numbers of view per rotation
+protocol.viewCount = 1000                   # total number of views in scan
+protocol.startViewId = 0                    # index of the first view in the scan
+protocol.stopViewId = protocol.startViewId+protocol.viewCount-1 # index of the last view in the scan
+protocol.airViewCount = 1                   # number of views averaged for air scan
+protocol.offsetViewCount = 1                # number of views averaged for offset scan
+protocol.rotationTime = 1.0                 # gantry rotation period (in seconds)
+protocol.rotationDirection = 1              # gantry rotation direction (1=CW, -1 CCW, seen from table foot-end)
+protocol.startAngle = 0                     # relative to vertical y-axis (n degrees)
+protocol.tableSpeed = 0                     # speed of table translation along positive z-axis (in mm/sec)
+protocol.startZ = 0                         # start z-position of table
+protocol.tiltAngle = 0                      # gantry tilt angle towards negative z-axis (in degrees)
+protocol.wobbleDistance = 0.0               # focalspot wobble distance
+protocol.focalspotOffset = [0, 0, 0]        # focalspot position offset
+
+# X-ray tube technique and filtration
+protocol.mA = 200                           # tube current (in mA)
+protocol.spectrumCallback = "Spectrum"      # name of function that reads and models the X-ray spectrum
+protocol.spectrumFilename = "tungsten_tar7.0_120_filt.dat" # name of the spectrum file
+protocol.spectrumScaling = 1                # scaling factor, works for both mono- and poly-chromatic spectra
+protocol.spectrumUnit_mm = 0;               # Is the spectrum file in units of photons/sec/mm^2/<current>?
+protocol.spectrumUnit_mA = 1;               # Is the spectrum file in units of photons/sec/<area>/mA?
+protocol.bowtie = "large.txt"               # name of the bowtie file (or [])
+protocol.filterCallback = "Xray_Filter"     # name of function to compute additional filtration
+protocol.flatFilter = ['Al',3.0]            # additional filtration - materials and thicknesses (in mm)
+protocol.dutyRatio = 1.0                    # tube ON time fraction (for pulsed tubes)
+
+# Pre- and Post-log processes
+protocol.maxPrep = 9                        # set the upper limit of prep, can be used for simple low signal correction

gecatsim/examples/cfg/Recon_Sample_2d.cfg

@@ -0,0 +1,6 @@
+# Reconstruction
+recon.fov = 500.0           # diameter of the reconstruction field-of-view (in mm)
+recon.imageSize = 512       # number of columns and rows to be reconstructed (square)
+recon.sliceCount = 1        # number of slices to reconstruct
+recon.sliceThickness = 1.0  # reconstruction inter-slice interval (in mm)
+recon.centerOffset = [0.0,0.0,0.0] # reconstruction offset relative to center of rotation (in mm)

gecatsim/examples/cfg/Recon_Sample_Helical.cfg

@@ -0,0 +1,8 @@
+# Reconstruction
+recon.fov = 500.0           # diameter of the reconstruction field-of-view (in mm)
+recon.imageSize = 512       # number of columns and rows to be reconstructed (square)
+recon.sliceCount = 50        # number of slices to reconstruct
+recon.sliceThickness = 0.568  # reconstruction inter-slice interval (in mm)
+recon.reconType = 'helical_equiAngle'       # Name of the recon function to call
+recon.kernelType = 'standard'
+recon.centerOffset = [0.0,0.0,0.0] # reconstruction offset relative to center of rotation (in mm)

gecatsim/examples/cfg/Scanner_PCCT.cfg

@@ -0,0 +1,30 @@
+# Scanner geometry
+scanner.detectorCallback = "Detector_ThirdgenCurved" # name of function that defines the detector shape and model
+scanner.sid = 540.0                         # source-to-iso distance (in mm)
+scanner.sdd = 950.0                         # source-to-detector distance (in mm)
+scanner.detectorColsPerMod = 100            # number of detector columns per module
+scanner.detectorRowsPerMod = 2              # number of detector rows per module
+scanner.detectorColOffset = -10.75          # detector column offset relative to centered position (in detector columns)
+scanner.detectorRowOffset = 0.0             # detector row offset relative to centered position (in detector rows)
+scanner.detectorColSize = 0.25              # detector column pitch or size (in mm)
+scanner.detectorRowSize = 0.25              # detector row pitch or size (in mm)
+scanner.detectorColCount = 3600             # total number of detector columns
+scanner.detectorRowCount = scanner.detectorRowsPerMod     # total number of detector rows
+scanner.detectorPrefilter = []              # detector filter
+
+# X-ray tube
+scanner.focalspotCallback = "SetFocalspot"  # name of function that defines the focal spot shape and model
+scanner.focalspotShape = "Uniform"          # Parameterize the model
+scanner.targetAngle = 10.5                  # target angle relative to scanner XY-plane (in degrees)
+scanner.focalspotWidth = 1.0
+scanner.focalspotLength = 1.0
+
+# Detector
+scanner.detectorMaterial = "CZT"            # detector sensor material
+scanner.detectorDepth = 1.6                 # detector sensor depth (in mm)
+scanner.detectorColFillFraction = 0.88      # active fraction of each detector cell in the column direction
+scanner.detectorRowFillFraction = 0.88      # active fraction of each detector cell in the row direction
+scanner.detectionCallback = "Detection_PC"  # name of function that defines the detection process (conversion from X-rays to detector signal)
+scanner.detectionResponseFilename = 'PC_spectral_response_CZT0.25x0.25x1.6.mat'   # name of the response data file
+scanner.detectorBinThreshold = [24, 30, 40, 60, 80, 100, 160]  # energy thresholds (keV), n bins has n+1 thresholds; the first and last are the min and max energy thresholds.
+scanner.detectorSumBins = 0                 # 1: sum all bins (gray scale output), data dim [view row col]; 0: output multiple bins [view row col bin]

gecatsim/examples/cfg/Scanner_Sample_generic.cfg

@@ -0,0 +1,29 @@
+# Scanner geometry
+scanner.detectorCallback = "Detector_ThirdgenCurved" # name of function that defines the detector shape and model
+scanner.sid = 540.0                         # source-to-iso distance (in mm)
+scanner.sdd = 950.0                         # source-to-detector distance (in mm)
+scanner.detectorColsPerMod = 1              # number of detector columns per module
+scanner.detectorRowsPerMod = 16             # number of detector rows per module
+scanner.detectorColOffset = 0.25             # detector column offset relative to centered position (in detector columns)
+scanner.detectorRowOffset = 0.0             # detector row offset relative to centered position (in detector rows)
+scanner.detectorColSize = 1.0               # detector column pitch or size (in mm)
+scanner.detectorRowSize = 1.0               # detector row pitch or size (in mm)
+scanner.detectorColCount = 900              # total number of detector columns
+scanner.detectorRowCount = scanner.detectorRowsPerMod     # total number of detector rows
+scanner.detectorPrefilter = ['graphite', 1.0]  # detector filter
+
+# X-ray tube
+scanner.focalspotCallback = "SetFocalspot"  # name of function that defines the focal spot shape and model
+scanner.focalspotShape = "Uniform"          # Parameterize the model
+scanner.targetAngle = 7.0                   # target angle relative to scanner XY-plane (in degrees)
+scanner.focalspotWidth = 1.0
+scanner.focalspotLength = 1.0
+
+# Detector
+scanner.detectorMaterial = "Lumex"          # detector sensor material
+scanner.detectorDepth = 3.0                 # detector sensor depth (in mm)
+scanner.detectionCallback = "Detection_EI"  # name of function that defines the detection process (conversion from X-rays to detector signal)
+scanner.detectionGain = 15.0                # factor to convert energy to electrons (electrons / keV)
+scanner.detectorColFillFraction = 0.9       # active fraction of each detector cell in the column direction
+scanner.detectorRowFillFraction = 0.9       # active fraction of each detector cell in the row direction
+scanner.eNoise = 5000.0                     # standard deviation of Gaussian electronic noise (in electrons)

gecatsim/examples/vct_examples/Phantom_Sample_Analytic.cfg

@@ -0,0 +1,6 @@
+# Phantom
+phantom.callback = "Phantom_Analytic"           # name of function that reads and models phantom
+phantom.projectorCallback = "C_Projector_Analytic" # name of function that performs projection through phantom
+phantom.filename = 'CTDI_16cm_WaterAirPEBoneChambers.ppm'       # phantom filename
+phantom.centerOffset = [0.0, 0.0, 0.0]      # offset of phantom center relative to origin (in mm)
+phantom.scale = 1.0                         # re-scale the size of phantom

gecatsim/examples/vct_examples/Physics_Sample.cfg

@@ -0,0 +1,41 @@
+# Geometric and energy sampling
+physics.energyCount = 20
+physics.monochromatic = -1
+physics.colSampleCount = 2
+physics.rowSampleCount = 2
+physics.srcXSampleCount = 2
+physics.srcYSampleCount = 2
+physics.viewSampleCount = 2
+
+# Flags to determine what has to be recalculated each view
+physics.recalcDet = 0
+physics.recalcSrc = 0
+physics.recalcRayAngle = 0
+physics.recalcSpec = 0
+physics.recalcFilt = 0
+physics.recalcFlux = 0
+physics.recalcPht = 0
+physics.recalcDet = 0
+
+# Noise on/off settings
+physics.enableQuantumNoise = 1
+physics.enableElectronicNoise = 1
+
+# Internal physics models
+physics.rayAngleCallback = "Detector_RayAngles_2D"
+physics.fluxCallback = "Detection_Flux"
+physics.scatterCallback = ""
+physics.prefilterCallback = "Detection_prefilter"
+physics.crosstalkCallback = "CalcCrossTalk"
+physics.col_crosstalk = 0.01
+physics.row_crosstalk = 0.01
+physics.lagCallback = "Detection_Lag"
+physics.lag_taus = [1., 10.]
+physics.lag_alphas = [0.9, 0.1]
+physics.opticalCrosstalkCallback = "CalcOptCrossTalk"
+physics.col_crosstalk_opt = 0.01
+physics.row_crosstalk_opt = 0.01
+physics.DASCallback = "Detection_DAS"
+
+# I/O preferences
+physics.outputCallback = "WriteRawView"

gecatsim/examples/vct_examples/Protocol_Sample_axial.cfg

@@ -0,0 +1,33 @@
+protocol.scanTypes = [1, 1, 1, 1]           # flags for airscan, offset scan, phantom scan, prep
+
+# Table and gantry
+protocol.scanTrajectory = "Gantry_Helical"  # name of the function that defines the scanning trajectory and model
+protocol.viewsPerRotation = 1000            # total numbers of view per rotation
+protocol.viewCount = 1000                   # total number of views in scan
+protocol.startViewId = 0                    # index of the first view in the scan
+protocol.stopViewId = protocol.startViewId+protocol.viewCount-1 # index of the last view in the scan
+protocol.airViewCount = 1                   # number of views averaged for air scan
+protocol.offsetViewCount = 1                # number of views averaged for offset scan
+protocol.rotationTime = 1.0                 # gantry rotation period (in seconds)
+protocol.rotationDirection = 1              # gantry rotation direction (1=CW, -1 CCW, seen from table foot-end)
+protocol.startAngle = 0                     # relative to vertical y-axis (n degrees)
+protocol.tableSpeed = 0                     # speed of table translation along positive z-axis (in mm/sec)
+protocol.startZ = 0                         # start z-position of table
+protocol.tiltAngle = 0                      # gantry tilt angle towards negative z-axis (in degrees)
+protocol.wobbleDistance = 0.0               # focalspot wobble distance
+protocol.focalspotOffset = [0, 0, 0]        # focalspot position offset
+
+# X-ray tube technique and filtration
+protocol.mA = 200                           # tube current (in mA)
+protocol.spectrumCallback = "Spectrum"      # name of function that reads and models the X-ray spectrum
+protocol.spectrumFilename = "xcist_kVp80_tar7_bin1.dat" # name of the spectrum file
+protocol.spectrumScaling = 1                # scaling factor, works for both mono- and poly-chromatic spectra
+protocol.spectrumUnit_mm = 1;               # Is the spectrum file in units of photons/sec/mm^2/<current>?
+protocol.spectrumUnit_mA = 1;               # Is the spectrum file in units of photons/sec/<area>/mA?
+protocol.bowtie = "large.txt"               # name of the bowtie file (or [])
+protocol.filterCallback = "Xray_Filter"     # name of function to compute additional filtration
+protocol.flatFilter = ['Al',3.0]            # additional filtration - materials and thicknesses (in mm)
+protocol.dutyRatio = 1.0                    # tube ON time fraction (for pulsed tubes)
+
+# Pre- and Post-log processes
+protocol.maxPrep = 9                        # set the upper limit of prep, can be used for simple low signal correction

gecatsim/examples/vct_examples/Recon_Sample_2d.cfg

@@ -0,0 +1,6 @@
+# Reconstruction
+recon.fov = 500.0           # diameter of the reconstruction field-of-view (in mm)
+recon.imageSize = 512       # number of columns and rows to be reconstructed (square)
+recon.sliceCount = 1        # number of slices to reconstruct
+recon.sliceThickness = 1.0  # reconstruction inter-slice interval (in mm)
+recon.centerOffset = [0.0,0.0,0.0] # reconstruction offset relative to center of rotation (in mm)

gecatsim/examples/vct_examples/Scanner_Sample_generic.cfg

@@ -0,0 +1,31 @@
+# Scanner geometry
+scanner.detectorCallback = "Detector_ThirdgenCurved" # name of function that defines the detector shape and model
+scanner.sid = 540.0                         # source-to-iso distance (in mm)
+scanner.sdd = 950.0                         # source-to-detector distance (in mm)
+scanner.detectorColsPerMod = 1              # number of detector columns per module
+scanner.detectorRowsPerMod = 16             # number of detector rows per module
+scanner.detectorColOffset = 0.25             # detector column offset relative to centered position (in detector columns)
+scanner.detectorRowOffset = 0.0             # detector row offset relative to centered position (in detector rows)
+scanner.detectorColSize = 1.0               # detector column pitch or size (in mm)
+scanner.detectorRowSize = 1.0               # detector row pitch or size (in mm)
+scanner.detectorColCount = 900              # total number of detector columns
+scanner.detectorRowCount = scanner.detectorRowsPerMod     # total number of detector rows
+scanner.detectorPrefilter = ['graphite', 1.0]  # detector filter
+
+# X-ray tube
+scanner.focalspotCallback = "SetFocalspot"  # name of function that defines the focal spot shape and model
+scanner.focalspotData = "vct_small_fs.npz"
+scanner.targetAngle = 7.0                   # target angle relative to scanner XY-plane (in degrees)
+scanner.focalspotLength = 0.5646668803041932
+scanner.focalspotWidth = 0.8974411689872975
+scanner.focalspotWidthThreshold =0.5
+scanner.focalspotLengthThreshold =0.5
+
+# Detector
+scanner.detectorMaterial = "Lumex"          # detector sensor material
+scanner.detectorDepth = 3.0                 # detector sensor depth (in mm)
+scanner.detectionCallback = "Detection_EI"  # name of function that defines the detection process (conversion from X-rays to detector signal)
+scanner.detectionGain = 0.0671               # factor to convert energy to electrons (electrons / keV)
+scanner.detectorColFillFraction = 0.9       # active fraction of each detector cell in the column direction
+scanner.detectorRowFillFraction = 0.9       # active fraction of each detector cell in the row direction
+scanner.eNoise = 15.35                     # standard deviation of Gaussian electronic noise (in electrons)

gecatsim/examples/vct_examples/read.md

@@ -0,0 +1 @@
+note that Gain factor is kVp dependent: 80kVp-0.0671, 100kVp-0.0673, 120kVp-0.0684, 140kVp-0.0709

gecatsim/focal_spot/readme.md

@@ -0,0 +1,31 @@
+This folder contains focal spot images in `.npz` format for XCIST. Currently it only has VCT focal spot images.
+
+The `.npz` format should have following keys:
+* 'data': 2d focal spot images
+* 'pixsize_x': pixel size in mm in the X direction
+* 'pixsize_z': pixel size in mm in the Z direction
+* [optional]
+
+
+Example usages of the new focal spot model are:
+1. Uniform:
+`scanner.focalspotCallback = "SetFocalspot"
+scanner.focalspotShape = "Uniform"
+scanner.focalspotWidth = 1.0
+scanner.focalspotLength = 1.0`
+2. Gaussian
+`scanner.focalspotCallback = "SetFocalspot"
+scanner.focalspotShape = "Gaussian"
+scanner.focalspotWidth = 1.0
+scanner.focalspotLength = 1.0
+scanner.focalspotWidthThreshold = 0.5
+scanner.focalspotLengthThreshold = 0.5`
+3. Customized focal spot image
+`scanner.focalspotCallback = "SetFocalspot"
+scanner.focalspotData = "xxx/xxx/xxx.npz"
+scanner.focalspotWidth = 1.0
+scanner.focalspotLength = 1.0
+scanner.focalspotPixSizeX = 0.04
+scanner.focalspotPixSizeZ = 0.04
+#scanner.focalspotWidthThreshold = 0.5
+scanner.focalspotLengthThreshold = 0.5`

gecatsim/material/Ac

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Ac (actinium)
+
+# Number of elements:
+1
+
+# Density:
+10.070000
+
+# Atomic number(s) and corresponding mass fraction(s):
+89    1.000000    # actinium

gecatsim/material/Ag

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Ag (silver)
+
+# Number of elements:
+1
+
+# Density:
+10.50
+
+# Atomic number(s) and corresponding mass fraction(s):
+47    1.00    # silver

gecatsim/material/Al

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Al (aluminum)
+
+# Number of elements:
+1
+
+# Density:
+2.699
+
+# Atomic number(s) and corresponding mass fraction(s):
+13    1.000    # aluminum

gecatsim/material/Ar

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Ar (argon)
+
+# Number of elements:
+1
+
+# Density:
+0.001662
+
+# Atomic number(s) and corresponding mass fraction(s):
+18    1.000000    # argon

gecatsim/material/As

@@ -0,0 +1,11 @@
+# Elemental material: 
+# As (arsenic)
+
+# Number of elements:
+1
+
+# Density:
+5.73
+
+# Atomic number(s) and corresponding mass fraction(s):
+33    1.00    # arsenic

gecatsim/material/At

@@ -0,0 +1,11 @@
+# Elemental material: 
+# At (astatine)
+
+# Number of elements:
+1
+
+# Density:
+10.00
+
+# Atomic number(s) and corresponding mass fraction(s):
+85    1.00    # astatine

gecatsim/material/Au

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Au (gold)
+
+# Number of elements:
+1
+
+# Density:
+19.32
+
+# Atomic number(s) and corresponding mass fraction(s):
+79    1.00    # gold

gecatsim/material/B

@@ -0,0 +1,11 @@
+# Elemental material: 
+# B (boron)
+
+# Number of elements:
+1
+
+# Density:
+2.37
+
+# Atomic number(s) and corresponding mass fraction(s):
+5    1.00    # boron

gecatsim/material/Ba

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Ba (barium)
+
+# Number of elements:
+1
+
+# Density:
+3.50
+
+# Atomic number(s) and corresponding mass fraction(s):
+56    1.00    # barium

gecatsim/material/Be

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Be (beryllium)
+
+# Number of elements:
+1
+
+# Density:
+1.848
+
+# Atomic number(s) and corresponding mass fraction(s):
+4    1.000    # beryllium

gecatsim/material/Bi

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Bi (bismuth)
+
+# Number of elements:
+1
+
+# Density:
+9.747
+
+# Atomic number(s) and corresponding mass fraction(s):
+83    1.000    # bismuth

gecatsim/material/Br

@@ -0,0 +1,11 @@
+# Elemental material: 
+# Br (bromine)
+
+# Number of elements:
+1
+
+# Density:
+0.007072
+
+# Atomic number(s) and corresponding mass fraction(s):
+35    1.000000    # bromine

gecatsim/material/C

@@ -0,0 +1,11 @@
+# Elemental material: 
+# C (carbon)
+
+# Number of elements:
+1
+
+# Density:
+1.700000
+
+# Atomic number(s) and corresponding mass fraction(s):
+6    1.000000    # carbon

gecatsim/material/CIRS_adipose_adult2

@@ -0,0 +1,20 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: May 6, 2014
+# Material composition source: same as density
+
+# Material name:
+# CIRS_adipose
+
+# Number of elements:
+6
+
+# Density:
+0.962000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.100560    # hydrogen
+6    0.698419    # carbon
+7    0.016310    # nitrogen
+8    0.166400    # oxygen
+12    0.016510    # magnesium
+17    0.001801    # chlorine

gecatsim/material/CIRS_adipose_adult3

@@ -0,0 +1,20 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: September 30, 2015
+# Material composition source: same as density
+
+# Material name:
+# CIRS_adipose_adult3
+
+# Number of elements:
+6
+
+# Density:
+0.940000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.103072    # hydrogen
+6    0.726108    # carbon
+7    0.016912    # nitrogen
+8    0.151306    # oxygen
+17    0.001801    # chlorine
+20    0.000801    # calcium

gecatsim/material/CIRS_bone_10yo

@@ -0,0 +1,21 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: August 17, 2012
+# Material composition source: same as density
+
+# Material name:
+# CIRS_bone_10yo
+
+# Number of elements:
+7
+
+# Density:
+1.545000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.050715    # hydrogen
+6    0.401620    # carbon
+7    0.010603    # nitrogen
+8    0.340702    # oxygen
+12    0.041312    # magnesium
+17    0.000500    # chlorine
+20    0.154546    # calcium

gecatsim/material/CIRS_bone_1yo

@@ -0,0 +1,21 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: August 17, 2012
+# Material composition source: same as density
+
+# Material name:
+# CIRS_bone_1yo
+
+# Number of elements:
+7
+
+# Density:
+1.450000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.057717    # hydrogen
+6    0.450635    # carbon
+7    0.012304    # nitrogen
+8    0.316095    # oxygen
+12    0.034010    # magnesium
+17    0.000600    # chlorine
+20    0.128639    # calcium

gecatsim/material/CIRS_bone_5yo

@@ -0,0 +1,21 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: August 17, 2012
+# Material composition source: same as density
+
+# Material name:
+# CIRS_bone_5yo
+
+# Number of elements:
+7
+
+# Density:
+1.518000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.052321    # hydrogen
+6    0.416467    # carbon
+7    0.011104    # nitrogen
+8    0.333233    # oxygen
+12    0.035414    # magnesium
+17    0.000500    # chlorine
+20    0.150960    # calcium

gecatsim/material/CIRS_bone_adult_and_15yo

@@ -0,0 +1,21 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: August 17, 2012
+# Material composition source: same as density
+
+# Material name:
+# CIRS_bone_adult_and_15yo
+
+# Number of elements:
+7
+
+# Density:
+1.600000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.048314    # hydrogen
+6    0.370411    # carbon
+7    0.009703    # nitrogen
+8    0.356707    # oxygen
+12    0.061919    # magnesium
+17    0.000500    # chlorine
+20    0.152446    # calcium

gecatsim/material/CIRS_bone_newborn

@@ -0,0 +1,21 @@
+# Molecular weights from https://www.nist.gov/pml/atomic-weights-and-isotopic-compositions-relative-atomic-masses.
+# Density source: CIRS, as used for GE CRC's phantoms, Rev: August 17, 2012
+# Material composition source: same as density
+
+# Material name:
+# CIRS_bone_newborn
+
+# Number of elements:
+7
+
+# Density:
+1.407000
+
+# Atomic number(s) and corresponding mass fraction(s):
+1    0.064732    # hydrogen
+6    0.456528    # carbon
+7    0.011106    # nitrogen
+8    0.306653    # oxygen
+12    0.069535    # magnesium
+17    0.000500    # chlorine
+20    0.090945    # calcium