Add status to r packet last interaction tracker

tardis/transport/montecarlo/packet_trackers.py

@@ -4,6 +4,7 @@
 import numpy as np
 import pandas as pd
 
+from tardis.transport.montecarlo.r_packet import PacketStatus
 
 boundary_interaction_dtype = np.dtype(
     [
@@ -211,6 +212,7 @@
 
 rpacket_last_interaction_tracker_spec = [
     ("index", int64),
+    ("status", int64),
     ("r", float64),
     ("nu", float64),
     ("energy", float64),
@@ -227,6 +229,10 @@
     ----------
         index : int
             Index position of each RPacket
+        status : ENUM
+            status of the RPacket whether it has been reabsorbed or emitted.
+            1: EMITTED
+            2: REABSORBED
         r : float
             Radius of the shell where the RPacket is present
         nu : float
@@ -244,6 +250,7 @@
         Initialize properties with default values
         """
         self.index = -1
+        self.status = PacketStatus.IN_PROCESS
         self.r = -1.0
         self.nu = 0.0
         self.energy = 0.0
@@ -255,6 +262,7 @@
         Track properties of RPacket and override the previous values
         """
         self.index = r_packet.index
+        self.status = r_packet.status
         self.r = r_packet.r
         self.nu = r_packet.nu
         self.energy = r_packet.energy