[Rivet] Rivet Jet Analysis

[Hendrik Hoeth]

Hi Sam,

Thus spake Andy Buckley (andy.buckley at ed.ac.uk):

> > For the CDF_2008_LEADINGJETS  <Nch> against pT(leading jet) plots, the
> > runs from each pT cut are needed to fill a single histogram. Should I
> > deal with this by chopping each histogram down to its valid region (say
> > 0-20 pT for the minimum bias run, then 20-100pT for the pTmin=10GeV run
> > etc) then stick all these together to make a complete distribution? Or
> > is there a smarter way of doing this.

Currently we do exactly that. There is a script that can do the
histogram slicing and merging for you: rivet-mergeruns. It's a bit
awkward to use, since it's not exactly aimed at end-users, but
nevertheless. Here is how I used it in a recent Pythia-8 production:

rivet-mergeruns CDF_2008_LEADINGJETS_CKIN0/Rivet.aida:1960:0 \
                CDF_2008_LEADINGJETS_CKIN10/Rivet.aida:1960:30 \
                CDF_2008_LEADINGJETS_CKIN20/Rivet.aida:1960:50 \
                CDF_2008_LEADINGJETS_CKIN50/Rivet.aida:1960:80 \
                CDF_2008_LEADINGJETS_CKIN100/Rivet.aida:1960:130 \
                CDF_2008_LEADINGJETS_CKIN150/Rivet.aida:1960:180 \
                > foo.dat
flat2aida foo.dat

You have to specify <aidafile>:<energy>:<mergingpT>. Note that you need
to set the histogram merging pT a fair bit higher than the generator
level cut!

HTH,

     Hendrik

-- 
It pays to be obvious, especially if you have
a reputation for subtlety.    -- Isaac Asimov