Xronos

lcurve (3 Series)


Enter multiples files for each time series Multiple files for each time series can be entered using an ascii file containing the list of input filenames. This list has the following format: filenames of the individual lightcurve or event file are listed one per line and the group of files that belong to different time series are separated by three slashes. For example to run lcurve with three time series, each with four files, the input file list, input_list.txt will have the following format
     file1_ser1
     file2_ser1
     file3_ser1
     file4_ser1
     ///
     file1_ser2
     file2_ser2
     file3_ser2
     file4_ser2
     ///
     file1_ser3
     file2_ser3
     file3_ser3
     file4_ser3
where fileN_setM are the individual lightcurve or event files.

Running lcurve to create hardness and color color diagram The example uses the file list created above. The output lightcurve is binned with 300 seconds and data are organized in an interval with 219 points. When run with 3 time series lcurve calculates 2 hardness ratio. The three 300 sec bin time series, the two hardness ratio and the sum of all the time series are assembled in one array and passed to "PLT >" for plotting. By default lcurve pre-prepared a plot by answering :

 Colour-Colour [1] ; Intensity vs Time (or Phase)[ # of series (3 or 4)];
Enter PLOT style number (default=1)[1] 1
however since all the calculated data are available in "PLT >" , any other plot (e.g. hardness vs time) can be obtained by using within "PLT >" the commands to turn on-off arrays. See the example below .
lcurve 1.0 (xronos5.18)

Number of time series for this task[1] 3
Ser. 1 filename +options (or @file of filenames +options)[@sax.txt] @filelist.txt
 Series 1 file   1:file1_ser1.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 11:37:01.600
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 14:28:13.627
 No. of Rows .......        10265        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .964785340       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 1 file   2:file2_ser1.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 14:43:57.629
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 18:24:13.662
 No. of Rows .......        13209        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909090880       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 1 file   3:file3_ser1.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 21:31:57.686
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 01:47:49.708
 No. of Rows .......        15345        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909918130       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 1 file   4:file4_ser1.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6369 02:05:09.709
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 05:47:01.715
 No. of Rows .......        13305        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .945537030       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 2 file   1:file1_ser2.lc

 Selected FITS extensions: 1 - RATE TABLE;

Source ............ EXO2030+40          Start Time (d) ....  6368 11:37:01.600
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 14:28:13.627
 No. of Rows .......        10265        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .964785340       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 2 file   2:file2_ser2.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 14:43:57.629
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 18:24:13.662
 No. of Rows .......        13209        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909090880       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 2 file   3:file3_ser2.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 21:31:57.686
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 01:47:49.708
 No. of Rows .......        15345        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909918130       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 2 file   4:file4_ser2.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6369 02:05:09.709
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 05:47:01.715
 No. of Rows .......        13305        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .945537030       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 3 file   1:file1_ser3.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 11:37:01.600
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 14:28:13.627
 No. of Rows .......        10265        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .964785340       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 3 file   2:file2_ser3.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 14:43:57.629
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6368 18:24:13.662
 No. of Rows .......        13209        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909090880       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 3 file   3:file3_ser3.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6368 21:31:57.686
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 01:47:49.708
 No. of Rows .......        15345        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .983090820       .909918130       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

 Series 3 file   4:file4_ser3.lc

 Selected FITS extensions: 1 - RATE TABLE;

 Source ............ EXO2030+40          Start Time (d) ....  6369 02:05:09.709
 FITS Extension ....  1 - `RATE      `   Stop Time (d) .....  6369 05:47:01.715
 No. of Rows .......        13305        Bin Time (s) ......    1.000
 Right Ascension ... 307.624755960       Internal time sys.. Converted to TJD
 Declination ....... 37.416656430        Experiment ........ EXOSAT   ME

 Corrections applied: Vignetting - Yes; Deadtime - Yes; Bkgd - Yes; Clock - No
             values: .945537030       .911577050       1.00000000

 Selected Columns:  1- Y-axis;  2- Y-error;

 File contains binned data.

Name of the window file ('-' for default window)[-]

 Expected Start ...  6368.48404629571  (days)      11:37: 1:600  (h:m:s:ms)
 Expected Stop ....  6369.24099207078  (days)       5:47: 1:715  (h:m:s:ms)

 Minimum Newbin Time    1.0000000      (s)
 for Maximum Newbin No..            65401

 Default Newbin Time is:   127.88446    (s) (to have 1 Intv. of     512 Newbins)
 Type INDEF to accept the default value

Newbin Time or negative rebinning[500] 300

 Newbin Time ......     300.00000      (s)
 Maximum Newbin No.               219

 Default Newbins per Interval are:       219
 (giving       1 Interval of       219 Newbins)
 Type INDEF to accept the default value

Number of Newbins/Interval[79856] 219
 Maximum of       1 Intvs. with       219 Newbins of       300.000     (s)
Name of output file[saxcygx-2.flc] myoutfile.flc
Do you want to plot your results?[yes]
Enter PGPLOT device[/xw]

By setting the plot style to 1 "PLT >" plots a colour-colour diagram, setting the plot style to 3 (for 3 time series) "PLT >" plots the hardness versus time. After the default plot is done, within "PLT >" is possible to plot any axis vs any other by using the command "color on" or "color off" and "xaxis N" or "yplot M". See the example below .
 Colour-Colour [1] ; Intensity vs Time (or Phase)[ # of series (3 or 4)];
Enter PLOT style number (default=1)[1] 1

     219 analysis results per interval

 99% completed
 Intv    1   Start  6368 11:39:31
     Ser.1     Avg  18.72        Chisq 0.1161E+06   Var  124.1     Newbs.    178
               Min  9.094          Max  53.43    expVar 0.1983      Bins  51868

     Ser.2     Avg  33.13        Chisq 0.2765E+06   Var  429.5     Newbs.    178
               Min  15.30          Max  92.51    expVar 0.2874      Bins  51868

     Ser.3     Avg  51.85        Chisq 0.3866E+06   Var  1014.     Newbs.    178
               Min  25.12          Max  145.1    expVar 0.4856      Bins  51868
PLT>
color-color example 


If, however, the Intensity vs. Time plot type (3 in this case) had been chosen, the following plot would result:
3 light curve example

It is possible to use the PLT prompt to customize the output plot.

The command "info" lists all the columns available which are:

 
 1st Time
 2nd lightcurve series 1 
 3rd lightcurve series 2 
 4th lightcurve series 3 
 5th hardness ser2/ser1
 6th hardness ser3/ser2
 7th sum 1+2+3
 8th only used for lcurve run with 1 time series
 
when lcurve is run with three time series

PLT> info

 Scales:
 Grp  Wind    Label     XData Min    XData Max    YData Min    YData Max
  1    -1  Time (s)   :  2371.    ,  6.7471E+04 :  2371.    ,  6.7471E+04
  2     2  Ser 1      :  2371.    ,  6.7471E+04 :  9.094    ,   53.43    
  3     3  Ser 2      :  2371.    ,  6.7471E+04 :  15.30    ,   92.51    
  4     4  Ser 3      :  2371.    ,  6.7471E+04 :  25.12    ,   145.1    
  5    -1  Ser 2/Ser  :  2371.    ,  6.7471E+04 :  1.475    ,   1.943    
  6    -1  Ser 3/Ser  :  2371.    ,  6.7471E+04 :  1.515    ,   1.678    
  7    -1  Ser 1+2+3  :  2371.    ,  6.7471E+04 :  50.23    ,   290.3    
  8    -1             :  2371.    ,  6.7471E+04 :  1.000    ,   1.000    

For example to plot on the x-axis time and on the y-axis the two hardness and the sum lightcurve (a stack plot with 3 windows), first turn off the color of the current plotted groups and turn on the new groups (using color index of 1 -- white)
PLT> color off 2 3 4
PLT> color 1 on 5 6 7 
PLT> xaxis 1
the following command will able the stack plot (Plot Vertical) and rescale the x-axis to include the min and max value of the time axis
PLT> p v
PLT> r x
plot example

At the PLT> prompt you will be able to use any of the "QDP" commands So for example to change the above plot, such that the lightcurve-sum is plotted in the first window here are the sequence of commands to use. Many other commands to manipulate plots are availble within the qdp help file.

Use WHEAD command to see the location (LOC) values of the current plot windows. The relevant plot commands from WINFO are:

WIN      5
LOC  0 0.6083333 1 0.9416667
LAB  NX OFF
LAB  T  EXO2030+40
WIN      6
LOC  0 0.3583333 1 0.6916667
LAB  NX OFF
WIN      7
LOC  0 0.1083333 1 0.4416667
In the current plot, groups 5 through 7 are ordered from top to bottom. In the new plot, we want group 7 on top where 5 is, 5 in the center where 6 is, and 6 on the bottom where 7 is. Also, the bottom window (7) has x-axis number labelling turned on. We want to turn it off for its new location and add axis labelling to the new bottom plot (6). The title is currently attached to the top window (5), which should be removed and added to the new top window (7).
PLT> WIN      5
PLT> LOC  0 0.3583333 1 0.6916667
PLT> LAB T
PLT> WIN      6
PLT> LOC  0 0.1083333 1 0.4416667
PLT> LAB NX ON
PLT> WIN      7
PLT> LOC  0 0.6083333 1 0.9416667
PLT> LAB NX OFF
PLT> LAB  T  EXO2030+40
PLT> p
plot example 

Xronos Home Page Xanadu Home Page Xanadu ftp site

Please send reports of errors to : xanprob@athena.gsfc.nasa.gov
HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public

Last modified: Friday, 26-Mar-2004 16:35:04 EST

HEASARC Staff Scientist Position - Applications are now being accepted for a Staff Scientist with significant experience and interest in the technical aspects of astrophysics research, to work in the High Energy Astrophysics Science Archive Research Center (HEASARC) at NASA Goddard Space Flight Center (GSFC) in Greenbelt, MD. Refer to the AAS Job register for full details.