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
|
#!/usr/bin/perl -w
use strict;
package GetTimeObs;
sub main
{
my ($fh,$site,$file,$whatState,$whatObservable)=@_;
my $labelDoubleOcc = "site nupNdown(gs) nupNdown(timevector) time";
my $labelDensity = "site nUp+nDown(gs) nup+ndown(timevector) time";
my $label = $labelDoubleOcc;
$label = $labelDensity if ($whatObservable eq "density");
defined($site) or die "$0: Undefined site\n";
defined($file) or die "$0: Undefined file\n";
defined($whatState) or die "$0: Must specify what state gs or time\n";
defined($whatObservable) or die "$0: Must specify what observable (density or nd)\n";
my $stateIndex = 2;
$stateIndex = 1 if ($whatState eq "gs");
my $sd = getSuperDensity($site,$file);
$sd = 1 if ($whatState eq "gs");
open(FILE,$file) or die "Cannot open file $file: $!\n";
while(<FILE>) {
if (/\Q${label}/i) {
last;
}
}
my $prevT = -1;
while(<FILE>) {
next if (/^VectorWithOffsets/);
last if (/^#/);
if (/SuperDensity.*=\(([^,]+),/) {
$sd = $1 unless ($whatState eq "gs");
next;
}
my @temp=split;
last unless $temp[0]=~/^(\d+\.?\d*|\.\d+)$/; # match valid number
if ($temp[0]==$site) {
my $val = $temp[$stateIndex];
$val =~ s/\(//;
$val =~ s/,.*\)//;
next if ($prevT == $temp[3]);
$prevT = $temp[3];
$val /= $sd;
print $fh "$temp[3] $val\n";
}
}
close(FILE);
}
sub getSuperDensity
{
my ($site,$file)=@_;
my $sd;
open(FILE,$file) or die "Cannot open file $file: $!\n";
while(<FILE>) {
if (/SuperDensity.*=\(([^,]+),/) {
$sd = $1;
last;
}
}
close(FILE);
defined $sd or die "SuperDensity is not defined\n";
return $sd;
}
1;
|
