Lab 1
 
time series analysis



AIM
To get familiar with some common timing analysis tools, such as power density spectra, autocorrelation functions and epoch folding. Simulate and investigate the appearence of a few different types of signals and study the influence of noise. Find, download and analyse archive data from RXTE all sky monitor, ASM. ASM is continously scanning the sky in soft X-rays. Each day most bright X-ray sources are measured a few times. These data are available online and can therefore be downloaded and analysed.

Background
Fourier transform, autocorrelation function, epoch folding, detrending, Monte Carlo simulation, chi square, error estimation.

Linux
If you are new to Linux/Unix take a look at this command summary. For an online unix course.

To start the DES7 program.
Open a second terminal window and cd to DES7 in both windows. Run the program in one and idl in the other by typing in window 1 and 2 respectively:
1.     des7
2.     idl
        desdraw
 
Before the lab you should take a look at the online des7 manual.

1. Run simulations and compute power density spectra
  In each case do the following sequence 34. Simulation, 3. Detrend 4. 16. Power spectrum
  Compute one power spectrum without detrending first. How and why does the power spectrum differ?

- Poisson noise
- Coherent oscillation (sinusoid)
   Try a few different amplitude - noise level combinations and two different data lengths. 
 - Amplitude modulated sinusoid
 - A few cases with different parameter settings of either a quasi-periodic oscillation
    (2nd order autogressive process) or shot noise.

2. Epoch folding
-  Simulate a sinusoid in channel2 (x2) (in simulation store errors, internal) and make an epoch folding analysis. (Option 24, then 31. DES-Folding in the submenue, but remember to set folding period [Note: Period not frequency])

   Plot chi-square and pulse profile (In "40. View ...."). Rerun with Monte-Carlo simulation option set to at
    least 10 to estimate errors for the period and amplitude obtained from the analysis.
    (Note: Disregard the messages written during execution (for testing purposes only).
     Compare with a period determination from a power density spectrum. 

3. Period analysis of RXTE/ASM data 
     Go to the ASM lightcurve page.
    Click on "data downloads" in the left column.
     Select a source and possibly a start and stop time. (You may e.g. select an accreting X-ray pulsar)
    Make sure the "No Header" option is set.
    Click on the button "Retrive Entire Light Curve Files via HTTP links"
    In the file list window that opens, click on the file you want to download with the right mouse button and select "Save link target as".
    Now use the tools above to analyse the periods that can study from these observations in your particular source.