[논문 리뷰] A simple alternative to the Crystal Ball function

S. Das|arXiv (Cornell University)|2016. 03. 06.

Particle physics theoretical and experimental studies참고 문헌 4인용 수 7

한 줄 요약

이 논문은 고에너지물리학에서 지수 꼬리가 있는 피크를 피팅하기 위한 크리스탈 볼 함수의 더 간단한 대안으로 GaussExp 및 ExpGaussExp 함수를 제안한다. 거듭제곱 꼬리 대신 지수 꼬리를 사용함으로써, 새로운 함수들은 피팅의 불안정성을 감소시키고 파라미터 수를 줄이며, 특히 대규모 데이터셋이나 높은 꼬리 파라미터를 가진 경우 시뮬레이션 및 실제 데이터에서 더 뛰어난 수렴성과 안정성을 보인다.

ABSTRACT

We present a simple alternative to the Crystal Ball function that has an exponential tail stitched to a Gaussian core. It has one parameter less than the Crystal Ball function and, where appropriate, offers more stable fits to peaks that continue into exponential tails. The function may also be extended with two exponential tails on each side of the Gaussian, and this has two parameters less than the corresponding double-shouldered Crystal Ball function. This function has been used to model background and signal processes in a recent Higgs pair production search and may be of versatile use in experimental physics and other fields.

연구 동기 및 목표

• 특히 장기간 지수 꼬리를 가진 분포를 피팅할 때 크리스탈 볼 함수의 불안정성과 높은 파라미터 민감성 문제를 해결하기 위해.
• 더 적은 수의 파라미터를 가진 연속 미분 가능(C¹) 함수를 개발하여 가우시안 코어와 지수 꼬리 간의 매끄러운 스티칭을 유지하기 위해.
• 크리스탈 볼 함수가 거대한 거듭제곱 꼬리 지수 n의 극단적 값으로 인해 실패하는 고통계량 데이터셋에서 피팅 수렴성과 파라미터 신뢰도를 향상시키기 위해.
• 헤이그스 보손 탐색과 같은 실험 입자물리학에서의 관측 질량 분포를 모델링하기 위한 강력한 대안을 제공하기 위해.
• 최근 CMS 힉스 쌍 생성 탐색에서 배경(가우시안)과 신호(지수 꼬리) 모델링에 응용된 바와 같이, 새로운 함수들이 실제 응용에 유용한지 입증하기 위해.

제안 방법

• GaussExp 함수는 가우시안 코어의 한쪽 측면에 지수 꼬리를 사용하며, 전이 지점에서 함수 값과 일阶 도함수의 연속성을 강제한다.
• 전이 지점은 크리스탈 볼 함수의 두 파라미터(α 및 n)를 대체하는 단일 감쇠 상수 k로 정의되며, 파라미터 수를 한 개 줄인다.
• 이중어깨형 ExpGaussExp 함수는 가우시안 코어의 양쪽에 지수 꼬리를 적용하며, 각각 kL 및 kH로 지정된 감쇠 상수를 사용하며, 이중어깨형 크리스탈 볼 함수 대비 파라미터 수를 두 개 줄인다.
• 함수는 전이 지점에서 값과 기울기가 일치하도록 구성되어 있어 C¹ 연속성을 보장한다.
• 피팅은 ROOT/MINUIT 프레임워크 내에서 최대우도법을 사용하며, 시뮬레이션 및 실제 데이터에서 크리스탈 볼 함수와의 비교를 수행한다.
• 메서드는 지수 감쇠 분포와 시뮬레이션된 H→γγ 피크에서 검증되었으며, 다양한 사건 수에서 수렴성과 파라미터 불확실성의 정도를 평가하였다.

실험 결과

연구 질문

• RQ1 크리스탈 볼 함수의 거듭제곱 꼬리를 지수 꼬리로 대체하면서도 C¹ 연속성을 유지하는 함수를 개발할 수 있는가?
• RQ2 제안된 GaussExp 함수는 파라미터 수가 적고, 특히 고통계량 영역에서 크리스탈 볼 함수보다 더 안정적인 피팅을 달성할 수 있는가?
• RQ3 ExpGaussExp 함수는 이중어깨형 크리스탈 볼 함수보다 신호 피크의 비대칭 꼬리를 더 신뢰성 있게 모델링할 수 있는가?
• RQ4 장기간 지수 꼬리가 있는 대규모 데이터셋을 피팅할 때, 새로운 함수들의 수렴성과 파라미터 불확실성은 어떻게 되는가?
• RQ5 GaussExp 및 ExpGaussExp 함수는 힉스 보손 질량 재구성과 같은 실제 실험 데이터에 성공적으로 적용될 수 있는가?

주요 결과

• GaussExp 함수는 테스트된 모든 경우(1,000~1,000,000 건의 사건)에서 수렴에 성공했지만, 크리스탈 볼 함수는 100,000 및 1,000,000 건의 사건 데이터셋에서 지수 n의 불안정성으로 인해 수렴 실패를 보였다.
• 100,000 건의 사건에서 크리스탈 볼 함수는 n = 127 ± 40을 반환하여 파라미터 제약이 매우 낮음을 보였고, 반면 GaussExp는 진짜 값과 일치하는 파라미터와 작은 불확실성을 반환하였다.
• H→γγ 신호 피크 피팅에서 크리스탈 볼 함수는 수렴에 실패했으며, nL = 100 ± 120 및 nH = 63 ± 131을 반환하여 매우 불안정하고 물리적으로 비합리적인 결과를 보였다.
• ExpGaussExp 함수는 χ²/ndf가 1에 매우 가까운 결과를 얻었고, 작은 파라미터 불확실성을 보고했지만, 크리스탈 볼 함수는 더 나쁜 χ²/ndf를 기록했고 수렴에 실패했다.
• GaussExp 및 ExpGaussExp 함수는 대규모 샘플과 고꼬리 파라미터 영역을 포함한 모든 테스트 데이터셋에서 뛰어난 수렴성과 안정성을 보였다.
• 이러한 신규 함수들은 최근 CMS 힉스 쌍 생성 탐색에서 성공적으로 활용되었으며, 배경(가우시안)과 신호(이중지수 꼬리)의 관측 질량 분포를 각각 모델링하는 데 사용되었다.