This PR comes from discussions with @Angela-White on slack about the most appropriate file to pass as --light_det_noise_filename to simulate_pixels.py for a spill level h5 file using a full DUNE ND hall geometry.
She advised that nothing existed but that module 1 information from could be extracted from this file and "multiplied" 35 times to make an NDLAr made up of 35 module 1s.
I am very much a non-expert here so including the logic I used to make the new "35 module 1s" file.
ln_2x2 = np.load("light_noise_2x2_4mod_July2023.npy")
n_chan = 96
module_number = 1
ln_module = ln_2x2[n_chan*module_number:n_chan*(module_number+1)]
n_modules_ndlar = 35
# Skip one count of module here, after 1 iteration of this
# loop we'll have a 2 module setup already.
for module in range(n_modules_ndlar-1):
if not module:
ln_ndlar = np.concatenate((ln_module, ln_module), axis=0)
continue
ln_ndlar = np.concatenate((ln_ndlar,ln_module), axis=0)
print(ln_2x2.shape)
print(ln_module.shape)
print(ln_ndlar.shape)
np.save('light_noise_ndlar_35xmodule1_August2023.npy', ln_ndlar)
This PR comes from discussions with @Angela-White on slack about the most appropriate file to pass as
--light_det_noise_filename
tosimulate_pixels.py
for a spill level h5 file using a full DUNE ND hall geometry.She advised that nothing existed but that module 1 information from could be extracted from this file and "multiplied" 35 times to make an NDLAr made up of 35 module 1s.
I am very much a non-expert here so including the logic I used to make the new "35 module 1s" file.