Summary

This is the fourth installment of the report of activities that I do as a participant in the participatory research project Citizen science particle physics project on Hive led by @lemouth and undertaken by a number of blockchain science enthusiasts. It describes the tasks that enable _Deciphering top quark production at CERN's Large Hadron Collider_, instructions for which can be found at this post. My first attempt is recorded here, in which I completed tasks 1, 2 and 3 completely, (the third with a doubt) and I report the bugs that I experienced with task 4 that was not completed. As for homework 1 and 2, I couldn't do it as I still have to learn more about it. I dedicate myself to describe the basic procedures according to the tutorial without going into technical details about the physical processes involved. This post may receive updates in the short term if I make other attempts. I thank @lemouth in advance for the respective review and advice. ---

Task 1: The event content

Following the instructions, after opening my madanalysis installation, I ran:

import /home/eniolw/CitizenScience/madanalysis/ANALYSIS_0/Output/SAF/_defaultset/lheEvents0_0/myeventstoo.lhe.gz as ttbar
ma5> set ttbar.xsection = 505.491

I used an absolute path to avoid any inconvenience. I then proceeded to repeat the procedure of the previous task with the new settings:

set main.lumi = 140
plot NAPID
submit
open

This yielded similar results, as expected.

Task 2: Lepton and photon multiplicity

I begin to simulate the experiment with **final-state leptons** by implementing the restriction or "cut" by typing the following commands in the MA5 command line interface:

define l = l+ l-
plot N(l) 5 0 5
select (l) PT > 20
plot N(l) 5 0 5
plot NAPID
resubmit
open

This yields new results: It can be seen that there is a small difference between the two graphs, confirming the expected results according to the tutorial. I now proceed to do the same for the **photons** case:

plot N(a) 5 0 5
select (a) PT > 20
plot N(a) 5 0 5
resubmit
open

Task 3: Selection of a subset of all simulated collisions

I proceed to make another cut to focus on the events contained in exactly one lepton in order to prepare the study of its properties. To do so, I execute the following instruction:

select N(l)==1

However, when I opened the HTML document again I could not find a table showing approximately the information expected in the tutorial. However, it did appear as soon as I tackled the next task, so I don't know if it was some command I ignored the first time. Here is the table obtained:

Task 4: Lepton multiplicity after the cut

To study the properties of the lepton, such as the transverse momentum and ABSETA (the absolute value of its pseudo-rapidity) I run the following in the MadAnalysis command line:

plot PT(l) 25 0 200 [logY]
plot ABSETA(l) 25 0 2.5 [logY]
resubmit
open

At the moment that's all I have. For the next attempt, I'll take a closer look at the exception raised, though if @lemouth or anyone else in the community is already familiar with the bug, I'd appreciate the feedback. I hope to report my second attempt with the solution as soon as I can.

Some references

• Towards a citizen science particle physics research project on Hive
• Citizen science on Hive - Deciphering top quark production at CERN’s Large Hadron Collider

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
• "Citizen Science" Project Activity Report 3 - Detector Effects and Event Reconstruction

NOTES

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

Resumen

Esta es la cuarta entrega del reporte de actividades que hago como participante 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 ciencia en la blockchain. En él se describen las tareas que permiten _descifrar la producción de quarks top en el Gran Colisionador de Hadrones del CERN_, cuyas instrucciones se encuentran en este post. Se registra aquí mi primer intento, en el cual realicé las tareas 1, 2 y 3 completamente, (la tercera con una duda) y reporto los bugs que experimenté con la tarea 4 que no fue completada. En cuanto a las tareas para la casa 1 y 2, no pude realizarlas, ya que aún debo aprender más sobre el asunto. Me dedico a describir los procedimientos básicos según el tutorial sin entrar en pormenores técnicos sobre los procesos físicos involucrados. Este post puede recibir actualizaciones en el corto plazo si realizo otros intentos. Agradezco de antemano a @lemouth la revisión y la asesoría respectivas. ---

Tarea 1: El contenido del evento

Siguiendo las instrucciones, tras abrir mi instalación de madanalysis, ejecuté:

import /home/eniolw/CitizenScience/madanalysis/ANALYSIS_0/Output/SAF/_defaultset/lheEvents0_0/myeventstoo.lhe.gz as ttbar
ma5> set ttbar.xsection = 505.491

Usé una ruta absoluta para evitar algún inconveniente. Seguidamente procedí a repetir el procedimiento de la tarea anterior con los nuevos ajustes:

set main.lumi = 140
plot NAPID
submit
open

Esto arrojó resultados similares, tal como se esperaba.

Tarea 2: Multiplicidad de leptones y fotones

Comienzo a simular el experimento con **leptones de estado final** al implementar la restricción o "corte" escribiendo los siguientes comandos en la interfaz de línea de comandos del MA5:

define l = l+ l-
plot N(l) 5 0 5
select (l) PT > 20
plot N(l) 5 0 5
plot NAPID
resubmit
open

Con lo cual se obtienen nuevos resultados: Se observa que hay una pequeña diferencia entre las dos gráficas, confirmando los resultados esperados según el tutorial. Ahora procedo a hacer lo propio para el caso de los **fotones**:

plot N(a) 5 0 5
select (a) PT > 20
plot N(a) 5 0 5
resubmit
open

Tarea 3: Selección de un subconjunto de todas las colisiones simuladas

Procedo a hacer otro corte para enfocarme en los eventos contenidos en exactamente un leptón a fin de preparar el estudio de sus propiedades. Para ello ejecuto la siguiente instrucción:

select N(l)==1

Sin embargo, al abrir el documento HTML de nuevo no pude encontrar una tabla que mostrara aproximadamente la información esperada en el tutorial. Sin embargo, sí apareció tan pronto como abordé la siguiente tarea, por lo que no sé si se trató de algún comando que ignoré la primera vez. He aquí la tabla obtenida:

Tarea 4: Propiedades del leptón

Para estudiar las propiedades del leptón, tales como el impulso transversal y el ABSETA (el valor absoluto de su pseudo-rapidez) ejecuto en la línea de comandos de MadAnalysis lo siguiente:

plot PT(l) 25 0 200 [logY]
plot ABSETA(l) 25 0 2.5 [logY]
resubmit
open

: De momento es todo lo que tengo. Para el próximo intento, miraré más de cerca a la excepción surgida, aunque si @lemouth o alguien más de la comunidad ya está familiarizado con el bug, agradecería el feedback. Espero reportar mi segundo intento con la solución tan pronto como pueda.

Algunas referencias:

• Towards a citizen science particle physics research project on Hive
• Citizen science on Hive - Deciphering top quark production at CERN’s Large Hadron Collider

Mis reportes previos

• "Citizen Science" Project Activity Report 1
• "Citizen Science" Project Activity Report 2 - Simulating top quark production at CERN’s Large Hadron Collider
• Reporte de actividades del proyecto "Ciencia Ciudadana" 3 - Efectos del detector y reconstrucción del evento

Notas

• A menos que se indique lo contrario, las imágenes de esta publicación son del autor.