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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
|
package Demeter::Plugins::SSRLA; # -*- cperl -*-
use Moose;
extends 'Demeter::Plugins::FileType';
use Demeter::Constants qw{$EPSILON3};
has '+is_binary' => (default => 0);
has '+description' => (default => "ASCII data from the SSRL XAFS Data Collector");
has '+version' => (default => 0.2);
my %special = (chr(169) => '(c)',
chr(176) => 'deg',
);
my $re = '[' . join('|', keys(%special)) . ']';
sub is {
my ($self) = @_;
my $null = chr(0);
open D, $self->file or $self->Croak("could not open " . $self->file . " as data (SSRL ASCII)\n");
my $line = <D>;
my $is_ssrl = ($line =~ m{^\s*SSRL\s+EXAFS Data Collector});
$line = <D>;
my $is_ascii = ($line !~ m{$null});
close D;
return 1 if ($is_ssrl and $is_ascii);
return 0;
};
sub fix {
my ($self) = @_;
my $new = File::Spec->catfile($self->stash_folder, $self->filename);
($new = File::Spec->catfile($self->stash_folder, "toss")) if (length($new) > 127);
open D, $self->file or die "could not open " . $self->file . " as data (fix in SSRLA)\n";
open N, ">".$new or die "could not write to $new (fix in SSRLA)\n";
my @labels;
my @offsets;
my ($header, $labels) = (1, 0);
while (<D>) {
chomp;
if ($_ =~ /^\s*Data:/) {
(($header, $labels) = (0,1));
next;
};
if (($_ =~ /^\s*$/) and $labels) {
(($header, $labels) = (0,0));
@labels = ($labels[2], $labels[1], $labels[0], @labels[3..$#labels]);
print N "# ", "-"x30, $/;
print N "# ", join(" ", @labels), $/;
next;
};
if ($labels) {
my $this = $_;
$this =~ s/\s+$//;
$this =~ s/\s+/_/g;
push @labels, $this;
} elsif ($header) { # comment header
if ($_ =~ /^\s*Offsets/) {
print N "# ", $_, $/;
my $line = <D>;
@offsets = split(" ", $line);
@offsets = ($offsets[2], $offsets[1], $offsets[0], @offsets[3..$#offsets]);
print N "# ", join(" ", @offsets), $/;
} else {
my $text = $_;
$text =~ s{($re)}{$special{$1}}g; # files from ROBL contain © and ° characters
print N "# ", $text, $/;
};
} else { # data columns
my @line = split(" ", $_);
next if $#line < 0;
next if $line[2] < $EPSILON3;
@line = ($line[2], $line[1], $line[0], @line[3..$#line]);
my $nn = $#line+1;
my $pattern = "%.4f " x $nn . $/;
#@line = map{ $line[$_] - $offsets[$_] } (0 .. $#line);
printf N $pattern, @line;
};
};
close N;
close D;
$self->fixed($new);
return $new;
}
sub suggest {
my ($self, $which) = @_;
$which ||= 'transmission';
if ($which eq 'transmission') {
return (energy => '$1',
numerator => '$4',
denominator => '$5',
ln => 1,);
} else {
return (energy => '$1',
numerator => '$6',
denominator => '$4',
ln => 0,);
};
};
__PACKAGE__->meta->make_immutable;
1;
=head1 NAME
Demeter::Plugin::SSRLA - SSRL XAFS Data Collector 1.3 ASCII filetype plugin
=head1 VERSION
This documentation refers to Demeter version 0.9.26.
=head1 SYNOPSIS
This plugin directly reads the files written by the SSRL XAFS Data Collector.
=head1 SSRL files
This plugin comments out the header lines, constructs a column label
line out of the Data: section, moves the first column (real time
clock) to the third column, and swaps the requested and acheived
energy columns.
=head1 ROBL files
The ROBL beamline at ESRF uses the same file format as the SSRL Data
Collector 1.1 except that it includes a couple of high-ASCII
characters. These characters may give a Wx::TextCtrl some trouble, so
they get stripped from the file and replaced with ASCII look-alikes.
© ==> (c)
° ==> deg
=head1 VERSION HISTORY
=over 4
=item 0.2
Allow variable lines after the C<Data:> header. Remove trailing lines
of all zeros from the file.
=back
=head1 AUTHOR
Bruce Ravel, L<http://bruceravel.github.io/home>
http://bruceravel.github.io/demeter
=cut
|
