adding more debugging prints

Author: jessicapilling

xga/sourcetools/_common.py

@@ -258,7 +258,9 @@ def _setup_inv_abel_dens_onion_temp(sources: Union[GalaxyCluster, ClusterSample]
     sources, outer_rads, has_glob_temp = _setup_global(sources, outer_radius, global_radius, abund_table, group_spec,
                                                        spec_min_counts, spec_min_sn, over_sample, num_cores, psf_bins,
                                                        stacked_spectra, src_telescopes)
-    
+    print("has_global_temp")
+    print(has_glob_temp)
+
     rads_dict = {str(sources[r_ind]): r for r_ind, r in enumerate(outer_rads)}
 
     # This checks and sets up a predictable structure for the models needed for this measurement.

xga/sourcetools/temperature.py

@@ -814,6 +814,8 @@ def onion_deproj_temp_prof(sources: Union[GalaxyCluster, ClusterSample], outer_r
                                           psf_corr, psf_model, psf_bins, psf_algo, psf_iter, group_spec, min_counts,
                                           min_sn, over_sample, one_rmf, freeze_met, abund_table, temp_lo_en, temp_hi_en,
                                           num_cores, stacked_spectra=stacked_spectra, telescope=telescope)
+        print("ann_rads")
+        print(ann_rads)
     elif annulus_method == "growth":
         raise NotImplementedError("This method isn't implemented yet")
 
@@ -837,9 +839,12 @@ def onion_deproj_temp_prof(sources: Union[GalaxyCluster, ClusterSample], outer_r
                 # The normalisation profile(s) from the fit that produced the projected temperature profile.
                 apec_norm_prof = src.get_apec_norm_profiles(cur_rads, group_spec, min_counts,
                                                             min_sn, over_sample, telescope=tel)
+                print("profiles here")
+
             except NoProductAvailableError:
                 # warn("{s} doesn't have a matching projected temperature profile, skipping.")
                 all_3d_temp_profs[tel][src_ind] = None
+                print("no profiles!")
                 continue
 
             # We need to check if a matching 3D temperature profile has already been generated, as then we
@@ -855,10 +860,14 @@ def onion_deproj_temp_prof(sources: Union[GalaxyCluster, ClusterSample], outer_r
 
             obs_id = 'combined'
             inst = 'combined'
+            print(tel)
+            print("len(proj_temp)")
+            print(len(proj_temp))
             # There are reasons that a projected temperature profile can be considered unusable, so we must check. Also
             #  make sure to only use those profiles that have a minimum of 4 annuli. The len operator retrieves the number
             #  of radial data points a profile has
             if proj_temp.usable and len(proj_temp) > 3:
+                print("In here")
                 # Also make an Emission Measure profile, used for weighting the contributions from different
                 #  shells to annuli
                 em_prof = apec_norm_prof.emission_measure_profile(src.redshift, src.cosmo, abund_table,