Register to reply

Mass invarint of ttbar pair by Pythia

by Safinaz
Tags: invarint, mass, pair, pythia, ttbar
Share this thread:
Sep9-13, 02:06 PM
P: 57
Hi all,

I'm not so expert in C++, but I try to make a a simple program to calculate
and plot the ttbar pair mass invariant distributtion by pythia in pp collition
(pp -> ttbar), taking the events from .lhe file in pythia exampls.

Now about the program attached here , i made the steps as following but
nothing shown in the Histogram .. may be there is some thing wrong in the loop ..
the steps :

- booked the Histogram
- start the event loop
- defined t , tbar
-start the particle loop
- make a contour to find the pairs
- calculate the invM.

So any help about that ? to know how to make such program
if any one knows about pythia ..

Thank you,
Attached Files
File Type: txt Mttbar.txt (1.8 KB, 11 views)
Phys.Org News Partner Physics news on
Engineers develop new sensor to detect tiny individual nanoparticles
Tiny particles have big potential in debate over nuclear proliferation
Ray tracing and beyond
Sep9-13, 03:52 PM
P: 11,928
for (int iEvent = 0; iEvent<100; ++iEvent)
I would expect {}-brackets for this loop.
If that is not the issue, add debug output in the loop (This should always be the first step if you cannot spot the issue). Do you get something meaningful from pythia? Do you get something meaningful for the histogram entry?

The indentation is really messy.

if (N + n > 1) {
This does not check for pairs, two b without a bbar pass this test as well (okay, that is rare).

What happened to the W boson from the top decay?
Sep10-13, 10:07 AM
P: 57

The ask about W boson, i think you mean we should to take into account the decay of the first generation
particles, so i defined before the particle loop a variable "iTop " and inside the loop
" if (pythia.event[i].id() == 6) iTop = i;" , i think that should work,

If : if (N+n >1) { --} dosn't account for the pair , what could be the probable loop ? or how to calculate the mass of these pairs ?


Sep10-13, 12:08 PM
P: 11,928
Mass invarint of ttbar pair by Pythia

Your topic is "invariant mass of ttbar pair", but you calculate a mass based on b-quarks. As far as I can see, you are missing the energy and momentum of the two W.

"if (N+n >1)" -> "if (N >0 && n>0)"
Just check if both a b and a bbar are present in the event.
Sep10-13, 01:08 PM
P: 57

Any way thank you for your comments, i may use another tool to do this calculation until i know better in c++.


Register to reply

Related Discussions
Geant4 Or Pythia ? Math & Science Software 2
Pair Production and Mass Energy Equivalence Quantum Physics 14
About pythia High Energy, Nuclear, Particle Physics 4
Pythia and differential cross sections High Energy, Nuclear, Particle Physics 5
Partial and Covariant derivatives in invarint actions General Math 0