Command-Driven walk through

Move to the partial command-driven user interface by typing:

[16]Xronos> cui c
Current user interface is :            command driven

Input filename(s) are entered by using the data command, as in the partial question/answer user interface. The same input files as in the previous section are used in the examples below. To check the specified input file type:

[17]Xronos> data/show
 
Series A :  file  1 : C2030
 
Series B :  file  1 : A2030
 
Series C :  file  1 : B2030

lc1

In order to use file a2030.rbf (SERIES B) and repeat the first lc1 example you should type:

[18]Xronos> lc1 b

The difference in this case is that, after the input filenames, all default values are taken automatically and lc1 is executed. The following summary is output:

Maximum of      1 Intvs. with       512 Newbins of       1.45349     (s)
 
Intv    1   Start  302 15:59: 6
    Ser.1     Avg  42.83        Chisq  3685.       Var  329.1     Newbs.   511
              Min -7.916          Max  117.1    expVar  49.47      Bins    743
            Light curve ready !
Frme    1   Written to output file : xronos.qlc1

before the light curve is plotted and the XRONOS prompt is issued again.

In the command-driven user interface the execution of a XRONOS application is controlled by QUALIFIERS. To display the list of allowed qualifiers for lc1 type:

[19]Xronos> lc1/?
 
 Allowed qualifiers are:
 
/CMDQDPFILE    /FULLINFO      /HEADER        /IMAGE
/LPARM         /NBINS         /NINT          /NOECHO
/NOPLOT        /OUTFILE       /PLT           /TNEWBIN
/TREND         /TRENDMODE     /TUNIT         /WINDOW

The /FULLINFO qualifier will run lc1 displaying the same information as in the question/answer user interfaces; type:

[20]Xronos> lc1/fullinfo b

The /HEADER qualifier will display the same information, but without executing the lc1 command; this is especially useful for checking the lc1 set up in the command-driven user interface:

[21]Xronos> lc1/header b
>>>> Program lc1          9-JUL-1992 13:25:55         Xronos vers. 3.00 Jul 92
 
     Enter up to  50 input filenames and options for series 1 (or rtn)
lc1> Ser. 1 filename   1 +options[.rbf]=> A2030
 
Source ............EXO 2030+40    +   Start Time ........1985 302  15:59: 4
Bin Time (s) ......   1.000           Stop  Time ........1985 302  16:11:20
No. of Points .....         743       Experiment ........   EXOSAT ME
Instrum., OBC rec..     AR , E5       Detectors on Source       5-  8
PHA Bin range .....       2- 16       Detectors on Bkgd..       1-  4
 
lc1> Ser. 1 filename   2 +options[.rbf]=>
lc1> Window Options ?          {y/n}[n]=>                              ! WINDOW
 
   0 Time,    0 Phase,    0 Intensity,    2 Exposure Windows
 
Expected Start ...   302.66604129903  (days)      15:59: 5:968  (h:m:s:ms)
Expected Stop ....   302.67464446032  (days)      16:11:29:281  (h:m:s:ms)
Minimum Newbin Time    1.0000000      (s)
Maximum Newbin No..              744
Default Newbin Time    1.4534868      (s) (to have 1 Intv. of    512 Newbins)
 
lc1> Enter Newbin Time (s) or neg.rebinning
     of Min. Newbin Time [   1.453    ]=>                        ! TNEWBIN,NINT
 
Newbin Time ......     1.4534868      (s)
Maximum Newbin No.               512    (or      1 Intvs. of       512 Newbins)
 
lc1> No. of Newbins/Intv.  (max: 65536)
     {value or neg.pow.of 2}[      512]=>                               ! NBINS
 
Maximum of      1 Intvs. with       512 Newbins of       1.45349     (s)
 
lc1> Order of polynomial trend removal
           {0=none,1=1st,....,4=4th}[0]=>                               ! TREND
lc1> Time axis units
       {0=s from start,1=s,2=h,3=d} [1]=>                               ! TUNIT
lc1> Enter output filename
     [xronos.qlc1                     ]=>                       ! OUTFILE,IMAGE
 
lc1> Plot options ?            {n/y}[n]=> y
lc1> Plot results ?            {y/n}[y]=> y                            ! NOPLOT
lc1> Interactive or automatic? {i/a}[i]=> a                               ! PLT
lc1> Pgplot device   [/vt125          ]=>
lc1> Specify QDP command file? {y/n}[n]=>                          ! CMDQDPFILE

The /LPARM (list parameters) qualifier does the same thing, except that the input file(s) information is not displayed. Note that the qualifiers controlling the replies to individual questions are indicated on the right hand side after the exclamation mark.

Now run lc1 on the same file, this time by using a NEWBIN time of 2 s and INTERVALS of 100 newbins (i.e. 200 s); add the /LPARM to check that the set up is correct:

[22]Xronos> lc1/tnewbin=2./nbins=100/lparm b

If you are satisfied with the setup remove the /lparm qualifier and execute lc1 (note that qualifiers can be abbreviated as long as the remain unambiguous):

[23]Xronos> lc1/tn=2./nb=100/lparm b

In order to leave QDP/PLT in the interactive plotting mode (PLT> prompt) after each plot and allow customisation, add the /PLT qualifier:

[24]Xronos> lc1/tn=2./nb=100/plt b

psd

Now calculate the power spectrum of file c2030.rbf (SERIES A) with application psd

[25]Xronos> psd/fullinfo a

Repeat by using a NEWBIN time of 3 s and INTERVALS of 64 : first check that the setup is correct:

[26]Xronos> psd/tn=3./nb=64/lparm a

then remove /LPARM and run psd (note that a, the default, can be omitted)

[27]Xronos> psd/tn=3./nb=64

The spectra corresponding to each of the 4 intervals are averaged to produce the power spectrum of a FRAME (i.e. 1 FRAME of 4 INTERVALS); this is written to file xronos.qpsd and plotted automatically with QDP.

Repeat by producing and plotting interactively a power spectrum for each interval (i.e. 4 FRAMES, each of 1 INTERVAL):

[28]Xronos> psd/tn=3./nb=64/plt/eachpsd

Repeat by producing an average power spectrum for each group of two intervals (i.e. 2 FRAMES each of 2 INTERVALS):

[29]Xronos> psd/tn=3./nb=64/plt/navpsd=2

Now obtain again a single power spectrum from the 1 s resolved light curve (i.e. NEWBINS of 1 s) in a single 1024 newbin INTERVAL (i.e. 512 Fourier frequencies) and decrease artificially the Fourier resolution by a factor of four by rebinning the power spectrum; first check the setup:

[30]Xronos> psd/rebinpsd=4/lparm

and then run psd:

[31]Xronos> psd/rebinpsd=4

Now rebin the power spectrum in a differential way such that the frequency bins are nearly equispaced in a log-frequency plot. To do this specify the (negative) step of the geometric series by typing:

[32]Xronos> psd/rebinpsd=-1.2

Alternatively you might specify the dedicated qualifier /logrebinpsd, which will use a geometric series step of about 1.19 (producing a spacing of about 1/4 of octave):

[33]Xronos> psd/logrebinpsd

lc2

Now plot the soft and hard light curves of the 42 s pulsar EXO2030+375 and their ratio with lc2: the corresponding files are a2030.rbf (SERIES B) and b2030.rbf (SERIES C).

[34]Xronos> lc2 b c
 
Maximum of      1 Intvs. with       512 Newbins of       1.45349     (s)
 
 
Intv    1   Start  302 15:59: 6
    Ser.1     Avg  42.83        Chisq  3685.       Var  329.1     Newbs.   511
              Min -7.916          Max  117.1    expVar  49.47      Bins    743
 
    Ser.2     Avg  35.58        Chisq  3182.       Var  283.1     Newbs.   511
              Min-0.1958          Max  106.8    expVar  48.55      Bins    743
            Light curves ready !
Frme    1   Written to output file : xronos.qlc2

In order to leave QDP/PLT in the interactive plotting mode (PLT> prompt) and plot the hardness ratio of the two light curves vs. the sum of the two light curves, type:

[35]Xronos> lc2/hdns/plt b c

In order to plot in an automatic fashion the sum of the light curves and their ratio vs. time type:

[36]Xronos> lc2/sumandhdns b c

Next: Commands Up: Walk through Previous: Partial Question/Answer walk

Lorella Angelini
Thu 12 Oct 16:35:19 1995