3DapReader.cpp File Reference

#include "3DapReader.h"
#include <vector>
#include <iostream>
#include <cstdlib>
#include <fstream>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_fit.h>
#include <gsl/gsl_multifit.h>
#include "common/basics.h"

Include dependency graph for 3DapReader.cpp:

Classes
struct	SINGLE_HIT
struct	VOLTAGE_DATA
struct	THREEDAP_DATA
struct	THREEDAP_EXPERIMENT

Functions
float	evaluate_polyBase (float x, float y, unsigned int baseNum, unsigned int order)
unsigned int	polyFitData (const vector< float > &obsX, const vector< float > &obsY, const vector< float > &obsZ, vector< float > &coeff, unsigned int polyOrder)
unsigned int	readOps (const char *file, THREEDAP_EXPERIMENT &data, unsigned int &badLine, unsigned int &progress, ATOMIC_BOOL &wantAbort, bool strictMode=false)
void	computeMToC (THREEDAP_EXPERIMENT &experiment)
void	cylinderReconstruct (const THREEDAP_EXPERIMENT &experiment, float ionVolume, float scaleRadius, std::vector< IonHit > &ions)
unsigned int	reconstruct3DapData (const char *file, const THREEDAP_RECON_PARAMS &params, unsigned int &progress, ATOMIC_BOOL &wantAbort, std::vector< IonHit > &ionData)
	Reconstruct an "OPS" formatted file (from a 3Dap) to an IonHit vector. More...

Variables
const unsigned int	NUM_XY =2
const float	CHARGE_ELECTRON = 1.602176e-19
const float	AMU_PER_KG = 6.0221409e26
const char *	OPS_ENUM_ERRSTRINGS []

Function Documentation

computeMToC()

void computeMToC

(

THREEDAP_EXPERIMENT &

experiment

)

References THREEDAP_DATA::alpha, ASSERT, VOLTAGE_DATA::beta, THREEDAP_EXPERIMENT::eventData, THREEDAP_DATA::flightPath, VOLTAGE_DATA::nextHitGroupOffset, THREEDAP_EXPERIMENT::params, VOLTAGE_DATA::pulseVolt, THREEDAP_DATA::tZero, VOLTAGE_DATA::voltage, and THREEDAP_EXPERIMENT::voltageData.

Referenced by reconstruct3DapData().

Here is the caller graph for this function:

cylinderReconstruct()

void cylinderReconstruct	(	const THREEDAP_EXPERIMENT &	experiment,
		float	ionVolume,
		float	scaleRadius,
		std::vector< IonHit > &	ions
	)

References ASSERT, THREEDAP_DATA::detectorChannels, THREEDAP_EXPERIMENT::eventData, THREEDAP_EXPERIMENT::params, IonHit::setMassToCharge(), and IonHit::setPos().

Referenced by reconstruct3DapData().

Here is the caller graph for this function:

evaluate_polyBase()

float evaluate_polyBase	(	float	x,
		float	y,
		unsigned int	baseNum,
		unsigned int	order
	)

polyFitData()

unsigned int polyFitData	(	const vector< float > &	obsX,
		const vector< float > &	obsY,
		const vector< float > &	obsZ,
		vector< float > &	coeff,
		unsigned int	polyOrder
	)

readOps()

unsigned int readOps	(	const char *	file,
		THREEDAP_EXPERIMENT &	data,
		unsigned int &	badLine,
		unsigned int &	progress,
		ATOMIC_BOOL &	wantAbort,
		bool	strictMode = false
	)

References THREEDAP_DATA::alpha, ASSERT, VOLTAGE_DATA::beta, THREEDAP_DATA::beta, COMPILE_ASSERT, THREEDAP_DATA::detectorChannels, THREEDAP_DATA::detectorRadius, THREEDAP_EXPERIMENT::eventData, THREEDAP_EXPERIMENT::eventPulseNumber, THREEDAP_DATA::flightPath, getFilesize(), VOLTAGE_DATA::nextHitGroupOffset, NUM_XY, OPS_ENUM_ERRSTRINGS, OPSREADER_ABORT_ERR, OPSREADER_CHANNELS_DATA_ERR, OPSREADER_ENUM_END, OPSREADER_FORMAT_CHANNELS_ERR, OPSREADER_FORMAT_CLINE_ERR, OPSREADER_FORMAT_DETECTORLINE_ERR, OPSREADER_FORMAT_DOUBLEDASH, OPSREADER_FORMAT_DUPLICATE_DETECTORSIZE, OPSREADER_FORMAT_DUPLICATE_SYSDATA, OPSREADER_FORMAT_LINE_DASH_ERR, OPSREADER_FORMAT_LINE_VOLTAGE_DATA_ERR, OPSREADER_FORMAT_LINE_VOLTAGE_DATACOUNT_ERR, OPSREADER_FORMAT_LINE_VOLTAGE_NOBETA, OPSREADER_FORMAT_LINETYPE_ERR, OPSREADER_FORMAT_SLINE_EVENTCOUNT_ERR, OPSREADER_FORMAT_SLINE_EVENTDATA_ERR, OPSREADER_FORMAT_SLINE_FORMAT_ERR, OPSREADER_FORMAT_SLINE_PREFIX_ERR, OPSREADER_FORMAT_TRAILING_DASH_ERR, OPSREADER_OPEN_ERR, OPSREADER_READ_ERR, THREEDAP_EXPERIMENT::params, VOLTAGE_DATA::pulseVolt, splitStrsRef(), stream_cast(), stripWhite(), THREEDEP_ARRAYSIZE, SINGLE_HIT::tof, THREEDAP_DATA::tZero, VOLTAGE_DATA::voltage, THREEDAP_EXPERIMENT::voltageData, SINGLE_HIT::x, and SINGLE_HIT::y.

Referenced by reconstruct3DapData().

Here is the caller graph for this function:

reconstruct3DapData()

unsigned int reconstruct3DapData	(	const char *	file,
		const THREEDAP_RECON_PARAMS &	params,
		unsigned int &	progress,
		ATOMIC_BOOL &	wantAbort,
		std::vector< IonHit > &	ionData
	)

Reconstruct an "OPS" formatted file (from a 3Dap) to an IonHit vector.

References computeMToC(), cylinderReconstruct(), THREEDAP_RECON_PARAMS::detectionEfficiency, THREEDAP_RECON_PARAMS::ionVolume, THREEDAP_RECON_PARAMS::radiusCylinder, and readOps().

Referenced by DataLoadFilter::refresh().

Here is the caller graph for this function:

Variable Documentation

AMU_PER_KG

const float AMU_PER_KG = 6.0221409e26

CHARGE_ELECTRON

const float CHARGE_ELECTRON = 1.602176e-19

NUM_XY

const unsigned int NUM_XY =2

Referenced by readOps().

OPS_ENUM_ERRSTRINGS

const char* OPS_ENUM_ERRSTRINGS[]

Initial value:

=
{
        "\"C\" line (exp. parameters) not in expected format",
        "\"I\"  line (Detector parameters) not in expected format",
        "\"-\" line (Time-of-flights) not in expected format",
        "\"V\" line (voltage data) not in expected format",
        "Error interpreting \"V\" line (voltage data) data",
        "Missing beta value from \"V\" line - needs to be in voltage, or in system header (\"C\" line",
        "Incorrect number of \"V\" line (voltage data) data entries",
        "Unknown linetype in data file",
        "\"P\" line (detector channels) not in expected format",
        "Unable to interpret channels line",
        "Incorrect number of events in \"S\" (hit position) line",
        "Unable to interpret event data on \"S\" (hit position) line",
        "Unable to parse \"S\" (hit position) line event count (eg 8 in S8 ...)",
        "\"S\" (hit position) line data not preceded by TOF data, as it should have been",
        "Duplicate system data/experiment setup (\"C\") entry found -- there can only be one",
        "Duplicate detector (\"I\") data entry found -- there can only be one",
        "Trailing \"-\" line found -- should have be followed by a \"S\" line, but wasn't.",
        "Duplicate\"-\" line found -- should have be followed by a \"S\" line, but wasn't.",
        "Unable to open file",
        "unable to read file, after opening",
        "Abort requested"
}

Referenced by readOps(), and DataLoadFilter::refresh().