La validación empleada hasta el momento ha sido revisada nuevamente y pareciera tener errores.
Uno esperaría ver un porcentaje considerable de partículas mermadas para la geometría empleada. SInembargo se encontró que al ser partículas tan energéticas, para 45 metros de recorrido, más del 80% logran atravesarla.
La mala interpretación previa se debe a errores en el código de S3 (última parte del análisis de datos generado por LAGO, S3 es el análisis de las partículas interaccionando con el detector de agua y su electrónica ya a nivel de suelo) e incompatibilidades con versiones que tiene este código en las versiones actuales tanto de geant4 y de root.
La manera en la que se analizan estas partículas en esta fase de la simulación, es mediante un código en Geant4 que recibe como input lo producido por ARTI. Este código separa dato por dato (de los miles de secundarios que hay) el momento, el tipo de partícula y otras cantidades físicas que luego serán utilizadas para la interacción con el detector, el agua la electrónica, etc.
Justo es en este código que produce error, a manera de ayuda y con ganas de avanzar pronto, usé este código para el proyecto de la validación de las cajas que estoy empleando pero esto parece que ha sido la falla. Usé el código puesto que ya está hecha la función que introduce e itera por partícula toda la información mencionada antes sin tener que yo volver a escribir un código similar con mis objetivos ya que es básicamente la misma información la que se debe extraer de los archivos finales de ARTI, los archivos.shw.
Como detecté este error a nivel de estadística en el proceso de tabular las partículas y revisar la relación de cambio entre las iniciales y las finales, es decir ya en parte del proceso final, tuve que volver a escribir este código para solventar el error.
En el proceso de escribir un nuevo código desde cero, me ecnontré con algunos bugs que hay que solucionar y que he estado corrigiendo. Nuevamente, este código ya existe pero está implementado para software viejo con dependencias hardcodeadas que no realizan el conteo esperado que busco.
Una vez implementada esta función de input de datos correctamente, ya todo debería corregirse y el conteo debería estar más a pegado a lo que se espera, al final del día, son partículas que atraviesan cierto volumen de aire y van perdiendo energía en el proceso.
Como método en paralelo de validación y para no perder el estado actual del proyecto, he estado introduciendo manualmente partículas con sus respectivos momentos y hasta el tiempo en el que estoy escribiendo este issue, he encontrado que para la geometría empleada previamente (45 metros de longitud de la caja) atraviesan mucho más del 80% de partículas introducidas, nuevamente, esta estadística NO es confiable puesto que varía mucho en cada iteración que realizo por ser una muestra pequeña del total de partículas ingresadas.
Para ver cambios en realidad, he optado por cambiar la geometría por partícula ingresada, es decir, ingreso cierta partícula con su respectivo momento y hago el análisis de G4 para los 45 metros iniciales, luego aumento la longitu de la caja en valores de 50 metros, 100 m, 500 m y 1 km. Aquí se observa que, evidentemente las partículas no tan energéticas, se van quedando en el camino y no logran el recorrido (sobretodo para las distancias de 500 m y de 1 km.)
La validación empleada hasta el momento ha sido revisada nuevamente y pareciera tener errores.
Uno esperaría ver un porcentaje considerable de partículas mermadas para la geometría empleada. SInembargo se encontró que al ser partículas tan energéticas, para 45 metros de recorrido, más del 80% logran atravesarla.
La mala interpretación previa se debe a errores en el código de S3 (última parte del análisis de datos generado por LAGO, S3 es el análisis de las partículas interaccionando con el detector de agua y su electrónica ya a nivel de suelo) e incompatibilidades con versiones que tiene este código en las versiones actuales tanto de geant4 y de root.
La manera en la que se analizan estas partículas en esta fase de la simulación, es mediante un código en Geant4 que recibe como input lo producido por ARTI. Este código separa dato por dato (de los miles de secundarios que hay) el momento, el tipo de partícula y otras cantidades físicas que luego serán utilizadas para la interacción con el detector, el agua la electrónica, etc.
Justo es en este código que produce error, a manera de ayuda y con ganas de avanzar pronto, usé este código para el proyecto de la validación de las cajas que estoy empleando pero esto parece que ha sido la falla. Usé el código puesto que ya está hecha la función que introduce e itera por partícula toda la información mencionada antes sin tener que yo volver a escribir un código similar con mis objetivos ya que es básicamente la misma información la que se debe extraer de los archivos finales de ARTI, los archivos.shw.
Como detecté este error a nivel de estadística en el proceso de tabular las partículas y revisar la relación de cambio entre las iniciales y las finales, es decir ya en parte del proceso final, tuve que volver a escribir este código para solventar el error.
En el proceso de escribir un nuevo código desde cero, me ecnontré con algunos bugs que hay que solucionar y que he estado corrigiendo. Nuevamente, este código ya existe pero está implementado para software viejo con dependencias hardcodeadas que no realizan el conteo esperado que busco.
Una vez implementada esta función de input de datos correctamente, ya todo debería corregirse y el conteo debería estar más a pegado a lo que se espera, al final del día, son partículas que atraviesan cierto volumen de aire y van perdiendo energía en el proceso.
Como método en paralelo de validación y para no perder el estado actual del proyecto, he estado introduciendo manualmente partículas con sus respectivos momentos y hasta el tiempo en el que estoy escribiendo este issue, he encontrado que para la geometría empleada previamente (45 metros de longitud de la caja) atraviesan mucho más del 80% de partículas introducidas, nuevamente, esta estadística NO es confiable puesto que varía mucho en cada iteración que realizo por ser una muestra pequeña del total de partículas ingresadas.
Para ver cambios en realidad, he optado por cambiar la geometría por partícula ingresada, es decir, ingreso cierta partícula con su respectivo momento y hago el análisis de G4 para los 45 metros iniciales, luego aumento la longitu de la caja en valores de 50 metros, 100 m, 500 m y 1 km. Aquí se observa que, evidentemente las partículas no tan energéticas, se van quedando en el camino y no logran el recorrido (sobretodo para las distancias de 500 m y de 1 km.)
@caredg