|
Discussion Forums: general  Start New Thread
By: Juan A. Garzon
 No attachment found
 Simulaciones [ reply ]
2016-11-16 18:03
|
Hola Corsarios.
Espero que todo os vaya bien.
Os pongo al dia de las ultimas novedades con objeto de ir retomando
el trabajo e intentar salir del aparente letargo y avanzar mas deprisa.
Hago un resumen rapido.
1. Los que aun no esteis dados de alta en la forja del CESGA, os ruego
que lo hagais lo antes posible para canalizar toda nuestra informacion a
traves del foro CORSAIR.
2. Se ha incorporado a nuestro grupo de investigacion en Santiago un
nuevo estudiante de doctorado, Yanis Fontenla, una de cuyas tareas va a
ser hacer las simulaciones para Santiago de Compostela y, preparar
trabajo para que sea repetido en cada uno de vuestros centros. En cierto
modo, actuara como coordinador de todo vuestro trabajo.
3. Yanis ha pedido cuenta en el CESGA para correr alli nuestras
simulaciones. Para evitar los problemas que podais tener localmente en
vuestros centros, si es necesario, contemplo la posibilidad de que casi
todas o todas las simulaciones se puedan correr en el CESGA y que
despues, en cada centro, hagais vuestros analisis, parametrizaciones, etc.
4. Nosotros ya hemos hecho nuestra primera tarea. Es muy simple pero,
ademas de que me ha permitido familiarizarme con Corsika, pues no tenia
ni idea, me esta permitiendo entender un poco mas lo que observamos y lo
que esperamos observar. Podeis ver nuestros resultados en la hoja de
calculo CORSAIR_3361 a la que creo que todos teneis acceso. Todavia los
resultados son preliminares y estan todavia a falta de normalizacion
pero ya permiten ver un poco los resutados. Os hago un resumen:
Resumen
1. Hemos generado sistematicamente simulaciones inducidas por protones y
nucleos de helio, carbono y hierro, ademas de gammas, para Santiago de
Compostela, con energias entre 10^10 y 10^14 eV, por nucleon. Hemos
generado dos puntos por decada y, en algun caso, cuatro. Quizas lo
ampliemos hasta 10^15eV, cuando acabemos de instalar CORSIKA en el CESGA
porque en un PC normal, los tiempos de calculo se hacen demasiado altos.
Hemos abierto una hoja de calculo para cada uno de los primarios.
2. Los numeros 3361 en el nombre de la hoja se refieren en que hemos
puesto los umbrales de corte en la generacion de particula en los
valores: 0.3, 0.3, 0.06 y 0.1 , en GeV, respectivamente para hadrones,
muones, electrones y gammas. Son valores un poco altos pero de esa forma
son bastante validos para compararlos con los que resultados que en un
futuro pueda proporcionar un detector si tenemos en cuenta que puede
estar dentro de une edificio o que los electrones puedan ser detectados
como una traza visible.
3. Hemos hecho las generaciones para diversos angulos de incidencia del
primario para ver como afecta el espesor de la atmosfera a la incidencia
en tierra. Hemos elegido intervalos constantes en coseno del angulo
cenital Theta, excepto el intervalo vertical. Los intervalos han sido:
0-12.8º, 12.8º-22.3º, 22.3º-29.0º, 29.0º-34.4º. El primero tiene una
anchura de 0.025 en cosTheta y el resto tienen anchuras de 0.05 en cosTheta.
4. Nos hemos olvidado de las salidas root y solamente hemos mirado los
contajes de las principales particulas y sus incertidumbres (estas
empiezan por s): protones, neutrones, gammas, electrones y muones. Esos
contajes son los que se muestran en cada una de las columnas. Sobre
fondo amarillo estan los valores sacados de Corsika y a la derecha del
todo, estan las sumas de nucleones, gamas, electrones y positrones,
muones de ambos signos, cocientes, etc.
5. A la derecha del todo, se muestran algunos valores diferenciales
relativos al numero de particulas de distintas especies que llegan por
intervalo de energia. Por ejemplo, a traves de esos valores ya podemos
deducir que la energia mas probable del proton primario que ha generado
un muon que observamos con nuestro detector TRAGALDABAS en Santiago
proviene de protones de energia alrededor de 42 GeV (Columna NTm*FxInt:
numero total de muones por proton multiplicado por el flujo total de
protones, integrado.
6. Ahora, poco a poco, vamos a repetirlo para otros angulos y otros
primarios
7. Observad que los errores son relativamente grandes. Por eso hemos
decidido generar un numero grande de eventos que podria ser mas grande.
Como en principio queremos ver pautas generales entre los distintos
lugares, primarios, energias, angulos, etc. nos contentamos con dichos
valores. Aunque, si podeis, no estaria mal aumentar el numero de eventos
generados.
Proximamente voy a viajar a Mexico y Tunja por lo que ya podre
hablar con algunos de vosotros para explicaros detalles, pensar en
mejoras, etc.
No dudeis en preguntarme todo lo que querais y no veais claro.
Adelante!
Hans
|
|
