1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
/*
* beam-parameters.c
*
* Beam parameters
*
* (c) 2006-2010 Thomas White <taw@physics.org>
*
* Part of CrystFEL - crystallography with a FEL
*
*/
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include <stdio.h>
#include <stdlib.h>
#include "beam-parameters.h"
#include "utils.h"
struct beam_params *get_beam_parameters(const char *filename)
{
FILE *fh;
struct beam_params *b;
char *rval;
int reject;
fh = fopen(filename, "r");
if ( fh == NULL ) return NULL;
b = calloc(1, sizeof(struct beam_params));
if ( b == NULL ) return NULL;
b->fluence = -1.0;
b->beam_radius = -1.0;
b->photon_energy = -1.0;
b->dqe = -1.0;
b->adu_per_photon = -1.0;
b->water_radius = -1.0;
do {
int n1;
char line[1024];
int i;
char **bits;
rval = fgets(line, 1023, fh);
if ( rval == NULL ) break;
chomp(line);
n1 = assplode(line, " \t", &bits, ASSPLODE_NONE);
if ( n1 < 3 ) {
for ( i=0; i<n1; i++ ) free(bits[i]);
free(bits);
continue;
}
if ( bits[1][0] != '=' ) {
for ( i=0; i<n1; i++ ) free(bits[i]);
free(bits);
continue;
}
if ( strcmp(bits[0], "beam/fluence") == 0 ) {
b->fluence = atof(bits[0]);
} else if ( strcmp(bits[0], "beam/radius") == 0 ) {
b->beam_radius = atof(bits[0]);
} else if ( strcmp(bits[0], "beam/photon_energy") == 0 ) {
b->photon_energy = atof(bits[0]);
} else if ( strcmp(bits[0], "detector/dqe") == 0 ) {
b->dqe = atof(bits[0]);
} else if ( strcmp(bits[0], "detector/adu_per_photon") == 0 ) {
b->adu_per_photon = atof(bits[0]);
} else if ( strcmp(bits[0], "jet/radius") == 0 ) {
b->water_radius = atof(bits[0]);
} else {
ERROR("Unrecognised field '%s'\n", bits[0]);
}
for ( i=0; i<n1; i++ ) free(bits[i]);
free(bits);
} while ( rval != NULL );
fclose(fh);
reject = 0;
if ( b->fluence < 0.0 ) {
ERROR("Invalid or unspecified value for 'beam/fluence'.\n");
reject = 1;
}
if ( b->beam_radius < 0.0 ) {
ERROR("Invalid or unspecified value for 'beam/radius'.\n");
reject = 1;
}
if ( b->photon_energy < 0.0 ) {
ERROR("Invalid or unspecified value for"
" 'beam/photon_energy'.\n");
reject = 1;
}
if ( b->dqe < 0.0 ) {
ERROR("Invalid or unspecified value for 'detector/dqe'.\n");
reject = 1;
}
if ( b->adu_per_photon < 0.0 ) {
ERROR("Invalid or unspecified value for"
" 'detector/adu_per_photon'.\n");
reject = 1;
}
if ( b->water_radius < 0.0 ) {
ERROR("Invalid or unspecified value for 'jet/radius'.\n");
reject = 1;
}
if ( reject ) {
ERROR("Please fix the above problems with the beam"
" parameters file and try again.\n");
return NULL;
}
return b;
}
|
