diff --git a/gstlal-inspiral/python/stats/inspiral_intrinsics.py b/gstlal-inspiral/python/stats/inspiral_intrinsics.py
index ab886f21f912f20b32cf04fec3a5e6e96018a4fa..bce7b8a4deb2de4456b21c49d562f57c726f600e 100644
--- a/gstlal-inspiral/python/stats/inspiral_intrinsics.py
+++ b/gstlal-inspiral/python/stats/inspiral_intrinsics.py
@@ -115,10 +115,11 @@ class SourcePopulationModel(object):
 				coefficients = model['coefficients'].value
 				snr_bp = model['SNR'].value
                                 try:
-                                        model_ids = model['event_id'].value
+                                        model_ids = model['template_id'].value
                                 except KeyError:
                                         # FIXME: assume sequential order if model['event_id'] doesn't exist
-                                        model_ids = numpy.arange(numpy.shape(model['coefficients'].value)[-1])
+                                        #model_ids = numpy.arange(numpy.shape(model['coefficients'].value)[-1])
+                                        model_ids = model['event_id'].value
 			# PPoly can construct an array of polynomials by just
 			# feeding it the coefficients array all in one go, but then
 			# it insists on evaluating all of them at once.  we don't