GETSignalFunctions.hh

Go to the documentation of this file.

//======================================================================
/*! \file GETSignalFunctions.hh
 *
 *  Include file for signal analysis functions.
 */
//======================================================================

#include "GETChannel.hh"

//----------------------------------------------------------
/** @name Functions for fitting the signals */
//@{

// global variable for integration step.
extern u_int  getIntegPoints;

// convolution of a square distribution with a Gauss function
extern double GETSquareCurrent     ( const double x[], const double p[] );
extern double GETSquareCurrentGaus ( const double x[], const double p[] );
extern double GETSquareCurrentGSC  ( const double x[], const double p[] );

// same as the previous one, but the time information is the dist. center
extern double GETSquareCurrent2     ( const double x[], const double p[] );
extern double GETSquareCurrent2Gaus ( const double x[], const double p[] );
extern double GETSquareCurrent2GSC  ( const double x[], const double p[] );

//@}


//----------------------------------------------------------
/** @name Gauss function multiplied by sin(x)/x */
//@{

//  -- normalisation
extern TGraph * getGSCNorm;
extern void     GETSetNormalizationGSC ( );

//  -- function
extern double GETGausSinC ( const double x[], const double p[] ); // inline
extern double GETGausSinC ( double t, double amp, double ctr, double sig, double w );

//  -- convolution with the GSC function
extern u_int  getGSCConvPoints;
extern double getGSCConvMax;

extern double GETConvolutionGSC ( const double x[], const double pg[],
                                  double(*fct)(const double[], const double[]),
                                  const double pf[] );
//@}

//======================================================================
//  Inline functions

/* Product of a Gaus function with sinc function: G(x) * sin(x/w) / (x/w).
 *  - p[0] amplitude
 *  - p[1] center
 *  - p[2] sigma
 *  - p[3] sin cardinal dilatation factor
 */
inline double GETGausSinC ( const double x[], const double p[] )
  { return ( GETGausSinC (x[0], p[0], p[1], p[2], p[3]) ); }