Mass invarint of ttbar pair by Pythia

  • Thread starter Safinaz
  • Start date
  • #1
219
2

Main Question or Discussion Point

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
 

Attachments

Answers and Replies

  • #2
34,044
9,892
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?
 
  • #3
219
2
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.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.
 
  • #4
34,044
9,892
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.
 
  • #5
219
2
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.
 

Related Threads on Mass invarint of ttbar pair by Pythia

  • Last Post
Replies
4
Views
3K
Replies
5
Views
7K
Replies
5
Views
975
Replies
12
Views
1K
Replies
5
Views
1K
Replies
15
Views
2K
Replies
5
Views
4K
Replies
1
Views
2K
Top