diff --git a/pixel_mask_analysis.py b/pixel_mask_analysis.py
index 3201432..5669ba7 100755
--- a/pixel_mask_analysis.py
+++ b/pixel_mask_analysis.py
@@ -67,6 +67,12 @@ clargs = parser.parse_args()
 with h5py.File(clargs.fname) as f:
     image = f["pixel_mask"][()]
 
+
+data = image[image != 0].flatten()
+plt.hist(data, bins=np.arange(min(data), max(data)+2)-0.5)
+plt.show()
+
+
 values = get_bit_values(image)
 
 n_plots = len(values) + 2