common.h

Go to the documentation of this file.

/*
        peakfinder
        Copyright (C) 2005 Hal Finkel

        This program is free software; you can redistribute it and/or modify
        it under the terms of the GNU General Public License as published by
        the Free Software Foundation; either version 2 of the License, or
        (at your option) any later version.

        This program is distributed in the hope that it will be useful,
        but WITHOUT ANY WARRANTY; without even the implied warranty of
        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
        GNU General Public License for more details.

        You should have received a copy of the GNU General Public License
        along with this program; if not, write to the Free Software
        Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
*/

#ifndef PF_COMMON_H
#define PF_COMMON_H

#define NCHAN 4096

#define MAXCALPTS 1024

#define MAXUNIT 256

struct pf_peak_bg {
        double a; 
        double b; 
        double c; 
        double d; 
        double a_err; 
        double b_err; 
        double c_err; 
        double d_err; 
        double chisq; 
        double chisq_dof; 
        int conv; 
};

struct pf_peak_fit {
        double a; 
        double k; 
        double x0; 
        double m; 
        double a_err; 
        double k_err; 
        double x0_err; 
        double m_err; 
        double chisq; 
        double chisq_dof; 
        int conv; 
};

struct pf_peak {
        int num; 
        int lbin; 
        int cbin; 
        int rbin; 
        struct pf_peak_bg bg; 
        struct pf_peak_fit fit; 
        double center; 
        double width; 
        double lbound; 
        double rbound; 
        double bglbound; 
        double bgrbound; 
        double area; 
        double height; 
};

struct pf_cal_fit {
        char unit[MAXUNIT]; 
        double a; 
        double b; 
        double c; 
        double a_err; 
        double b_err; 
        double c_err; 
        double chisq; 
        double chisq_dof; 
        int conv; 
        int valid; 
};

struct pf_cal_pt {
        int num; 
        double pt; 
        double ri; 
};

struct pf_cal_pts {
        char unit[MAXUNIT]; 
        int npts; 
        struct pf_cal_pt pts[MAXCALPTS]; 
        struct pf_cal_pt* pts_by_value[MAXCALPTS]; 
        struct pf_cal_pt* pts_by_ri[MAXCALPTS]; 
};

struct pf_data {
        double data_raw[NCHAN]; 
        double data_norm[NCHAN]; 
        double max_raw; 
        double min_raw; 
        double max_norm; 
        double min_norm; 
        int total_counts; 
        double total_chisq; 
        int data_bg[NCHAN]; 
        int npeaks; 
        struct pf_peak peaks[NCHAN]; 
        struct pf_peak* peaks_by_number[NCHAN]; 
        struct pf_peak* peaks_by_height[NCHAN]; 
        struct pf_peak* peaks_by_area[NCHAN]; 
};

struct pf_plot_params {
        double xstart; 
        double xend; 
        int xstartcal; 
        int xendcal; 
        int peak; 
        char* fn; 
        char* dev; 
        double rot; 
        int norm; 
        int rescale; 
        int cal; 
        int marked; 
        int annot; 
};

enum pf_peak_sort {
        sort_number, 
        sort_height, 
        sort_area 
};

int pf_parse_command(char* command);

int pf_loaddata(char* fn, struct pf_data* d);

double pf_chanwindowavg(double* data, int ds, int ci, int ac);

void pf_plotdata(struct pf_data* d, struct pf_plot_params* p, struct pf_cal_fit* cal);

void pf_plotpeak(struct pf_data* d, struct pf_plot_params* p, struct pf_cal_fit* cal);

int pf_printdata(struct pf_data* d, char* fn, int norm);

int pf_printpeaks(struct pf_data* d, char* fn, enum pf_peak_sort srt);

int pf_printpeak(struct pf_data* d, struct pf_cal_fit* cal, int peak, char* fn);

int pf_printcal(struct pf_data* d, struct pf_cal_fit* f, char* fn);

void pf_fitpeaks(struct pf_data* d);

void pf_prunepeaks(struct pf_data* d);

double pf_bgfunc(struct pf_data* d, int peak, double x);

double pf_peakfunc(struct pf_data* d, int peak, double x);

double pf_fitfunc(struct pf_data* d, int peak, double x);

int pf_calibrate(struct pf_data* d, struct pf_cal_fit* f, char* fn);

double pf_calfunc(struct pf_cal_fit* f, double ch);

#endif

---