From f9228669718715cfbbfe3d925217b4e29bf7688d Mon Sep 17 00:00:00 2001 From: Mariana Vivas Date: Tue, 5 Nov 2024 21:16:47 -0400 Subject: [PATCH] app en espanol --- app.py | 18 +- docs/__init__.py | 0 docs/english/__init__.py | 0 .../english/analyses/03_analyses.py | 0 docs/english/analyses/__init__.py | 0 docs/spanish/__init__.py | 0 docs/spanish/analyses/03_analyses.py | 585 ++++++++++++++++++ docs/spanish/analyses/__init__.py | 0 .../__pycache__/03_analyses.cpython-311.pyc | Bin 0 -> 40576 bytes .../__pycache__/__init__.cpython-311.pyc | Bin 0 -> 201 bytes 10 files changed, 599 insertions(+), 4 deletions(-) create mode 100644 docs/__init__.py create mode 100644 docs/english/__init__.py rename 03_analyses.py => docs/english/analyses/03_analyses.py (100%) create mode 100644 docs/english/analyses/__init__.py create mode 100644 docs/spanish/__init__.py create mode 100644 docs/spanish/analyses/03_analyses.py create mode 100644 docs/spanish/analyses/__init__.py create mode 100644 docs/spanish/analyses/__pycache__/03_analyses.cpython-311.pyc create mode 100644 docs/spanish/analyses/__pycache__/__init__.cpython-311.pyc diff --git a/app.py b/app.py index b52f218..3756362 100644 --- a/app.py +++ b/app.py @@ -1,10 +1,18 @@ import streamlit as st from streamlit_option_menu import option_menu +from st_social_media_links import SocialMediaIcons + import importlib -from utils import load_sidebar_tabs import json +import importlib.util import os -from st_social_media_links import SocialMediaIcons + +def load_module_from_path(module_name, file_path): + spec = importlib.util.spec_from_file_location(module_name, file_path) + module = importlib.util.module_from_spec(spec) + spec.loader.exec_module(module) + return module + st.set_page_config( page_title="ATLAS for Teachers", @@ -113,8 +121,10 @@ def reset_language(): module.run(selected_language) elif selected_tab == tabs[3]: - module = importlib.import_module("03_analyses") - module.run(selected_language) + # Dynamically import 03_analysis.py from the language-specific folder + module_path = f"docs.{selected_language.lower()}.analyses.03_analyses" + module = importlib.import_module(module_path) + module.run(selected_language) elif selected_tab == tabs[4]: module = importlib.import_module("04_extrapython") diff --git a/docs/__init__.py b/docs/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/docs/english/__init__.py b/docs/english/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/03_analyses.py b/docs/english/analyses/03_analyses.py similarity index 100% rename from 03_analyses.py rename to docs/english/analyses/03_analyses.py diff --git a/docs/english/analyses/__init__.py b/docs/english/analyses/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/docs/spanish/__init__.py b/docs/spanish/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/docs/spanish/analyses/03_analyses.py b/docs/spanish/analyses/03_analyses.py new file mode 100644 index 0000000..ab88d70 --- /dev/null +++ b/docs/spanish/analyses/03_analyses.py @@ -0,0 +1,585 @@ +import streamlit as st +from streamlit_theme import st_theme +import numpy as np +import uproot # for reading .root files +import awkward as ak # for handling complex and nested data structures efficiently +import pandas as pd +from utils_analysis import * +from PIL import Image +import base64 +import os +import json +import random +import time + + +def run(selected_language): + # Initialize everything needed + # Initialize flags to keep track of each step + if 'data_loaded' not in st.session_state: + st.session_state.data_loaded = False + + if 'nlepton_cut_applied' not in st.session_state: + st.session_state.nlepton_cut_applied = False + + if 'leptontype_cut_applied' not in st.session_state: + st.session_state.leptontype_cut_applied = False + + if 'leptoncharge_cut_applied' not in st.session_state: + st.session_state.leptoncharge_cut_applied = False + + if 'leptonpt_cut_applied' not in st.session_state: + st.session_state.leptonpt_cut_applied = False + + if 'invariant_mass_calculated' not in st.session_state: + st.session_state.invariant_mass_calculated = False + + if 'mc_loaded' not in st.session_state: + st.session_state.mc_loaded = False + + # Initialize a special session state variable for the selectbox + # This is so that the first cut option resets when using the reset analysis button + if 'n_leptons_selection' not in st.session_state: + st.session_state['n_leptons_selection'] = "--" + + # Initilize the expanders for the quizzes + # First time they are expanded, when reseting they are not + if 'expand_all' not in st.session_state: + st.session_state['expand_all'] = True + + # Initialize session state for hint on lepton cut + if 'show_hint' not in st.session_state: + st.session_state['show_hint'] = False + + # Define a function to toggle the hint + def toggle_hint(): + st.session_state['show_hint'] = not st.session_state['show_hint'] + + if 'is_higgs' not in st.session_state: + st.session_state['is_higgs'] = False + + if 'is_z' not in st.session_state: + st.session_state['is_z'] = False + + # Get the current theme using st_theme + theme = st_theme() + + # This is the json with all the info for all analyses + # Build the path to the JSON file + json_file_path = os.path.join('analyses', 'event_counts.json') + # Open and load the JSON file + with open(json_file_path, 'r') as json_file: + analyses = json.load(json_file) + + + ################ HERE IS WERE THE APP STARTS ################ + st.title("¡Descubre el Bosón Z y el Bosón de Higgs tú mismo!") + + # Introduction + st.markdown(""" + Con esta aplicación interactiva podrás descubrir los bosones Z y de Higgs. + ¡Comprendiendo y seleccionando eventos descubrirás las partículas tú mismo! + """) + + st.markdown("## ¿Cuántos datos quieres utilizar?") + st.markdown("""Comience su análisis eligiendo la cantidad de datos con los que desea trabajar. Utilice el control deslizante a continuación para seleccionar la **luminosidad integrada**, que es una medida de la cantidad de datos que ha recopilado el detector ATLAS. + +Cuantos más datos analice, más posibilidades tendrá de detectar eventos raros como el bosón de Higgs. Pero tenga en cuenta que *más datos también pueden significar más tiempo de procesamiento*.""") + # Create a slider for luminosity + lumi = st.slider( + 'Seleccionar luminosidad (fb$^{-1}$):', + min_value=12, + max_value=36, + step=12, + value=12 + ) + + if st.button("Carga los datos"): + # Reset the steps, so that people cannot break it clicking again + # Reset info for the events + if st.session_state.nlepton_cut_applied: + # Reset flags + st.session_state.nlepton_cut_applied = False + st.session_state.leptontype_cut_applied = False + st.session_state.leptoncharge_cut_applied = False + st.session_state.invariant_mass_calculated = False + st.session_state.mc_loaded = False + st.session_state.is_higgs = False + + # Delete the widget keys from session_state + for key in ['n_leptons_selection', 'flavor_selection', 'charge_pair_selection']: + if key in st.session_state: + del st.session_state[key] + + # Reading the data + random_sleep = random.randint(1, lumi) + # Display a spinner with the loading message + with st.spinner("Cargando datos... Por favor espere."): + # Simulate a time-consuming process with a random sleep + time.sleep(random_sleep) + + st.session_state.data_loaded = True + st.toast('¡Datos cargados exitosamente!', icon='📈') + + if st.session_state.data_loaded: + st.info(f" Número inicial de eventos: {analyses[f'{lumi}']['nEvents']}") + + with st.expander("🔍 Quiz", expanded=st.session_state['expand_all']): + st.markdown("##### ⁉️ Entendiendo la luminosidad") + st.markdown(f"Seleccionaste una luminosidad de **{lumi} fb⁻¹**. Pero, ¿qué representa realmente la luminosidad integrada en un experimento de física de partículas?") + + possible = ['La cantidad de tiempo que el detector está activo', + 'La cantidad de colisiones por segundo', + 'Una medida de la cantidad total de colisiones durante un período de tiempo', + 'La energía a la que ocurren las colisiones'] + answer = st.radio("Elige tu respuesta:", possible, index=None) + + if answer == possible[2]: + st.success("¡Correcto! La luminosidad representa el número total de colisiones en el conjunto de datos. Cuanto mayor sea la luminosidad, mayores serán las posibilidades de observar eventos raros como el bosón de Higgs") + elif answer: + st.error("Incorrecto. ¡Inténtelo de nuevo o lea más sobre luminosidad!") + + # Using a selectbox to let users choose between amounts of leptons + st.markdown("## Número de leptones en el estado final") + st.markdown("En los colisionadores de partículas, cuando se produce una partícula, puede decaer inmediatamente en otras partículas, que se detectan y analizan. Al identificar todas las partículas en el estado final, podemos inferir qué partículas se crearon inicialmente durante la colisión. Un ejemplo es el número de leptones en el estado final, ya que diferentes procesos producen diferentes cantidades de leptones.") + st.markdown("Below is a Feynman diagram showing a typical process that results in a final state with two leptons:") + # Diagram for Z decay + image_zdecay = f"images/Z_decay_{theme['base']}.png" + # Encode the image in base64 + with open(image_zdecay, "rb") as image_file: + encoded_string = base64.b64encode(image_file.read()).decode() + # Display the image centered and resizable using HTML with CSS + st.markdown( + f""" +
+ Centered Image +
El bosón Z decae a dos leptones.
+
+ """, + unsafe_allow_html=True + ) + + st.markdown("Decaimientos más complejas pueden implicar más leptones en el estado final") + # Diagram for H decay + image_zdecay = f"images/higgs4l_decay_{theme['base']}.png" + # Encode the image in base64 + with open(image_zdecay, "rb") as image_file: + encoded_string = base64.b64encode(image_file.read()).decode() + # Display the image centered and resizable using HTML with CSS + st.markdown( + f""" +
+ Centered Image +
El bosón de Higgs decae a bosones Z y, posteriormente, a leptones.
+
+ """, + unsafe_allow_html=True + ) + + st.markdown("A continuación se puede ver el recuento de leptones en todo el conjunto de datos. Notamos que, en general, es más común tener menos leptones en un evento") + # Get the appropriate plot file based on the theme + st.image(f"images/lepton_plot_{theme['base']}_{lumi}.png") + + st.markdown("Estudia los diagramas y los datos, y selecciona cuántos leptones esperas observar en tu estado final dependiendo del análisis que estés haciendo: encontrar el bosón Z o el bosón de Higgs.") + + st.warning(""" + Estás a punto de seleccionar la cantidad de leptones que deseas en tus datos. Sin embargo, es importante tener en cuenta que hemos agregado criterios adicionales para garantizar la calidad de estos leptones: + + - **Aislamiento:** cada leptón debe estar aislado, lo que significa que no está agrupado con otras partículas. Esto garantiza que nos estamos centrando en los leptones que probablemente se originaron directamente de la partícula que nos interesa, en lugar de interacciones de fondo. + - **Niveles de identificación:** los leptones deben cumplir con criterios de identificación específicos para confirmar su tipo con alta confianza. La reconstrucción de partículas es compleja, por lo que tenemos diferentes niveles de identificación para mayor certeza. Por ejemplo, los muones deben pasar un nivel de identificación medio, mientras que los electrones solo necesitan pasar un nivel de identificación flexible, ya que los electrones son más fáciles de detectar. + - **Condiciones de trigger:** los triggers son criterios establecidos para capturar eventos con ciertas características, lo que permite guardar solo los eventos que queremos analizar. Aquí, utilizamos triggers de electrones y muones para seleccionar eventos con señales significativas, refinando el conjunto de datos para aumentar las posibilidades de observar partículas que decaen a leptones o muones. + + Estos criterios ayudan a "limpiar" los datos, mejorando nuestras posibilidades de observar partículas como los bosones Z y Higgs. + """) + + # Define the options + n_leptons_options = ("--", 2, 3, 4) + + # Create the selectbox + n_leptons = st.selectbox( + '¿Cuántos leptones esperas en el estado final?', + n_leptons_options, + index=0, # Default index for "--" + key='n_leptons_selection' + ) + + # Access the selected value from session_state + n_leptons = st.session_state['n_leptons_selection'] + + if n_leptons == 2: + st.success("""Has elegido un estado final con **2 leptones**. + Esto sugiere que te interesa un proceso en el que una única partícula intermedia deca a un par de leptones. + Los pares de leptones son comunes en muchas interacciones de partículas, especialmente cuando se consideran intermediarios neutros""") + elif n_leptons == 4: + st.success("""Estás viendo **4 leptones** en el estado final. + Esto a menudo indica una cadena de decaimientos, donde múltiples partículas intermedias decae a pares de leptones. + Estos escenarios son interesantes para estudiar interacciones complejas""") + elif n_leptons != '--': + st.warning("""Tener un número impar de leptones es inusual en procesos de desintegración simples, ya que los leptones suelen producirse en pares debido a las leyes de conservación. + Sin embargo, esto podría sugerir que estás explorando modos de desintegración más exóticos. ¿Estás buscando materia oscura?""") + + # Number of leptons button + # We define a variable to avoid the page breaking when clicked more than once + if st.button("Aplicar selección por número de leptones"): + if st.session_state.nlepton_cut_applied: + st.toast("Ya aplicaste una selección. Para reiniciar el análisis ve al final de la página.", icon='❌') + elif n_leptons != '--': + random_sleep = random.randint(1, round(lumi/3)) + # Display a spinner with the loading message + with st.spinner("Seleccionando eventos... Por favor espere."): + # Simulate a time-consuming process with a random sleep + time.sleep(random_sleep) + st.session_state.nlepton_cut_applied = True + st.toast("Selección aplicada exitosamente.", icon='✂️') + else: + st.error("Seleccione un número válido de leptones.") + + # Step 2: Dynamically generate selection for lepton flavors + if st.session_state.nlepton_cut_applied: + st.info(f"Eventos luego de la selección: {analyses[f'{lumi}'][f'{n_leptons}leptons']['nEvents']}") + + with st.expander("🔍 Quiz", expanded=st.session_state['expand_all']): + st.markdown("##### ⁉️ Estado final de Leptones") + st.markdown(f""" + ¿Qué puede indicar el número de leptones? + """) + + possible_lepton = ['Dos leptones indican una desintegración simple de una partícula neutra, como un bosón Z', + 'Cuatro leptones indican una cadena de desintegración más compleja, que posiblemente involucre a un bosón de Higgs', + 'Un número impar de leptones sugiere un proceso exótico', + 'Todas las anteriores'] + answer_lepton = st.radio("Elige tu respuesta:", possible_lepton, index=None, key="lepton_selection_quiz") + + if answer_lepton == possible_lepton[3]: + st.success("¡Correcto! Cada opción representa una posibilidad en función de tu selección.") + elif answer_lepton: + st.error("Incorrecto. Intenta pensar en cómo se producen los leptones en pares.") + + st.markdown("## Garanticemos la conservación") + st.markdown("En las interacciones entre partículas, ciertas propiedades siempre se conservan, como la *carga* y el *sabor* leptónico. Comprender estas leyes de conservación ayuda a reducir las posibilidades de las partículas que están involucradas en el estado final.") + st.markdown("En tu análisis, puedes ver el 'sabor' de los leptones (es decir, si son electrones o muones) y su carga (positiva o negativa). Las partículas con carga opuesta se llaman antipartículas, un electron positivo se llama positron y un muón positivo, anti-muón. El plot siguiente muestra la distribución del sabor de los leptones, con una barra para leptones cargados positivamente y otra para leptones cargados negativamente. Esto ayuda a identificar si el estado final obedece reglas de conservación.") + + # Display the pre-generated plot based on the theme + st.image(f"images/lepton_barplot_{theme['base']}_{lumi}.png") + + st.markdown("Con esto en mente, realicemos la siguiente selección. Si no estás seguro, vuelve a ver los diagramas de Feynmann que están arriba. Puede que en ellos encuentres información que te pueda ayudar.") + + flavor_options = ["--", 'Mismo', 'Diferente'] + flavor = st.selectbox(f'¿Deberían los pares de leptones tener el mismo o diferente sabor?', flavor_options, key=f"flavor_selection") + + if flavor == 'Mismo': + flavor = 'Same' #! Como arreglar esto para los diferentes idiomas + st.success(""" + Seleccionar dos leptones del mismo sabor significa que estás considerando escenarios donde las propiedades de los leptones son idénticas, + como dos electrones o dos muones. Pares del mismo sabor se presentan en escenarios donde los decaimientos producto de la interacción entre + partículas respeta la conservación de sabor. + """) + + elif flavor!= '--': + flavor = 'Different' + st.warning(""" + Escoger leptones de diferentes sabores indica que estás examinando una situación donde los leptones no son idénticos, como un electrón y un muón. + Mientras esto puede ocurrir en algunos procesos, es menos común en decaimientos simples debido a las leyes de conservación de sabor. + """) + + # Apply lepton type cut based on flavor selection + if st.button("Aplica selección de tipo de leptones"): + if st.session_state.leptontype_cut_applied: + st.toast("Ya aplicaste una selección. Para reiniciar el análisis ve al final de la página.", icon='❌') + elif flavor != '--': + # Display a spinner with the loading message + random_sleep = random.randint(1, round(lumi/3)) + with st.spinner("Seleccionando eventos... Por favor espere."): + # Simulate a time-consuming process with a random sleep + time.sleep(random_sleep) + + # Display the cut result + st.session_state.leptontype_cut_applied = True + st.toast("Selección aplicada exitosamente.", icon='✂️') + + else: + st.error("Selecciona una opción para el sabor de los leptones") + + # Step 3: Dynamically generate selection for lepton charges + if st.session_state.leptontype_cut_applied: + st.info(f"Eventos luego de la selección: {analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}']['nEvents']}") + + # Offer options for charge pairing: Same charge or Opposite charge + charge_pair_options = ["--",'Igual', 'Opuesta'] + charge = st.selectbox('¿Deberían los pares de leptones tener la carga igual o contraria?', charge_pair_options) + + # Define the condition for the charge mask based on the selection + if charge == 'Igual': + charge = 'Same' + st.warning("""Pares de leptones con la misma carga son inusuales porque usualmente la carga se conserva en las interacciones. + Sin embargo, algunos procesos más exótico, o identificación erronea puede resultar en pares con igual cargaa. + """) + elif charge != '--': + charge = 'Opposite' + st.success("""Haz escogido leptones con carga opuesta. La mayoría de las interacciones conservan la carga, + por lo que es usual ver producción conjunta de un lepton y su antiparticulaso, resultando en cargas opuestas. + """) + + # Apply lepton type cut based on flavor selection + if st.button("Apply lepton charge selection"): + if st.session_state.leptoncharge_cut_applied: + st.toast("Ya aplicaste una selección. Para reiniciar el análisis ve al final de la página.", icon='❌') + elif charge != '--': + # Display a spinner with the loading message + random_sleep = random.randint(1, round(lumi/3)) + with st.spinner("Seleccionando eventos... Por favor espere."): + # Simulate a time-consuming process with a random sleep + time.sleep(random_sleep) + st.session_state.leptoncharge_cut_applied = True + + # Provide feedback to the user + st.toast("Selección aplicada exitosamente.", icon='✂️') + else: + st.error("Selecciona una opción para la carga de los leptones.") + + if st.session_state.leptoncharge_cut_applied: + st.info(f"Eventos luego de la selección: {analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['nEvents']}") + + + with st.expander("🔍 Quiz", expanded=st.session_state['expand_all']): + st.markdown("##### ⁉️ Entendiendo la Selección de Carga") + st.markdown(f""" + Seleccionaste pares de leptones con **carga {charge}**. ¿Por qué es la carga de los leptones importante en física de partículas? + """) + possible_charge = [ + 'La conservación de la carga es un principio fundamental y los productos de decaimientos deberían conservar la carga', + 'Leptones con la misma carga son más comunes en decaimientos de partículas como el decaimiento del bosón Z', + 'La carga no juega un rol importante en las interacciones de particulas', + 'Pares de leptones de la misma carga indican colisiones altamente energéticas' + ] + + answer_charge = st.radio("Elige tu respuesta:", + possible_charge, + index=None, key="charge_selection_quiz") + + # Checking for the correct answer and giving feedback + if answer_charge == possible_charge[0]: + st.success("¡Correcto! La conservación de la carga es un fundamental en la física de partículas, y leptones de cargas opuestas son típicos en decaimientos como los del bosón Z.") + elif answer_charge: + st.error("Incorrecto. Considera cómo la conservación de la carga funciona en los decaimientos. Leptone de carga opuesta son esperados en muchos decaimientos.") + + + if n_leptons==2 and flavor=='Same' and charge=='Opposite': + st.session_state.is_z = True + elif n_leptons==4 and flavor=='Same' and charge=='Opposite': + st.session_state.is_higgs = True + + # Step 4: Cuts on leptons pT only for Higgs + if st.session_state.leptoncharge_cut_applied and st.session_state.is_higgs: + st.markdown("## Cortes en el p$_T$ de los Leptones ") + st.markdown("""En nuestra búsqueda por el **bosón de Higgs**, nos valemos de **cortes** que nos ayudan a encontrarlo en una vasta cantidad de datos. Una forma en la que podemos hacer esto es enfocandonos en una variable llamada **momento transverso** (p$_T$), que representa el momento de las partículas perpendiculas al eje del haz. + +Cada lepton en nuestros datos tiene un valor de p$_T$, que hemos guardado en orden descendiente como **principal** (mayor p$_T$), **secundario**, y así sucesivamente. Al examinar estos valores, podemos aplicar cortes para sólo quedarnos con los datos con las características más parecidas a aquellas del bosón de Higgs, mientras eliminamos datos que es menos probable que lo contengan. + +Una forma de aislar al Higgs es aplicar cortes que remuevan regiones donde el **fondo** (datos no relacionados al Higgs, pero con características similares) tienden a dominar, dejando más eventos que concuerda con una huella como la del Higgs. Esto significa establecer limites inferiores en el p$_T$ de los leptones, ya que mayores valores de p$_T$ values más probablemente capturen los eventos en los que estamos interesados. + +Los gráficos siguientes muestran la distribución de p$_T$ del primer, segundo y tercer leptón en señal simulada (el Higgs que estamos buscando) y simulación del fondo. Aplicando cortes a los valores bajos de p$_T$ puede ayudar a reducir el fondo y mejorar la visibilidad de potenciales eventos del Higgs. +""") + + # Display initial image + if not st.session_state['show_hint']: + st.image(f"images/lepton_pt_{theme['base']}.png", caption="Distribución de pT de los tres leptones más energéticos en cada evento.") + else: + st.image(f"images/lepton_pt_{theme['base']}_lines.png", caption="istribución de pT de los tres leptones más energéticos en cada evento con posibles cortes.") + + st.markdown("Con esto en mente, consideremos los mejores cortes inferiores del p$_T$ que ayudarían a filtrar eventos de fondo mientras retenemos eventos que posiblmente son candidatos del Higgs.") + with st.expander("🔍 Quiz", expanded=True): + st.markdown("##### ⁉️ Selecciona Cortes Apropiados para el p$_T$ de los leptones.") + + st.markdown(""" + En base a lo discutido, ¿cuáles de los siguientes valores para cortes nos ayudarían a enfocarnos en eventos como los del Higgs, reduciendo el fondo? + """) + + cut_options = [ + "20, 15, 10 ", + "25, 20, 15", + "10, 5, 2", + "Ningún corte es necesario" + ] + + answer_cut = st.radio("Selecciona la mejor opción para los cortes en p$_T$:", cut_options, index=None, key="pt_cut_quiz") + + st.button("¿Quieres una pista? Haz click aquí y revisa la imágen de arriba" if not st.session_state['show_hint'] else "Esconder pista", on_click=toggle_hint) + # Checking for correct answer and providing feedback + if answer_cut == cut_options[0]: + st.success("¡Correcto! Establecer cortes en 20, 15, y 10 en los valores de p$_T$ de los leptones puede ayudar a enfocarnos en eventos similares a los del Higgs, filtrando eventos de fondo.") + elif answer_cut: + st.error("No realmente. Enfocate en los valores que efectivamente eliminarían fonda, todavía manteniendo la mayoría de la señal.") + + + st.markdown("Ahora, hagamos los cortes:") + + if st.button("Corte en el p$_T$ de los leptones"): + if st.session_state.leptonpt_cut_applied: + st.toast("Ya aplicaste una selección. Para reiniciar el análisis ve al final de la página.", icon='❌') + else: + st.session_state.leptonpt_cut_applied = True + # Display a spinner with the loading message + random_sleep = random.randint(1, round(lumi/3)) + with st.spinner("Seleccionando eventos... Por favor espere."): + # Simulate a time-consuming process with a random sleep + time.sleep(random_sleep) + st.session_state.leptoncharge_cut_applied = True + + # Provide feedback to the user + st.toast("Selección aplicada exitosamente.", icon='✂️') + + if st.session_state.leptonpt_cut_applied: + st.info(f"Eventos luego de la selección: {analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['ptLeptons']['nEvents']}") + + # Steep 5: invariant mass plot + if (st.session_state.leptoncharge_cut_applied and not st.session_state.is_higgs) or (st.session_state.is_higgs and st.session_state.leptonpt_cut_applied): + st.markdown("## Descubriendo Particulas con la Masa Invariante") + st.markdown("La *masa invariante* es una herramienta importante en la física de partículas. Nos permite reconstruir la masa de las partículas que se producen en las colisiones, incluso cuando no las observamos directamente. Al analizar la energía y momento de los leptones en el estado final, podemos calcular su *masa invariante* combinada.") + st.markdown("Al graficarla, la distribución de la masa invariante usualmente muestra picos donde las partículas como el bosón Z o el bosón de Higgs aparecen. Estos picos revelan la masa característica de la partícula, permitiéndonos \"verla\", incluso cuando ya decayó para el tiempo en el que estamos analizando los datos.") + st.markdown("Al calcular y graficar la masa invariante, podrás observar estos picos y, potencialmente, ¡descubrir partículas por ti mismo!") + + with st.expander("🔍 Quiz", expanded=st.session_state['expand_all']): + st.markdown("##### ⁉️ Masa Invariante") + st.markdown(""" + La masa invariante es una cantidad importante en la física de partículas. ¿Qué nos puede decir sobre una partícula? + """) + possible_mass = ["La energía de la partícula", + "El momento de la partícula", + "La masa en reposo de la partícula que produjo los leptones", + "El tipo de partícula que decayó"] + + answer_mass = st.radio("Elige tu respuesta:", + possible_mass, + index=None, key="invariant_mass_selection_quiz") + + if answer_mass == possible_mass[2]: + st.success("¡Correcto! La masa invariante nos dice la masa de la partícula que decayó en leptones.") + elif answer_mass: + st.error("Incorrecto. Recuerda, la masa invariante está relacionada con la masa en reposo de la partícula.") + + if st.button("Get invariant mass"): + st.session_state.invariant_mass_calculated = True + # Step 6: Discussion + if st.session_state.invariant_mass_calculated: + if st.session_state.is_higgs: + st.image('analyses/'+analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['ptLeptons'][f"plot_data_only_{theme['base']}"]) + else: + st.image('analyses/'+analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}'][f"plot_data_only_{theme['base']}"]) + + if not st.session_state.is_higgs: + with st.expander("🔍 Quiz", expanded=st.session_state['expand_all']): + st.markdown("##### ⁉️ Interpretando el Gráfico de la Masa Invariante") + st.markdown(""" + Generaste un gráfico de la masa invariante. Si un pico aparece alrededor de 91 GeV, a qué partícula puede que corresponda? + """) + + possible_final = ["Bosón de Higgs", "Bosón Z", "Fotón", "Quark top"] + answer_final = st.radio("Elige tu respuesta:", + possible_final, + index=None, key="invariant_mass_quiz") + + if answer_final == possible_final[1]: + st.success("¡Correcto! Un pico alrededor de 91 GeV tipicamente corresponde al bosón Z.") + elif answer_final: + st.error("Incorrecto. Un pico alrededor de 91 GeV usualmente indican la presencia de un bosón Z, ya que 91 GeV es su masa.") + + if answer_final == possible_final[1]: + if st.session_state.is_z: + st.balloons() + st.markdown("### Discusión") + st.markdown("Llegaste al final del análisis. Una vez que estés feliz con el resultado espera por la discusión, o resetea el análisis para intentar uno nuevo.") + + else: + st.markdown("## ¿Cómo sabemos que hallamos el Higgs?") + + st.markdown(""" + Para determinar si hemos observado al bosón de Higgs, vamos a comparar nuestros datos del detectos con simulaciones de procesos de fondo conocidos y simulaciones de la señal. Esta aproximación nos permite ver si en nuestros datos se presenta un pico donde esperamos que aparezca el bosóns Higgs. + """) + + # Step 1: Show data only + st.markdown("### Paso 1: Observa los Datos") + st.markdown("Empecemos observando únicamente los datos. Mira cuidadosamente: ¿notas alguna característica especial? Sin información adicional, puede ser difícil decir si alguno de los picos son debido a procesos de fondo o de señal.") + higgs_data_only = 'analyses/'+analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['ptLeptons'][f"plot_data_only_{theme['base']}"] + st.image(higgs_data_only, caption="Únicamente datos del detector") + + # Quiz question for the data-only plot + quiz_data = st.radio( + "¿Qué observas en el gráfico con únicamente datos del detector?", + options=["Un pico evidente", "Algunas fluctuaciones, pero es díficil decir", "Ninguna característica en particular"], + index=None + ) + if quiz_data == "Algunas fluctuaciones, pero es díficil decir": + st.success("¡Exacto! Sin contexto, es difícil decir cuál fluctuacion es señal.") + elif quiz_data == "Un pico evidente": + st.warning("¿Ya encontraste el Higgs? Continuemos para ver si ese es el caso.") + elif quiz_data == "Ninguna característica en particular": + st.info("¡Vamos a aclarar con algunas simulaciones!") + + if quiz_data: + # Step 2: Show data with background simulation + st.markdown("### Paso 2: Añadimos Simulaciones de Fondo") + st.markdown("Ahora, hemos añadido simulaciones de procesos de fondo (sin ninguna señal). Esto nos muestra lo que esperamos ver de otras partículas e interacciones en la ausencia del bosón de Higgs. Observa detalladamente: ¿hay algún pico que aparezca en adición al fondo?") + higgs_data_bkg = 'analyses/'+analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['ptLeptons'][f"plot_data_backgrounds_{theme['base']}"] + st.image(higgs_data_bkg, caption="Datos del detector y simulación del fondo") + + # Quiz question for data with background + quiz_background = st.radio( + "Con las simulaciones del fondo, ¿notas un pico adicional?", + options=["Sí, parece que hay un pico de más", "No, todo parece fondo", "No estoy seguro"], + index=None + ) + if quiz_background == "Sí, parece que hay un pico de más": + st.success("Buena observación! Puede que estemos viendo algo más además del fondo.") + + if quiz_background: + # Step 3: Show data with background and simulated Higgs signal + st.markdown("### Paso 3: Agregamos la Simulación de la Señal del Higgs") + st.markdown(""" + Finalmente, hemos agregado la simulación de la señal del bosón de Higgs para ver qué tan bien acuerda con los datos. Si los datos se alinean con el fondo más la simulación de la señal del Higgs, tenemos una fuerte evidencia de la existencia del bosón de Higgs. ¿Puedes ver un pico evidente donde esperamos el Higgs? + """) + higgs_data_bkg_sig = 'analyses/'+analyses[f'{lumi}'][f'{n_leptons}leptons'][f'flavor{flavor}'][f'charge{charge}']['ptLeptons'][f"plot_data_backgrounds_signal_{theme['base']}"] + st.image(higgs_data_bkg_sig, caption="Datos con fondo y simulación de la señal del Higgs") + + # Final quiz question + quiz_signal = st.radio( + "Al agregar la simulación de la señal del Higgs, ¿qué podemos concluir?", + options=["Hay un pico que concuerda con la señal del Higgs", "Todavía no es claro", "El fondo concuerda mejor con los datos"], + index=None + ) + if quiz_signal == "Hay un pico que concuerda con la señal del Higgs": + st.success("¡Eso es correcto! La alineación del pico con la simulación de la señal del Higgs provee evidencia de que estamos observando el bosón de Higgs.") + + if quiz_signal: + if quiz_signal == "Hay un pico que concuerda con la señal del Higgs": + st.balloons() + st.markdown("---") + st.markdown("### Discusión") + st.markdown("Llegaste al final del análisis. Una vez que estés feliz con el resultado espera por la discusión, o resetea el análisis para intentar uno nuevo.") + + + # Reset button to start the analysis again + st.markdown('---') + st.write("""Si quieres volver a seleccionar eventos, haz click en el botón `Reset analysis`. + No te preocupes, ¡no necesitas volver a cargar los datos! Vas a empezar nuevamente en la selección del número de leptones.""") + if st.button("Reset analysis"): + # Reset flags + st.session_state.nlepton_cut_applied = False + st.session_state.leptontype_cut_applied = False + st.session_state.leptoncharge_cut_applied = False + st.session_state.leptonpt_cut_applied = False + st.session_state.invariant_mass_calculated = False + st.session_state.mc_loaded = False + st.session_state.expand_all = False + st.session_state.is_higgs = False + st.session_state.is_z = False + + # Delete the widget keys from session_state + for key in ['n_leptons_selection', 'flavor_selection', 'charge_pair_selection']: + if key in st.session_state: + del st.session_state[key] + + st.rerun() + st.toast("El análisis ha sido reseteado.") \ No newline at end of file diff --git a/docs/spanish/analyses/__init__.py b/docs/spanish/analyses/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/docs/spanish/analyses/__pycache__/03_analyses.cpython-311.pyc b/docs/spanish/analyses/__pycache__/03_analyses.cpython-311.pyc new file mode 100644 index 0000000000000000000000000000000000000000..7f72173abda35fecc1a26d90fd615422097c77c4 GIT binary patch literal 40576 zcmd^odvILWdEf5hwRkMBARYurfLuH%5fHF^f-jK-K@xmFMT(M0k6)r%_LMqxwb-EQbX5q$~a9ei;a^%GLyT^b-J??P8oO9jAxoDSd*>e zX*&IV-?`5PmJit;yA64H_U=8e?|kp?JLfKcFA%8m;LrbA);spApZ9qF3w<&F_>7w$ zPT}T@9@!&%g1Za>*Q+g8JN44dlummIlw(@KjWG9%5|The;xM^euADs>X~jZG3N2Si`)0C zuhYJCqxHpew%WYER%6~dKjnIP-n-!V_bjfTE=e!6TTGS64et5@&va?J-EMKV&b*&0 zeHRk;o+V||q#N+p-QVj;RCjq+cv7Vod>2bz@?7wGJui8@_?xzRdOT-?=BRR`d(?t< z^;Wy4$86X9nC<32wB0;w-WIu4UVwkC?=O5868N4ah12D(6Aumyf*E3QDDB+a%77 zaY4AioL3IX9q)I#R|6Q$YJA$(?5ribmrIl%xyG|3S?-!Dn=1bukF_w5+-=|gj^}M( z;W}O8w3!pjOXh|_uC_;e|0y@zj9sT}+cNv~Lua+yEhX6Hp0`U1+mCtV-YyIS$Xau_aMncgRF%@`YQ5NMy{xeH zPoM>D3@k4k)it}_ZpCA^TUppnUNz;DSKBxNz6ahenXYzOgKBwAA#JB?(nIFj!Id56 z{Z!4hPV?UR>9$sGD6VxmTblRo9@BxJ3O?6P1+c<(Zn>3a;&%)#&%lyFGWe**}=KL)WZMVdlX{fNBymh+P1=Ukl+ig>| z@^-oM{o!|EL*KLP=yVZ}Pp-7(XS!~xPVRq%X4z2)OJV@Mc0Qulu35bnpjXv7)^G)S z$zJ#F!d|X5+uq3@cRS@T=hpZR-(!!1}pjmOOpXq=CmXRW9$( zz>jsq<(Ao^mi~ad*U!PyTXv@imJZ5?=D<=h>vGt=o}xa1HS1-c@@uY{ZJ*$|LZ4u& z4m$D3+di9%e<8etqqEjLw>R;8VLNyS$KV|VUF&yx2b6huVU~2)Pt{Ka<^47m9wCn} z7S0R*puRv3ZU5kYI@@cvWKpBA26__nRglj+Zxn>t1Hr!_0|=ZCJDpG<$Aoob-y%GM>1 z(x~A6bKrir{hW8hn^NZvg?L3PMgOA|lC$kR;H|N-PM{!f9DCb+$L0)Q`3(lXE z^HWnzPbLel*?(3!i1aMH%vDnm^SnZ?`qW zBfN~sLO9`N%rEL?WC~k5qLD?gc(O5D3v@uzxSv-vTeZ}b*hH5jK52K1Gv6^zebI63 z4B(SxhRx~m>BxK8oa$4gB2kJpBlP$|7lmw}=(R#z6=^p$tKFP%$pM?DsSejzw%;^8 z%Wp!qMSlJLPr0=_vgM)29^nm)@afqjxcb|x{-s%~aq_)G+K?f6{EXf|d+4b4?4K!Y z2fO#rAw?}NXl*I_<*61((Z84_MXj#(mMxl_8(fs}$(Isd`B$dA?|-%s=Tohi?eQh^ zR}0qy$!skunQs=helq>;TCa`I&lQg9YG>g^{xx?yAwV%*Nzg<-Jtr35r zWB#9`<$fe3hFvuM>09GpxAoSeywoDFS{hlYFVdpX;;ZgQwSfPzrAQl#m(V8c!ip^s zm))FRG^_g?kJ;|?kJ;`wAKK2LSLa^aC)`PS>!GvSSouO>JLuhi3B9|(g+*KMJ~Ood z`ss_a^zOo`g{XP_J(Mej^Fl{1EUF`atFZM?pd)SF_uKAKaR#KAv|GEKbJOUoRe!lu z`Sl`g8ZAzD`yI9q@H>x@f14dY&i106a(lJVuW-hm8Gz)*0sMY zFUZ6$o?0w_@CeTRH-(tMzSm;;_vh@j#N-(H#;jfo(Ch0)FL&O<_KyB-VXvvh@EZQW zg#(+9zERk^NV`8MZ1*#v4agIVUYUraI(5niIcvcJNz; zn3(DSeoKqs?T-sv7istJ3)}rn@cSo^#BZkyK3i5g<^S;5`29BU`{A7U?V9S6zw-$E zepHBwsV?PtN57v;zuqbw-L+m@E`%D}@)Xv+P z@9(PfBuB-jQ!e#G@1GWMU+urbncSiFhn3=bOc|^yEa; z^}IzqkDQT|F?UU)sF_H++PB6M=Ora3CDJ3~QYtkPk>U|GE~PO=<#@z^0U(Ts-l>sP zG#-tlW68vief|CUs_4ZDDIt$Y@%YhvW$H}w!pNCeA{{ZYL7&8c#JE4)f#Qor4)w`i z+YPKK+GZn0AVl@}4k#q6J-)mzmCgrJ(Nqc`jil0m{HvZHqRR8*k!8=PqDpGnIImAy zwmlh9$D@gKYFRoeMb4nA|3XqZ2cReVrF2|M^(QByiGITUGC3JZElXi(vDBGmQbLMP zVv1$USB;qW;fcw7bvijV7LW3ZGcAgm2ycuhfxy(}u%l)zU`d5gyy3ai7uNeW-aMiA z?bep={+v(uw_I76^)J!*Qs`u7-XBY0z_GDZz7%(vyx*KLQy)DK=tm++6>A8eO(hd~ zucEF&r;py*9!*8mQ6(CR#zR|^saxMjgpP+M9pB_=XcvtdO5gfcXgro0Pj;wZL)WSh z{%lJoLQ&u@B!RgiQY1$GLB^;efz;=v&_q&JZoQKV$zqzA5{f5Nq0uB@h^BZZbNnz5 z_|ZGtlH(IfG$F_EF9}GPZ6cr%xFs3lm0$#GJQ^Gn2u(;z`ql@;flf$yOqCp?uCCBW z-`=L)dMA;l?h=)+sBECqzb#so-l zDHH+M$7Bg-Pou_A5g>wA_Ns~}(Nt6lr4?yZIx8vR(DO7r#sh}YEUhHt)P!aPqCyxA z4`T_{z#bMNITj6sNul1}xH=w7BvU*o;de}tWU05ek0+!hB&5)IR7O3)G$)7z+8HUN zL?g+GSX?4x5+G>^Wl{+Z9o;u{BwSUsO$Ag~$T(At(L{U%`1{1O2|zK51!4J^Hw`L* zz%;2CN>VHuDUyPLBFS-{2{b3>u*e-c6jhSce@qHR6QPI-7L|B8y^gul()eiX)>jju z2{kICSt>S`h@A$jDFT{wEIK}sq){i7WF(rBh<52@Z#dI^g!BzF4shr|-_p~g-KSpb zUoqX?v*8DH*8d0Otn`#ZZN(BJ=cTwB#jSMFxGhbkqZ9cuex9k@2JI*@rJ03#N|ai3 z&IeD&rSnO})^YhJLy=BMF~{Q{kmk>H@hHU3JV-bk4jlpsPSX$&rwK4+SS7y^;>i@i za>qb0aN)%mh$E3Qigu`_XfpG=zj+C-_>{)LMJh!dp#!(RHBP9FC1R17L?kvTvLTeO zNDPxMNU4?F_pQ%`4y&;Yy)-2Baa2~HMoky}34Q<5zw+O2{(5K_lbP}e&T`-~(`6$W zXkspORMXVGy=sb}Iz9URPyg{pf7IKHwF0_6NaL&Ot*?R|A?+!EM4Tzbd1b`E&S5Nu zCq_}BaV`QyF>1<`KYi10|U-_cuE6_$Cd0V}xmL6TS?8WJ$6N=0`_ zc}EDtI}z(ZK9G?x7-ufXF}qM-XJJDaXu-N6p)W$?(j;SsVRlRbr9m~&6tJjJzPtvv z3IPI=qo{>9VKQpyp{he2#Hw~DBF3u1p^x6#ok&ATq@!^LEuluwCqv0l9BP<#Nh(Pe z&H+}3O6e7_xG5|ELij0wRIHWI3G$qdLHnt1LvDwKhBpa-N?H~=I!sR=^cdTz6e}cI zh1w_8Y&Yo>nhgCSQdEJElGcXfAQ&H0Ct)~T=1k~#3}`6W1dJl9WY~{zXeb_vK{2Kc zrI=32=v`QC3IwPe?denHDFo70K?-09#p$oC}nNa!y zF%r4asSTOB7`;iET6TPd=r=i%FCCRq(M)N0A~BZtDWjPY_BhBjqCXqu*m;aK8INx2 zOh+%K`(a4O5*tET6I`TcJ2$&pbC1~gSSY1LHg!V2qzydXG7SDK&9i28pZY(*x1Av= zp5D~CjYo;fA%4BHGc;>38qk>iLQGDd*$^67zT{a*DX9Anjq)h9Pl)pmi|;M6WtJ69$p_Ef0mdbscZ-V zA%^d^$6x|5hnQ!Y4DzxYy*h3XS#5;$|G$NE2%T`+P1HL)ZW4VMYbBOcSfKkr851et zC!&ga6$BUp;hp{~X+FYx=seUL^Z+amxpu6W&1net&{XUfqIaQ9IRMj1stLBQkD88= zNtEKaw@G~KTNr|#!|jD*Y^{LWlZ!xJPBPPA5L8S&0_~I@$$LlC4*~daDy>382xm@c z3h>ON?c6~(IUI3lb|ZMOS3#~C`q<_RFsU36nU(`e8=;{Q%*a589Z|3g0=MI=^CGYgG&@{&yloB|K-)&iPlerxF@xce3?~$^Tb%Sd zMK)t7z%Gck$h(0qHsdJRTL&38R9V`rA3o$zEEP9g`whLl=!Q6l79k>vQTTGu*h)yE zW-=g#6{4kC(E@)K32_8;hsr*@So(yN0_;!~JStHj$_a-tz9YNYb8Axjc+=vJ( z4HKHOBnygf95NmnQzbUT3GQ4WC3A{Pz;nofXKsLg`4?7pP6wX`dWa3>l~P! zG^xcvI(3&PNA)WACyo*c>{GB3P(oFS`6nO005Q<6GZ%C z+Ok1~V#*$$La|f6*v0E+!>q=-HoRY9mxC2?nv8ALT#SsaReLN0_%naza| zCC62gc_3z7MRtb*28f1JCEHMd1zBYx0V5279T$=|gn9+N9iYM+OF@yS$PT8I<-o7n zP9di->gPq`vA1`%wE(E|*#s)SD#k5k2oox(gaiTMjL-*58pQ`tL%1>>={~sXNFqQJ zeJc*tZWwPXzsl=o43nvWGyum~x*KN?2L@{^ zpn}P{KVs6D0>e(*eVGe`ahvS38!?B9@k zhs+m(EY@q|@SoYz!+&$Ff-@e`i?_Z3jh##(JMwLVjYicJoIn6Nj=c`-+=P-y{VUR- zN}6HFh>;9UONN?^v58iP?5Ic3{3SExZ>Ih1elWb`=qXW=HDFKH_Il5QS>i^NO@o4S zZoNbPQCKD2^!>NrR7f*rdXJDXD9roMoht1>2ctQ`E%qYdAv&|AVHMu8=lfTG2^mdZ zur)sE5o?|2Dc=tiIo2CiL+HtAQ{o3TIwq>EEIsC=%&Op&l(F03%%3_nA_m^a*#smV zc3@nWI{fI{hdB$w@sp5JM&BG%7<6?+uJbqs+wNpo(VwCrl5aPU2u`O}y0cS41;O)V z!vegRA&5$+myqj&ezHZvIfGM?xaVSiG^JJ~v)H(qE;9BJlW`TPCZVZ~837tcuz@0Wa&1FGLjQQ zH#kzIj56C@AU$weO<2w?B0z^QFp-_f%H4p2vfT&^$yr3_$gOW+i-h84GvDP>Jf@lv ze2ARoP7y*y*tLQHxfG3xUmW99xqYz0JSG`^ zj7oWAlWRaDHpvCPN?8400uP{urF;-ykYuII)Wp*a|I;k&PRUuy3LH@NA%bM< zcd=2BfCiuh1>>X;CBy(ilOQF4o&eP0|61EqCKiOtGs&@av9>0fxUCU$n@nLxpNK4C zg?=bpb5IePbz~ceL!c^5P8GH`)9g)kWTw7I{ZOKe z+$iMak(@KU62lcH)`h^3*eGG_JJw+>`Lg}AA)l|@ZYIdpZKTA$y&cg#f;O-fjSOlD zmqh42+#5NZxY*FgR*4m53;GW#T^OgwtE|9>vYMtUJ%YT2`rD8ym*p}2cgviTc6JIA z@f05$FgVMK)NHg+2Eg=*VS62(p>RW44LDR32=3C5>_5u10VHD7IRRvgMv$381vFsY zcp%aO?6xqP3g9KZEOsbCDdHW=Z#AI95KiaWP~+e*f>W~HSCTNHgy#g%A^hM-d&s>uXG z*O+G(Y(!Wb24QiMB?+(E^eF@Z2phKUntMI_NncSk%mT~GTWp5HV1N`-92>E~zUk9b z%w?p5F&AZs(a1#_!^6kJM|$j!KsCF0!fA#&U{_UWP>~BG(wlA`iH|KF>~Wc11?$mb zN;h9FP8C?fz8l$=c)sGGupLh%>cmKFsbB=)e{-@9m_-yYsH_KkY9-j^&@4Rx<|7^? z`85X+!5YU?z&2t8;t>8Sw7xZy!+{_=kPuUxVg<}|liU>!6%tsa#?08J56iGEg|J{z zg%WcxP^@3U$Jp5uUY1aTv><^3Mu#*v5PAQ>3ASZX^#PExOUjVgCppDol9dG0ZOl2_ z$vH^cy=UPsSh6L|j8%r9Z*E;+=RtvF5xmlVitrG{QZYYewRk4D00J)Il$gFmjL$>C zP#K6Ya@IUCw4dOk)Z948v=8CLAhtu8jYHx%i(qxy2UG_h!POQY5F#L?Cf>sNaxoTX zpdNP6#E=<7thN#RW;V@2mn~ToiV7+_q>steE8NzQDgkIZg7VR~DLKqr-oS4$P@TIT z;NT+XOlCoZU`)k4`{vNy)=WfIidSNY8zy2nW`GkVye|OZ6kAY&)RaZUYFtW;^4G$+ zg9VI(6Xw8boMtvO53e0Hh2m>QWM(?6V01#F>4stmaeT|s2GH1hQ<)vKPF5{Y`DYOw zbH0rZxjFWX=8}I}7tJJ?RNW8256x0QL`nAo%%IaehXG8|rh`yElnolY^;Px*tRq;g zo{48gnD;eIhMyE@B-}^ntV>f6bKqyL|_(Y5I|N8X`kaVfUnSOGID*-1XAE5Rogq^3p?ArEhEvcn`2 zo`u~rPazh?(IQINVJ+Z^iSCi3-Hbl8wU|)-Iu3ykBZ4=gj?r7+O2KBpnWY31y7l%J z#Qb!wk?cCo+=(q9I!1=$K9~u4PqH^u%dRk`)U7>JY(W4J@=dS@K>JXn2c5?>>6i+R zc#v<+7!gNL%#$$Kg%wdNxuB>&jU&0xm^7C)EDR>LPejay74-IwCyDcLkP%tA6jXF7 z2@sbuCVE7&!;v|_V$Jrho5Ibt-YI=K9 z*vyBbLLP-^gmXx#TOS|_rv0z;LIn9llHnST@Pu+OwY^0mzQKclM7cJkghzk@XiE$Ah>*Y^2$4x`T3QO&;AB!CRH5lnVjl+1 zEblNHhk+@#V7QM^zPyRc%m%Ok3__+c9Js2gvRMej@{Th$p5w-O3ye^(IAT}qg5gLS zqA^6Q&_ie>JZ5Ho-p3_m3qnGRD$zJ%KZq?!pE)$JN2fXoEwhjSw-zUQf#M$K10pTl z0-R$;f*rYh9ERbvFK3V)`*RozAt(0lcvO;`GmJZYh0=#&Arvy=jQHG+-QbY(J|e{t zn4ENgI7r5^1Lj9hmpU-jqLo}X&jm8S5=IP;C8>v|FbM0~Os6AADj$RdUd(Mk%th`6 zp_zapg0h>54RoLn^xB6fV>o*#_O7kyp5i`qplujhX+a*u0v^ktqo$M!630Hl^>AWI zQ_zw`Oouw!heSv;jsxhx=TZ}KF1XADJg1JhshxHayZIpzQ4qBZ*+2{k_)R3?L^-}$ zs1jbT$O(YrW-bkr^KbxF>?~4xEp#4dSdD|*JVFv2O$UU5YZI3?y{j_y?i>o<{!sV8 z6m6e{sG~F*0m-Qglg@^#ZrNkhgHp&6L-R5VhH<8=!Z=ixsTzsLAQp6tPCbcPep-_- zW0NNjrA;b@)J<6ZL(H&^w5U4`1aaVvIgwWdz>GYrpOe7}QG&gR?*suj-9iTl;}lpq zaB2|4QauXJ+>5gfCf;cf0i|$bD8?e&1C=HmKvN9WT>#zXGu!U3X^V%Tys*Pex`P#~ z1*$ln=Cy$&Pck8_4#CFnee1l}9MW(b|tro&H!LSqbv6q&5aa~Vu)nP=gq zlCdRY%6gr2gQ0_jngN^^-%+GO`Cu3&m??ve3%5El)hn0xg;osU-||qVY9(&OgN%O# zD(I`a082Z7GoO?MV4z61^Fa+VUWLI6q=RW~W&3$C#M8!t05p?(L%!O04T^URwqh3A zvoT&HfdZWX=SqViiWnovDV`&91QR)lbwU0(p%a7qjT3EzN|6oB)KW?Y`>lxYs5_m; zrRU=al;Su~nRHMV? zR8HK_Q7!)V}oE*azmzL8)*RoqnKdiD<{(X*fmaNR=~&I4lfC>NQp@g zSrH3(gH6YN9Qq0EHeZH|s{4>h*$3AUZ^cpa&v;2LiVb`5nhxF{!=YsILr`rWt0TNX zd(8PRLpXGRq}e!gV;!jCZ4=-C^N?J^=K%095ho}F>gaxK%S!`wp7s{P!DRa+5#9l_ivH)YO zkH8h?rYG`C_BSj|9I2InWb8iR%sXU}$imSaL;G+owIp_mCFgfYvCJ|$F zB6^5RowKMiNz}z*^IPAr^ttg8i{N3FAYhp2Wxf{xY!e;We1z3EAdNtvh44uWmxRay zK4k*&JA2j%D6*60Iwz^#@Eq7!!8*J{aRf;S6qthnDZER^MDo8wrQaTLS+4yBa`T&L z;Ic;hX2E8P5Nx9e>>gWwjGc8dOu}gbN&&8TF@!u{8}Px93kqa&2y_5r@jDBfqlZer zp)gz!+d1o>6$3D)K{P~!cmZ-2OR%rV;(At)&lDrILyXU-EHSwaD|a!tBj0W*8}UYN zLEhsf4C1sTbD)>dgQG#DOB`yY$b76&j7(El9mVL*G{tjt>>Hs4%W46Ok>hpi9iYJs zb3_1QPFfs8;mo|9(X_2p3CL8YlDM8J|Ea#8c8PNlOrm<_v_+1+J0BCp@ZX^ z!o}EWga!sV3c2!kF@TF7Xv1qKu>w*I9B`Ylg)uh#1;#r@r|>x+=^(D)RFZF)dPH{-7i?G^isU>>^L8O*^#7h zB~y7AXEx8lq)y};3OI(H8hWe3$tTa7z|RJn_BycTCY)9qMmRRM5)G@bVn8QrJZwP6 zMuM5>16zmXgH23j)F*|9D@+$&)|I*R_z#x?a!b(414He$wM?Q z9TLWwTtxh|WiB#4h6V?$HbHY~ z#&}!Qlqn;V!nC0E5YkxAG7(0<;FXk-$8g~-(@qR>NP;e2u_1I&yoAVxb2}kG{af(M z@Hh^>Gc^oklG+f5P6gZBb~dy>20BDkN;#S@>25&ml)&pOl*hqlC92&CEM_8#!Uy@4 z3)jZK_3i{rE^lWc*AMj!7ZE4Jfpa!sF)^^PCRo1kMks)^_BdygiCY`Q;1`W={WhJV z6$_e0d_^_m9K0cDV1c7~u%xB@h(DZYnLWy2PsxvirtBjcqSN-#2mQeD5TgcihGSs7 zr%Wuv4iX7jNaptr%q2gdxA`fwn-iPRRi=Zk;5QyBlKj$IijOQSbl#57V}10_@I|B% zlO2q7fKD&xAtUei5cG6h%2JlDk95dunioiRDSQ2Ul_h4#=P<3#7i-lKi2n`_%CC6y0BaU_? zo;zSVlDx@lhDHpL#2vuV&T*;zY@S()Q+;?>Lf~Cs04JQt1?FtL>0gT+E|v!Ome)b+ z#elP3vk4Wx(gJaZijY*i3@$XiWkb;2$U>Sj3bjOfQEZvv6~vFeeMUkyI?nMvn_9OT z00QY<3Ch-y-X~-X=9_J2VN{Bo8&mN67;JVs?hfXi6(Xk|j(sb%>6kBng2qzeYUgOnwI1e7de<4;ePL~6DINHYVkReG=pVkt? zDmkzgqpbp?@ta=M8zeJ%O`wIJh;%3oz=Ts1g(i#q7W%UT21}ffN3WaFKFJH_qL)$q!!Zx|J9fdnM_kg%Dg73g>r9ELSD zb9$tfPGd^SF|d*|+!hRj^T@z?aqh;hH0gP9KG9?qQ&XB-29{u%$eryF5@Pf1b8|PP z@TK^ON{3PelhcZslMu&06ybuD#9c%1YM4yy`8>Ll^=l7S?2ny?xEfQ4{a0YXU65XL z!NmT6$w@XB-31tp=tYNnIV%KD3Roj`nan&jYzU@f5Ru%Yg<3RM$jwzL=7ceT13>+q zKa(K%!OqXKB=QK4U`qy%^ym3#0SC~+AgHVr9nJ*kSis8ROi6!#zj_D2flc9*)9ZjC zoRTJ8O1Lmj##>otN_CkTNANrk$$Jy8JV%Zw=LIJ5Gd8b;-N7yLO7xly{A>Ki!~_}L zkKQ4-6l#dxXmd;ozlB~zwD8vvIz~@wvmQWLQRV!?cA8kLtNr=P(KV~dhmfQBQY230d^ub^`c}p}G5LxMNT9?1;@7hH z1l?%zB8DWJr7R%e{G;@0JoYg1Wf%CpvmiY!Q=<5Nsh;k9z33H{N06q%RxpkLwr~Uf zucGs{gp3j7bfY7Ts7mX*LfH+V2|uC3$3gQo0t+Kzv3}$-cz-2djvep|z(oziJCfnN zM7~OVqmlEqVtq1xZVM-0#`5*{uY(zmTE5bHgkatn+4u>ud_YWKeyufrbjS4Z=uLrq zgY$QH%o@r9IEE8Vtl^BW)$}WACKy%5RCtyTSa49fAzy7&@t0bwj)Td3jZrDS)t#YV z&KiL)EbAI#Hx)@v1Xyg09^10i&Y- zZbjz@oqF#Et%9!Hzy_mY<=u*TZ+GeodbJ9=a!Y!Riq09o$KURRYvmfQ(ptk=|1zp= ze0!z7@ENUwu3S$)*IsGM`k$h;HaJIJYE;tisLhmn%l3F@JdP60=e<2oZ|%`4=rU1r zWkHsZ6Zpb&deaiEg05WGQZY?))*pIM?g=#g%APm(yfOU1?`dv%>xD1A{P~x!^=jdv zY|~b~Y3myYKCTR0N?eI#D?9bdPMUM+wb-48`ERYi@=~_pX}#g;Oa7S>Jmo4b?op%r zRM^kvy8F?Oel*wP2UVW3iqGuQ0>jzT9eU{wt#k*##G}14o|1b1Jx@tl<$WsfwX9OK z-m3-IW=q%UrR%iPb$3e3sjaVmU&%d>uWWzGeJ=3?_gH34?VIa!{%~$yBVB>ITy1cs zy1b_Ro(IK!D&E-sAn5T|eP&Xt@61+o=@nh~JgsG0ynv&oR1lZxr1S;CTwE^6H8cyl zH}q-^y?0s`=q*dHH)t*WKf>%Si{5(unsmJ}J3p+?4_~VJI8gWIajj)jHn3R_Y{m$c zOSt&BCiv!_%Tw8!ZoQ^ktLdI8Gi$gSXsKKRtQf8Beg$_-p*dRlJxeXoyW_8V)34R9 zxLL0G2ebY`-9Lz%k1HzmiusJVSG{FBz0|5g3m&@Rzq$QptG;~Ohr9LV2k-z_)_+*{ zAJ+VbX*DOVEc`&xmkoYcqc7X91rB5b2lT*!8Bd^+&_v}O&~LBSGL)?W4!3GGTlv9^ z$J2VzcT*Oh&hnne6PWbf=Mv9xKT}%LvJ<1(<)vMM{a@4;t-A5zhex%hJ=vx`defdu zJMPvsU1_=2adqMCx^Atm`%ZA)TU8jQEeYoG)K7r9sKgVi%+=2aUk3mw+WNKDV|mV{&c0`g6b>SqJqUR!tKY2F8bO6Ex_^u2-$Gcg253!b zFt?{t1$R@4Cn&j6Ys_T*uj&5RH2-UN=Wo2;bp8DM?b`f};vxpCzub4NG#gl=2f$k1 z%01re3+}Xc>Fxd3&;EY$YsqYQP!A7g+qdZLTmDn~&fD!fv+cX}_T85bKr5iN=p#xf z20VdgD0uE{!YMWUM!(X;Mm_-bI;?i4BV#zUr{uB z+WRzrD=5}9|E+zwnoYTy_T0RNTyV)uWmyv_gyKFGmnxxko963{-P)23dgHT~D(_Y| zUh#f*^QFxs;(MT+T1EiBFYpTBrxH)>vfNWU|4PSaPhC1?wT^hHwG5V{l3H6Ox9+|j zSfm9O8ST9-r!cxVFoMyk1S_J1fgyvTH4SB}x9ZhfN%2FU)z#~Di}kuyGoI=Iz6wcC zXlh0azLpKXrUzfc4X#}Mf-Cav`c7DqJBycGef_%h&88dan?qke|Mm7iXwNPl)E5t4 zuK5_KIj)8FWt;ZvP5S|8eLoi;H@3XB=gL&Ju}^R8(;EAr$Bi1U#_H<(u{h^5yR&^K z^!V&n<(_pLzu?n@Aw9U-Si;5+*J;6h+2B4sxDPkDJ{D^m&IWhr!5vzVtg*j}wfi&o zJcDJ2P33;}#=4v7o4fT@JN~>!+xfD%^i?OdQ)Bw7Gg|O0uB`u@?mwsb&pk}n*HuP9 zEGk*GYb}G>nk{x-x(lpy^PqO* z1?^;7Yg4mrs@|qv+JCpM?ONRpuO3=?yKbdcx01EGwlI+mCiP%a<11Ich&TX12n0yk z;HVxP)qo`&^1A4=}DMwI=IdtNYh#{BRUY(Hur_-3?HSYmZBJ&QCT0x1qQC=#)TU~$o zh1=C_T6J5lY2l?EhJ=;2wA}NQ23qe^F=H1sp4!IC@@G$8I+<$I;Si`LTRUuY0KqD3!e>nHX4$xG#EP&02X`ud1sKB5>35=*0ivTeMa zCe@s4T5QpGc56%5uFh5u=+y%tXytPjeb>Ry3E=BbmA>WmhoNhzEY!8RE!)4a6~;(! z+ypFxEp}<`TRyy~H6F?~9?}~R;btaacNc>ohQLq=BCKJ>^I|0{g1H<22-Z-L`%pG` zSPvf7f`>n@4siIf`krS|z~u>pPTV42!MJhh+5TJ+x~3ufF6Xye}Er=H`*I>o3_oA z;FG^_LH4YP>{%1Oq9`xxnm6Myi_3KnYCRoYU-oM~n{!Lo=E5trbq8S)I}VeXM~R*& zF9&9-JdGXK4&E5jyVu-qSfe$pxzpJE7W@ZKh(bb76NQ9S($dT_*VGAhUw__5ZFiw& zkxELf={zkOCACFa?rChlB7J`2<&D>NT;K7f1J@4RIDa#9d+ly*?e1LPa)V=fw8oxX z)1qtM&%b>6a?Pn_IrxXxKk9U}B2Vb&qg%-@=}a$AUZN z8&^+m@L+>!nwQ`6tcPj-=5zYMus(1=d^24BjamOD-M>lmZ!(NWMH385MP=E2DrTs7 z(C7)ygNJjiTCeMa?yVg65x4Eqn#>rXX(R4&0WUkf%qxMTgfNuEB}bV{%<^N=cdHR! zG;P)RLf`$~IJX?#GEG&|*_!{rlwO%=Kz5eQ{ z+jY-qbW5_9i%1KAEG6r#1pH`)h{H!1ZyNkh?)h*2(0=S2rQJ;7UeTGG6YKQ zF=FmJkmN+z?0aUh!a+o{2g?l~WqH=WLiexG{41cD$Vpz8tzNHJucwIAHq4+*%|i|I zUCRv}cmmUGtIc0D)eX&Uc*u|)Z`p{Ctuhhjl<)BjGukTNZA8o0yhHS1e4`JLdG&UB zHtajZUk=lk=X@vl%ggiy^N13RQAzO(m3V>@YP+%;|04s)29;5^`A<6a4SVzrNBGP4 zUOJ(_6w_Z)pohvffA7L1dS3MIq5&>?PlC2Whlc7ea zr0L``tMYlSJjRvBxbm1?`2ulp<%iQ{0oS?$xV%0Yl~H+}?-7$QKZc^*Vgj&CR z;#$k~qt|=%j@36SwShh2(mVEQ`%Y@lpP+KzNxkEgHZq}|Nl-a5p?AEhDbsq#>mc+N zpF~u9&d2KKMc*l+&!u?s@RHfF za1f(Q$P;Y4f;@-5c=_$R&piF^^BhRY3DRmJ*TNq0&0~10L=xPkJRtBY|VDPX1i9iy{J`(V2ieEvo-7Vnsr*uIzusd zTBpdg2ByihqQu}RbBlVeFTEk(Sg-eP`mj~od_-J&?@{gf)7mMS%I8n(y<^&$N$u4O zRGyjCdo$=gC7kDH4CnbQfqTM7$tjQI8zW9T;X6a$cKf)FKJf6Fb3~cN=X{qA5@fqB zk7(foTFt?1%|X59Aa1bS)`eI7x!}*%}w;?2ez zF};7whr6|*W8%{LU(jBRY0?;#FQW3Sc5X^bXQ(_krT0&x_ouvj2vOm^guHz|=7Xrf z;W^){R2B6pjyh^W?b&bdRDF;vfa<@pNLz7I3!cgbPwBx^xB+6jmV9|>uJKT=EvfO9 zTd)ZD?3xFBcFhAmyXIZ4qbQW*g7xjpww+b9Wm$?|J6e10TfD zTIgi9@s!?pN^3mD91gMXpkTO(g5f5LgHe=wdIzw*(0DVZ5A6E$-P-O`;?f63v{zo$ zVo54rc~u`Er!y@=v}xf67~kjvYh1(r=uQ!-tzqx`sBatba@vRe4g4aZZwv8v+BZZG zjXL_E&f7@yujJM7D84Z=Ar5iJ}%6yw|}N*(!63Itf?JgsW!4 zRrBS&z*P%y_0sM6FKP2%!lu<#a$>tDiN#wp1TjiHL3ww6C%o;ZDPL~kQ(EuI+Y3)> z3s2r@TXc2rZ|%Rj|H^*qTzzdg+tRDI^lB}%19r7q>t3C08PHn>w3Y#WfU%am#-mU1 z=umI$QhrMF4d$*#<_FnHYqruaToxWUmkMXnveDK(teryhDhxXDQRufv2?drhW!8Nl4 zP&4SMS{z*aZeP>RK*^s4N_W;*{8^(H9OXDzi#!S{F~~gm-|Zo zrmytn4CwW#8U9bT0GGraYw$xz8Zgi0s7EB@hWFUtR2 zwsAU8@{fVi>H3O)Z1kf1&QR~r#-VL1OTM#hW$DniH5K34w%*I1yBblkr>*2Wd)i8e z_Jk_Fv!}~TpFNE@fPz>P{wCAAY(d13PUi1CZpdW@28VzDs^`X>Mv>NwI!OYGQH*P)fHvzbLynBfmg5u_Pz4SU0~QHBUDsu_RGH zB|o`XzqlYVFS9s9AE-H}vN*L^KR!M)FS8^*Uaz3?7l%!5eoARhs$CH)&@zy7iur-W R2WCb_#t#fIqKFwN1^`=%H01yQ literal 0 HcmV?d00001