[Rivet] [PATCH] filter nan's in make-plots

Chris Pollard cpollard at cern.ch
Sun Aug 13 19:30:07 BST 2017


Hi,

Is this the behavior we want in make-plots? I'm not sure if it's better to
throw an error in NaN-valued bins or to just plot nothing, but by default I
prefer the former.

Chris

On Wed, Jul 12, 2017 at 8:09 AM, Dmitry Kalinkin <dmitry.kalinkin at gmail.com>
wrote:

> Hi,
>
> There is a bug in make-plots which prevents it from generating correct
> latex when I try to compare TProfile1D’s with occasional empty bins. Patch
> below fixes the problem for me. Could you please integrate it.
>
> Thanks,
>
> Best,
> Dmitry
>
> diff --git a/bin/make-plots b/bin/make-plots
> index 4cbdfacf..ecaf4dbe 100755
> --- a/bin/make-plots
> +++ b/bin/make-plots
> @@ -1810,6 +1810,9 @@ class Histogram(DrawableObject, Described):
>                  if self.getErrorBands():
>                      self.description['SmoothLine'] = 0
>                      for b in self.data:
> +                        if isnan(b.val) or isnan(b.err[0]) or
> isnan(b.err[1]):
> +                            seen_nan = True
> +                            continue
>                          out += ('\\psframe[dimen=inner,linewi
> dth=0pt,linestyle=none,fillstyle=solid,fillcolor=%s,opacity=%s]' %
> (self.getErrorBandColor(),self.getErrorBandOpacity()))
>                          out += ('(' + coors.strphys2frameX(b.xmin) + ', '
> \
>                                      + coors.strphys2frameY(b.val -
> b.err[0]) + ')(' \
>
> _______________________________________________
> 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/20170813/91337f7b/attachment.html>


More information about the Rivet mailing list