TGo4RollingGraph.cxx

Go to the documentation of this file.

// $Id$
//-----------------------------------------------------------------------
//       The GSI Online Offline Object Oriented (Go4) Project
//         Experiment Data Processing at EE department, GSI
//-----------------------------------------------------------------------
// Copyright (C) 2000- GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
//                     Planckstr. 1, 64291 Darmstadt, Germany
// Contact:            http://go4.gsi.de
//-----------------------------------------------------------------------
// This software can be used under the license agreements as stated
// in Go4License.txt file which is part of the distribution.
//-----------------------------------------------------------------------
 
#include "TGo4RollingGraph.h"
 
#include "TMath.h"

TGo4RollingGraph::TGo4RollingGraph()
   : TGraphErrors(), fiNumPoints(0), fiUpdateInterval(1), fiCounter(0), fdSum(0.), fdSumAveraged(0.)
{
}

TGo4RollingGraph::TGo4RollingGraph(Int_t NumPoints, Int_t UpdateInterval)
   : TGraphErrors(), fiNumPoints(NumPoints), fiUpdateInterval(UpdateInterval), fiCounter(0), fdSum(0.),
     fdSumAveraged(0.)
{
}

TGo4RollingGraph::~TGo4RollingGraph()
{
}

void TGo4RollingGraph::Fill(Double_t value, Double_t xerror, Double_t yerror)
{
   if (fiNumPoints)
   {
      Int_t NumAverage = fiCounter % fiUpdateInterval + 1;
      // Create a new bin and reset internal value for averaging if counter is a multiple of fiUpdateInterval
      if (NumAverage == 1)
      {
         NextBin();
         fdSum = 0.;
      }
 
      // Add value for averaging
      fdSum += value;
      fdSumAveraged = fdSum / static_cast<Double_t>(NumAverage);
 
      // Draw new averaged value
      SetPoint(fNpoints-1, static_cast<Double_t>(fiCounter / fiUpdateInterval), fdSumAveraged);
      if (xerror < 0)
         xerror = 0.5;
      if (yerror < 0)
         yerror = fdSumAveraged / TMath::Sqrt(NumAverage);
      SetPointError(fNpoints-1, xerror, yerror);
   }
   fiCounter++;
}

void TGo4RollingGraph::NextBin()
{
   if (fiNumPoints)
   {
      if (fiNumPoints <= fNpoints)
      {
         // Move points to the left
         memmove (&fX[0],  &fX[1],  (fNpoints-1)*sizeof(Double_t));
         memmove (&fEX[0], &fEX[1], (fNpoints-1)*sizeof(Double_t));
         memmove (&fY[0],  &fY[1],  (fNpoints-1)*sizeof(Double_t));
         memmove (&fEY[0], &fEY[1], (fNpoints-1)*sizeof(Double_t));
      }
      else
      {
         // Create a new point
         Set(fNpoints+1);
      }
   }
}

void TGo4RollingGraph::Clear(Option_t *)
{
   for (Int_t i = 0; i < fiNumPoints; ++i)
      RemovePoint(i);
 
   Set(0);
   fdSum = 0.;
}
 
ClassImp(TGo4RollingGraph);