Monte Carlo from Scratch
From the lab to the computer¶

Monte Carlo simulation of radiation transport is a scientific instrument by its own right. Cultivating the skill and the art is no less than gaining mastery in instrumentation. Monte Carlo is never a plug-and-play software tool. Sadly, Monte Carlo is often abused as a black box. Throughout the years I have also been ‘offered’ co-authorship of to-be papers which were to become publishable if I ran Monte Carlo for them. No, Monte Carlo cannot make an otherwise unpublishable work publishable.
A piece of study begins with a study or experimental design, considering the choice of samples and the control of variables. Letting all variables lose is not going to lead us to anything conclusive.
Embarking on a Monte Carlo project, I would think of the experimental setup the similar way as I would do when setting up a physical experiment in the radiation laboratory. The difference is that now:
we need not worry about ALARA or ALARP, radiation exposure or radiation protection; we can do cleaner experiments, we shall be able to control controlled variables much better; we can do experiments which would otherwise be
unethical e.g. shooting test beams into human bodies;
impractical e.g. knock down shielding walls and rebuild new ones iteratively as part of the design process;
impossible e.g.
detecting particles directly and perfectly without perturbing the radiation field and without - requiring the particles to interact and deposit energy;
finding out about the particle’s ancestry at the time we detect the particle; we call this technique ‘latching’.
Experiments¶

Irradiation condition¶

Emission and transmission¶

Shielding¶

A source can be very simple¶

A source can be very complex¶

Photon cross-sections¶

A simple estimation¶

FLUKA input decks¶

Geometry and materials¶

The input file¶

A simple simulation¶

100% escaped¶

Carbon instead of vacuum¶

80% transmission¶

More histories¶

70% transmission¶

Lower cutoffs¶

Carbon beam instead of carbon shell¶

Beam-start position¶

Pion+¶

Alpha¶

Helium-3¶

Deuterium¶

Tritium¶

Gamma¶

Electron¶

Positron¶

Proton¶

Proton¶

Proton¶

Counts of different particles¶

So many neutron¶

Pion-¶

Carbon-12¶

Nine protons¶

The 9 protons¶

2nd-generation protons¶

Two gammas¶

Pion+¶

The pion+¶

The other pion+¶

Particles no longer tracked¶

A history out of many¶

FLUKA¶

Fortran coding is optional¶

Housebrands¶
