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

Timestamp:

Jul 7, 2000, 10:41:11 AM (25 years ago)

Author:

ansari

Message:

cmv 7/7/2000

hisprof.cc,h:

fonction IsOk()
GetMean -> GetValue (mauvais nom)
float Error2() -> double Error2()
nouveau: GetError(TVector<r_8>& v)
HProf::PrintF() sur-ecriture de Histo::PrintF
protection dans createur par copie dans alloc

SumY... pour le cas ou H.bins==0

protection dans UpdateHisto pour HProf cree par defaut (bins==0)
Dans WriteSelf UpdateHisto appele que si necessaire (IsOk()?)

histos2.cc : le Print() dit si les erreurs ont ete demandees
histos.cc,h:

protection dans createur par copie dans alloc

pour le cas ou H.bins==0

protection dans Zero() pour Histo cree par defaut (bins==0)
protection dans operator=() pour Histo cree par defaut (bins==0)
le Print() dit si les erreurs ont ete demandees

cmv 7/7/2000

File:

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

trunk/SophyaLib/HiStats/hisprof.cc

@@ -47 +47 @@
 HProf::HProf(const HProf& H)
 : Histo(H)
-, SumY(new double[H.bins]), SumY2(new double[H.bins]), SumW(new double[H.bins])
+, SumY((H.bins>0)  ? new double[H.bins] : NULL)
+, SumY2((H.bins>0) ? new double[H.bins] : NULL)
+, SumW((H.bins>0)  ? new double[H.bins] : NULL)
 , Ok(H.Ok), YMin(H.YMin), YMax(H.YMax), Opt(H.Opt)
 {
-  memcpy(SumY,  H.SumY,  bins*sizeof(double));
-  memcpy(SumY2, H.SumY2, bins*sizeof(double));
-  memcpy(SumW,  H.SumW,  bins*sizeof(double));
+  if(bins>0) {
+    memcpy(SumY,  H.SumY,  bins*sizeof(double));
+    memcpy(SumY2, H.SumY2, bins*sizeof(double));
+    memcpy(SumW,  H.SumW,  bins*sizeof(double));
+  }
   END_CONSTRUCTOR
 }
@@ -125 +129 @@
 void HProf::UpdateHisto() const
 {
-
 float m,e2;
+if(bins<=0) return;
 for(int i=0;i<bins;i++) {
   if(SumW[i]<=0.) {
@@ -185 +189 @@
 /********* Methode *********/
 /*!
-  Operateur H = H1
+  Operateur HProf H = H1
 */
 HProf& HProf::operator = (const HProf& h)
@@ -209 +213 @@
 /********* Methode *********/
 /*!
-  Operateur H += H1
+  Operateur HProf H += H1
 
   Attention dans cette addition il n'y a pas de gestion
@@ -273 +277 @@
 char strg[256];
 
-dobj->UpdateHisto();
+if(!(dobj->IsOk())) dobj->UpdateHisto();
 
 // Ecriture entete pour identifier facilement