[Rivet] Comparing histograms and ratios of certain histograms (2D scatter plots) with data

Andy Buckley andy.buckley at cern.ch
Mon Jul 4 16:57:46 BST 2016


Hi Raghav,

I'm not really sure what the question is: the ratio "binning" is that of 
the numerator & denominator histograms passed to divide(). [For now, at 
least, the num/denom binnings have to match exactly, although one could 
work harder and do bin-merging to the compatible subset.]

If not filling it from a call to divide(), you can also call 
bookScatter2D() with an argument to pre-populate the scatter with points 
whose x-values and x-errors match the binning in the provided reference 
data. Check the Doxygen documentation for the details.

Andy


On 01/07/16 13:57, Raghav Kunnawalkam Elayavalli wrote:
> Hi Rivet Experts,
>
> I want to write an analysis compares cross section of different
> objects and their ratios with data. Say I have 4 cross sections, A,
> B, C and D defined as Histo1D and I make a set of ratios which are
> B/A, C/A and D/A and i have data points for A, B, C, D, B/A, C/A and
> D/A. Now in rivet, these ratios will be done during the finalize
> method where i create scatter2D plots for the ratios BRatio, CRatio
> and DRatio. divide(A, B, BRatio); Now when i write this analysis, i
> plan on taking the bins for the histogram from the data yoda file A =
> bookHisto1D("name", refData(1, 1, 1)) something like that yeah. My
> question is how do i declare the ratio plots in the analysis? can
> they also take the binning from the data? and another technical
> question is all these A,B,C,D etc... are all arrays of histograms. so
> when i declare them, im doing something like this inside the loops.
> stringstream jetpt, raa; jetpt << "Jet_pT_radbin" << r << "_centbin"
> << i; _h_jetpT[r][i] = bookHisto1D(jetpt, _ptedges);
>
> Cheers Raghav _______________________________________________ Rivet
> mailing list Rivet at projects.hepforge.org
> https://www.hepforge.org/lists/listinfo/rivet
>


-- 
Dr Andy Buckley, Lecturer / Royal Society University Research Fellow
Particle Physics Expt Group, University of Glasgow


More information about the Rivet mailing list