Top Banner
SPIRE-P Mapmaking, Data Processing and Future Developments Matt Griffin
40

SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Aug 23, 2020

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE-P Mapmaking, Data Processing and Future

Developments

Matt Griffin

Page 2: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Broadband Photometer Flux

Calibration

Matt Griffin

Chris North, Bernhard Schulz, Alex Amaral-Rogers, George Bendo, Jamie Bock, Alex Conley, Darren Dowell,

Marc Ferlet, Jason Glenn, Tanya Lim, Chris Pearson, Michael Pohlen, Bruce Sibthorpe, Locke Spencer,

Bruce Swinyard, IVan Valtchanov

Page 3: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Point Source Calibration

Page 4: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Measured Flux Density for a Point Source

•  Detector absorbed power ∝ flux density, weighted by aperture efficiency and SRF, integrated across the band

ν

η(ν)

S(ν)

F(ν)

S(ν) = source spectrum

F(ν) = instrument SRF

η(ν) = aperture efficiency

Page 5: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Primary Calibrator for SPIRE: Neptune Brightness Temperature Spectrum

Moreno model 2012 (“ESA-3” or “ESA-4”)

Page 6: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Neptune Spectrum and SPIRE Bands

Page 7: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Source Spectrum •  Spectrum characterised by value at some frequency

and a shape function

•  Examples:

Power law with index α :

Grey body (T, β) :

Page 8: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Monochromatic Flux Density at νo •  Measured SRF-weighted flux density is converted to

monochromatic flux density at standard frequencyνo

•  A source spectrum must be assumed

•  Power law case

•  Standard: adopt αo = –1, i.e., ν S(ν) = constant

•  Pipeline produces

Page 9: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Colour Correction •  Pipeline output = monochromatic

flux density for αs = –1

•  Colour correction factor KColP must be applied

•  Depends on “actual” source spectrum E.g., for grey body assumption:

Log(ν)

Log(Sν)

Log(ν)

Log(Sν)

Page 10: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Extended Source Calibration

Page 11: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Extended Source •  Small region at off-axis angles (θ, φ)

- Brightness - Size - Flux density - Normalised beam response

•  Measured flux density for the small region:

θ

φ

Page 12: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Extended Source •  Integrate over beam for total measured flux density

•  Same as point source case except S(ν ) replaced by sky intensity integrated over the beam

Page 13: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

•  Spectral shape: e.g. power law, grey body . . .

•  Spatial shape: e.g., Gaussian, power law, uniform (g = 1) . . .

•  Beam profile is also circularly symmetric:

•  Deriving peak surface brightness from measured flux density:

•  KMonE converts flux density (Jy/beam) to surface brightness (Jy/pixel or MJy/sr)

Some Simplifying Assumptions •  Source characterised by spectral and spatial shape functions

•  Spatial distribution is circularly symmetric

Page 14: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Case of Uniform Extended Emission: g(θ, θo) = 1 •  Integral over the beam ≡ frequency-dependent beam solid

angle:

•  This depends on how the beam solid angle varies across the band - Not measured directly - Must be modeled from band-averaged profile as measured on a point source with a known spectrum

Jy/beam → surface brightness

Page 15: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Conversion from Point Source Pipeline to Extended Emission Pipeline

•  KColE converts to peak intensity of an extended source with brightness profile g(θ, θo) and spectrum f(ν, ν o)

•  Extended pipeline surface brightness, which assumes a fully extended source with spectral index –1 :

•  Conversion to actual source peak intensity:

Page 16: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Beam Profile and Aperture Efficiency

Page 17: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Conversion from a Point Source Pipeline to Extended Emission Case

•  Treatment above is very general - Beam properties (aperture efficiency and

beam profile) are allowed to vary across the band

- Source spectral and spatial distributions can also vary

•  Details depend on detector and optical system

architecture

•  Two common cases - Absorber-coupling with 0.5λ/D square pixels - Antenna-coupling with 2λ/D circular feedhorns

Page 18: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Angle (degrees)

Absorber-Coupled Detectors

•  Array of square pixels

•  Pixel size ≡ 0.5λ/D for instantaneous full

spatial sampling

•  Wide pixel field of view ( ~ π sr)

•  Cold stop limits field of view

•  Near top-hat illumination of the telescope

•  Pixel couples to Airy disk

Cold enclosure

D

Illumination pattern

Page 19: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

2λ/D

Feedhorn-Coupled Detectors

• Single-mode feedhorns

• Tapered (~ Gaussian) illumination of the telescope

• Maximum aperture efficiency

when horn diameter ≡ 2λ/D • Throughput is λ2

Illumination pattern

Page 20: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Aperture Efficiency vs. Detector Size

0 0.5 1 1.5 2 2.5 30

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

Horn diameter or pixel side ( F*Lambda)

Aper

ture

effic

iency

G

Ape

rtur

e Ef

ficie

ncy

Size (units of λ/D)

Absorber-coupled pixel Feedhorn

Page 21: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE SRFs and Aperture Efficiencies

Page 22: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Primary Illumination Edge Taper Variation across the Band

2λ/D Feedhorn

Page 23: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Beam Profile for 2λ/D Feedhorn with λ /Δλ = 3

λ/Δλ = 3

2λ/D Feedhorn

Page 24: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Corresponding Beam Solid Angle

λ/Δλ = 3

Ω ∝ ν-1.75

2λ/D Feedhorn

Page 25: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Earlier SPIRE Method of Dealing With This Effect: λ2 Weighting of the SRF

P = Point source (no weighting) E = Extended source (λ2 weighting)

P

P P

E E E

Frequency (GHz)

Tran

smis

sion

( ) ( )2

oo ⎟

⎞⎜⎝

⎛=νν

νν RR

Page 26: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Beam Solid Angle and Aperture Efficiency for 0.5λ/D Pixel

η(νo) = 0.158

Page 27: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

η(νo) = 0.158

0.5λ/D Absorber-Coupled Pixel: Beam Solid Angle and Aperture Efficiency

Page 28: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Ω(ν) η(ν) = constant

η(νo) = 0.158

0.5λ/D Absorber-Coupled Pixel: Beam Solid Angle and Aperture Efficiency

Page 29: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Examples

Page 30: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Feedhorn Coupling, Square Passband with νo at Centre

λ/Δλ =

10

5

3

Point source colour correction factor vs. source

spectral index

Page 31: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Different Choice of Nominal Frequency

λ/Δλ = 3 Centre

- 3%

+ 3%

Page 32: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Colour Correction (Power-Law SED)

Page 33: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Colour Correction (Grey Body SED)

Page 34: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Beams

Circularised profiles

• Broadband beams measured on Neptune (α = 1.3 – 1.5)

• Averaged over all detectors

ΩMeas = (450, 795, 1665) sq. arcsec. with 4% uncertainty

Page 35: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Beams • Modelled 250-µm beam at νeff and band edges

Assumption is that measured broadband beam = monochromatic beam at some frequency νeff

Page 36: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Beams • Effective Beam Solid Angle vs. Source Spectral Index

Page 37: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Beams • Effective Beam Solid Angle vs. Grey Body Temperature

Page 38: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Calibration of Partly Extended Emission

Conversion factor from extended source pipeline to total flux density for a Gaussian source brightness profile

Page 39: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

SPIRE Calibration of Partly Extended Emission

Conversion factor from extended source pipeline to peak surface brightness for a Gaussian source brightness profile

Page 40: SPIRE-P Mapmaking, Data Processing and Future Developmentsherschel.esac.esa.int/CalibrationWorkshop5/... · George Bendo, Jamie Bock, Alex Conley, Darren Dowell, Marc Ferlet, Jason

Conclusions •  General scheme for calibration of broadband photometric

measurements of point-like, semi- or fully-extended emission

•  Can be applied to absorber-coupled or antenna-coupled systems

•  Calibration implementation for SPIRE: see - Talks by George Bendo and Bernhard Schulz - Poster by Chris North - Updated SPIRE Observers’ Manual

•  New method produces brightness values ~ (7, 7, 12)% lower than previous method for a ν3 uniformly extended source

•  Future work: - Develop practical scheme for dealing with semi-extended emission - Reduce beam solid angle ΩMeas uncertainties using Neptune shadow map to remove zodiacal light and background point sources