QOnlineScope.hh

Go to the documentation of this file.

#ifndef _QONLINE_SCOPE_HH_
#define _QONLINE_SCOPE_HH_

// ROOT Includes
#include "QCuore.hh"
#include "TH1F.h"

#define DEFAULT_BUFFER_SIZE 60000 
#define DEFAULT_HISTO_SIZE 10000  
#define NUMBER_OF_CHANNELS 4      
// shared memory structure
struct QOnlineScopeBuffer_t 
{
   unsigned long fId[NUMBER_OF_CHANNELS];    
   unsigned long fSize;                      
   unsigned long fRead[NUMBER_OF_CHANNELS];  
   unsigned long fWrite[NUMBER_OF_CHANNELS]; 
   long fHistoSize[NUMBER_OF_CHANNELS];      
   float fData0[DEFAULT_BUFFER_SIZE];        
   float fData1[DEFAULT_BUFFER_SIZE];        
   float fData2[DEFAULT_BUFFER_SIZE];        
   float fData3[DEFAULT_BUFFER_SIZE];        
};

class QOnlineScope {
public:
   QOnlineScope(bool attach);
   
   virtual ~QOnlineScope();
   
   inline void Add(float data, int index=0);
   
   inline unsigned long GetLg(int index=0) const {return p_buffer->fId[index];}
   
   inline void SetLg(unsigned long lg, int index=0);

   // Update histogram
   inline TH1F* UpdateHisto(int index=0);

   // Get mean value and sigma of the histogram (on Y axis!)
   void GetMeanSigma(double &mean, double &rms, int index=0);
   
   // Getters for getting Mean or Sigma individually
   const double GetMean(int index=0);
   const double GetSigma(int index=0);
   
private:
   
   int  fHandle;
   int  fShmid;                        // shared memory id
   bool fOwn;                          // true if owner of shared memory
   QOnlineScopeBuffer_t *p_buffer;     // shared memory buffer pointer
   TH1F *fHisto[NUMBER_OF_CHANNELS];   // root histogram with scope plot
   float *GetData(int index);          // return pointer to buffer corresponding to index
   
};


//_____________________________________________________________________________
inline void QOnlineScope::SetLg(unsigned long lg,int index) 
{
   //
  // set new logical channel 
  // all buffer is cleared
  //

  if (lg == GetLg(index)) return;
  if (fHisto[index]) delete fHisto[index];
  
  fHisto[index]=0;
  p_buffer->fId[index] = lg;
  p_buffer->fRead[index]=0;
  p_buffer->fWrite[index]=0;
  
  return;
}


//_____________________________________________________________________________
inline float*  QOnlineScope::GetData(int index) 
{
  //
  // Get the right buffer
  //
  
  switch(index) 
  {
  case 1:
    return p_buffer->fData1; 
    break;
  case 2: 
    return p_buffer->fData2; 
    break;
  case 3: 
    return p_buffer->fData3; 
    break;
  case 0: 
  default:
    return p_buffer->fData0; 
    break;
  }
}
                        

//_____________________________________________________________________________
inline void QOnlineScope::Add(float data, int index) 
{
  //
  // Add data to buffer 
  //
  
  float *pData=GetData(index);
  pData[p_buffer->fWrite[index]] = data;
  p_buffer->fWrite[index]++;
  if (p_buffer->fWrite[index] == p_buffer->fSize)
    p_buffer->fWrite[index]=0;
}

//_____________________________________________________________________________
inline TH1F *QOnlineScope::UpdateHisto(int index) 
{
  if (GetLg(index)==0) return (TH1F*)0;
  
  if (NULL == fHisto[index]) 
  {
    char title[200];
    char name[200];
    sprintf(title,"Scope");
    sprintf(name,"Real Time Waveform Channel=%lu",GetLg(index));
    fHisto[index] = new TH1F(title,name,p_buffer->fHistoSize[index],0.,10.);
  } 
  
  // get available number of words in buffer 
  int nw = p_buffer->fWrite[index] - p_buffer->fRead[index];
  
  if (nw<0) nw += p_buffer->fSize;
  
  // do nothing if there are not enough new samples
  if (nw<200) return fHisto[index];
  
  if (nw>200) nw=200;
    
  float *pData = GetData(index);
    
  // update histogram shifting data to left 
  if (nw > p_buffer->fHistoSize[index]) 
  {
    nw = p_buffer->fHistoSize[index];
    for(int i=0; i<p_buffer->fHistoSize[index]; i++) 
    {
      int j = (p_buffer->fRead[index]+i)%p_buffer->fSize;
      fHisto[index]->SetBinContent(i+1,pData[j]);
    }
  }
  else 
  {
    for(int i=1; i<=p_buffer->fHistoSize[index]-nw; i++) 
    {
      double val =  fHisto[index]->GetBinContent(i+nw);
      fHisto[index]->SetBinContent(i,val);
    }
    for(int i=0; i<nw; i++) 
    {
      int j = (p_buffer->fRead[index]+i)%p_buffer->fSize;
      fHisto[index]->SetBinContent(p_buffer->fHistoSize[index]-nw+i+1,pData[j]);
    }
  }     
    
  // release data 
  p_buffer->fRead[index] += nw;
  if (p_buffer->fRead[index] >= p_buffer->fSize)
    p_buffer->fRead[index] -= p_buffer->fSize;
    
  return fHisto[index];
}

#endif // _QONLINE_SCOPE_HH_

---