#include <calibpulse.h>

Public Member Functions

CalibPulse ()
virtual ~CalibPulse ()
void findPosition (unsigned short *data, int n, float *t1, float *t2, int fLen=4, int prec=-1, int iBegin=0)
float getIntegral (unsigned short *data, int n, float *t1, float *t2, float prec=0, int iBegin=0)
float getGain (unsigned short *data, int n, int nRamp=4)
void getFilteredData (float *t, float *y, int *n)
void getInterpolatedData (float *t, float *y, int *n)
void getPulseShape (float *t, float *y, int *n)
void displayHeader (FILE *fout)
void display (FILE *fout, int col=0, int row=0)
void displayIntegral (FILE *fout, int col=0, int row=0)
float getSNR ()
float getGradient ()
float getPedestal ()
float getPulseHeight ()
int isPulse (float minSNR=1.)
int isMuonPulse (float minGradient=100.)

Private Attributes

float tMin
float tMax
float tBegin
float tEnd
float tIntBegin
float tIntEnd
int filterLength
int len
int indexBegin
long dy1 [FLT_PAGE_SIZE]
long dy2 [FLT_PAGE_SIZE]
int nPed
float meanPed
float varPed
int nPulse
float meanPulse
float varPulse
float maxPulse
float gain
float gradient
float integral

Detailed Description

Analyse the ADC data to find a the calibration pulse and characterize it.

Todo:: Calculate the integral of a pulse!

Definition at line 27 of file calibpulse.h.

Constructor & Destructor Documentation

CalibPulse::CalibPulse ( )

virtual CalibPulse::~CalibPulse ( ) [virtual]

Member Function Documentation

void CalibPulse::display	(	FILE *	fout,
		int	col = `0`,
		int	row = `0`
	)

Display gain and statistical parameters.

void CalibPulse::displayHeader ( FILE * fout )

Display the header of the tabular with the results.

void CalibPulse::displayIntegral	(	FILE *	fout,
		int	col = `0`,
		int	row = `0`
	)

Dispaly integration of the pulse.

void CalibPulse::findPosition	(	unsigned short *	data,
		int	n,
		float *	t1,
		float *	t2,
		int	fLen = `4`,
		int	prec = `-1`,
		int	iBegin = `0`
	)

The function is intended to find the position of the calibration pulse in the ADC data.

The ADC data is expected to be an array of size n it spans 100us.

Parameters:

data	Pointer to the ADC data
n	Length of the ADC data array
t1	Calculated start of the pulse
t2	Calculated end of the pulse
fLen	Length of the edge filter used to determine the position of the peak
prec	Desired precision of the result. The parameter specifies the the number of samples used for the calculation: -1 interpolate between the samples, 1 use every sample, 2 use every second sample, i use every i.th sample.
iBegin	Index of the first element of the given data. Only necessary if the analysis# should not start at the beginning of the ADC data.

void CalibPulse::getFilteredData	(	float *	t,
		float *	y,
		int *	n
	)

Return the filtered data to display the filter result.

float CalibPulse::getGain	(	unsigned short *	data,
		int	n,
		int	nRamp = `4`
	)

Calculate the gain of the complete signal chain from the PMT to the digital interface.

The calculation will use the position of the pulse determined by the function findPosition - the pointer to the ADC data has to be the same in both calls!

Parameters:

data	Pointer to the ADC data
n	Length of the ADC data array
nRamp	Number of samples to leave out at the edges of the pulse

float CalibPulse::getGradient ( )

Calculate gradient of the pulse.

Alogrithm: (mean_(pulse) - mean(var)) / ( t_2 - t_1)

float CalibPulse::getIntegral	(	unsigned short *	data,
		int	n,
		float *	t1,
		float *	t2,
		float	prec = `0`,
		int	iBegin = `0`
	)

Calculate the integral charge of the pulse.

The maximum value and the pedestal obtained from getGain is used to determine the length of the pulse and the charge.

Parameters:

data	Pointer to the ADC data
n
t1	Lower limit
t2	Upper limit
prec	Minimal value for the integral border
iBegin	Index of the first element of the given data. Only necessary if the analysis# should not start at the beginning of the ADC data.
iBegin