Additional histograms can be created easily in VBFNLO by modifying the file utilities/histograms.F
.
InitHistograms
create a new histogram with the following command: call CreateHist(HIST_ID, "title",nobins,xmin,xmax)
HIST_ID
= integer identification of histogram (make sure to use a different number for each histogram) title
= title of histogram, will be used as plot title – e.g. dS/dmjj : dijet massnobins
= number of histogram binsxmin
= lower bound of histogramxmax
= upper bound of histogramHistogramEvent
call FillHist(HIST_ID, value, dw, NLO)
HIST_ID
= integer identification of histogram, same number as in CreateHistvalue
= value of your distribution variable (e.g. mjj)dw
= the variable dw (this is already calculated by the code – it is the weight*cross section)NLO
= the variable NLO (this is already filled by the code – flag determining the order of the calculation)
Information about the final state particles is stored in the arrays jets
, leptons
, invisible
(neutrinos) and photons
. The entries are as follows
jets(i,j)
: information about the jets is stored here, entries are $p_{T}$-orderedleptons(i,j)
: information about the charged leptons is stored hereinvisible(i,j)
: information about the invisible particles (i.e. neutrinos) is stored herephotons(i,j)
: information about the photons is stored here