[Rivet] Two-dimensional histograms

James Robinson james.robinson at cern.ch
Tue Sep 16 22:53:15 BST 2014


Dear Andy and David,

On 16 September 2014 20:58, Andy Buckley <andy.buckley at cern.ch> wrote:

> In fact, there is even an efficiency(pass, tot) function that takes 1D
> histos and returns a Scatter2D with appropriate binomial statistics
> treatment!
>
> Andy
>

​I think this is exactly what I want! The binomial error treatment is
important for this analysis as the efficiencies tend towards 1.0 in some
parts of phase space.​ What is the appropriate class called?


On 16/09/14 19:00, David Bjergaard wrote:
> > Hi James,
> >
> > Why can't you hold two histos:
> > hist_pass and hist_fail?
> >
> > Then the 1D (operator+= ) methods for adding and dividing would work
> (under the
> > assumptions imposed by dividing). Then you can do:
> > hist_denom+=hist_pass;
> > hist_denom+=hist_fail;
> > divide(hist_pass,hist_denom);
> >
> > I (personally) keep a std::map<std::string,Histo1DPtr> container indexed
> > by string for bookkeeping.
> >
> >     Dave
>

​This is also fine, but I was trying to minimise the number of extraneous
histograms that the code booked (that I would then need to remove).

As far as Andy's other questions/suggestions go. I'm using the ATLAS
packaged version of Rivet. I'm not sure which Rivet version this is, but
the athena version is 19.0.2.1. Perhaps you'd recommend a newer release,
Andy?

Thank you both for your helpful answers,
James
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20140916/8b2e60ce/attachment.html>


More information about the Rivet mailing list