In summary, the statistical results collected between 5:00 and 12:00 AM on April 2, 2025. For each injection, we calculate the frequency mean value 
(f_mean) and the peak height (P) of 194Pb ions for frames with t > 0.14 s, correct each injection’s frequency by subtracting f_mean, and only sum 
the 
spectra from injections with P > 100; moreover, by comparing the 194Pb82+ spectrum with that of 184Pt78+—whose isomer at 1.8 MeV (T1/2 = 1 ms) is 
clearly visible—the setup is demonstrated to be suitable for detecting an isomer with a lifetime of about 1 ms.
===================================================================
How to dod analysis and plot this result:
1. Add up injections and create accumulated spectrum.
(pyroot) gpuuser@appc256:/data.local/G22-00018_00203/scripts/OnlineDataAnalysisSystem/shift/20250331/data$ sh generate_data.sh /data.local/G22-
00018_00203/results/analyzers/RSA01 RSA01- 
'*.tiq_spectrogram.npz' 2025.04.02.05.00.00.000 2025.04.02.12.00.00.000 0.001 408.59e6 408.62e6 0.1 1

After run this command, you will create the accumulated spectrum under the folder with time stamp of start time and stop time.

2. Download the code "plot_fft.cxx" and plot the spectrum with following command:
...2025.04.02.05.00.00.000_2025.04.02.12.00.00.000$ root
   ------------------------------------------------------------------
  | Welcome to ROOT 6.32.10                        https://root.cern |
  | (c) 1995-2024, The ROOT Team; conception: R. Brun, F. Rademakers |
  | Built for linuxx8664gcc on Feb 12 2025, 01:50:41                 |
  | From tags/6-32-10@6-32-10                                        |
  | With                                                             |
  | Try '.help'/'.?', '.demo', '.license', '.credits', '.quit'/'.q'  |
   ------------------------------------------------------------------

root [0] .L plot_fft.cxx
root [1] plot_fft("2025.4.2 5:00~12:00",411.14e6,411.18e6,0,100)

The ntcap data record 410 detector continuously. The frame length is  26 ms/frame * 32 frames= 832 ms.

I examined the spectrum of 194Pb at the 211 harmonic. It was observed that the peak corresponding to 194Pb82+ was mixed with signals 
from some unknown ions. Consequently, we switched to the 212 harmonic, where the frequency of 194Pb82 was measured at 410.5385024.

Currently, we only use 37% of the injections in our data analysis. However, 63% of the injections have a height below 100, and these may include 
single ion injections. We need to visually inspect these injections to identify the single ion injections.

Summary:

1.Simulated spectra for 194Pb and 188Hf.
2.The results are attached.
3. For the 188Hf simulation, some nuclides in the fragments are missing from the AME database, causing errors in the LiseReader program.
To fix this, please download the attached reader.py file and replace your local file at lisereader-main\lisereader\reader.py.
After replacing the file, recompile the program. Once recompiled, you should be able to simulate the 188Hf spectrum without issues.

The `queqiao` program used for RSA data analysis is located at:
  
  (pyroot) gpuuser@appc256:/data.local/G22-00018_00203/scripts/OnlineDataAnalysisSystem_rsa/queqiao

It is currently running in the screen session named:

  OnlineDataAnalysisSystem_rsa

1. Login the server gpuuser@appc256

2. Reattach to the session, use:

  screen -rd OnlineDataAnalysisSystem_rsa

3. Run the program manually, activate the pyroot environment and execute the following command:

  queqiao ../parameters/queqiaosettings_194Pb_rsa.toml

4. Then press "start" botton, the `queqiao` program copies data from:

  /mnt/ntcap-data2/G22-00018_00203/data/analyzers/RSA01
  /mnt/ntcap-data2/G22-00018_00203/data/analyzers/RSA02

to:

  /data.local/G22-00018_00203/data/analyzers/RSA01
  /data.local/G22-00018_00203/data/analyzers/RSA02

This program is used to combine spectrum of many injections into one spectrum.

Example of usage:

(pyroot) gpuuser@appc256:/data.local/G22-00018_00203/scripts/OnlineDataAnalysisSystem/shift/20250331/data$ bash generate_data.sh /data.local/G22-00018_00203/results/analyzers/RSA01 RSA01- 
'*.tiq_spectrogram.npz' 2025.04.01.04.00.00.000 2025.04.01.19.00.00.000 0.001 408590000 408620000 0.3

Until now, there is no 184Hf72+ and 186Hf72+ ions. 

still can not see both 186Hf72+ and 184Hf71+

This shows the spectrum of 188Hf72+ during "April06_01:00_02:00". Until April 6 2:00, the Bp = 7.918262 Tm has not been changed. 
The ntcap spectrum is from:
/data.local/G22-00018_00203/results/ntcap/iq/IQ_2025-04-05_17-43-49/0000756.iq.tdms.root


analysis settings:

2^18 samples per frame + average over 32 frame.
26.21 ms per frame
838.86 ms for 32 averaged frames

This is a preliminary result during 4:00-19:00 April 1st 2025. 
I correct the frequency drift of 194Pb82+ ions injection by injection.

I have analyzed all injection data during 4:00-19:00 and calculated the frequency of 194Pb82+ for each injection. This is the frequency of 194Pb82+ as 
a function of time. You can clearly see that the frequency of 194Pb82+ was drifted from 408.603 to 408.600 MHz. The variation of frequency is about 
3kHz.

datafile = "/data.local/G22-00018_00203/results/analyzers/RSA01/RSA01-2025.04.01.14.13.07.068.tiq_spectrum.npz"
filep = "/data.local/G22-00018_00203/scripts/OnlineDataAnalysisSystem/rionid_exercises/exercise_194Pb_Exp_2025/New_TEline2022-ESR-194Pb-208Pb-nostripper_524MeV.lpp"
alphap = "1.399"
harmonics = "211"
refion = "194Pb+82"
circumference = "108.34"
mode = "Frequency"
value = "408602000"
reload_data = true
nions = "100"

Now ESR is set to the 188Hf72+ setting.