The computation of the errors for profile histograms is
wrong in some cases. The following macro shows the problem:
v/create xv(3) r 0.2 0.4 1.5
v/create yv(3) r 0.8 1.2 1.1
v/write xv,yv test.dat '2(1x,f5.2)'
nt/create 200 'Test ntuple' 2 ! ! x y
nt/read 200 test.dat '2(1x,f5.2)'
profile 201 'Test profile histogram 1' 2 0. 2. 0. 2.
nt/pro 201 200.y%x
v/create con1(2) r
v/create err1(2) r
h/get_vect/con 201 con1
h/get_vect/err 201 err1
v/print con1
v/print err1
The output should be:
CON1(1) = 1
CON1(2) = 1.1
ERR1(1) = 0.141421
ERR1(2) = 1.04881
and not:
CON1(1) = 1
CON1(2) = 1.1
ERR1(1) = 0.141421
ERR1(2) = 0.000119604
|