[Rivet] Histogram normalisation

Frank Siegert frank.siegert at durham.ac.uk
Thu Oct 22 22:16:12 BST 2009


Frank Siegert, Tuesday 20 October 2009:
> Andy Buckley, Tuesday 20 October 2009:
> > One further thing, which I'm not sure counts in your "dropped"
> > KFactor proropsal: I don't see how we can automate the finalize steps
> > without always getting ~half of the observables very wrong,
>
> Can you elaborate? Why is it wrong if the finalize for each histogram
> automatically does
>
>   scale(_hist, crossSection()/sumOfWeights())
>
> ?
> The finalize method will then only consist of setting Scale/Norm for
> the histograms where that is necessary, e.g.:
>
>   _hist1->SetScale(1.0/inclusive_xs);
>   _hist2->SetNorm(data_integral);

As an addendum to this: We don't even have to wait for YODA to start this 
transition. Hendrik has improved the make-plots script such that it can 
read a global normalisation setting. So it's already possible to use 
scale(hist, crossSection()/sumOfWeights()) in finalize() together with 
NormalizeToIntegral=1 in *.plot as a replacement of normalize(hist).

I have started to migrate some analyses (in particular the ones where I 
was desperately waiting for the ability to merge runs).

Once we switch more analyses Professor will have to be adapted to make use 
of the normalisations in *.plot in case there are any.

Cheers,
Frank



More information about the Rivet mailing list