Skip to content
This repository has been archived by the owner on Jun 5, 2024. It is now read-only.

S2 Photon timing downsampled for large time in instruction #425

Open
FaroutYLq opened this issue Dec 14, 2023 · 1 comment
Open

S2 Photon timing downsampled for large time in instruction #425

FaroutYLq opened this issue Dec 14, 2023 · 1 comment
Labels
bug Something isn't working

Comments

@FaroutYLq
Copy link
Contributor

When taking a instruction with huge timestamp (also in int64), you will have weird S2 waveforms, where the photons timing are downsampled. More discussion here. It is true to all S2s simulated.

Example for simulation using big time instruction

import cutax
st = cutax.xenonnt_sim_base(output_folder='/scratch/midway3/yuanlq/wfsim')
st.set_config({"fax_file" : "/project/lgrandi/yuanlq/wimp/s2sim/021952-9-se-all-20.csv"}) # feel free to check
se_peaklets = st.get_array('021952_se', 'peaklets')

for p in se_peaklets[:100]:
    if p['dt']==10:
        plt.plot(p['data'])
plt.xlabel('Sample [10ns]')

image

Example for simulation using small time instruction (what people usually do)

st = cutax.xenonnt_sim_base(output_folder='/scratch/midway3/yuanlq/wfsim')
st.set_config({"fax_file" : "/project/lgrandi/yuanlq/wimp/s2sim/jy_se_f6191.csv"}) # the first 100 rows of /project/lgrandi/jyangqi/se_sim_stuff/se_f6191.csv

se_peaklets = st.get_array('se_f6191', 'peaklets')
for p in se_peaklets[:100]:
    if p['dt']==10:
        plt.plot(p['data'])
plt.xlabel('Sample [10ns]')

image

The issue is suspected to happen when int64 were casted into float.

@ramirezdiego
Copy link
Collaborator

Solution in XENONnT/fuse#171.

Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

2 participants