It is almost a month ago that I released the fourth episode of our citizen science project on Hive. Time flies fast, I know…
In this last episode, there were some assignments to which I promised to deliver the answers shortly later. Unfortunately, offline family time, a conference, and now COVID delayed this. However, patience is a great virtue (and a bad excuse ;) ).
In the present post, I provide the solutions of the various exercises that I proposed. A new episode of our project is planned to be released hopefully very soon (I started to write it)… It will be dedicated to the simulation of a potential new physics signal at CERN’s Large Hadron Collider.
Let’s first start with a small recaps of the previous episodes, which could be useful for anyone who would like to embark in this project. It is still time! It suffices to start with the four episodes below and publish associated reports on chain (with the #citizenscience tag). They will be reviewed and commented out, and I am always available on chain and on the STEMsocial Discord server to help.
- Ep. 1 - Getting started and installation of the MG5aMC software dedicated to particle collider simulations. We got seven reports from agreste, eniolw, gentleshaid, mengene, metabs, servelle and travelingmercies. @metabs’ report has to be highlighted as it consists of an excellent documentation on how to get started with a virtual machine running on Windows.
- Ep. 2 - Using MG5aMC to generate 10,000 simulated LHC collisions relative to the production of a top-antitop pair at the LHC. We got eight reports from agreste, eniolw, gentleshaid, isnochys, mengene, metabs, servelle and travelingmercies.
- Ep. 3 - Installation of MadAnalysis5, a piece of software dedicated to detector effects, the reconstruction of the output of the complex simulations into something easier to grasp, and the analysis of the produced events. We got seven contributions from agreste, eniolw, gentleshaid, isnochys, metabs, servelle and travelingmercies.
- Ep. 4 - A short study of top-antitop production at CERN’s Large Hadron Collider. We got five contributions from agreste, eniolw, gentleshaid, servelle and travelingmercies.
Just a special mention to @isnochys, @metabs and @mengene. Are you still interested in this project? Are you stuck somewhere? Don’t hesitate to to reach me out if needed.
[Credits: geralt (Pixabay)]
As usual, before moving on with the main material for the exercises of this week, I acknowledge all participants to this project and supporters from our community: @agmoore, @agreste, @aiovo, @alexanderalexis, @amestyj, @darlingtonoperez, @eniolw, @firstborn.pob, @gentleshaid, @gtg, @isnochys, @ivarbjorn, @linlove, @mengene, @mintrawa, @robotics101, @servelle, @travelingmercies and @yaziris. Please let me know if you want to be added or removed from this list.
1 - Jet multiplicity
In order to obtain the distribution in the number of jets associated with our simulated events, we mimic the syntax introduced for leptons and photons. The only difference is related to the usage of the label j (for jets) instead of l (for leptons) or a (for photons).
We draw the plots twice, once without any restriction on the jets, and once after focusing only on jets with a transverse momentum larger than 25 GeV. This is achieved through the code:
ma5> plot N(j) 10 0 10 ma5> select (j) PT > 25 ma5> plot N(j) 10 0 10 ma5> submit ma5> open
The important feature to notice here is that I requested histograms containing 10 bins (ranging from 0 to 10). This allows to have a small enough number of events in the overflow bin (less than 0.1% in my case).
Controlling overflow and underflow bins is always a good practice, and this can easily be achieved through the HTML page generated by the code. We get (with the statistics available from the generated web page):
[Credits: @lemouth]
The same exercise can be done with b-jets now (making use of the the symbol b in the code).
ma5> plot N(b) 5 0 5 ma5> select (b) PT > 25 ma5> plot N(b) 5 0 5 ma5> submit ma5> open
Here, a histogram containing 5 bins is enough and there is no stress to have with the overflow bin (it is empty in my case). We obtain as a result:
[Credits: @lemouth]
As can be seen there are little variations between the two figures. The reason is simple. The b-jets originate in the considered signal from the decay of massive top quarks. They should therefore carry a significant fraction of the available (mass and kinetic) energy of the decaying top quarks. Therefore, the restrictions imposed on the b-jets (25 GeV is only 1/7 of the top quark mass) almost automatically selects all available b-jets. The two figures are therefore almost identical.
2 - Lepton multiplicity after selection
In the course of the analysis, we imposed the number of leptons to be equal to 1. Therefore, if we plot the lepton multiplicity before and after the cut, we should see the distribution reducing to a histogram including a single populated bin, centred on 1. Any event featuring 0 or at least 2 leptons should indeed get rejected by the selection.
This is achieved through the commands
ma5> plot N(l) 5 0 5 ma5> select N(l)==1 ma5> plot N(l) 5 0 5
that leads to the figure:
[Credits: @lemouth]
The behaviour observed is that expected. Everything is thus fine!
3 - Missing transverse energy spectrum
The final question asked in the previous episode of our project concerned the the missing transverse energy distribution. I mentioned somewhere in that blog that the MET keyword was associated with this quantity. It corresponds to the total amount of energy carried away by the invisible particles produced in a (simulated) collision.
Here, we considered the production of two top quarks, one of them decaying into 1 b-jet and two lighter jets, and the other decaying into 1 b-jet, 1 lepton and a neutrino. The neutrino is an invisible particle and therefore leaves the detector… undetected. We should thus have a significant amount of missing energy in the events.
This can be tested by typing in the command line interface of the code:
ma5> plot MET 50 0 350 [logY]
I have required a histogram of 50 bins ranging from 0 GeV to 350 GeV, the bounds allowing to have not too many events populating the overflow bin (0.5%) and a reasonable bin size related to the number of generated events (so that we are not too sensitive to statistical fluctuations). The results read:
[Credits: @lemouth]
As expected, we observed a peak at about 50 GeV, which is what could be expected from neutrinos typically produced in LHC collisions.
Summary: deciphering top pair production at the LHC
As promised, I finally released the solutions to the assignments I proposed in the fourth episode of our citizen science project on Hive. This post hence contains what should be done to solve the three questions raised in the assignment that I proposed 3 weeks ago.
I apologise again in the delay of releasing this. I just under-estimated how life could keep us busy… Please be ready for episode 5, that I hope to release very soon (I first need to catch up with a few French adaptations of the previous physics blogs before doing so).
In the meantime, feel free to come back to me if needed, or join the effort! It is never too late!
You are actually moving ahead with this great initiative in HIVE to show these exercises. One question, which software are you using to perform these plots ? I know that it isn't matlab or R for example.
!1UP
We are using MadAnalysis 5, a high-energy physics package dedicated to the analysis of simulated collisions (see Ep. 3). This code takes care of not only the analysis of the simulated collisions, but also of handling the plots.
For the figures shown in the present blog, MadAnalysis 5 relies on matplotlib. Note that ROOT can also be optionally used (this is what @travelingmercies uses).
PS: You can join us on this citizen science adventure on Hive if you want to try this program ;)
Greetings dear @lemouth, excellent work has been done by you and the participants of this project, who are always active in developing the assigned activities.
See you later, have a great week.
Thanks for passing by and the support.
Although we have lost a few participants (note that I am trying to recover them ;) ), we still have a nice number of them on board. This makes me quite happy :)
Hi @lemouth! I am alive, but I have not been logging in here and unavailable for some personal issues that more or less are getting solved now! So... now I have to use this weekend to chatch up with the previous and (hopefully) this assignment.
Glad to read you are fine. I was really wondering whether everything was alright. I hope your personal life is getting better, and I am looking forward to read something from you ;)
Thanks for the mention. I'll do this as soon as I can.
No problem. I am patiently waiting for you ;)
The rewards earned on this comment will go directly to the people(@lemouth) sharing the post on Twitter as long as they are registered with @poshtoken. Sign up at https://hiveposh.com.
Congratulations @lemouth! You have completed the following achievement on the Hive blockchain and have been rewarded with new badge(s):
Your next target is to reach 13500 replies.
You can view your badges on your board and compare yourself to others in the Ranking
If you no longer want to receive notifications, reply to this comment with the word
STOPCheck out the last post from @hivebuzz:
Thanks for your contribution to the STEMsocial community. Feel free to join us on discord to get to know the rest of us!
Please consider delegating to the @stemsocial account (85% of the curation rewards are returned).
You may also include @stemsocial as a beneficiary of the rewards of this post to get a stronger support.
Your content has been voted as a part of Encouragement program. Keep up the good work!
Use Ecency daily to boost your growth on platform!
Support Ecency
Vote for new Proposal
Delegate HP and earn more
You have received a 1UP from @gwajnberg!
@stem-curator, @vyb-curator, @pob-curator, @neoxag-curatorAnd they will bring !PIZZA 🍕
Learn more about our delegation service to earn daily rewards. Join the family on Discord.
Congratulations @lemouth! Your post has been a top performer on the Hive blockchain and you have been rewarded with the following badge:
You can view your badges on your board and compare yourself to others in the Ranking
If you no longer want to receive notifications, reply to this comment with the word
STOPCheck out the last post from @hivebuzz: