Fork me on GitHub

Opened 6 years ago

Last modified 6 years ago

#1355 new Bug

The stransverse mass MT2

Reported by: Najimuddin Khan Owned by:
Priority: minor Milestone:
Component: Delphes code Version: Delphes 3
Keywords: Cc:

Description

Hi,
Is it possible to calculate the transverse mass MT2 using the standard variables in the Delphes program? Any suggestion !!!!

Thanks and regards,
Najim

Change History (5)

comment:1 by Pavel Demin, 6 years ago

I think that it could be calculated using the Mt2() method of the TLorentzVector class. Most of the Delphes classes have a method called P4() that returns an object of the TLorentzVector class.

For example, the following code should work for a jet:

  Double_t mt2 = jet->P4().Mt2()
Last edited 6 years ago by Pavel Demin (previous) (diff)

comment:2 by Najimuddin Khan, 6 years ago

Hi,
I think it (the stransverse mass MT2) is not the same as in REFs:
Chris Lester, David Summers Phys.Lett.B463:99-103,1999
Alan Barr, Christopher Lester, Phil Stephens J.Phys.G29:2343-2363,2003

Am I missing anything?

comment:3 by Pavel Demin, 6 years ago

Am I missing anything?

What is missing is a detailed description of what you're looking for.

comment:4 by Najimuddin Khan, 6 years ago

Hi,
Thank you so much for the reply.
I have two charged lepton (coming from W+W- or some BSM particles but not from the Z) plus missing ET in the final state.
I want to calculate the stransverse mass MT2 variable(see equation (8) of http://iopscience.iop.org/article/10.1088/0954-3899/29/10/304/pdf). The authors have used some maximization and minimization procedure using different codes.

Is this code available inside the Delphes or I have to introduce it by hand?
Is it same as you suggested in your earlier reply?: Double_t mt2 = jet->P4().Mt2()

Thank you
NK

comment:5 by Pavel Demin, 6 years ago

I don't have access to the full text linked article.

Anyway, if the variable that you're looking for is calculated from the kinematic parameters of several selected objects (lepton, missing ET, etc), then I think that it should be calculated in the code that analyses the Delphes output.

So, the short answer is that the code isn't available inside Delphes and you have to introduce it by hand.

Note: See TracTickets for help on using tickets.