Hi,
Could anybody tell me why errors are identical in all
bins in TProfile when that profile is filled with small values?
The problem happens in ROOT 3.02.07, If I repeat the same thing
in ROOT 3.02.00, everything looks fine.
The code is :
==================test.C=====================
#include "TProfile.h"
#include "float.h"
#include "TRandom.h"
#include "TString.h"
#include "TCanvas.h"
#include "iostream.h"
void test(){
TRandom rd;
TProfile* hprof= new TProfile("test", "test", 100,-4., 4., -1.*FLT_MAX,
FLT_MAX, "");
Float_t px, py, pz;
for ( Int_t j=0; j<25; j++) {
rd.Rannor(px,py);
pz = px*px + py*py;
hprof->Fill(px,pz*1e-7); <======== errors are calc. correctly if
*1e-7 is removed from this line. It can be also calc. correctly if run in ROOT 3.02.00
}
hprof->Print("all");
TCanvas *c1 = new TCanvas("c1","Profile histogram example",200,10,700,500);
hprof->Draw();
}
========================================
Errors are shown the same for all bins in the following output:
----------------------------------------
*******************************************
* *
* W E L C O M E to R O O T *
* *
* Version 3.02/07 1 March 2002 *
* *
* You are welcome to visit our Web site *
* http://root.cern.ch *
* *
*******************************************
Compiled for linux with thread support.
CINT/ROOT C/C++ Interpreter version 5.15.25, Jan 6 2002
Type ? for help. Commands must be C++ statements.
Enclose multiple statements between { }.
Welcome to the ROOT tutorials
Type ".x demos.C" to get a toolbar from which to execute the demos
Type ".x demoshelp.C" to see the help window
*** Start at Date : Tue Jul 16 15:36:50 2002
QAInfo:You are using STAR_LEVEL : pro, ROOT_LEVEL : 3.02.07 and node :
rcas6015.rcf.bnl.gov
root4star [0] .x test.C++
Info in <ACLiC>: creating shared library
/direct/star+data10/GC/aihong/anaSimV1_short/V2_5Pct_V1Vary_New/job4/dir8Test/./test_C.so
In file included from
/direct/star+data10/GC/aihong/anaSimV1_short/V2_5Pct_V1Vary_New/job4/dir8Test/filexIyM7Z.h:29,
from
/direct/star+data10/GC/aihong/anaSimV1_short/V2_5Pct_V1Vary_New/job4/dir8Test/./filexIyM7Z.C:5:
/direct/star+data10/GC/aihong/anaSimV1_short/V2_5Pct_V1Vary_New/job4/dir8Test/./test.C:
In function
`void test ()':
/direct/star+data10/GC/aihong/anaSimV1_short/V2_5Pct_V1Vary_New/job4/dir8Test/./test.C:25:
warning: unused
variable `TCanvas *c1'
TH1.Print Name= test, Entries= 25, Total sum= 2.96441e-06
fSumw[0]=0, x=-4.04, error=0
//bla... bla...
fSumw[32]=2.96728e-07, x=-1.48, error=5.04825e-08
fSumw[33]=0, x=-1.4, error=0
fSumw[34]=0, x=-1.32, error=0
fSumw[35]=1.60949e-07, x=-1.24, error=5.04825e-08
fSumw[36]=0, x=-1.16, error=0
fSumw[37]=3.59341e-07, x=-1.08, error=5.04825e-08
fSumw[38]=0, x=-1, error=0
fSumw[39]=0, x=-0.92, error=0
fSumw[40]=0, x=-0.84, error=0
fSumw[41]=0, x=-0.76, error=0
fSumw[42]=5.56881e-08, x=-0.68, error=5.04825e-08
fSumw[43]=1.46839e-07, x=-0.6, error=5.04825e-08
fSumw[44]=2.11833e-07, x=-0.52, error=5.04825e-08
fSumw[45]=2.29327e-07, x=-0.44, error=5.04825e-08
fSumw[46]=2.87134e-07, x=-0.36, error=5.04825e-08
fSumw[47]=3.0389e-07, x=-0.28, error=5.04825e-08
fSumw[48]=0, x=-0.2, error=0
fSumw[49]=0, x=-0.12, error=0
fSumw[50]=0, x=-0.04, error=0
fSumw[51]=0, x=0.04, error=0
fSumw[52]=1.50715e-07, x=0.12, error=5.04825e-08
fSumw[53]=8.96748e-08, x=0.2, error=5.04825e-08
fSumw[54]=0, x=0.28, error=0
fSumw[55]=9.49502e-08, x=0.36, error=5.04825e-08
fSumw[56]=0, x=0.44, error=0
fSumw[57]=3.50281e-08, x=0.52, error=5.04825e-08
fSumw[58]=0, x=0.6, error=0
fSumw[59]=0, x=0.68, error=0
fSumw[60]=0, x=0.76, error=0
fSumw[61]=0, x=0.84, error=0
fSumw[62]=1.21321e-07, x=0.92, error=5.04825e-08
fSumw[63]=1.92461e-07, x=1, error=5.04825e-08
fSumw[64]=0, x=1.08, error=0
fSumw[65]=0, x=1.16, error=0
fSumw[66]=2.2853e-07, x=1.24, error=5.04825e-08
----------------------------------------------------
Thanks!
Aihong
This archive was generated by hypermail 2b29 : Sat Jan 04 2003 - 23:51:00 MET