Analysis Implementation Questions

[Gianco99]

Hello all,

I am implementing my first MadAnalysis script and have found some success in defining my selections and producing my desired distributions. However, I had a few questions on how to streamline my analysis workflow to simplify a few things.

One of the main variables I use in the analysis is the number of constituents of the FatJet object. I was able to find the piece of code in the DelphesTreeReader.cpp that reads and defines the fatjet properties. I had to modify the variable for the constituents to read as it is automatically set to 0, i.e. to jet->ntracks_ = part->NCharged + part->NNeutrals; Is it possible for this to be merged into the codebase (for both the jets and fatjets)? Was there a reason it was being manually set to 0?

In addition, I had to write the code implementation to plot the number of fatjet constituents by hand as there is no particle defined for the fatjet. For example, in the reco-level, we can simply use "j" to refer to the standard jet collection and write plot NTRACKS(j[1]) 40 0 60. How would I introduce a new particle definition "FatJ" such that, for example, plot NTRACKS(Fatj[1]) 40 0 60, can be used to find the number of particles in the fatjet collection read from the Delphes output?

After modifying the user.cpp script by-hand with some additional selections or plot for the fatjet properties, is there a way to recompile the PDF output? I was able to get the SAF histograms in the XML format and plot these using my own script I wrote by hand, but I was wondering if the PDF Latex step could be run again after recompiling the user.cpp script with some manual changes?

Finally, this may be a trivial question, but is there a way to make 2D plots? For example, if I want to plot NTracks vs Jet pT in 2D?

Thanks in advance for any insight!