EmcWaveform.cxx

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oriented Simulation Tool                     //
//---------------------------------------------------------------------------//
// Description:
// Author: Hemiao
// Created: Oct 25, 2004
// Modified:
// Comment:
//---------------------------------------------------------------------------//
// $Id: EmcWaveform.cxx
 
#include "EmcWaveform/EmcWaveform.h"
#include <math.h>
 
// Constructors
 
EmcWaveform::EmcWaveform() {
  array_size   = 60;
  m_tau        = 1000.;
  m_sampleTime = 50.;
  m_peakTime   = 4000.;
  m_timeOffset = 2500.;
  emcWave      = new double[array_size];
 
  for ( int i = 0; i < array_size; i++ ) emcWave[i] = 0;
}
 
EmcWaveform::EmcWaveform( int size, double tau, double sampleTime )
    : m_tau( tau ), m_sampleTime( sampleTime ) {
  if ( size > 0 )
  {
    array_size   = size;
    emcWave      = new double[array_size];
    double* init = emcWave + array_size;
    while ( init != emcWave ) *--init = 0.0;
  }
}
 
// Destructors
EmcWaveform::~EmcWaveform() {
  delete[] emcWave;
  emcWave = 0;
}
 
// Operators
EmcWaveform& EmcWaveform::operator*=( const double scale ) {
  for ( int i = 0; i < array_size; i++ ) emcWave[i] *= scale;
  return *this;
}
 
EmcWaveform& EmcWaveform::operator/=( const double scale ) {
  for ( int i = 0; i < array_size; i++ ) emcWave[i] /= scale;
  return *this;
}
 
EmcWaveform& EmcWaveform::operator+=( const EmcWaveform& assign ) {
  for ( int i = 0; i < array_size; i++ ) emcWave[i] += assign[i];
  return *this;
}
 
EmcWaveform& EmcWaveform::operator=( const EmcWaveform& assign ) {
  if ( this != &assign )
  {
    if ( emcWave != 0 ) delete[] emcWave;
    emcWave    = new double[assign.array_size];
    array_size = assign.array_size;
    for ( int i = 0; i < array_size; i++ ) emcWave[i] = assign[i];
  }
  return *this;
}
 
double EmcWaveform::max( int& binOfMax ) const {
  double maxi = emcWave[0];
  binOfMax    = 0;
  for ( int i = 1; i < array_size; i++ )
  {
    if ( emcWave[i] > maxi )
    {
      maxi     = emcWave[i];
      binOfMax = i;
    }
  }
  return maxi;
}
 
double EmcWaveform::getADCTrg( int time ) {
  double adc = emcWave[0];
  for ( int i = 0; i < array_size; i++ )
  {
    if ( time >= i * m_sampleTime && time < ( i + 1 ) * m_sampleTime ) adc = emcWave[i];
  }
  return adc;
}
 
int EmcWaveform::frontEdgeTrg( double thres ) {
  int binOfMax = -1;
  for ( int i = 0; i < array_size; i++ )
  {
    if ( emcWave[i] > thres )
    {
      binOfMax = i;
      break;
    }
  }
  return binOfMax;
}
 
void EmcWaveform::makeWaveformTrg( double energy, double time ) {
  double time0    = time * m_sampleTime + m_timeOffset - m_peakTime; // start time
  double tempTime = 0;                                               // time for each bin
 
  double peak =
      m_peakTime * m_peakTime * m_peakTime * m_peakTime * exp( -m_peakTime / m_tau ) / 24;
 
  for ( int i = 0; i < array_size; i++ )
  {
    tempTime = i * m_sampleTime + m_timeOffset - time0;
    if ( tempTime > 0 )
      emcWave[i] += energy * tempTime * tempTime * tempTime * tempTime *
                    exp( -tempTime / m_tau ) / ( 24 * peak );
  }
}
 
void EmcWaveform::makeWaveform( double energy, double time ) {
  double time0    = time * m_sampleTime + m_timeOffset - m_peakTime; // start time
  double tempTime = 0;                                               // time for each bin
 
  double peak =
      m_peakTime * m_peakTime * m_peakTime * m_peakTime * exp( -m_peakTime / m_tau ) / 24;
 
  for ( int i = 0; i < array_size; i++ )
  {
    tempTime = i * m_sampleTime + m_timeOffset - time0;
    if ( tempTime > 0 )
      emcWave[i] += energy * tempTime * tempTime * tempTime * tempTime *
                    exp( -tempTime / m_tau ) / ( 24 * peak );
  }
}
 
void EmcWaveform::print() {
  cout << "New Wave!" << endl;
  for ( int i = 0; i < array_size; i++ ) cout << emcWave[i] << "\t";
  cout << endl;
}