ferella / polepp

Automatically exported from code.google.com/p/polepp
0 stars 0 forks source link

polelim example results in GSL ERROR #4

Closed GoogleCodeExporter closed 8 years ago

GoogleCodeExporter commented 8 years ago
I tried running the polelim example in the README file, which resulted in the 
following error:

$ ./polelim --effdist 0 --bkgdist 2 --bkgmeas 2 --bkgsigma 0.5 --nobs 4
Tools::calcIntRange() -> Unknown pdf type = 999

================ P O L E ==================
 1.0 - conf. level  : 0.1
 N observed         : 4
----------------------------------------------
 Coverage friendly  : No
 True signal        : 1
----------------------------------------------
 Efficiency meas    : 1
 Efficiency sigma   : 0.2
 Efficiency dist    : Undefined
 Efficiency scale   : 1
----------------------------------------------
 Background meas    : 2
 Background sigma   : 0.5
 Background dist    : Gauss
 Background scale   : 1
----------------------------------------------
 Bkg-Eff correlation: 0
----------------------------------------------
 GSL int. N(calls)  : 100
----------------------------------------------
 Int. eff. min      : 0
           max      : 5
----------------------------------------------
 Int. bkg. min      : 0
           max      : 5
----------------------------------------------
 Binary search thr. : 0.0001
 1-CL threshold     : 0.0001
 Min prob in belt   : 0.001
----------------------------------------------
 *Test hyp. min     : 0
 *Test hyp. max     : 35
 *Test hyp. step    : 0.01
----------------------------------------------
 Step mu_best       : 0.01
 Max N, mu_best     : 20
----------------------------------------------
 Method             : FHC2
----------------------------------------------
 Verbosity          : 0
----------------------------------------------
 Parameters prefixed with a * above are not   
 relevant for limit calculations.
==============================================
gsl: vegas.c:128: ERROR: xu must be greater than xl
Default GSL error handler invoked.
Abort trap

For some of the other examples I get slightly different output, in some cases 
because default values have changed, 

For one example I get a slightly different upper limit:
./polelim --nobs 30
result in  [ 17.812500, 47.708691 ]
whereas the README file states [ 17.81, 47.66 ]

Original issue reported on code.google.com by Deil.Chr...@googlemail.com on 1 Jul 2011 at 3:56

GoogleCodeExporter commented 8 years ago
Hello,
1. could you provide me with the exact output from ./polelim --nobs 30?

2. the first call fails since you have effdist 0 which is 'NONE'. As far as I 
remember this should be equivalent to the case when the efficiency is exactly 
known quantity.
It seems however that this functionality is broken. Need to look into that.

3. do you have details on system/compiler etc?

Thanks for your comments!

Original comment by teg...@gmail.com on 4 Jul 2011 at 8:13

GoogleCodeExporter commented 8 years ago

Original comment by teg...@gmail.com on 5 Jul 2011 at 7:04

GoogleCodeExporter commented 8 years ago
Fixed the issue with crashing code when the efficiency is set to be constant 
(--effdist 0).
Still need to trace down the cause for the difference in limits as observed.

Original comment by teg...@gmail.com on 8 Jul 2011 at 9:36