NA62 beam studies (1/11/22) -------------------------- -------------------------------------------------------------------------------------------------------------- SPILL - poisson plots 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 Fig 6e kspill7.kumac as kspill6.kumac - minimal plots Fig 6f kspill8.kumac as kspill7.kumac - minimal plots +prob of gaussian -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- AUTOCORR versions Fig 11 general test plot -autocorr2.kumac Fine foldtime bins (0.00005) Fig 11a autocorr3.kumac as autocode2 updated - 06/01/23 Fig 11b autocorr4.kumac - with statistics - FINAL 08/01/23 Fig 11c autocorr5.kumac - add autocorr for spill - FINAL 09/01/23 Fig 11d autocorr6.kumac - no Hurst - FINAL 09/01/23 -STANDARD Fig 11e autocorr7.kumac - no Hurst - FINAL 15/01/23 - STANDARD DATA 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 - minimal version -------------------------------------------------------------------------------------------------------------- -------------------------------------------------------------------------------------------------------------- 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 Fig Ea fftplot1dev.f to cft1.kumac - FFT stand-alone - cft version -file 1 Fig Eb fftplot2dev.f to cft2.kumac - FFT stand-alone - cft version -file 2 Fig Ec fftplot3dev.f to cft3.kumac - FFT 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 12a simulation, corrtest3.kumac - check autocorr , FFT , Hurst R/S Fig 12b simulation, corrtest4.kumac - check autocorr , FFT , Hurst + fit to log-log 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