"""Define bins for dark current distribution calculations"""
import numpy as np

x0 = 2.
x1 = 20000.
dx = 0.05
log_bins = np.arange(np.log10(x0), np.log10(x1), dx)
bins = 10 ** log_bins
bin_centers = 10 ** ((log_bins[1:] + log_bins[:-1]) / 2.0)

@@ -78,15 +78,90 @@
 np%0A%0A
+# Even log spacing from exactly 2 to 2000, then additional bins up to 11000%0A
 x0 = 2
-.
 %0Ax1 
@@ -170,42 +170,56 @@
 2000
-0.%0Adx = 0.05%0Alog_bins = np.arange(
+%0Ax2 = 11000%0Adx = 0.05%0A%0Alx0 = np.log10(x0)%0Alx1 = 
 np.l
@@ -228,11 +228,16 @@
 10(x
-0),
+1)%0Alx2 =
  np.
@@ -243,18 +243,112 @@
 .log10(x
-1)
+2)%0A%0An_bins = np.round((lx1 - lx0) / dx)%0Adx = (lx1 - lx0) / n_bins%0A%0Alog_bins = np.arange(lx0, lx2
 , dx)%0Abi

Update dark bins