[Rivet] numerical difference between data and mc aida

Albert Knutsson albert.knutsson at cern.ch
Wed Apr 3 10:39:19 BST 2013


Dear Rivet experts,

we have a numerical problem with an aida file.

The data aida file have points looking like this:

   <dataPoint>
     <measurement value="5.000000e-02" errorPlus="1.666667e-02" 
errorMinus="1.666667e-02"/>
     <measurement value="1.186100e+00" errorPlus="1.535834e-01" 
errorMinus="1.535834e-01"/>
   </dataPoint>
   <dataPoint>
     <measurement value="8.333333e-02" errorPlus="1.666667e-02" 
errorMinus="1.666667e-02"/>
     <measurement value="1.763658e+00" errorPlus="2.291718e-01" 
errorMinus="2.291718e-01"/>
   </dataPoint>


The corresponding points in the MC aida produced by rivet, looks like this:

    <dataPoint>
     <measurement value="5.000001e-02" errorPlus="1.666666e-02" 
errorMinus="1.666666e-02"/>
     <measurement value="1.281114e+00" errorPlus="2.092222e-02" 
errorMinus="2.092222e-02"/>
   </dataPoint>
   <dataPoint>
     <measurement value="8.333334e-02" errorPlus="1.666667e-02" 
errorMinus="1.666667e-02"/>
     <measurement value="1.924731e+00" errorPlus="2.560152e-02" 
errorMinus="2.560152e-02"/>
   </dataPoint>


The consequence is that rivet-mkhtml complains about
"+++ Error in Histogram.divide(): Binning of histograms differs"

The ratio plots are still there, but there is a discontinuity in the 
distribution at the second point above.

Do you have a solution to this?

Thanks in advance!

Cheers,
Sunil and Albert


More information about the Rivet mailing list