get_cube_shielding.m is a script that:
- generates 3D polygons in the shape of cubes of different sizes,
- generates coordinates for samples at the center of each side of the cube,
- uses Balco's code (
generate_cosmic_rays.m) to calculate the shielding correction for each sample - and calculates the approximations that describe the shielding correction vs. the size of the cube.
All the scripts that simulate the cosmic radiation in the samples are extracted from the supplementary material of the following publication:
Greg Balco (2014) Simple computer code for estimating cosmic-ray shielding by oddly shaped objects. Quaternary Geochronology, Volume 22, Pages 175-182.
Side shielding correction:
y=0.5+0.5*exp(-x*0.8)
Bottom shielding correction:
y=exp(-x*1.1)
where x is the side of the dice in metres
After repeating this experiment for different densities (d) and particle attenuation lengths (L):
Side shielding correction:
y=0.5+0.5*exp(-z*d/(2*L*0.8))
Bottom shielding correction:
y=exp(-z*d/(L*1.16))
where z is the side of the dice in cm
Ángel Rodés, 2022