Fitting a new detector

From RHESSI Wiki

Revision as of 19:00, 8 July 2009

Fitting a Detector Using Results from Previous Fits

Introduction

This section goes through the procedure for fitting a spectrum file created for a single detector using a script created after fitting another detector for the same event. In order to create single detector spectrum and SRM files, the SSW routine hsi_spectrum_sep_det_files.pro can be used.

Prerequisites for this Procedure

• A single detector spectrum and srm files for the detector to be fit. It should meet the following requirements:
  • Created over the same time period and energy range as the files for the detector already fit
  • Created with the same time bins and energy bins as the files for the detector already fit
• An OSPEX script and results file for the detector already fit

Using Previous Fits

The first step is editing the script created for the previous detector. Using your preferred IDL editor (the IDLDE is a common choice) open the script. The lines in the script specifying the spectrum and SRM files need to be edited to point to the new files to be used. For example, below are the original lines from a script created after fitting a flare on 10-Nov-2004 using the spectrum file for detector four:

       obj-> set, $                                                                                             
        spex_specfile= '/Users/andygopie/RHESSI_Flares/10_Nov_2004/hsi_spectrum_20041110_000000_d4.fits'
       obj-> set, $                                                                                             
        spex_drmfile= '/Users/andygopie/RHESSI_Flares/10_Nov_2004/hsi_srm_20041110_000000_d4.fits'

Note these lines give the path to the local location of the original spectrum and SRM files for detector four. In your script they will point to the location of the files on your local drive. The actual syntax of the path will depend on your operating system. To analyze the a different detector change the path to point to the spectrum file and SRM file for the new detector:

        obj-> set, $                                                                                             
         spex_specfile= '/Users/andygopie/RHESSI_Flares/10_Nov_2004/hsi_spectrum_20041110_000000_d5.fits'
        obj-> set, $                                                                                             
         spex_drmfile= '/Users/andygopie/RHESSI_Flares/10_Nov_2004/hsi_srm_20041110_000000_d5.fits' 

These lines now set the spectrum and SRM files to detector five.