Changeset 2604 in Sophya for trunk/SophyaLib/HiStats

Timestamp:

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

Author:

cmv

Message:

nouvelles methode dans HistoErr cmv 3/9/04

Location:

trunk/SophyaLib/HiStats

Files:

• histerr.cc (modified) (2 diffs)
• histerr.h (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;
 }
 

trunk/SophyaLib/HiStats/histerr.h

@@ -26 +26 @@
   //! Addition du contenu de l'histo pour le bin numBin poids w et l'erreur e
   void AddBin(int_4 numBin, r_8 w=1., r_8 e=1.);
+  //! remplissage contenu de l'histo pour le bin numBin poids w et l'erreur e
+  void SetBin(int_4 numBin, r_8 w=1., r_8 e=1., uint_4 nb=1);
+  //! remplissage nombre d'entrees pour le bin numBin
+  void SetNentB(int_4 numBin, uint_4 nb=1);
 
   //! Retourne le nombre d'entree dans le bin
@@ -34 +38 @@
   void GetNBin(TVector<int_4>& v) const;
   void PutNBin(TVector<int_4>& v);
+
+  // Compute the correlation histogram
+  void ToCorrel(void);
 
 protected: