small GUI tweak

bc_qualityMetrics_pipeline.m

@@ -63,13 +63,19 @@
 %% example: get the quality metrics for one unit
 % this is an example to get the quality metric for the unit with the
 % original kilosort and phy label of xx (0-indexed), which corresponds to
-% the unit with qMetric.clusterID == xx + 1. This is *NOT NECESSARILY* the
+% the unit with qMetric.clusterID == xx + 1, and to
+% qMetric.phy_clusterID == xx . This is *NOT NECESSARILY* the
 % (xx + 1)th row of the structure qMetric - some of the  clusters that kilosort
 % outputs are empty, because they were dropped in the last stages of the
 % algorithm. These empty clusters are not included in the qMetric structure
+% there are two ways to do this: 
+% 1:
 original_id_we_want_to_load = 0;
-id_we_want_to_load_1_indexed = original_id_we_want_to_load + 1;
+id_we_want_to_load_1_indexed = original_id_we_want_to_load + 1; 
 number_of_spikes_for_this_cluster = qMetric.nSpikes(qMetric.clusterID == id_we_want_to_load_1_indexed);
+% or 2:
+original_id_we_want_to_load = 0;
+number_of_spikes_for_this_cluster = qMetric.nSpikes(qMetric.phy_clusterID == original_id_we_want_to_load);
 
 
 

qualityMetrics/bc_runAllQualityMetrics.m

@@ -136,7 +136,9 @@
 
     % get this unit's attributes 
     thisUnit = uniqueTemplates(iUnit);
-    qMetric.clusterID(iUnit) = thisUnit;
+    qMetric.phy_clusterID(iUnit) = thisUnit - 1; % this is the cluster ID as it appears in phy
+    qMetric.clusterID(iUnit) = thisUnit; % this is the cluster ID as it appears in phy, 1-indexed (adding 1)
+
     theseSpikeTimes = spikeTimes_seconds(spikeTemplates == thisUnit);
     theseAmplis = templateAmplitudes(spikeTemplates == thisUnit);
 

visualizationTools/bc_unitQualityGUI.m

@@ -442,8 +442,7 @@ function updateUnit(unitQualityGuiHandle, memMapData, ephysData, rawWaveforms, i
 else
     set(guiData.isiTitle, 'String', '\color[rgb]{0 .5 0}ISI');
 end
-set(guiData.isiLegend, 'String', [num2str(qMetric.fractionRPVs_estimatedTauR(iCluster)*100), ' % r.p.v.', newline, ...
-    'estimated refractory period (ms): ', num2str(qMetric.RPV_tauR_estimate(iCluster))])
+set(guiData.isiLegend, 'String', [num2str(qMetric.fractionRPVs_estimatedTauR(iCluster)*100), ' % r.p.v.'])
 
 %% 6. plot isolation distance
 if param.computeDistanceMetrics