README.md



<<<<<<< README.md

Proyecto de Simulación de Interacción de Neutrones con Suelo
Este repositorio contiene dos carpetas principales que abordan las simulaciones realizadas para la interacción de neutrones con suelo húmedo.

1. Simulación de Flujo de Neutrones Monoenergeticos
La carpeta Simular_flujo_de_neutrones contiene archivos para simular la interacción de neutrones con suelo húmedo, contándolos a través de 6 anillos concéntricos de 2 metros, con radios de 10, 25, 45, 95, 145 y 195 metros. El sistema físico consta de un bloque de suelo de 800,000 m² dividido en tres partes. Durante las simulaciones, se varía únicamente la humedad de la primera capa. En la ubicación Simular_flujo_de_neutrones/src/, se encuentra el archivo PrimaryGeneratorAction.cc, que define la inyección de neutrones, su dirección y energía.
Para las simulaciones, se inyectaron 10,000 neutrones primarios en ángulos cenitales (θ) que variaban desde 0 grados hasta 45 grados, con incrementos de 5 grados. Asimismo, los neutrones son monoenergéticos en el rango de 1 eV a 10 MeV, aumentando de a 1 década.
La salida del código es un archivo .root que contiene 6 histogramas, uno para cada contador. Cada barra en estos histogramas representa el número de partículas que llegan a cada contador. El orden de las partículas en los histogramas es el siguiente:

Electrones y positrones
Gammas
Protones
Neutrones

Es posible realizar un seguimiento por paso de los siguientes parámetros: "ID", "X", "Y", "Z".
:  ID de la partícula y posición de la partícula.
Para imprimir lo anterior en pantalla, es necesario descomentar las líneas 130 a 133 en el archivo PrimaryGeneratorAction.cc. Al realizar esta acción, se mostrará en la pantalla el ID del neutrón y su posición en cada paso de la simulación.
Lo anterior se realiza descomentando las lineas 130 - 133, realizando esto se imprime en la pantalla el ID del neutron y su posicion en cada step.

2. Simulación de un Flujo de Neutrones
En la carpeta Simular_neutrones_monoenergeticos, se encuentran los archivos necesarios para simular la interacción de neutrones con suelo húmedo. El sistema físico consiste en un bloque de suelo de 400,000 m² y un contador en forma de cubo con arista de 2 metros, elevado a 3 metros sobre el suelo. Las simulaciones implican cambios en la humedad del bloque de suelo, y la modificación del flujo de neutrones se logra mediante el archivo salidasimu.shw.bz2 en la carpeta flujos.
La salida del codigo en un archivo .root donde se generan 2 histogramas.

Numero de neutrones que llegan al contador
Energia de los neutrones que llegan al contador

Es posible realizar un seguimiento por paso de los siguientes parámetros: "N", "ID", "X", "Y", "Z", "PX", "PY", "PZ", "THETA", "PHY".
: ID de la partícula, posición de la partícula, momento de la partícula, y ángulos cenital y azimutal de la partícula.
Para imprimir lo anterior en pantalla, es necesario descomentar las líneas 167 a 176 en el archivo PrimaryGeneratorAction.cc. Al realizar esta acción, se mostrará en la pantalla el numero del evento, el ID del neutrón, su posición, su momento y los angulos cenital y asumutal en cada paso de la simulación.

Archivos Adicionales

En la carpeta flujos, el archivo salidasimu.shw.bz2 contiene el flujo de neutrones sobre Bucaramanga (956 m s.n.m) se debe rehubicar en la carpeta Simular_flujo_de_neutrones.
En la carpeta flujos, el script flujo_artificial.py genera un flujo artificial de neutrones con las características deseadas.
En la carpeta flujos, el script auto_simulaciones.sh automatiza las simulaciones de un flujo monoenergético de neutrones. Si se desea utilizar, es necesario cambiar la energía de los neutrones antes de compilar este código.