Changeset 2312 in Sophya for trunk/SophyaLib/HiStats/hisprof.cc

Timestamp:

Jan 21, 2003, 7:36:01 PM (23 years ago)

Author:

cmv

Message:

bugs sur calcul sigma et erreur sur moyenne cmv 21/01/2003

File:

• trunk/SophyaLib/HiStats/hisprof.cc (modified) (3 diffs)

trunk/SophyaLib/HiStats/hisprof.cc

@@ -131 +131 @@
   \verbatim
     - Pour le bin (j):
-    H(j) = sum(y), E(j) = sum(y^2), L(j) = sum(w)
-    ->  s(j) = sqrt(E(j)/L(j) - (H(j)/L(j))^2)  dispersion
-    ->  e(j) = s(j)/sqrt(L(j)) erreur sur la moyenne
+    H(j) = sum(y*w), E(j) = sum(y^2*w), L(j) = sum(w)
+    ->  h(j) = H(j)/L(j)   moyenne dans le bin
+    ->  s(j) = sqrt(E(j)/L(j) - (H(j)/L(j))^2)  dispersion dans le bin
+    ->  e(j) = s(j)/sqrt(L(j)) erreur sur la moyenne dans le bin
     spread=true:  opt=0 : dispersion des donnees dans le bin = s(j)
     spread=false: opt=1 : erreur sur la moyenne du bin = e(j)
@@ -189 +190 @@
   else {
     Ok = false;
-    y *= w;
-    SumY[numBin]  += y;
-    SumY2[numBin] += y*y;
+    SumY[numBin]  += y*w;
+    SumY2[numBin] += y*y*w;
     SumW[numBin]  += w;
     nHist += w;
@@ -209 +209 @@
   else {
     Ok = false;
-    y *= w;
-    SumY[numBin]  += y;
-    SumY2[numBin] += y*y;
+    SumY[numBin]  += y*w;
+    SumY2[numBin] += y*y*w;
     SumW[numBin]  += w;
     nHist += w;