|
[Rivet] Comparing histograms and ratios of certain histograms (2D scatter plots) with dataAndy Buckley andy.buckley at cern.chMon 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 |