Analysis codes

---

calibrations for bin_dgs in GEBSort_nogeb

GEBSort_nogeb is the program that can analyze data from DGS, DFMA and GRETINA.

You can get if from here:

cd workdir (e.g., /home/gtuser/gebsort)
  svn co https://svn.anl.gov/repos/gs_analysis/GEBSort .
or
  wget http://www.phy.anl.gov/gretina/GEBSort/AAAtar.tgz
  tar -zxvf AAAtar.tgz

For DGS data, enable bin_dgs in the GEBSort.chat file. To find the PZ values to use, sort some data from a 207Bi source. Then extract the pz spectra in .spe format with the get_pz.cc script

  GEBSort_nogeb ....
  rootn.exe
  .x get_pz.cc

Now run:

  dgs_pz 350 100 dgs_pz.cal 0.993

where 350 100 are the M and K values you find in the runxx.save file. Specify the values in 10 nsec units. In this case I saw these lines in the .save file:

 caput GLBL:DIG:d_window 0.06   
 caput GLBL:DIG:k_window 0.20     
 caput GLBL:DIG:m_window 3.50
 caput GLBL:DIG:k0_window 0.80
 caput GLBL:DIG:d3_window 0.20
 caput GLBL:DIG:raw_data_window 0.32
 caput VME01:SDIG1:k0_window0 0.0
 caput VME01:SDIG1:k0_window1 0.0
 etc

for the K value: sum up all the K and D values, in this case: 0.06+0.20+0.80+0.20 = 1.26 us or 126 in 10 nsec units. Notice that what is considered the K value also includes the D values (per SZ 6/25/18) as well as a D2 which is fixed at 0.15 (per JTA 6/26/18) and not in the listing above because the user cannot set it. Thus, in total, K in this example is 1.41 us or 141 in 10 ns units. The M value is 3.50 us or 350 in 10 nsec units. The 0.993 is a modification factor that needs to be determined by looking at energy vs baseline spectra.

  **Notice that M is also in the bin_dgs.c code** , so if you change the M value, you must first change it in the c code before following the instructions on this page.

Now, remove the energy calibration file:

 rm dgs_ehi.cal

so that the calibrations defaults to 0 and 1 for offset and gain and sort again using the new pz values that were extracted above. When you resort, the PZ values in dgs_pz.cal are read in and used. Extract the new clean, uncalibrated, ehi spectra as

  .x get_ecln.cc

and run the calibration program

  dgs_ecal dgs_ehi.cal 207Bi

you can also use "88Y", "60Co" for the source. The calibration will be 1keV/ch. Next when you run GEBSort_nogeb, both the new PZ values in dgs_pz.cal and the gain and offset values in dgs_ehi.cal are read in and used.

The energy processing in bin_dgs follows algorithms that were developed by Shoufei Zhu.

TODO:

  put M value in chat file to avoid recompile
  specify PZ cal file in chatfile
  specify ehi cal file in chatfile