Using the RHESSI Software from the IDL Command Line

The RHESSI data analysis software is based on the concept of object-oriented programming.  The software uses object classes to stepwise reconstruct images, spectra and light curves from raw data. Each step of the construction is associated with an object class. 

An object class is a software construct that describes a data type and the methods (i.e., procedures and functions) to access and handle it. Objects make it easier to control the large collection of data products and processes that are needed to produce a final product; the final and all intermediate products, the settings used to create them, and the current state of each level (processed or in need of updating) are stored within a single entity that a user can refer to by a simple variable name, e.g. 'o' .  

On this page:

External Links:

RHESSI Object Basics

RHESSI users normally interact with the three high-level object classes - those for images, spectra, and quicklook data, shown in the following table. 

Object ClassObject CreationObject Purpose
hsi_imageo=hsi_image()Generate single images or cubes of images at multiple times and/or energies
hsi_spectrumo=hsi_spectrum()Generate spectra and lightcurves for arbitrary time and energy binning
hsi_obs_summaryo=hsi_obs_summary()    Accumulate quicklook products from the observing summary files


Users control the data generation process by setting options (called Control Parameters) such as the time interval or energy binning, and by calling methods to get or display the data.  Informational parameters (called Info Parameters) that provide information about the data product can be retrieved. These methods are shown in the following table.

MethodSyntaxPurpose
seto->set,param=valueSet a control parameter
getvalue=o->get(/param)Retrieve the value of a control or info parameter
getdatadata=o->getdata()  Generate the data (image, spectrum, whatever)
plot or plotmano->plot or o->plotmanPlot the data in a simple window or a plotman interface

 

Each object keeps track of its state.  For example, if you change a control parameter and then request a plot, the object knows that it needs to redo some part of its chain of processing in order to display the updated final product.  You don't need to call getdata first.  Each object has additional methods that are not typically called by users, or are object-specific, and are documented in the description for that object.

All of the control and info parameters available to SET and GET in the RHESSI objects are listed in the RHESSI Object Parameter Tables.  All control parameters are initialized to default values when the object is created.  You only need to set the parameters that you want to be different from the default.  To see the default value of any parameter, use the GET method, e.g. to see the default image algorithm, type

print, o -> get(/image_algorithm)

A simple example follows, showing how to use the basic methods above to reconstruct and display a RHESSI image.  We create an image object, set a time interval, energy band, and image reconstruction algorithm, and request the image to be reconstructed and plotted.  We also retrieve the image array.

o = hsi_image()
o -> set, im_time_interval= ['12-Aug-2002 02:19:19.776', '12-Aug-2002 02:19:36.224']
o -> set, im_energy_binning= [6., 12.]
o -> set, image_algorithm= 'clean'
o->plotman
image = o -> getdata()

See the RHESSI Snippets for more examples of using the objects from the CLI.

RHESSI Object Chain

In fact, each of the high-level objects is comprised of a chain of objects (>100!) that handle every operation needed to go from the raw data to the final product (view an object tree diagram).  Normally users set control parameters in the high-level object, and it passes the information down the chain until the particular object that handles that option is reached.  Similarly, info parameters can be requested from the high-level object, which will  request the information from the appropriate underlying object.  Each level object generates its own data product which is used in the next higher level object, all contributing to the final product of the top level object.

All of the RHESSI object classes are similar in structure, i.e. they all have control and info parameters, and common methods for setting and retrieving them (GET and SET), as well as for requesting and displaying the final data product (GETDATA, PLOT, and PLOTMAN, if plotting is appropriate).

For many operations, users need not be concerned with the chain of underlying objects.  However, it is possible to

  • extract data from one of the lower-level objects,
  • extract a lower-level object and work with it directly, or to
  • create a stand-alone instance of any of the underlying objects.  

For example, in order to reconstruct an image, the image object needs pointing and roll angle information which it gets from an underlying aspect solution object.  If you wanted to retrieve the aspect data yourself, you could type (assuming o is your image object):

aspect = o -> getdata(class='hsi_aspect_solution')

or, the equivalent

aspect_obj = o -> get(/obj, class='hsi_aspect_solution')
aspect_data = aspect_obj -> getdata()

or to create a stand-alone aspect solution object:

aspect_obj = hsi_aspect_solution()
aspect_obj -> set, obs_time_interval=['12-Aug-2002 02:19:19.776', '12-Aug-2002 02:19:36.224'])
aspect_data = aspect_obj -> getdata()

Special GET calls

In addition to simply getting the value of a control or info parameter, you can also do the following:

params = o -> get()     ; returns a structure containing all control and info parameters
control = o -> get(/control)  ; returns a structure containing all control parameters
info = o -> get(/info)   ; returns a structure containing all info parameters
obj = o -> get(/obj, class='hsi_xxx')  ; returns the lower-level object called hsi_xxx

Defaults

The RHESSI objects are initialized to default values for all control parameters. You can override the defaults in a number of ways - in the object creation command, using the set command, running scripts, setting up a text file of customized defaults, and more.

On rare occasions, the software developers change some of the object initialization default values. You can force the use of the defaults in effect before that change.

Please see RHESSI OBJECT PARAMETER DEFAULTS for more details about these options.