This isn't working for me. I replaced the following:

```
ecen = get_bin_centers(self.ebins)                                       
czcen = get_bin_centers(self.czbins)                                     
for ie,egy in enumerate(ecen):                                           
    ebin_width = (self.ebins[ie+1] - self.ebins[ie])                     
    for icz,cz in enumerate(czcen):                                      
        czbin_width = (self.czbins[icz+1] - self.czbins[icz])            
        aeff_hist[ie][icz] /= (ebin_width*czbin_width*2.0*np.pi)   
 ```
with what I think you'd like me to do:
```
ebin_sizes = get_bin_sizes(ebins)
czbin_sizes = 2.0*np.pi*get_bin_sizes(czbins)
bin_sizes = np.meshgrid(ebin_sizes,czbin_sizes)
aeff_hist /= np.abs(bin_sizes[0]*bin_sizes[1])
```

and I got this error: _ValueError: operands could not be broadcast together with shapes (38,40) (40,38) (38,40)_
Do you know why this is?