SNOwGLoBES / snowglobes

SNOwGLoBES code: SuperNova Observatories with GLoBES
GNU General Public License v2.0
19 stars 19 forks source link

channels_argon_400bins.dat #24

Open mcolomermolla opened 2 months ago

mcolomermolla commented 2 months ago

This channel does not run:

running rc = RateCalculator(base_dir="/home/marta/snowglobes_git/") rates = rc.run(flux, detector='ar40kt', material='argon_400bins') one gets as error: AssertionError: Matrix shape (200, 200) is inconsistent with (bins_true-1, bins_smeared-1)=(400, 400)

mcolomermolla commented 2 months ago

We should figure out how to mix these different binning for different channels (ar_coh VS non coh channels)

JostMigenda commented 2 months ago

Hmm … is there an issue in the channels_argon_400bins.dat file? I.e. are you saying it shouldn't mix channels with 200 and 400 bins?

Or is this an issue with SNEWPY, which should be updated to support such mixed channels? (In that case, it would be better to open an issue on the SNEWPY repo.)

mcolomermolla commented 2 months ago

Now the only dependency on snowglobes from snewpy is those channel/eff/smear files... so yes, it is rather from snewpy, it can't handle the mix of different bins in different channels. But I am not sure the error may come already from before (matrix incompatibility between xsec or eff and smear binning), so more "snowglobes".

JostMigenda commented 1 month ago

I’ve opened an issue on the SNEWPY repo: https://github.com/SNEWS2/snewpy/issues/334