Probar con simulaciones restringiendo ángulo polar

[caredg]

Para simular correctamente lo que sucedería con un sismo, el efecto debería venir no de todo el alrededor polar hacia el detector sino solamente desde una región en particular. Por esto nos sería necesario restringir el ángulo $\phi$ en las simulaciones de Corsika, junto con la restricción para $\theta$.

Esto no está disponible en la simulación que maneja ARTI y tampo es algo que Corsika pregunta al momento de correr de manera interactiva. Entonces, se hace necesario correr:

• con las las banderas respectivas para $\theta$, $E$ y $\phi$
• probar qué es lo que hace la inclusión de esa restricción

Para hacer esto se podría proceder en dos fases:

• La primera consistiría en generar los scripts de simulación de la manera usual, sin modificar ARTI ni ninguno de sus scripts, pero luego, a mano, modificar el datacard correspondiente al ejemplo de simulación que se vaya a correr. Luego, estudiar el resultado de las simulaciones y entender si realmente funciona esta restricción polar (debería ser el caso).

• Una vez que nos aseguramos que existe el efecto y tiene sentido, se debería modificar los scripts necesarios de ARTI para añadir la opción de correr con ese ángulo polar, i.e., añadir en todo el andamiaje, la opción de incluir unos flags nuevos que den la opción de restringir este ángulo. Sin embargo, antes de intentar hacer esto:

  • también se debería probar correr con solamente flags para energía y ángulos, es decir, no correr interactivamente y asegurarse de que todo eso funciona bien

En resumen, los tasks que hay que hacer:

• [x] Correr la simulación con 20 o 30 min (o menos dependiendo de las otras pruebas) pero sin utilizar el modo interactivo, sino directamente utilizando los flags que soporta actualmente ARTI
• [x] Correr lo mismo de arriba pero modificando a mano el datacard correspondiente para restringir también el ángulo polar
• [ ] #34

[pablobuitron]

Añadido script para cambiar con la mano el ánglo polar aquí

[caredg]

Parecería ser que la restricción de ángulo polar funciona. Al menos los archivos logs de texto sí muestran que la restricción funciona a nivel de corsika. Por ejemplo: aquí y acá. También parecería que se puede ver un efecto en el gráfico: acá (simulación con restricción de ángulo polar de 0 a 10 grados).

@pablobuitron tiene la idea de restringir el ángulo de -180 a -170, e.g., para ver si se consigue el efecto alterno. Es una buena idea.

Al final, creo que podemos declarar, luego de hacer el test anterior, que la restricción polar sí funciona y mirar el último task de implementar esta restricción en todo el andamiaje de ARTI (o buscar alternativas).

[pablobuitron]

Esta simulación es para una lluvia de 30 minutos, de ángulos polares restringidos de -180° a -170°:

[caredg]

Ok, no creo que se puede ver realmente. Hay que mirar otros gráficos, inventar o pensar en otras variables a explorar...