Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng,...

14
Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan Ann Arbor, MI 48109-1040 Looking for Needles in Haystacks - Using the Swift/XRT and the Fermi/GBM to uncover GRB γ-rays in the Fermi/LAT detector F in d in g

Transcript of Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng,...

Page 1: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay

Physics DepartmentUniversity of Michigan

Ann Arbor, MI 48109-1040

Looking for Needles in Haystacks-

Using the Swift/XRT and the Fermi/GBM to uncover GRB γ-rays in the Fermi/LAT detector

Finding

Page 2: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

GRBs and High Energy γ–rays

• Max Eγ > 20 GeV (but limited by γγ→ e+e-)

• Both short and long bursts emit >100 MeV photons

• Prediction prior to launch: ~200 LAT GRB detections/year

• Fermi/GBM detections: ~250/year

• Fermi/LAT detections: ~9/year

• LAT fluence ~ GBM fluence

• GBM localization ~ 10°

• LAT localization ~ 1°

• Swift XRT localization ~ 1˝

Challenge: With better & faster localizations, more bursts could be optically detected and better characterized.

“Haystack I” (XRT): arXiv:1010.1436v2 [astro-ph.HE], ApJ 725 (2010)“Haystack II” (GBM): arXiv:1010.1588v2 [astro-ph.HE], ApJLett (2011)

Page 3: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Using XRT localizations to find LAT GRB γ-rays

Matched filter technique: w ~ p(signal) / p(background)

Modified matched filter weight for each GRB event: iw

2

22 ( )2 Ew e

, , , E t cw E t c w w w w

( ) ( )/ ; 1 3(3) (3)

GRB backc

GRB back

r i r iw ir r

0

0 bb

t

b cb

tc t tt

wtc t t tt

1/30 1 2

0 1 2

3 w w ww w w

2

14 ( )

back GRB

Eth

EwE E

Page 4: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Data Selection:

DATA_QUAL = 1

IN_SAA = 0

100 MeV < Eγ < 300 GeV

gbGRB > 10°

Θzenith < 105°

ΘGRB-bore < 66°

θγ < 10.5° relative to XRT localization

0 < t < 47.5 s relative to GRB trigger

41 fields selected:

6 with previous LAT GRB detections35 with no prior evidence for LAT GRBs

Circular 10.5° tiling of the LAT FoVfor defining random fields. Dotted circle at 68° to boresight.

Page 5: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

LAT photon event class counts for XRT-localized fields

transie

ntso

urcediffu

se

Page 6: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Cumulative distributions of modified matched filter weights, blue line represents 35 XRT-localized fields, red and green are for random fields.

Cumulative distributions of LAT photon number, blue line represents 35 XRT-localized fields

Independent confirmation: “diffuse” class photon rate for 35 XRT fields is >2σ higher than rate for random background fields.

Page 7: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

LAT vs GBM fluences for XRT-localized GRBs. The two green points have been identified by our analysis.

Æ LAT fluences may not be statistically robust! ¨

Page 8: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Using GBM localizations to find LAT GRB γ-rays

Data Selection: Same as before EXCEPT:

No gbGRB cut

ΘGBM-bore < 52°

θγ < 16.0° relative to GBM localization

GBM fluence > 5.0 μerg/cm2

22 fields selected with no prior evidence for LAT GRBs

ij i j ijQ w w D

iw

ijD1/3( )i j ijijk k jk ikR w w w D D D

GRB photon cluster finder:

a) Compute pair score, ,

where traces E, t and event class correlations

and traces pair angular correlation

b) Rank photon triplet weights,

c) Select triplet with largest as best GRB candidateiw

Page 9: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Cumulative distributions of modified matched filter weights, blue line represents 22 GBM-localized fields, red and green are for random fields.

Cumulative distributions of LAT photon number, blue line represents 22 GBM-localized fields

Page 10: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

N N

solid blue circle is Earth limb, outer green circle is LAT FoV

green circle is 16° locus from GBM estimated direction

Our localization is 8.7° from GBM (von Kienlinet al.) and 0.5° from IPN (Guirec et al.)

GRB 090228A

Page 11: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

solid blue circle is Earth limb, outer green circle is LAT FoV

green circle is 16° locus from GBM estimated direction

N NGRB 081006A

GBM fluence < 5.0

Page 12: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Conclusions:

1. High energy, low fluence GRBs can be identified using matched filter techniques.

2. The number of additional detectable LAT GRBs is somewhere between 25% and 50% of the current rate.

3. The real-time identification of such events is probably possible within the constraints of the Fermi flight software. Since the Fermi mission is not likely to be emulated soon, this enhancement deserves some immediate consideration from the GRB community.

4. The number of low fluence events is surprisingly small relative to the majority of the GRBs identified by the LAT to date. This dearth is extremely difficult to understand.

Page 13: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan

Estimating LAT GRB fluences - an unusual statistical problem

Estimating the population mean from a sample mean:

;sx x2

2s

xx n

However, for LAT GRB fluence estimates:

; 2 3th th

dn E cdE E E

1 ;2 thE E

2E

Conclusion: Fluence can not be robustly computed by summing photon energies. Instead, find median energy and multiply by number of detected photons. Even for finite Emax, theseconsiderations apply as long as .1

max min( / ) 500E En g

Page 14: Looking for Needles in Haystacks Using the Swift/XRT and ... · Carl W. Akerlof, Weikang Zheng, Shashi Bhushan Pandey & Timothy A. McKay Physics Department University of Michigan