가중 히스토그램 비교를 위한 CHICOM 프로그램

초록

본 논문은 가중 이벤트와 비가중 이벤트를 포함한 히스토그램을 비교하기 위한 Fortran‑77 기반 CHICOM 코드를 소개한다. 정규화 가중치와 비정규화 가중치 두 경우에 대해 χ² 형태의 검정 통계량을 계산하며, 하나의 비가중 히스토그램과 하나의 가중 히스토그램, 혹은 두 개의 가중 히스토그램을 동시에 비교할 수 있다. 코드 구조, 입력 형식, 사용 예시 및 검증 결과가 상세히 제시된다.

상세 분석

CHICOM은 실험 물리학 및 천문학 등에서 흔히 사용되는 히스토그램 기반 데이터 비교 문제를 해결하기 위해 설계된 독립 실행형 프로그램이다. 기존의 χ² 검정은 각 구간에 들어간 사건 수가 정수이며 동일한 가중치를 가정한다는 제한이 있었지만, 현대 시뮬레이션에서는 이벤트마다 가중치가 부여되는 경우가 빈번하다. 저자들은 이러한 상황을 두 가지 클래스로 구분한다. 첫 번째는 “정규화 가중치”로, 전체 가중치 합이 1이 되도록 스케일링된 경우이며, 두 번째는 “비정규화 가중치”로, 실제 발생 확률과 직접적인 비례 관계를 유지하지만 총합이 1이 아닐 수 있다.

두 경우 모두 기대값과 분산을 가중치의 합과 제곱합을 이용해 추정한다. 구체적으로, i번째 구간에 대한 기대값 μ_i는 전체 가중치 합 W와 구간별 가중치 합 w_i를 이용해 μ_i = N·(w_i/W) 로 정의되고, 분산 σ_i²는 (N·w_i/W)·(1 - w_i/W) + (N·s_i²/W²) 형태로 계산된다. 여기서 N은 전체 사건 수, s_i²는 구간별 가중치 제곱합이다. 이러한 식을 기반으로 χ² 통계량 χ² = Σ ( (n_i - μ_i)² / σ_i² ) 를 도출한다.

CHICOM은 입력 파일에 히스토그램 빈도와 가중치 정보를 순차적으로 기록하도록 요구한다. 파일 포맷은 Fortran‑77의 기본 입출력 방식을 따르며, 각 구간에 대해 (빈도, 가중치 합, 가중치 제곱합) 세 개의 실수를 제공한다. 프로그램은 옵션 플래그를 통해 (1) 비가중 vs 가중, (2) 두 가중 히스토그램 비교, (3) 정규화 vs 비정규화 가중치 선택을 제어한다. 출력은 χ² 값, 자유도, p‑값 및 각 구간별 기여도(잔차) 등을 포함한다.

검증 절차에서는 Monte‑Carlo 시뮬레이션을 이용해 가중치가 임의로 부여된 데이터셋을 생성하고, CHICOM이 제공하는 p‑값이 이론적 균등 분포를 따르는지를 확인하였다. 결과는 95 % 신뢰 구간 내에서 일관성을 보였으며, 특히 비정규화 가중치 경우에도 기대값과 분산 추정이 정확함을 입증했다. 또한, 실제 물리 실험(예: 입자 검출기 효율 측정) 데이터에 적용한 사례가 제시되어, 기존 χ² 검정 대비 가중치 효과를 정량적으로 반영함으로써 오차 감소와 검정력 향상을 확인하였다.

프로그램의 제한점으로는 (a) 가중치가 매우 큰 경우 수치적 불안정성이 발생할 수 있음, (b) 구간 수가 적어 자유도가 낮을 때 검정 결과가 과도하게 보수적일 수 있음, (c) Fortran‑77 기반이므로 현대 환경에서 컴파일 및 연동이 다소 번거로울 수 있다는 점을 언급한다. 향후 개선 방향으로는 고정밀 실수형 지원, 다중 스레드 병렬 처리, Python·C++ 인터페이스 제공 등이 제시된다.