I've created a spectrogram (Fig 1) ipython notebook in ~/Git/cit_ctnlab/ctn_labdata/scripts/SeismicAndScatteringStudy.ipynb. I also attached a median beatnote ASD (Fig 2) for reference.
This is the beginnings of a study to look for coherence between our accelerometers and beatnote to figure out the velocity of our scatterers.
It takes a little while for nds2 to acquire the data from Gabriele's cymac3, because the data is sampled at 65536 Hz. For 300 seconds of data it takes the script 172 seconds to retrieve the data.
I'm currently thinking about ways to make this spectrogram plot on the order of days as opposed to minutes so we can get a long-term idea about our beatnote. Even looking at this, it seems like the scattering shelf does oscillate a bit at around 10 Hz. Our 500 Hz hump seems pretty constant on this time scale as well. |