Routine for analyzing BPM turn by turn data

&bpm_data
 numbers = 19887 19891 19895

 tune_x = 0.56
 tune_y = 0.63
/         

• When the turn by turn data are acquired for all ~100 BPMs, the data are written to two consecutive files, 50 BPMs in each.
• It is only necessary to list the first of the two files, as the cbpm_anal program will automatically look for the consecutive file.
• The above input will instruct the program to analyze the data in BPM files (RD-19887 RD-19888), (RD-19891 RD-19892), (RD-19895 RD-19896).
• And to look in the fft in the vicinity of 0.56 and 0.63 for the horizontal and vertical tunes.

The program will create subdirectories RD_19887, RD_19891, and RD_19895. In each subdirectory is written a file for each BPM and each bunch.
For example:

• RD_19887_bunch_001_BPM_015 is the data for bunch 1 at BPM 15.
• Each file lists x and y position on each turn and x and y fft.
• Each also identifies average positions and peaks in the spectra, as well as other information from the RD file.
• There is a summary file (RD_19887_bunch_001.dat) that lists average and rms position at each BPM.

Code