Hey @scott.coughlin There are three notebooks that I was able to track down but not yet the function that actually generates the histograms…I think that’s on LLO
If you want to get an idea of the method for getting the VCO values, then at CIT there’s this one:
The only real difference is that usually upsample the VCO using a fit with IMC-F to get a “fast” channel. maybe that’s mentioned in there.
There are two methods for estimating whistle rates that I’ve come up with…both based on histograms. The code for the histograms is fairly simple…it just creates the fast VCO and finds its value at each omicron trigger and histograms that value…it then takes those bins and also fills them with the timeseries for the “fast” VCO over the same time span.
One method tries to fit the latter histogram to the former suing a simple linear fit…iteratively removing outlier bins (which should be the bins associated with whistles) until you stop improving the fit. Despite the name…that notebook is here:
The other method, uses emcee to estimate a whistle rate in each histogram bin and simultaneously fit the background essentially the same idea as the previous one…just with a 51-parameter emcee run. it takes longer and doesn’t work any better…so I probably wont’ use it. anyway…that notebook is here:
/home/meyers/whistles_nestle_sampler/whistles/tests/whistle_test_emcee_sampler-Copy1.ipynb
Like I said, I’ll try to get the simple function I use to get the VCO value off of LLO and send that to you as well
Hey @scott.coughlin There are three notebooks that I was able to track down but not yet the function that actually generates the histograms…I think that’s on LLO
If you want to get an idea of the method for getting the VCO values, then at CIT there’s this one:
/home/meyers/forPeople/me/Thesis/detchar/whistle_examples.ipynb
The only real difference is that usually upsample the VCO using a fit with IMC-F to get a “fast” channel. maybe that’s mentioned in there.
There are two methods for estimating whistle rates that I’ve come up with…both based on histograms. The code for the histograms is fairly simple…it just creates the fast VCO and finds its value at each omicron trigger and histograms that value…it then takes those bins and also fills them with the timeseries for the “fast” VCO over the same time span.
One method tries to fit the latter histogram to the former suing a simple linear fit…iteratively removing outlier bins (which should be the bins associated with whistles) until you stop improving the fit. Despite the name…that notebook is here:
/home/meyers/whistles_nestle_sampler/whistles/tests/whistle_test_emcee_sampler.ipynb
The other method, uses
emcee
to estimate a whistle rate in each histogram bin and simultaneously fit the background essentially the same idea as the previous one…just with a 51-parameteremcee
run. it takes longer and doesn’t work any better…so I probably wont’ use it. anyway…that notebook is here:/home/meyers/whistles_nestle_sampler/whistles/tests/whistle_test_emcee_sampler-Copy1.ipynb
Like I said, I’ll try to get the simple function I use to get the VCO value off of LLO and send that to you as well