[논문 리뷰] Gammapy - A prototype for the CTA science tools

Christoph Deil, R. Zanin외 29명|arXiv (Cornell University)|2017. 01. 01.

Gamma-ray bursts and supernovae참고 문헌 9인용 수 11

한 줄 요약

Gammapy는 Numpy, Scipy 및 Astropy를 기반으로 한 오픈소스 파이썬 패키지로, 고에너지 감마선 데이터 분석을 위한 것으로, 체렌코프 전파망원경 배열(Cerenkov Telescope Array, CTA) 과학 도구의 프로토타입 역할을 한다. 이는 이벤트 리스트와 기구 응답 함수(instrument response functions, IRFs)로부터 천체 이미지, 스펙트럼 및 빛의 세기 변화도를 생성할 수 있으며, H.E.S.S., Fermi-LAT 및 시뮬레이션된 CTA 데이터에 대한 검증된 응용 사례를 통해 CTA의 개방적이고 협업 기반의 과학 프레임워크에 적합함을 입증한다.

ABSTRACT

Gammapy is a Python package for high-level gamma-ray data analysis built on Numpy, Scipy and Astropy. It enables us to analyze gamma-ray data and to create sky images, spectra and lightcurves, from event lists and instrument response information, and to determine the position, morphology and spectra of gamma-ray sources. So far Gammapy has mostly been used to analyze data from H.E.S.S. and Fermi-LAT, and is now being used for the simulation and analysis of observations from the Cherenkov Telescope Array (CTA). We have proposed Gammapy as a prototype for the CTA science tools. This contribution gives an overview of the Gammapy package and project and shows an analysis application example with simulated CTA data.

연구 동기 및 목표

• 향후 체렌코프 전파망원경 배열(Cherenkov Telescope Array, CTA) 관측소를 지원하는 고성능, 유연한 오픈소스 소프트웨어 프레임워크를 개발하는 것.
• Astropy 및 Numpy와 같은 기존 천문학 소프트웨어 스택을 활용하여 과학적 파이썬 기반의 공동 개발 기반의 분석 파이프라인을 구축하는 것.
• 현재 및 향후 감마선 실험 간의 데이터 포맷 및 분석 방법의 상호운용성과 표준화를 보장하는 것.
• CTA를 위한 확장 가능하고 유지보수 용이한 과학 도구 스택의 프로토타입을 제공하여 1D 및 3D 분석 방법을 모두 지원하는 것.
• Astropy와 같은 널리 채택된 도구를 기반으로 하여 접근성 있고 잘 문서화된 소프트웨어를 제공함으로써 감마선 천문학 공동체 간의 협업을 촉진하는 것.

제안 방법

• 고성능, 가독성 및 유지보수 용이성을 확보하기 위해 과학적 파이썬 스택(Numpy, Scipy, Astropy)을 기반으로 한 코드 기반을 구축한다.
• 이벤트 리스트와 기구 응답 함수(instrument response functions, IRFs)는 FITS 형식을 사용하며, 천체 이미지 및 스펙트럼 생성을 위한 에너지 및 공간의 빈(bin) 처리를 지원한다.
• 링 모델 또는 데이터 기반 템플릿을 통해 배경 추정을 수행하고, 이벤트 에너지 및 방향 정보를 이용해 노출도 지도를 계산한다.
• 2D 및 3D(cube) 분석 방법을 지원하며, 후자는 원천 방출의 동시에 공간적 및 스펙트럼적 모델링을 가능하게 한다.
• 모델링 및 피팅을 위한 Sherpa, 시각화를 위한 Matplotlib, 그리고 영역(region), 재샘플링(reproject), Healpy를 활용한 공간 및 하늘 좌표 연산을 위한 옵션 도구를 통합한다.
• 이벤트 리스트, 천체 이미지, 배경 2D, 원천 스펙트럼 모델 등의 클래스를 포함한 모듈러하고 객체 지향적인 설계를 통해 재사용 가능하고 조합 가능한 분석 워크플로우를 제공한다.

실험 결과

연구 질문

• RQ1 Astropy와 Numpy를 기반으로 한 파이썬 기반 오픈소스 패키지가 CTA 과학 도구의 타당하고 확장 가능한 프로토타입으로 기능할 수 있는가?
• RQ2 Gammapy는 실제 원천 탐지, 형태 분석, 스펙트럼 추출 및 빛의 세기 변화도 생성과 같은 고수준 감마선 분석 작업을 얼마나 잘 지원하는가?
• RQ3 Gammapy는 현재의 IACTs(H.E.S.S., MAGIC 등) 및 Fermi-LAT에서 사용 중인 기존 데이터 포맷과 파이프라인과 얼마나 잘 통합되는가?
• RQ4 시뮬레이션된 CTA 데이터에 대해 Gammapy의 3D(cube) 분석 접근법이 전통적인 1D 스펙트럼 분석 대비 얼마나 효과적인가?
• RQ5 Gammapy의 아키텍처와 개발 모델은 CTA의 30년 운영 수명 동안의 장기적 진화와 공동 기여를 지원할 수 있는가?

주요 결과

• Gammapy는 이벤트 리스트와 기구 응답 함수(instrument response functions, IRFs)를 사용하여 천체 이미지, 스펙트럼 및 빛의 세기 변화도 생성을 포함한 고수준 감마선 데이터 분석을 성공적으로 지원한다.
• 이 패키지는 H.E.S.S. 및 Fermi-LAT의 실제 데이터를 대상으로 검증되었으며, 현재 CTA 데이터의 시뮬레이션 및 분석에 적극적으로 사용되고 있다.
• Gammapy는 2D 및 3D 분석 방법을 모두 지원하며, 3D 큐브 분석은 복잡한 원천 영역에서 원천 특성 분석 향상에 유망한 가능성을 보이고 있다.
• Astropy와 Numpy의 사용 덕분에 큰 규모의 이벤트 리스트 및 천체 이미지에서 효율적인 계산이 가능하며, Numpy와 Astropy의 C 확장 기능을 통해 성능이 향상된다.
• Gammapy는 PyPI, conda 및 여러 리눅스 패키지 매니저를 통해 배포되어 과학 공동체의 광범위한 접근성과 설치 용이성을 확보하고 있다.
• 이 프로젝트는 약 30명의 감마선 천문학자로부터 기여를 유치하여 강력한 공동체 수용과 협업 기반 개발 잠재력을 입증하고 있다.