Changeset 2604 in Sophya for trunk/SophyaLib/HiStats/histerr.cc

Timestamp:

Sep 3, 2004, 6:24:38 PM (21 years ago)

Author:

cmv

Message:

nouvelles methode dans HistoErr cmv 3/9/04

File:

• trunk/SophyaLib/HiStats/histerr.cc (modified) (2 diffs)

trunk/SophyaLib/HiStats/histerr.cc

@@ -83 +83 @@
 }
 
+/********* Methode *********/
+/*!
+ Remplissage du contenu de l'histo pour le bin numBin poids w et l'erreur e
+*/
+void HistoErr::SetBin(int_4 numBin, r_8 w, r_8 e, uint_4 nb)
+{
+Histo::SetBin(numBin,w);
+Histo::SetErr2(numBin,e*e);
+SetNentB(numBin,nb);
+return;
+}
+
+/********* Methode *********/
+/*!
+ Remplissage nombre d'entrees pour le bin numBin
+*/
+void HistoErr::SetNentB(int_4 numBin, uint_4 nb)
+{
+ if(numBin>=0 && numBin<mBins) mNData[numBin] = nb;
+return;
+}
+
 /*!
   Remplissage d'un tableau avec les nombres d'entrees dans le bin de l'histo
@@ -99 +121 @@
 void HistoErr::PutNBin(TVector<int_4> &v)
 {
-uint_4 n = (v.NElts()<(uint_4) mBins) ? v.NElts(): (uint_4) mBins;
-if(n>0) for(uint_4 i=0;i<n;i++) mNData[i] = v(i);
-return;
+int_4 n = (v.NElts()<mBins) ? v.NElts(): mBins;
+if(n>0) for(int_4 i=0;i<n;i++) mNData[i] = v(i);
+return;
+}
+
+
+/********* Methode *********/
+/*!
+  Compute the correlation histogram.
+  Each bin content is divided by the number of entries in that bin.
+  Each squared error is divided by the number of entries in that bin.
+  The number of entries by bin is set to 1
+*/
+void HistoErr::ToCorrel(void)
+{
+ if(mBins<1) return;
+ for(int_4 i=0;i<mBins;i++) {
+   if(mNData[i]<1) continue;
+   mData[i] /= (r_8) mNData[i];
+   mErr2[i] /= (r_8) mNData[i];
+   mNData[i] = 1;
+ }
+ return;
 }