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

Timestamp:

Jul 7, 2000, 11:44:16 AM (25 years ago)

Author:

ansari

Message:

cmv 7/7/2000

• Introduction methode HRebin pour les HProf
• Passage de HRebin de Histo et HProf en virtual

File:

• trunk/SophyaLib/HiStats/hisprof.cc (modified) (1 diff)

trunk/SophyaLib/HiStats/hisprof.cc

@@ -234 +234 @@
 }
 
+/********* Methode *********/
+/*!
+  Pour rebinner l'histogramme de profile sur nbinew bins
+*/
+void HProf::HRebin(int nbinew)
+{
+  if( bins <= 0 ) return; // createur par default
+  if( nbinew <= 0 ) return;
+
+  // memorisation de l'HProf original
+  HProf H(*this);
+
+  // Rebin de la partie Histo
+  int binold = bins;
+  Histo::HRebin(nbinew);
+
+  // Le nombre de bins est il plus grand ??
+  if( bins > binold ) {
+    delete [] SumY;  SumY  = new double[nbinew]; memset(SumY, 0,bins*sizeof(double));
+    delete [] SumY2; SumY2 = new double[nbinew]; memset(SumY2,0,bins*sizeof(double));
+    delete [] SumW;  SumW  = new double[nbinew]; memset(SumW, 0,bins*sizeof(double));
+  }
+
+  // Remplissage des parties propres au HPprof
+  for(int i=0;i<bins;i++) {
+    float xmi = BinLowEdge(i);
+    float xma = BinHighEdge(i);
+    int imi =  H.FindBin(xmi);
+    if( imi < 0 ) imi = 0;
+    int ima =  H.FindBin(xma);
+    if( ima >= H.bins ) ima = H.bins-1;
+    double wY = 0., wY2 = 0., wW = 0.;
+    if( ima == imi ) {
+      wY  = H.SumY[imi]  * binWidth/H.binWidth;
+      wY2 = H.SumY2[imi] * binWidth/H.binWidth;
+      wW  = H.SumW[imi]  * binWidth/H.binWidth;
+    } else {
+      wY  += H.SumY[imi]  * (H.BinHighEdge(imi)-xmi)/H.binWidth;
+      wY  += H.SumY[ima]  * (xma -H.BinLowEdge(ima))/H.binWidth;
+      wY2 += H.SumY2[imi] * (H.BinHighEdge(imi)-xmi)/H.binWidth;
+      wY2 += H.SumY2[ima] * (xma -H.BinLowEdge(ima))/H.binWidth;
+      wW  += H.SumW[imi]  * (H.BinHighEdge(imi)-xmi)/H.binWidth;
+      wW  += H.SumW[ima]  * (xma -H.BinLowEdge(ima))/H.binWidth;
+      if(ima>imi+1) for(int ii=imi+1;ii<ima;ii++)
+        {wY += H.SumY[ii]; wY2 += H.SumY2[ii]; wW += H.SumW[ii];}
+    }
+    SumY[i] += wY; SumY2[i] += wY2; SumW[i] += wW;
+  }
+  // On synchronise les tableaux Sum?? et l'Histogramme
+  UpdateHisto();
+}
+
 ///////////////////////////////////////////////////////////
 // --------------------------------------------------------