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