Pulsar Workshop MAGIC pulsar workshop Padova Feb 2010
Pulsar Workshop MAGIC pulsar workshop, Padova, Feb 2010 1
Pulsar Timing Analysis MAGIC pulsar workshop, Padova, Feb 2010 2
/hadron separation (I) MAGIC n Different kind primary particles different showers different images n /hadron separation based on image parameter distributions • -images are smaller and point to camera center • Hadron showers are broader are randomly oriented MAGIC pulsar workshop, Padova, Feb 2010 3
/hadron separation (II) MAGIC n After applying /hadron cuts based on image shape, exploit shower direction MAGIC pulsar workshop, Padova, Feb 2010 4
Timing analysis MAGIC Goal: Find the periodic signal of the pulsar, hidden in the arrival times of the atmospheric showers n The timing analysis involves 4 steps: • Barycenter correction • Obtain the Light curve • Application of Uniformity test • Upper limits calculation All these steps have been implemented in a dedicated software, for the pulsar Analysis in MAGIC pulsar workshop, Padova, Feb 2010 5
Timing analysis (I) MAGIC Barycenter correction n n Remove the effect of the earth movement on the arrival times t. UTC. Transform the measured arrival times to the Solar System Barycenter: Total Barycenter correction MAGIC pulsar workshop, Padova, Feb 2010 6
Timing analysis (II) MAGIC Ligth curve n If F is the known rotational frequency of the pulsar at time T 0, the number of revolutions in dt=t-T 0 is: Taylor n Integrating, and taking the fractional part, we get the rotational phase : where t is the barycenter time MAGIC pulsar workshop, Padova, Feb 2010 7
Timing analysis (III) n MAGIC Ephemeris are usually taken from radio observations but affected by irregularities in pulsar rotation: Timing noise Glitches Crab pulsar Crab is “glitching” once every 3 -6 years ! Need to have contemporaneous ephemeris We use monthly ephemeris by Jodrell Bank MAGIC pulsar workshop, Padova, Feb 2010 8
Timing analysis (III) MAGIC Uniformity tests n If no periodicity events will be uniformly distributed ² - test – Applied to k-bins histograms – Powerful test for narrow pulses – Probability follows a 2 of k-1 d. o. f Zm 2 -test • Independent of k H test • Find optimal m • Based on trigonometric moments of pulse profile m=number of harmonics • Powerful against a large variety of light curves • Sensitive to narrower pulses as m increases 9 MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC n Developed a full pulsar analysis toolkit for MAGIC n You can do all you need, both for optical & gamma-ray data n Input files can be ascii, raw or root files: – – – Extract Central Pixel signal Do barycenter correction Search best ephemeris for each arrival time Calc pulsar phase & create light curves Apply periodicity tests And moreover. . v v v 10 v Apply cuts Calc. collection areas Calc. upper limits Calc. spectrum MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC n Code at MARS/mtemp/mucm/mpulsar n More that 30000 lines of code n Manual at: MARS/mtemp/mucm/mpulsar/docmpulsar. pdf 11 MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC How to compile it – Compile the classes adding the following lines to the Makefile: SUBDIRS = manalysis mastro … mtemp/mucm/classes � mtemp/mucm/pulsar – Compile the Search. Period program 12 PROGRAMS = callisto celestina … Search. Period MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC n Main Classes 13 MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC How to use it n n You can use the classes writing your macro. E. g. : mtemp/mucm/mpulsar/macros Or using the compiled program Search. Period (similar to callisto, . . ) • It has a basic GUI (search. tcl) 14 MAGIC pulsar workshop, Padova,
MAGIC Timing Analysis Toolkit MAGIC Tests n Toolkit tested on: – Crab optical signal – EGRET data – XMM data n Results checked against TEMPO 15 MAGIC pulsar workshop, Padova,
- Slides: 15