[Rivet] Normalization?

[Salvatore Rappoccio]

Hi, Andy,

Aha, that worked great. Thanks!

Cheers,
Sal


On Thu, Mar 15, 2018 at 5:13 PM, Andy Buckley <andy.buckley at cern.ch> wrote:

> Hi Sal,
>
> I suspect we're going to iterate a little until I understand exactly what
> you want!
>
> First, I assume the N_PT_BINS thing is unrelated to the bins in your
> histogram, i.e. there's a histo per pT bin?
>
> So for each histo you effectively call normalize(histo, 1000), which will
> normalize that histogram to have an area of 1000 (including the overflows).
> This is working, since you can see the "Area" info printout in the file
> says 1000.
>
> So what do you want? The normalization before normalize() was called? You
> can get this by rescaling by 1/h.annotation("ScaledBy") -- that attribute
> records how much scaling has been done, so you can merge parallel runs of
> normalised observables together. [And we've got a more sophisticated riff
> on that idea in the pipeline... it'll emerge into daylight some day, I
> promise.]
>
> Does that help? I feel I probably missed the point, so enlighten me!
>
> Andy
>
> *Dr Andy Buckley, Lecturer / Royal Society University Research Fellow*
> Particle Physics Experiment Group, University of Glasgow
> On Mar 15 2018, at 9:01 pm, Salvatore Rappoccio <rappoccio at gmail.com>
> wrote:
>
>
> Hi, Folks,
>
> I'm a bit confused about the normalization of Yoda files from RIVET. I
> have this at the end of my module, where I normalize a bunch of histograms
> with unequal bin widths:
>
>     /// Normalise histograms etc., after the run
>     void finalize() {
>       const double normalizationVal = 1000;
>       for (size_t i = 0; i < N_PT_BINS_dj; ++i) {
>         normalize(_h_ungroomedJetMass_dj[i], normalizationVal);
>         normalize(_h_sdJetMass_dj[i], normalizationVal);
>       }
>     }
>
> Then I get some yoda output as below (*). How can I get the absolute
> normalization before any normalization from Rivet? Is that somehow in
> "ScaledBy" or perhaps "ScaledBy/Area" ??
>
> Cheers,
> Sal
>
>
>
> BEGIN YODA_HISTO1D /CMS_SMP_16_010/d09-x01-y01
> Path=/CMS_SMP_16_010/d09-x01-y01
> ScaledBy=2.1193569211303091e-06
> Title=
> Type=Histo1D
> XLabel=
> YLabel=
> # Mean: 1.518686e+02
> # Area: 1.000000e+03
> # ID ID sumw sumw2 sumwx sumwx2 numEntries
> Total   Total   1.000000e+03 9.342226e+01 1.518686e+05 3.082456e+07 10718
> Underflow Underflow 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0
> Overflow Overflow 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0
> # xlow xhigh sumw sumw2 sumwx sumwx2 numEntries
> 0.000000e+00 1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
> 0.000000e+00 0
> 1.000000e+00 5.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
> 0.000000e+00 0
> 5.000000e+00 1.000000e+01 0.000000e+00 0.000000e+00 0.000000e+00
> 0.000000e+00 0
> 1.000000e+01 2.000000e+01 9.119381e-02 8.316311e-03 1.729448e+00
> 3.279817e+01 1
> 2.000000e+01 4.000000e+01 1.060794e+01 9.968554e-01 3.627352e+02
> 1.265246e+04 113
> 4.000000e+01 6.000000e+01 6.498255e+01 6.082119e+00 3.357938e+03
> 1.753501e+05 695
> 6.000000e+01 8.000000e+01 1.255057e+02 1.173247e+01 8.858834e+03
> 6.293726e+05 1344
> 8.000000e+01 1.000000e+02 1.435704e+02 1.341602e+01 1.288541e+04
> 1.161209e+06 1538
> 1.000000e+02 1.500000e+02 2.662677e+02 2.487217e+01 3.269421e+04
> 4.068761e+06 2854
> 1.500000e+02 2.000000e+02 1.520766e+02 1.419840e+01 2.636244e+04
> 4.601090e+06 1631
> 2.000000e+02 2.500000e+02 9.095334e+01 8.469162e+00 2.028071e+04
> 4.540426e+06 978
> 2.500000e+02 3.000000e+02 6.430813e+01 6.001638e+00 1.763051e+04
> 4.846869e+06 690
> 3.000000e+02 3.500000e+02 4.010750e+01 3.760009e+00 1.294571e+04
> 4.186287e+06 429
> 3.500000e+02 4.000000e+02 2.480605e+01 2.329297e+00 9.231820e+03
> 3.440384e+06 265
> 4.000000e+02 4.500000e+02 1.311072e+01 1.220962e+00 5.534208e+03
> 2.338449e+06 141
> 4.500000e+02 5.000000e+02 3.152555e+00 2.925908e-01 1.479963e+03
> 6.954321e+05 34
> 5.000000e+02 5.500000e+02 3.687638e-01 3.400377e-02 1.891971e+02
> 9.710031e+04 4
> 5.500000e+02 6.000000e+02 9.075043e-02 8.235644e-03 5.316258e+01
> 3.114319e+04 1
> _______________________________________________
> Rivet mailing list
> Rivet at projects.hepforge.org
> https://www.hepforge.org/lists/listinfo/rivet
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20180315/9ca857a9/attachment.html>