Re: flair 1.0-1 and plotting 1-D projections of usrbin distributions - Is this a bug?

Hello everyone,

I set this issue aside while tending to other things but have had to come back to it.
I definitely see an odd behaviour when plotting 1-D projections using flair 1.0-1, even when I use plotting styles not
involving error bars, i.e. "points", "lines", "linespoints", etc.

Wed Nov 14 11:17:13 2012 Plotting: TR24_Mo100_DoseEq_Z
Starting gnuplot
set mouse
show version
Gnuplot version = 4.2
show terminal
Default terminal = x11
bin(x,m,s)=floor((x-m)/s)*s+m
logbin(x,m,s)=10**(floor((log10(x)-m)/s)*s+m)
rnd(x)=abs(int((x*1103515245+12345)/65536)%32768)
set terminal x11 3 enhanced raise title 'TR24.flair-3: TR24 - Mo100 - Dose Eq (uSv/h) - for a 2 meter slice in y around
the beam'
reset

set title 'TR24 - Mo100 - Dose Eq (uSv/h) - for a 2 meter slice in y around the beam'
set xlabel 'Z (cm)'
set ylabel 'Dose Eq (uSv/h)'
set grid
set key default
set cblabel ''
set style line 1 lt -1 lw 1
set logscale y
plot 'TR24_Mo100_DoseEq_Z.dat' us 1:($3*1.12347e+13) w points lt 1 lw 1 lc rgb 'black' pt 1 ps 1 notitle <=======

1) I no longer see the command used to generate TR24_Mo100_DoseEq_Z.dat (the data file used as input to gnuplot)
2) In the above gnuplopt command, the errorbars are plotted as a function of z and not the quantity scored because of
the $3.

3) If I choose to plot with errorbars, this is the command that is passed to gnuplot within flair:
plot 'TR24_Mo100_DoseEq_Z.dat' us 1:($3*1.12347e+13):(($3*$4/100)*1.12347e+13) w errorbars lt 1 lw 1 lc rgb 'black' pt 1
ps 1 notitle

And that's why I get the error: "line 0: x range is invalid"

When I use gnuplot directly (not via flair) with the corrected commands (i.e. the right column number used to plot the
quantity and error bars), things look fine.

plot 'TR24_Mo100_DoseEq_Z.dat' us 1:($2*1.12347e+13) w points lt 1 lw 1 lc rgb 'black' pt 1 ps 1 notitle
plot 'TR24_Mo100_DoseEq_Z.dat' us 1:($2*1.12347e+13):(($2*$3/100)*1.12347e+13) w errorbars lt 1 lw 1 lc rgb 'black' pt 1
ps 1 notitle

Does anyone know why this is happening within flair?

Thank you in advance,
Mina

On 12-11-01 01:08 PM, Mina Nozar wrote:
> Hello everyone,
>
> I am having a problem plotting 1-D projections of 2-D plots generated from a USRBIN estimator with errorbars, so
> basically all options under 'Plotting Style' that involves errorbars. I am not sure why this is.
>
> The error message I get is: " <e< line 0: x range is invalid "
>
>
> This is what the '.dat' file looks like:
> -970. 3.60040433E-16 75.3390884
> -960. 4.82170911E-16 72.745163
> -950. 4.82170911E-16 70.335434
> -940. 4.82170911E-16 100.
> -930. 4.82170911E-16 72.6674881
> -920. 4.82170911E-16 77.4434357
> -910. 4.82170911E-16 88.6815872
> -900. 4.82170911E-16 82.6831665
> -890. 4.82170911E-16 70.9556732
> -880. 4.98189045E-16 67.284584
> -870. 2.63280159E-15 22.167038
> -860. 9.28354294E-15 13.8853865
> -850. 2.48675507E-14 11.0314779
> -840. 5.83361687E-14 8.37456894
> -830. 1.20444832E-13 7.85978127
> -820. 3.04576839E-13 5.7727437
> -810. 8.30334716E-13 4.78959179
> -800. 2.10400078E-12 3.74105716
> -790. 4.89483714E-12 2.79912043
> -780. 1.29408264E-11 2.53790736
> -770. 3.35692758E-11 1.88130307
> -760. 8.68509709E-11 1.61146486
> -750. 2.13434118E-10 1.17526901
> -740. 5.33403655E-10 0.886002541
> -730. 1.33893463E-09 0.795381248
> -720. 3.36977357E-09 0.614626765
> -710. 8.44234993E-09 0.562169731
> -700. 2.11343867E-08 0.484024316
> etc...
>
> So, nothing odd about it. Is anyone else having this issue?
>
> Thanks a lot and best wishes,
> Mina
>
>

-- 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Mina Nozar, PhD
Scientific Computing Analyst
TRIUMF
4004 Wesbrook Mall
Vancouver, BC, V6T 2A3
http://www.triumf.ca
nozarm_at_triumf.ca
(604)222-1047 ext 6005