Analysis codes: Difference between revisions

From GammaSphere DAQ
Jump to navigation Jump to search
No edit summary
No edit summary
Line 42: Line 42:
   ''' **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.
   ''' **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 sort again (resort) using the new pz values that were extracted above.  
Now,
 
  remove the energy calibration file: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.
When you resort, the PZ values in dgs_pz.cal are read in and used.
Extract the new clean ehi spectra as  
Extract the new clean, uncalibrated, ehi spectra as  


   .x get_ecln.cc
   .x get_ecln.cc
Line 56: Line 61:


'''The energy processing in bin_dgs follows algorithms that were developed by Shoufei Zhu.'''
'''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

Revision as of 21:41, June 12, 2018


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 this command

  .x get_pz.cc

Now run:

  dgs_pz 350 100 dgs_pz.cal

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 values, in this case: 0.20+0.80 = 1.00 us or 100 in 10 nsec units. The M value is 3.50 us or 350 in 10 nsec units.

  **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: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