Analysis codes: Difference between revisions
Tlauritsen (talk | contribs) No edit summary |
Tlauritsen (talk | contribs) No edit summary |
||
Line 31: | Line 31: | ||
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. | 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. | 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 sort again (resort) using the new pz values that were extracted above. | Now sort again (resort) using the new pz values that were extracted above. |
Revision as of 20:34, June 12, 2018
GEBSort_nogeb
GEBSort_nogeb is the program that can analyze data from DGS, DFMA and GRETINA.
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 sort again (resort) 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 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.