fix tests

src/chromo/common.py

@@ -594,15 +594,17 @@ def __call__(self, nevents):
                     event = self._event_class(self)
                     # boost into frame requested by user
                     self.kinematics.apply_boost(event, self._frame)
-                    if self._apply_decay_handler and (
-                        not np.all(
+                    if self._apply_decay_handler:
+                        status1_pids = event.pid[event.status == 1]
+                        may_decay = np.isin(status1_pids, all_decaying_pids)
+
+                        if not np.all(
                             np.isin(
-                                event.pid[event.status == 1],
+                                status1_pids[may_decay],
                                 self._final_state_particles,
                             )
-                        )
-                    ):
-                        self._decay_handler(event)
+                        ):
+                            self._decay_handler(event)
                     yield event
                     continue
                 nretries += 1

src/chromo/models/dpmjetIII.py

@@ -111,6 +111,7 @@ def __init__(self, evt_kin, *, seed=None):
         # Set PYTHIA decay flags to follow all changes to MDCY
         self._lib.pydat1.mstj[21 - 1] = 1
         self._lib.pydat1.mstj[22 - 1] = 2
+        self._lib.dtfrpa.ovwtdc = False
         self._set_final_state_particles()
 
     def _cross_section(self, kin=None):

tests/test_setstable.py

@@ -10,6 +10,7 @@
 
 decay_list = [
     lp.pi_plus.pdgid,
+    lp.pi_minus.pdgid,
     lp.K_plus.pdgid,
     lp.pi_0.pdgid,
     lp.K_S_0.pdgid,