- yuxip's home page
- Posts
- 2015
- 2014
- December (2)
- November (1)
- October (2)
- September (6)
- August (2)
- July (1)
- June (3)
- May (4)
- April (1)
- March (2)
- 2013
- December (1)
- November (1)
- October (3)
- September (3)
- August (2)
- July (1)
- June (2)
- May (1)
- April (3)
- March (1)
- February (1)
- January (1)
- 2012
- 2011
- My blog
- Post new blog entry
- All blogs
estimate gain factors by adc spectrum
1. dataset: 6 runs on day080, 6 runs on day095
I combined runs that belong to the same day, and plot adc spectrum for each one of the 18 cells in group 1. Then I compared the slopes of a exponential fit to the region above threshold:
www.star.bnl.gov/protected/spin/yuxip/adcCellslopes.pdf
Fit was from ADC channel 600 to 2000. There are two frames for each cells, corresponding to the above two datasets respectively.
2. Following is the math that shows slope1 / slope2 = gain1 / gain2:
Assume the real energy spectrum obeys exponential distribution;
N = N0*exp(-k*E), with the factor k unknown, N is the number of events observed in energy bin E. Using E = ADC*gain as in our convention, then
N = N0*exp(-k*ADC*gain).
Plotting N as a function of ADC, then the slope above threshold (shoulder) region is equal to k*gain. Compare the slopes of two datasets, we can show that
slope1 / slope2 = gain1 / gain2. The ratio gain1 / gain2 can be equivalently calculated by the number of triggers above a certain threshold ADCth, as a fraction of the
total number of events, namely
ln(N/N0) = -k*ADCth*gain.
By fixing the same ADCth for each dataset, and compute ln(N/N0) in which N the the number of events above ADCth, we can show that
ln(N1/N10) / ln(N2/N20) = gain1 / gain2
3. Here is a comparion of the gains. www.star.bnl.gov/protected/spin/yuxip/adcCellslopes_1_095080_6.pdf
For each cell, there are 4 types of gain factors.
gcorrOld: gain correction factors calibrated for run12080001
gcorrAdc: calculated by the above method using slopes. ( from adc spectra of 6uns on day080 and 6 runs on day095)
gpredicted: by following the mean ADC variation of LED events. (with 6 run from day095 combined, then compared to the baserun 12080001 )
gpredicted = nomianl gcor for baserun*( 1st pink point / mean of the rest 6 pink points), where the 1st pink point is the mean LED adc of
baserun
gdynamic: following the mean ADC variation of LED events, but taking into account of the change of pulser intensity, which was described by the blue
dots. gdynamic = nominal gcor for baserun* (mean of 6 blue points / mean of 6 pink points) , 6 pink points are the ones from day095
My conclusion is that for those cells which have enough number of triggers above threshold, therefore a relatively good exponential fit, the gain factors calculated
from slopes of adc spectra are more close to gdynamic.
ps. Here is the LED channel map: www.star.bnl.gov/protected/spin/yuxip/ledmap/
and the related macros: www.star.bnl.gov/protected/spin/yuxip/macro/. This can be explained later, if necessary
- yuxip's blog
- Login or register to post comments