Summary

This article is the third activity report in the context of the participatory research project Citizen science particle physics project on Hive led by @lemouth and undertaken by a number of blockchain enthusiasts. In it we will see the realization of the project tasks concerning _detector effects and event reconstruction_. The instructions laid out in this article have been followed. In the first attempt there were some exceptions related to the generation of PDF and DVI files. Configuration differences were also noticed when comparing with the tutorial in the mentioned article. An output HTML document could be obtained, but without the expected histogram. I ask @lemouth to please review my procedures in order to correct them and make the corresponding updates. ---

Task 1: Getting started with MadAnalysis 5

In order to simulate the particle physics detector such as the one at CERN, it is necessary to install MadAnalysis5. To do this I ran git clone https://github.com/MadAnalysis/madanalysis5 to download the source code, which can be run immediately. However, when executing it, the program throws an exception since the "python" command is not understood by the system. This tells me that the software needs to access Python using that alias. Since my GNU/Linux distribution installation only uses the python3 alias, I had to do some configuration tweaking following the instructions at this post to create that alias and have python call to the Python 3 program. Once that was done, MadAnalysis5 ran normally. Then a series of compilations begins, which is normal according to the instructions. In the end I get this output: Comparing to @lemouth's tutorial I can see there are differences. On my system the SciPy, pyhf and matplotlib packages are disabled. Since the instructions indicate that matplotlib is crucial, I decided to install it, a little bit late though, that is, after the other procedures. In the end I realized I had done something wrong and found out what this was, so I installed matplotlib for the whole system.

Task 2: Additional packages

I installed zlib and fastjet using install zlib and, after restarting MadAnalysis5, I used install fastjet (funny names by the way). Around this time I realized I would need matplotlib and installed it. I started MadAnalysis5 again and the log already showed matplotlib as enabled. However, PySci is still disabled since I installed MadAnalysis5. I hope that won't be a problem.

Task 3: Top-antitop simulations

I'm still getting ready to do some physics. I copied the path to the file tag_1_pythia8_events.hepmc.gz generated with the past activity and then restarted Ma5 with ./bin/ma5 -R madanalysis/input/ATLAS_default.ma5. I ran the other instructions (import , set main.outputfile=myeventstoo.lhe.gz and submit) . This took a while to run, but in the end I got this output: As we can see, some things did not go well. It did generate a file of approximately 7mb as expected based on the instructions, but neither the PDF nor the DVI file was generated. I thought this might be because I didn't have SciPy installed, so I installed it by following this documentation. However, repeating the process did not fix the problem. I checked the suggested log and there seems to be a problem with Latex. As indicated in the log, there is a problem with the pdftexcmds.sty package. According to what I googled, that package comes separately in the most recent versions of Tex Live or something like that, so it would have to be installed separately. However, I didn't know how to fix this easily. Apparently I need to do a full TexLive install, but if I understand correctly, that would take about 5 GB of extra space. Did I understand correctly? Is all this really necessary? Why does MadAnalysis5 indicate that the latex and pdflatex packages are enabled and yet the log indicates something else (see image 5)? This is a bit confusing :S Similarly, the log for the DVI files gave this:

! LaTeX Error: File `pdftexcmds.sty' not found.

Type X to quit or  to proceed,
or enter new name. (Default extension: sty)

Enter file name: 
! Emergency stop.
 

l.108 \RequirePackage{pdftexcmds}[2009/04/10]
                                             ^^M
No pages of output.
Transcript written on main.log.

Which is basically the same exception.

Task 3: Some physics

I decided to go ahead and see if this could work after all. I ran the commands:

import ANALYSIS_0/Output/SAF/_defaultset/lheEvents0_0/myeventstoo.lhe.gz as ttbar
set ttbar.xsection = 505.491
plot NAPID
submit

And after a short wait I got this output: Again the same exceptions appeared regarding the generation of PDF and DVI files. An HTML document was generated and using the Open command, as suggested by the output, opened the web browser displaying this page: The _histos and cuts_ part was empty, a histogram plot was not generated as expected. At the moment this is what I have. I hope @lemouth you can help me fix these issues and get the histogram generated. Thanks in advance.

UPDATE 2022/04/26:

Finally I have obtained the expected outputs. Getting to this was quite a "journey" that had many stops and definitely deserves its own post. Essentially I had to get another computer where I could fully install TexLive (2022), reinstall and repeat all the tasks of this project from scratch and read a lot of documentation and support forums to overcome every exception that appeared. But it was worth it! Thanks to @lemouth for his feedback.

Some References

• Towards a citizen science particle physics research project on Hive
• Citizen science on Hive - detector effects and event reconstruction
• Setting python3 as default on Linux
• Documentation of SciPy - Installation
• How to Install LaTEX on your Computer

My previous reports

• "Citizen Science" Project Activity Report 1
• "Citizen Science" Project Activity Report 2 - Simulating top quark production at CERN’s Large Hadron Collider

NOTES

• Unless otherwise stated, the images in this post are the author's.

Resumen

Este artículo es el tercer reporte de actividades en el contexto del proyecto de investigación participativa Citizen science particle physics project on Hive dirigido por @lemouth y emprendido por una serie de entusiastas de la blockchain. En él veremos la realización de las tareas de proyecto en lo referido a los _efectos de los detectores y reconstrucción de eventos_. Se han seguido las instrucciones planteadas en este artículo. En el primer intento se produjeron unas excepciones relacionadas con la generación de archivos PDF y DVI. También se apreciaron diferencias de configuración al comparar con el tutorial en el artículo mencionado. Se pudo obtener un documento HTML de salida, pero sin el histograma esperado. Agradezco a @lemouth la revisión de mis procedimientos a fin de corregirlos y hacer las actualizaciones correspondientes. ---

Tarea 1: Empezar a utilizar MadAnalysis 5

Para poder hacer la simulación del detector de física de partículas tal como el del CERN es necesario instalar MadAnalysis5. Para ello ejecuté git clone https://github.com/MadAnalysis/madanalysis5 para descargar el código fuente, que puede ser ejecutado inmediatamente. Sin embargo, al ejecutarlo el programa produce excepción puesto que la orden "python" no es entendida por el sistema. Esto me indica que el software necesita acceder a Python usando ese alias. Dado que la instalación de mi distribución GNU/Linux solo usa el alias python3, tuve que hacer unos ajuste de configuración siguiendo las instrucciones en este post para crear ese alias y hacer que python invoque al programa Python 3. Una vez hecho eso, MadAnalysis5 se ejecutó normalmente. Acto seguido empieza una serie de compilaciones, lo cual es normal de acuerdo a las instrucciones. Al final obtengo esta salida: Al comparar con el tutorial de @lemouth, puedo ver que hay diferencias. En mi sistema, los paquetes SciPy, pyhf y matplotlib están deshabilitados. Dado que las instrucciones indican que matplotlib es crucial, decidí instalarlo, aunque tardíamente, es decir, después de los demás procedimientos. Al final me di cuenta que había hecho algo mal y descubrí que era esto, por lo que instalé matplotlib para todo el sistema.

Tarea 2: Paquetes adicionales

Instalé zlib y fastjet usando install zlib y, tras el reinicio de MadAnalysis5, install fastjet (curiosos nombres por cierto). Alrededor de este momento me di cuenta que necesitaría matplotlib y lo instalé. Inicié nuevamente MadAnalysis5 y el log ya me mostraba a matplotlib como enabled. Sin embargo, PySci sigue deshabilitado desde que hice la isntalación de MadAnalysis5. Espero que eso no vaya a ser un problema.

Tarea 3: Simulaciones de top-antitop

Sigo preparándome para hacer algo de física. Copié la ruta al archivo tag_1_pythia8_events.hepmc.gz generado con la actividad pasada y acto seguido reinicié a Ma5 con ./bin/ma5 -R madanalysis/input/ATLAS_default.ma5. Ejecuté las demás instrucciones (import , set main.outputfile=myeventstoo.lhe.gz y submit). Esto tardó algo en ejecutarse, pero al final obtuve esta salida: Como podemos ver, algunas cosas no salieron bien. Sí se generó un archivo de aproximadamente 7 mb tal como se esperaba según las instrucciones, pero no se generó el PDF ni el archivo DVI. Pensé que esto podía tener que ver con que no tenía SciPy instalado, así que lo instalé siguiendo esta documentación. Sin embargo, al repetir el proceso, esto no solucionó el problema. Verifiqué el log sugerido y al parecer hay un problema con Latex. Como se indica en el log, hay un problema con el paquete pdftexcmds.sty. Según googleé, ese paquete viene separado en las versiones más recientes de Tex Live o algo por el estilo, por lo que habría que instalarlo separadamente. Sin embargo, no supe cómo solucionar esto de manera sencilla. Al parecer debo hacer una instalación completa de TexLive, pero si entendí bien, ello necesitaría como 5 GB de espacio adicional. ¿Entendí bien? ¿Es realmente todo esto necesario? ¿Por qué MadAnalysis5 indica que los paquetes latex y pdflatex están habilitados y sin embargo el log indica otra cosa (ver imagen 5)? Esto es un poco confuso :S Similarmente, el log para los archivos DVI arrojó esto:

! LaTeX Error: File `pdftexcmds.sty' not found.

Type X to quit or  to proceed,
or enter new name. (Default extension: sty)

Enter file name: 
! Emergency stop.
 

l.108 \RequirePackage{pdftexcmds}[2009/04/10]
                                             ^^M
No pages of output.
Transcript written on main.log.

Que es básicamente la misma excepción.

Tarea 4: Algo de física

Decidí seguir adelante para ver si esto podía funcionar después de todo. Ejecuté los comandos:

import ANALYSIS_0/Output/SAF/_defaultset/lheEvents0_0/myeventstoo.lhe.gz as ttbar
set ttbar.xsection = 505.491
plot NAPID
submit

Y tras una breve espera obtuve esta salida: