Which detectors can I use to analyze this flare?

From RHESSI Wiki

(Difference between revisions)
Jump to: navigation, search
m (Added a fourth plot and renumbered figures)
m (Added a fourth plot, renumbered figures, and edited text)
Line 12: Line 12:
== Introduction ==
== Introduction ==
   
   
-
This is the most common question we are asked by people trying to analyze RHESSI imaging and/or spectroscopy data.  The answer is not always straightforward and it is becoming more difficult as the individual detectors respond differently to radiation damage over time and to the four detector anneals completed so far (Nov 2007, Mar 2010, Jan 2012, and Jun 2014).
+
This is the most common question we are asked by people trying to analyze RHESSI imaging and/or spectroscopy data.  The answer is not always straightforward and it is becoming more difficult as the individual detectors respond differently to radiation damage and the four detector anneals completed so far (Nov. 2007, Mar. 2010, Jan. 2012, and 26 June to 13 August 2014).
 +
 
 +
[[Image:Hsi_sepdet_spectrum_20140903_133928to134028.png|300px|thumb|left|text-top|'''Figure 1''': Total count and photon flux spectra for individual detector front segments during the M-class flare on 03 September 2014 at 13:40 UT.]]
    
    
-
We now have quicklook front-segment spectral plots for times throughout the mission that can help answer this question.  An example for a recent M-class flare is shown in Figure 1. These plots should help you decide which detectors are OK to use for specific applications such as light curves, spectroscopy, and imaging for a specific time/flare.The plots are available in Browser via the 'Detector Spectra' check-box under the 'Flare Quicklook' option (click the magnifying glass to expand).  They show the count-flux spectrum and the corresponding 'semi-calibrated' photon spectrum for the nine detector front segments for a one minute interval centered on the peak time of each contiguous event in the RHESSI flare catalog. Use the flare-seeking buttons on Browser to step through the data by flare. These plots should help you decide which detectors are OK to use for specific applications such as light curves, spectroscopy, and imaging for a specific time/flare.
+
We now have quicklook front-segment spectral plots for times throughout the mission that can help answer this question.  An example for the recent M-class flare on 03 September is shown in Figure 1. These plots can help you decide which detectors are OK to use for the analysis of data a specific flare for various applications such as light curves, spectroscopy, and imaging. The plots are available in Browser via the 'Detector Spectra' check-box under the 'Flare Quicklook' option (click the magnifying glass to expand).  They show the count-flux spectrum and the corresponding 'semi-calibrated' photon spectrum for the nine detector front segments for a one minute interval centered on the peak time of each contiguous event in the RHESSI flare catalog. Use the flare-seeking buttons on Browser to step through the data by flare. These plots should help you decide which detectors are OK to use for specific applications such as light curves, spectroscopy, and imaging for a specific time/flare.
-
 
+
-
[[Image:Hsi_sepdet_spectrum_20140903_133928to134028.png|300px|thumb|left|text-top|'''Figure 1''': Total count and photon flux spectra for individual detector front segments during a flare on 03 September 2014 at 13:40 UT.]]
+
== Description of Plots ==
== Description of Plots ==
 +
Each plot for a given flare is divided into two sections - an upper panel showing the count-flux spectrum for each of the nine detector front segments, and a lower panel with the corresponding 'semi-calibrated' photon flux spectra.
 +
 +
The count flux spectra are in units of counts cm<sup>-2</sup> s<sup>-1</sup> keV<sup>-1</sup> corrected for decimation and live time but with no background subtraction.  The energy bin widths are 0.3 keV below 10 keV and increase logarithmically to 3 keV at 100 keV to give a total of 138 bins.
-
The upper panel shows the count flux spectrum of each front segment in counts cm<sup>-2</sup> s<sup>-1</sup> keV<sup>-1</sup> corrected for decimation and live time with no background subtraction.  Energy bins are 0.3 keV wide below 10 keV, increasing logarithmically to 100 keV.
+
The 'semi-calibrated' photon flux spectra are in units of photons cm<sup>-2</sup> s<sup>-1</sup> keV<sup>-1</sup>.  They are computed by multiplying the corresponding count flux spectrum by the diagonal elements of the detector response matrix for the attenuator state at that time  
 +
[http://sprg.ssl.berkeley.edu/~tohban/nuggets/?page=article&article_id=27]. Since no background was subtracted, this is equivalent to assuming that all the counts were from solar photons and that the off-diagonal elements are negligible. Photon spectra determined in this way approximate the actual solar spectrum only in the energy range where the solar counts are significantly above background and the off-diagonal elements are small ((typically ~10 to 100 keV for large flares). These photon plots serve primarily to show the relative sensitivity of the front segments and how well the detector response matrices correct for the known sensitivities of the different detectors.  More accurate spectral analysis with background subtraction and using all elements of the detector response matrix for each detector requires the use of OSPEX
 +
[http://http://hesperia.gsfc.nasa.gov/rhessi2/home/software/spectroscopy/spectral-analysis-software/ link OSPEX].
-
The lower panel shows the 'semi-calibrated' photon flux spectrum of each front segment in photons cm<sup>-2</sup> s<sup>-1</sup> keV<sup>-1</sup>, computed by multiplying the corresponding count flux spectrum by the diagonal elements of the detector response matrix
+
Interpreting these plots depends on the attenuator state at the time the spectrum was obtained.  (Attenuators are aluminum disks that can be automatically brought into the field of view of the detectors to reduce the effects of pulse pile-up during the largest flares <ref>The RHESSI Spectrometer, Smith et al. 2002, Sol. Phys. 210: 33-60)</ref>.)  The attenuator-state dependent cross-hatched energy range on each plot indicates the lower energy range that generally should not be used for further data analysisThe upper energy limit (not shown on these plots) depends on the flare intensity and the energy at which the solar flux dips below the background level.  
-
[http://sprg.ssl.berkeley.edu/~tohban/nuggets/?page=article&article_id=27]. Since no background was subtracted, this is equivalent to assuming that all the counts were from solar photons and that the off-diagonal elements are negligiblePhoton spectra determined in this way approximate the actual solar spectrum only in the energy range where the solar counts are significantly above background and the off-diagonal elements are small. These photon plots serve primarily to show the relative sensitivity of the front segments and how well the detector response matrices correct for the known sensitivities of the different detectors.
+
-
The spectra differ considerably in the three attenuator states.  Attenuators are aluminum disks that can be brought into the field of view of the detectors to prevent saturation during the largest flares. The cross-hatched energy range on the plots indicates the portion of the spectrum that should not be used depending on the attenuator state.
+
Three examples are given below to illustrate how to interpret the plots in the three possible attenuator states - A0 when no attenuators are in place above the detectors, A1 when the thin attenuators are above the detectors, and A3 when both the thick and thin attenuators are above the detectors. (Note that the A2 state with just the thick attenuators above the detectors is not used.)
*'''Attenuator state A0 (no attenuators)'''  An estimate of the photon spectrum can be made down to ~3 keV, the approximate energy of the electronic lower threshold level for each detector and the effective cutoff from the absorption of material in front of the detectors including the thermal blankets and the beryllium windows on the cryostat.  Counts recorded below ~3 keV are electronic noise and should not be used in the determination of the photon spectrum.
*'''Attenuator state A0 (no attenuators)'''  An estimate of the photon spectrum can be made down to ~3 keV, the approximate energy of the electronic lower threshold level for each detector and the effective cutoff from the absorption of material in front of the detectors including the thermal blankets and the beryllium windows on the cryostat.  Counts recorded below ~3 keV are electronic noise and should not be used in the determination of the photon spectrum.
*'''Attenuator state A1 (thin attenuator)'''  The solar count flux spectrum peaks at ~11 keV and falls dramatically below 6 keV due to the overlying material, which now includes the thin attenuator. Counts below ~6 keV should not be used.
*'''Attenuator state A1 (thin attenuator)'''  The solar count flux spectrum peaks at ~11 keV and falls dramatically below 6 keV due to the overlying material, which now includes the thin attenuator. Counts below ~6 keV should not be used.

Revision as of 22:37, 8 September 2014


Nugget
Number: 23?
1st Author: Brian Dennis
2nd Author: Kim Tolbert
Published: September, 2014
Next Nugget: TBD
Previous Nugget: [1]
List all



Contents

Introduction

This is the most common question we are asked by people trying to analyze RHESSI imaging and/or spectroscopy data. The answer is not always straightforward and it is becoming more difficult as the individual detectors respond differently to radiation damage and the four detector anneals completed so far (Nov. 2007, Mar. 2010, Jan. 2012, and 26 June to 13 August 2014).

Figure 1: Total count and photon flux spectra for individual detector front segments during the M-class flare on 03 September 2014 at 13:40 UT.

We now have quicklook front-segment spectral plots for times throughout the mission that can help answer this question. An example for the recent M-class flare on 03 September is shown in Figure 1. These plots can help you decide which detectors are OK to use for the analysis of data a specific flare for various applications such as light curves, spectroscopy, and imaging. The plots are available in Browser via the 'Detector Spectra' check-box under the 'Flare Quicklook' option (click the magnifying glass to expand). They show the count-flux spectrum and the corresponding 'semi-calibrated' photon spectrum for the nine detector front segments for a one minute interval centered on the peak time of each contiguous event in the RHESSI flare catalog. Use the flare-seeking buttons on Browser to step through the data by flare. These plots should help you decide which detectors are OK to use for specific applications such as light curves, spectroscopy, and imaging for a specific time/flare.

Description of Plots

Each plot for a given flare is divided into two sections - an upper panel showing the count-flux spectrum for each of the nine detector front segments, and a lower panel with the corresponding 'semi-calibrated' photon flux spectra.

The count flux spectra are in units of counts cm-2 s-1 keV-1 corrected for decimation and live time but with no background subtraction. The energy bin widths are 0.3 keV below 10 keV and increase logarithmically to 3 keV at 100 keV to give a total of 138 bins.

The 'semi-calibrated' photon flux spectra are in units of photons cm-2 s-1 keV-1. They are computed by multiplying the corresponding count flux spectrum by the diagonal elements of the detector response matrix for the attenuator state at that time [2]. Since no background was subtracted, this is equivalent to assuming that all the counts were from solar photons and that the off-diagonal elements are negligible. Photon spectra determined in this way approximate the actual solar spectrum only in the energy range where the solar counts are significantly above background and the off-diagonal elements are small ((typically ~10 to 100 keV for large flares). These photon plots serve primarily to show the relative sensitivity of the front segments and how well the detector response matrices correct for the known sensitivities of the different detectors. More accurate spectral analysis with background subtraction and using all elements of the detector response matrix for each detector requires the use of OSPEX link OSPEX.

Interpreting these plots depends on the attenuator state at the time the spectrum was obtained. (Attenuators are aluminum disks that can be automatically brought into the field of view of the detectors to reduce the effects of pulse pile-up during the largest flares <ref>The RHESSI Spectrometer, Smith et al. 2002, Sol. Phys. 210: 33-60)</ref>.) The attenuator-state dependent cross-hatched energy range on each plot indicates the lower energy range that generally should not be used for further data analysis. The upper energy limit (not shown on these plots) depends on the flare intensity and the energy at which the solar flux dips below the background level.

Three examples are given below to illustrate how to interpret the plots in the three possible attenuator states - A0 when no attenuators are in place above the detectors, A1 when the thin attenuators are above the detectors, and A3 when both the thick and thin attenuators are above the detectors. (Note that the A2 state with just the thick attenuators above the detectors is not used.)

Note the date and time of plot creation in the lower right corner. The detector response matrices may be updated as we learn more about the variation of the detector responses with time. In cases where the photon spectrum is incorrect for a particular detector, you could check whether an improved response matrix has been generated since these quicklook plots were created by remaking these plots yourself in IDL.

Figure 2: Total count and photon flux spectra for individual detector front segments during a flare on 02 May 2014 at 00:54:36 UT.

Example 1 in the A0 Attenuator State

An example of these plots is shown in Figure 2 for a small flare in the A0 attenuator state.

The count flux plot shows immediately which detectors have recorded an increase in count rate from the flare. The solar signal is the peak at around 7 keV that falls off at lower energies down to 3 keV and extends up to ~25 keV before merging with the background spectrum at higher energies. Here is a breakdown of what you can determine about each of the nine detectors from these two plots:


Figure 3: Similar to Figure 1 during a C-class flare on 03 September 2014 at 02:52 UT.

Example 2 in the A1 Attenuator State

An example of a flare in the A1 attenuator state after the June 2014 anneal is shown in Figure 3.

After this anneal, only five detectors were segmented. Consequently, the spectrum of solar photons can only be seen for those five detectors - 1F, 3F, 6F and 8F, and 9F. The Fe-line complex peak at 6.7 keV can be seen in the spectra from the detectors with the best energy resolution. The background spectrum dominates at energies above ~25 keV in this example.


Figure 4: Similar to Figure 1 during an X-class flare on 20 January 2005 at 06:50:40 UT.

Example 3 in the A3 Attenuator State

An example of an intense flare in the A3 attenuator state on 20 January 2005 is shown in Figure 4.

The close alignment of the photon spectra for all detectors indicates that they can all be used above ~6 keV although there are still small differences between detectors as indicated below. As with the example in attenuator state A1, counts below ~6 keV are almost all K-escape events and should not be used.

Personal tools
Namespaces
Variants
Actions
Navigation
Toolbox