show_pulses.cxx File Reference

#include <stdio.h>
#include <stdlib.h>
#include "midas.h"
#include "../analyzer/compress/mucap_compress.h"
#include "TH1.h"
#include "TCanvas.h"

Go to the source code of this file.

Macros
#define	_FILE_OFFSET_BITS   64
#define	PRINT_SAMPLES   1

Functions
void	process_file (char *input_filename, int addr, int fadc)
int	process_event (void *pevent, int addr, int fadc)
bool	get_input_event (FILE *fp, char *pevent, int event_number)
int	processIsland (int t0, int nsamples, int *samples, int channel)

Variables
TCanvas *	c = 0

Macro Definition Documentation

#define _FILE_OFFSET_BITS   64

Definition at line 1 of file show_pulses.cxx.

#define PRINT_SAMPLES   1

Definition at line 22 of file show_pulses.cxx.

Function Documentation

bool get_input_event	(	FILE *	fp,
		char *	pevent,
		int	event_number
	)

Definition at line 29 of file show_pulses.cxx.

References printf(), and SUCCESS.

{
        int event_size, size_read;
        EVENT_HEADER *header;

        size_read = fread(pevent, 1, sizeof(EVENT_HEADER), fp);
        if (size_read != sizeof(EVENT_HEADER))
                return !SUCCESS;

        header = (EVENT_HEADER *) pevent;
        printf("block# %d\n", header->serial_number);

        if (event_number == 0 && header->trigger_mask != MIDAS_MAGIC) {
                fprintf(stderr, "Input is not a valid MIDAS file\n");
                exit(1);
        }

        event_size = header->data_size;

        if ((unsigned int) event_size >= MAX_EVENT_SIZE - sizeof(EVENT_HEADER)) {
                fprintf(stderr, "Invalid event size %d\n", event_size);
                exit(1);
        }

        size_read = fread(pevent + sizeof(EVENT_HEADER), 1, event_size, fp);
        if (size_read != event_size)
                return !SUCCESS;

        return SUCCESS;
}

int process_event	(	void *	pevent,
		int	addr,
		int	fadc
	)

void process_file	(	char *	input_filename,
		int	addr,
		int	fadc
	)

int processIsland	(	int	t0,
		int	nsamples,
		int *	samples,
		int	channel
	)

Definition at line 104 of file show_pulses.cxx.

References c, i, nsamples, printf(), sprintf(), and title.

{

        c->Clear();
        char title[80];
        sprintf(title, "Waveform: channel %d, t0=%d", channel, t0);

        double sum = 0;
        int min = 4097;
        int max = -1;

        double binWidth = 100./100.00;
        TH1D *waveform = new TH1D("waveform", title, nsamples, -0.5*binWidth, (nsamples-0.5)*binWidth);
        for(int i = 0; i < nsamples; i++) {
                waveform->SetBinContent(i+1, samples[i]);
                sum += samples[i];

                if(samples[i] < min) min = samples[i];
                if(samples[i] > max) max = samples[i];
        }
        waveform->Draw();

        double mean = sum/nsamples;
        double sumSqrDeviations = 0;

        for(int i = 0; i < nsamples; i++) {
                double deviation = samples[i] - mean;
                sumSqrDeviations += deviation*deviation;
        }
        double rms = sqrt(sumSqrDeviations/(nsamples-1));

        printf("Number of samples = %d, Min = %d, Max = %d, Mean = %lf, RMS = %lf\n", nsamples, min, max, mean, rms);

        //  if(max > 3500) {
        c->Update();
        c->Modified();

        char dummy[80];
        fgets(dummy, sizeof(dummy), stdin);
        if(dummy[0] == 'q') {
                return 1;
        }
//}

        delete waveform;
        //delete c;

return 0;
}

Variable Documentation

TCanvas* c = 0

Definition at line 24 of file show_pulses.cxx.