Mass invarint of ttbar pair by Pythia


by Safinaz
Tags: invarint, mass, pair, pythia, ttbar
Safinaz
Safinaz is offline
#1
Sep9-13, 02:06 PM
P: 36
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,
Safinaz
Attached Files
File Type: txt Mttbar.txt (1.8 KB, 11 views)
Phys.Org News Partner Physics news on Phys.org
Physicists design quantum switches which can be activated by single photons
'Dressed' laser aimed at clouds may be key to inducing rain, lightning
Higher-order nonlinear optical processes observed using the SACLA X-ray free-electron laser
mfb
mfb is offline
#2
Sep9-13, 03:52 PM
Mentor
P: 10,808
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?
Safinaz
Safinaz is offline
#3
Sep10-13, 10:07 AM
P: 36
Hi,

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 ?

S.

mfb
mfb is offline
#4
Sep10-13, 12:08 PM
Mentor
P: 10,808

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.
Safinaz
Safinaz is offline
#5
Sep10-13, 01:08 PM
P: 36
Hi,

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

Best,
S.


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