Accessing jet energy fractions?

Asked by Billy Ford on 2019-11-04

Hi,

Just wondering if it is possible to access jet constituency information in MA5, either in the normal mode or by editing certain elements of the SampleAnalyzer code?

In particular we're trying to apply some cuts to b-jets based on the energy fractions of neutral/charged hadrons and leptons, as well as plotting these energy fractions as histograms for our analysis.

Best regards,
Billy

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Benjamin Fuks
Solved:
2019-11-06
Last query:
2019-11-06
Last reply:
2019-11-06
Benjamin Fuks (fuks) said : #1

Hi Billy,

> Just wondering if it is possible to access jet constituency information in MA5, either in the normal mode or by editing certain elements of the SampleAnalyzer code?
This is possible in the expert mode via the constituents() method of the RecJetFormat class. I think this should solve your issue. Please let me know.

Cheers,

Benjamin

Billy Ford (billyford) said : #2

Thanks for the reply Benjamin,

Are there any good starting points for learning how to use the expert mode to put together an analysis, I've given the manuals a read and am still unsure on how exactly to build up an analysis from scratch?

Cheers,
Billy

Best Benjamin Fuks (fuks) said : #3

Hi Billy,

Please have a look to
  https://arxiv.org/abs/1808.00480 <https://arxiv.org/abs/1808.00480>
This document should include all the information you need.

Cheers,

Benjamin

> On 5 Nov 2019, at 18:53 , Billy Ford <email address hidden> wrote:
>
> Question #685609 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/685609
>
> Billy Ford posted a new comment:
> Thanks for the reply Benjamin,
>
> Are there any good starting points for learning how to use the expert
> mode to put together an analysis, I've given the manuals a read and am
> still unsure on how exactly to build up an analysis from scratch?
>
> Cheers,
> Billy
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Billy Ford (billyford) said : #4

Thanks Benjamin Fuks, that solved my question.

Billy Ford (billyford) said : #5

Hi,

Sorry to dig this issue up again, I've implemented a first analysis in the expert mode and am trying now to apply the constituent jet energy cuts.

I wonder if you could shed some light on the constituents() method. Does this contain truth information on the type of each particle in the given jet? For reference I'm only using fastjet (no delphes), so am unsure exactly on how to grab each particle type from the jet constituents to sum the energy contribution. The contributors we're interesting are electrons, muons, photons+neutral hadrons, charged hadrons, it seems that charge is not recognised, and pdg won't work for the hadron groups.

Best,
Billy

Benjamin Fuks (fuks) said : #6

Hi Billy,

Maybe an example is better in this case (otherwise, let's continue discussing). Please have a look to
  tools/SampleAnalyzer/Process/JetClustering/TauTagger.cpp
There is an example of how jet constituents could be used around line 76.

I hope this helps.

Cheers,

Benjamin

> On 2 Dec 2019, at 05:23 , Billy Ford <email address hidden> wrote:
>
> Question #685609 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/685609
>
> Billy Ford posted a new comment:
> Hi,
>
> Sorry to dig this issue up again, I've implemented a first analysis in
> the expert mode and am trying now to apply the constituent jet energy
> cuts.
>
> I wonder if you could shed some light on the constituents() method. Does
> this contain truth information on the type of each particle in the given
> jet? For reference I'm only using fastjet (no delphes), so am unsure
> exactly on how to grab each particle type from the jet constituents to
> sum the energy contribution. The contributors we're interesting are
> electrons, muons, photons+neutral hadrons, charged hadrons, it seems
> that charge is not recognised, and pdg won't work for the hadron groups.
>
> Best,
> Billy
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Billy Ford (billyford) said : #7

Benjamin,

This looks to be just what I'm after, thanks a lot for the continued assistance!

Best,
Billy