NA62 beam studies (1/11/22) -------------------------- -------------------------------------------------------------------------------------------------------------- SPILL - poisson plots Use kspill6 and kspillff Fig 6a general test plot - kspill4.kumac reads 3 files Uses heapsort. code bsed on kstamp - integer based.. kumac only Fig 6b kspill5.kumac reads 3 files as kspill4 + additions FINAL 060123 Fig 6c kspill5dev.kumac reads 3 files as kspill4 + additions DEV ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** Fig 6d kspill6.kumac new plots P1 for mean/variance + sample variance (20/2/23) 3/3/23 9 files gaussian for 1 ms limited range fit. prob gaussian. Best version 3/3/23. ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** Fig 6dran kspill6.kumac new plots P1 for mean/variance + random sample variance (20/2/23) Fig 6e kspill7.kumac as kspill6.kumac - minimal plots ************************************************************************************************************** Fig 6f kspill8.kumac as kspill7.kumac - minimal plots +prob of gaussian 4 files. Now outdated (3/3/23) use kspill6.kumac : 9 files + 1 ms gaussia + sample variance tests. ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** Fig 6g kspillff.kumac to kspillf.ps (pdf) for filter tests. 3/3/23 version: filter tests for 100 ms and 1 sec. Demonstrates effect of filtering and shows the low frequency component. ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** ************************************************************************************************************** Fig 6h kspillp.kumac , test periodogram 1 sec (kspillff derivative) Fig 6h kspillp2.kumac , test periodogram + correlogram (kspillff derivative) Fig 6h kspillp4.kumac , test periodogram * 4 (4 secs) (kspillff derivative) Displays kspillp4m.kumac outputs kspillp4.pdf , uses modified periodogram that dispays amplitude( equiv autocprr10m.kumac) -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- AUTOCORR versions Fig 11f autocorr7r.kumac - random time for test in autocorr for spill and in input for FFT. 110123-v2 has sin(wt) for spill and FFT. Gap for folded time. Fig 11g autocorr8.kumac - DEV - has Hurst with log-log fit. Fig 11G autocorr8a.kumac - DEV - has Hurst with log-log fit. Possible bug removed in autocorr8a.kumac - arrays zeroed. Fig 11f autocorr9.kumac - low frequency FFT ******************************************************************************************* ******************************************************************************************* Fig 11g autocorr10.kumac - FFT minimal version Updated to autocorr10m.kumac FFT with trend subtracted. (o/p to autocorr10.pdf) 04/03/23 ******************************************************************************************* ******************************************************************************************* *********************************************************************************************** Fig 11h autocorr11.kumac - FFT for 1 sec spill regions. i********************************************************************************************** -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- FFT tests Fig A fftplot1.kumac - FFT stand-alone Fig B fftplot2.kumac - FFT stand-alone - low frequency. Fig C fftplot3.kumac - FFT stand-alone - 3 files - low frequency Fig D fftplot32.kumac - FFT stand-alone - 3 files - low + high frequency Fig E fftplot1dev.f to cft.kumac - FFT stand-alone - cft version Shows three plots for three intervals in spill. 100 Hz signal absent in second selected region. ------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- Programs 05/02/23 CFT D702 tests Fig Ea fftplot1dev.f to cft1.kumac - FT stand-alone - cft version -file 1 Fig Eb fftplot2dev.f to cft2.kumac - FT stand-alone - cft version -file 2 Fig Ec fftplot3dev.f to cft3.kumac - FT stand-alone - cft version -file 3 Fig E1 fftplot41.f to cft41.kumac - FT stand-alone - cft version -file 1 Fig E2 fftplot42.f to cft42.kumac - FT stand-alone - cft version -file 2 Fig E3 fftplot43.f to cft43.kumac - FT stand-alone - cft version -file 3 Fig Ed fft1test.f to cft1test.kumac - FFT stand-alone - cft version -file 1 test version Fig Ed fft2test.f to cft2test.kumac - FFT stand-alone - cft version -file 2 test version Fig Ed fft1test.f to cft3test.kumac - FFT stand-alone - cft version -file 3 test version ------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- Simulation FFT + R/S Fig 12 simulation, corrtest2.kumac - check autocorr , FFT ************************************************************************************************************** ************************************************************************************************************** Fig 122 simulation, corrtest22.kumac - check autocorr , FFT, standard periodogram. Set to show output from corrtest22m.kumac. corrtest22m.kumac FFT + modified periodogram for f = 0 and amplitude. corrtest22m.kumac to corrtest22.ps mod to give same results as FFT (09/03/2023 ). Use this version for tests. Demonstrates: Peak at f = 0, = number of input events Peaks at f .ne. 0, = (0.5 * amplitude) of sin/cos terms. Set to show 100Hz ( k = 5 ) Fig 12dd simulation of data, corrtestdd.kumac. structure as corrtest22m.kumac Fig 12nn simulation of noise, corrtestnn.kumac. ************************************************************************************************************** ************************************************************************************************************** Fig 12a simulation, corrtest3.kumac - check autocorr , FFT , Hurst R/S Fig 12b simulation, corrtest4.kumac - check autocorr , FFT , Hurst + fit to log-log ***************************************************************************************************************************** Fig 12bran simulation, corrtest41.kumac - random only check FFT , Hurst + fit to log-log Comment: corrtest41.kumac is the only Hurst code that has been checked to agree with the literature. It does not give the expect slope for log R/S vs log time. ***************************************************************************************************************************** Fig 12c simulation, corrtest5.kumac - as corrtest4 , debug version Fig 12d simulation, corrtest9.kumac - random only - for study/debug Showing corrtest10.kumac for 5, 10, 15 Hz. Rerun corrtest9 for standard plots. -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- Fourier Transform - Python code Fig 6 Fourier transform of spill (3 - 4 sec) 12465/170 Fig 7 Fourier transform of spill (3- 5 sec) 12465/170 Fig 8 simulate 50, 100, 160 Hz spectrum (fftcode1.py ) 30/11/22 Fig 9 Fourier transform of simulation fig 8 (fftcode1.py) 30/11/22 -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- plots in ~/public_html/newplots