Gamma sampling documentation

docs/physics/tardisgamma/index.rst

@@ -5,3 +5,4 @@ TARDIS gamma :math:`\gamma`
 
 .. toctree::
     packetinitialization
+    packetsampling

docs/physics/tardisgamma/packetsampling.rst

@@ -0,0 +1,13 @@
+Packet Sampling
+===============
+
+
+In the simulation, we use spherical symmetry where the ejecta is divided into multiple spherical shells. To find the mass of each shell we multiply the mass fractions of each of the elements and isotopes with the density. We can set the timespace for the simulation
+to be between t\ :sub:`start`\  and t\ :sub:`end`\  using either a linear or a logarithmic scale.
+
+
+We use the `radioactivedecay <https://radioactivedecay.github.io>` python package to calculate the number of decays in each channel for each shell and at every timestep for each isotope.
+The composition is also updated after each timestep following :cite:`2024arXiv240308769G`.
+Then to calculate the decay energy of each channel, we multiply the total number of decays with the energy of each channel.
+We give weight to the decay energy of each channel of each isotope to sample the packets which takes into account the spatial distribution
+of radiocative isotopes present in different shells and the evolution of the composition over time.

docs/tardis.bib

@@ -348,3 +348,20 @@ @ARTICLE{Boyle2017
        adsurl = {https://ui.adsabs.harvard.edu/abs/2017A&A...599A..46B},
       adsnote = {Provided by the SAO/NASA Astrophysics Data System}
 }
+
+@ARTICLE{2024arXiv240308769G,
+       author = {{Guttman}, Or and {Shenhar}, Ben and {Sarkar}, Arnab and {Waxman}, Eli},
+        title = "{The thermalization of $\gamma$-rays in radioactive expanding ejecta: A simple model and its application for Kilonovae and Ia SNe}",
+      journal = {arXiv e-prints},
+     keywords = {Astrophysics - High Energy Astrophysical Phenomena},
+         year = 2024,
+        month = mar,
+          eid = {arXiv:2403.08769},
+        pages = {arXiv:2403.08769},
+          doi = {10.48550/arXiv.2403.08769},
+archivePrefix = {arXiv},
+       eprint = {2403.08769},
+ primaryClass = {astro-ph.HE},
+       adsurl = {https://ui.adsabs.harvard.edu/abs/2024arXiv240308769G},
+      adsnote = {Provided by the SAO/NASA Astrophysics Data System}
+}