Using Muon Rings for the Calibration of the Cherenkov Telescope Array: An Analytical Solution for the Dual-Mirror Telescope Using Vector Geometry

Resumen Técnico: Solución Analítica para la Calibración de Muones en Telescopios Cherenkov de Doble Espejo

Planteamiento del Problema
Los Telescopios de Imágenes de Cherenkov Atmosférico (IACT, por sus siglas en inglés) dependen de las imágenes de anillos de muones para la calibración precisa del rendimiento óptico y el monitoreo de la función de dispersión de punto (PSF). Aunque las soluciones analíticas para telescopios de un solo espejo con reflectores planos están bien establecidas (por ejemplo, Vacanti et al. 1994), una descripción analítica completa para telescopios de doble espejo ha sido esquiva. La complejidad surge de la necesidad de considerar el espejo secundario, su estructura de soporte (deflectores o baffles) y el agujero central en el espejo primario, particularmente para muones inclinados y aquellos con parámetros de impacto no nulos. Intentos previos de aplicar formalismos de un solo espejo a sistemas de doble espejo (como el telescopio ASTRI-Horn) han ignorado la sombra del espejo secundario o han utilizado factores de sombra promedio simplificados, lo que podría introducir errores significos en la calibración.

Metodología
Los autores derivan una solución analítica completa utilizando un formalismo de geometría vectorial. El enfoque consiste en:

1. Construcción Vectorial: Definir el vector de velocidad del muon y el vector de emisión de fotones Cherenkov mediante la fórmula de rotación de Rodrigues para dar cuenta de la inclinación del muon ($\nu$), el azimut ($\psi$) y el ángulo Cherenkov ($\theta_c$).
2. Propagación Geométrica: Resolver la intersección de las trayectorias de los fotones con el espejo primario (aproximado como una superficie parabólica $z \approx c \cdot \rho^2$) y, posteriormente, con el plano focal.
3. Expansión Simbólica: Utilizar el sistema de álgebra computacional SageMath para realizar la manipulación simbólica y las expansiones de Taylor de las expresiones derivadas. Las expansiones se llevan a cabo hasta primer o segundo orden en parámetros pequeños ($\nu$, $\theta_c$ y la curvatura del espejo $c$), asegurando que las desviaciones de la solución completa permanezcan por debajo del 1%.
4. Análisis de Sombreado: Derivar condiciones explícitas para el sombreado de la luz Cherenkov por parte del espejo secundario, sus deflectores y el agujero central del espejo primario. Esto implica resolver para las alturas de emisión máxima ($L_{2,max}$) y mínima ($L_{2,min}$) a lo largo de la trayectoria del muon donde los fotones son obstruidos.

Contribuciones Clave

• Primera Solución Analítica de Doble Espejo: El artículo presenta la primera descripción analítica completa de la detección de luz Cherenkov en IACT de doble espejo, considerando explícitamente el espejo secundario y la geometría de la cámara.
• Corrección por Curvatura del Espejo: Los autores derivan una corrección de primer orden para la altura de emisión máxima de los fotones Cherenkov derivada de la curvatura del espejo primario, un efecto que se ha pasado por alto en los análisis previos de espejos planos.
• Cuantificación de la Aberración de Coma: El formalismo cuantifica el impacto de la aberración de coma de tercer orden (Seidel) en los anillos de muones en telescopios parabólicos de un solo espejo, proporcionando una expresión analítica para el sesgo resultante en la reconstrucción del anillo.
• Formalismo de Sombreado de Deflectores: Se proporciona un tratamiento analítico detallado para los efectos de sombreado de los deflectores del espejo secundario, distinguiendo entre muones que atraviesan los deflectores, el espejo secundario o ninguno de los dos.
• Marco Unificado: Las expresiones derivadas forman un marco unificado para calcular la longitud total de la trayectoria del muon imbuida ($L_{max} - L_{2,max} + L_{2,min}$) bajo diversas obstrucciones geométricas.

Resultados

• Validación: El formalismo recupera con éxito las soluciones analíticas conocidas para casos de un solo espejo y espejo plano (Vacanti et al. 1994) y reproduce los comportos esperados para la aberración de coma en casos límite.
• Magnitud de las Desviaciones: Al compararse con algoritmos previos que tratan al espejo secundario simplemente como un gran agujero en el primario (ignorando la inclinación y la geometría específica de los deflectores), la nueva solución revela desviaciones de hasta un 40% para muones inclinados.
• Efectos de Sombreado: Para los Telescopios Schwarzschild-Couder (SCT) con deflectores, las fracciones de sombra pueden superar el 80%. La inclusión de la geometría de los deflectores altera significativamente la condición de sombra, particularmente para muones inclinados alejándose del espejo secundario.
• Sesgo de Coma: Si bien la aberración de coma introduce un sesgo en el radio del anillo reconstruido, el artículo señala que este efecto es generalmente pequeño (del orden del 1%) tras aplicar las correcciones estándar de escala de placa, aunque puede alcanzar varios puntos porcentuales para parámetros de impacto extremos.
• Corrección de Altura de Emisión: La corrección de primer orden a la altura máxima de emisión debido a la curvatura del espejo se encuentra en aproximadamente un 1% para los parámetros típicos del telescopio, un factor no despreciable que se había pasado por alto anteriormente.

Significancia
El artículo afirma que este trabajo permite una predicción significativamente más precisa de la distribución de luz esperada a lo largo de los anillos de muones en instrumentos de doble espejo. Al proporcionar una solución analítica rigurosa, los autores sostienen que el formalismo es directamente aplicable a la calibración basada en muones del Observatorio del Telescopio de la Red de Cherenkov (CTAO), específicamente para sus telescopios Schwarzschild-Couder. El trabajo demuestra que la complejidad de los sistemas de doble espejo puede tratarse analíticamente con suficiente precisión (con una desviación inferior al 1%) para cumplir con los estrictos requisitos de calibración de la próxima generación de observatorios de rayos gamma. Todas las fórmulas derivadas han sido implementadas en scripts de Python y se ponen a disposición pública para facilitar su aplicación inmediata en el análisis de datos de CTAO.