[Rivet] YODA/Histo1D - possible division by xMid ?

Antonin Maire (CERN) antonin.maire at cern.ch
Tue Jul 24 16:52:21 BST 2018


Good afternoon everybody,

   I wanted to know whether there is any possibility to play with the 
individual height of each bin in Histo1D ?
So far looking at the code :
- https://yoda.hepforge.org/trac/browser/include/YODA/HistoBin1D.h
- https://yoda.hepforge.org/trac/browser/include/YODA/Histo1D.h

I haven't found such an option, I believe.

The only way I found would be to read the bin height of the histo, play 
with it and then set up a from that a Scatter2D plot
Like done in the Histo1D::divide methods : 
https://yoda.hepforge.org/trac/browser/src/Histo1D.cc

Do you confirm this is the only way ?

What I want to do :
- I would like to rivetise an analysis which is in the form of the 
productio pT-spectrum : 1/(Nevt 2pi pT) d²N/dpTdy.
Usually the way to do that in Rivet is to fill on-the-fly the Histo1D 
with weights :
... ->fill(pT, weight * 1/pT);

What I would like to do is to stick to data (and its possible bias) and 
really do : 1/(Nevt 2pi pTmid) d²N/dpTdy,
i.e. divide each bin content (root-language) or bin height 
(YODA-language) by its HistoBin1D::xMid().

Regards,
    Antonin









-------------- next part --------------
A non-text attachment was scrubbed...
Name: antonin_maire.vcf
Type: text/x-vcard
Size: 4 bytes
Desc: not available
URL: <https://www.hepforge.org/lists-archive/rivet/attachments/20180724/44753642/attachment.vcf>


More information about the Rivet mailing list