tof_in.c

#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <jibal_units.h>
#ifdef WIN32
#include "win_compat.h"
#endif
#include "message.h"
#include "tof_in.h"


int tofin_file_parse_header(tofin_file *tofin, const char *header, const char *data) {
    char *end;
    switch (jibal_option_get_value(tofin_headers, header)) {
        case TOFIN_HEADER_BEAM:
            break;
        case TOFIN_HEADER_ENERGY:
            break;
        case TOFIN_HEADER_DETECTOR_ANGLE:
            break;
        case TOFIN_HEADER_TARGET_ANGLE:
            break;
        case TOFIN_HEADER_TOFLEN:
            tofin->toflen = strtod(data, &end);
            if(end == data) {
                tofe_list_msg(TOFE_LIST_ERROR, "ToF length could not be parsed.");
                return -1;
            }
            break;
        case TOFIN_HEADER_CARBON_FOIL_THICKNESS:
            tofin->foil_thickness = strtod(data, &end);
            if(end == data) {
                tofe_list_msg(TOFE_LIST_ERROR, "Carbon foil thickness could not be parsed.");
            }
            tofin->foil_thickness *= (C_UG/C_CM2);
            break;
        case TOFIN_HEADER_TARGET_DENSITY:
            break;
        case TOFIN_HEADER_TOF_CALIBRATION:
            if(sscanf(data, "%lf %lf", &tofin->tof_slope, &tofin->tof_offset) != 2) {
                tofe_list_msg(TOFE_LIST_ERROR, "ToF calibration could not be parsed.");
                return -1;
            }
            break;
        case TOFIN_HEADER_ANGLE_CALIBRATION:
            if(sscanf(data, "%lf %lf", &tofin->angle_slope, &tofin->angle_offset) != 2) {
                tofe_list_msg(TOFE_LIST_ERROR, "Angle calibration could not be parsed.");
                return -1;
            }
            break;
        case TOFIN_HEADER_DEPTH_STEPS:
            break;
        case TOFIN_HEADER_DEPTH_STEP_STOPPING:
            break;
        case TOFIN_HEADER_DEPTH_STEP_OUTPUT:
            break;
        case TOFIN_HEADER_DEPTHS_CONCENTRATION_SCALING:
            break;
        case TOFIN_HEADER_CROSS_SECTION:
            break;
        case TOFIN_HEADER_ITERATIONS:
            break;
        case TOFIN_HEADER_EFFICIENCY_DIRECTORY:
            free(tofin->efficiency_directory);
            tofin->efficiency_directory = strdup(data);
            break;
        case TOFIN_HEADER_NONE:
        default:
            tofe_list_msg(TOFE_LIST_WARNING, "Unknown header \"%s\"", header);
            break;
    }
    return 0;
}

tofin_file *tofin_file_load(const char *filename) {
    FILE *f = fopen(filename, "r");
    if(!f) {
        return NULL;
    }
    tofin_file *tofin = calloc(1, sizeof(tofin_file));
    char *line = NULL, *line_data;
    size_t line_size = 0;
    int lineno = 0;
    int error = 0;
    while(getline(&line, &line_size, f) > 0) {
        lineno++;
        if(*line == '#') /* Comments are allowed */
            continue;
        if(*line == '\r' || *line == '\n') { /* Empty line signals end of headers */
            break;
        }
        line[strcspn(line, "\r\n")] = 0; /* Strips all kinds of newlines! */
        line_data = line;
        strsep(&line_data, ":");
        while(*line_data == ' ') {line_data++;} /* Skip whitespace */
        if(tofin_file_parse_header(tofin, line, line_data)) {
            tofe_list_msg(TOFE_LIST_ERROR, "Parse error on line %i of file \"%s\".", lineno, filename);
            error++;
            break;
        }
    }
    if(error) {
        tofin_file_free(tofin);
        return NULL;
    }
    free(line);
    fclose(f);
    tofe_list_msg(TOFE_LIST_INFO, "Read %i lines from \"%s\".\n", lineno, filename);
    return tofin;
}

void tofin_file_free(tofin_file *tofin) {
    free(tofin->efficiency_directory);
    free(tofin);
}