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
|
#!/usr/bin/perl -w
use strict;
open(ROT, "< rotations.lst");
open(COORD, "< coords.dat");
open(OUT, "> cspad-feb2011.geom");
my $quadrant = 0;
printf(OUT "n_panels = 64\n\n");
my $x = 0;
my $y = 0;
my $p = 0;
my $global_x = 0.0;
my $global_y = 0.0;
while ( my $coord = <COORD> ) {
my $rot = <ROT>;
my $cx = 0.0;
my $cy = 0.0;
my $minx = $x*194;
my $miny = $y*185;
my $sc1;
my $sc2;
chomp $coord;
chomp $rot;
if ( $coord =~ /\(([0-9\.\-]+),\s([0-9\.\-]+)\)/ ) {
$sc1 = $1;
$sc2 = $2;
printf("%f %f\n", $sc1, $sc2);
} else {
printf("!!!\n");
}
my $a = 388.0; # Fast scan
my $b = 185.0; # Slow scan
# FIXME: Remove the $b for versions of Stephan's code after 20th Feb
$sc1 -= (1308.696-$b);
$sc2 -= (980.3862-$b);
my $sx = -$sc2;
my $sy = -$sc1;
printf(OUT "; Quadrant %i, asic %i\n", $quadrant, ($x%2)+2*$y);
printf(OUT "%i/min_fs = %i\n", $p, $minx);
printf(OUT "%i/min_ss = %i\n", $p, $miny);
printf(OUT "%i/max_fs = %i\n", $p, ($x+1)*194-1);
printf(OUT "%i/max_ss = %i\n", $p, ($y+1)*185-1);
printf(OUT "%i/badrow_direction = -\n", $p);
printf(OUT "%i/res = 9090.91\n", $p);
printf(OUT "%i/peak_sep = 6.0\n", $p);
printf(OUT "%i/clen = 77.0e-3\n", $p);
if ( $rot == "0" ) {
printf(OUT "%i/fs = -x\n", $p);
printf(OUT "%i/ss = -y\n", $p);
$cx = $sx;
$cy = $sy;
} elsif ( $rot == "90" ) {
printf(OUT "%i/fs = +y\n", $p);
printf(OUT "%i/ss = -x\n", $p);
$cx = $sx;
$cy = $sy - $a - 5.0;
} elsif ( $rot == "180" ) {
printf(OUT "%i/fs = +x\n", $p);
printf(OUT "%i/ss = +y\n", $p);
$cx = $sx - $a - 5.0;
$cy = $sy - $b;
} elsif ( $rot == "270" ) {
printf(OUT "%i/fs = -y\n", $p);
printf(OUT "%i/ss = +x\n", $p);
$cx = $sx - $b;
$cy = $sy;
}
printf(OUT "%i/corner_x = %5.2f\n", $p, $cx+$global_x);
printf(OUT "%i/corner_y = %5.2f\n", $p, $cy+$global_y);
printf(OUT "%i/no_index = 0\n", $p);
printf(OUT "\n");
#if ( ($y < 20) && ($x < 1) ) {
printf(STDERR "%f %f %f %f\n", $cx, $cy, $sc1, $sc2);
#}
$x++;
$p++;
printf(OUT "; Quadrant %i, asic %i\n", $quadrant, ($x%2)+2*$y);
printf(OUT "%i/min_fs = %i\n", $p, $x*194);
printf(OUT "%i/min_ss = %i\n", $p, $y*185);
printf(OUT "%i/max_fs = %i\n", $p, ($x+1)*194-1);
printf(OUT "%i/max_ss = %i\n", $p, ($y+1)*185-1);
printf(OUT "%i/badrow_direction = -\n", $p);
printf(OUT "%i/res = 9090.91\n", $p);
printf(OUT "%i/peak_sep = 6.0\n", $p);
printf(OUT "%i/clen = 77.0e-3\n", $p);
if ( $rot == "0" ) {
printf(OUT "%i/fs = -x\n", $p);
printf(OUT "%i/ss = -y\n", $p);
$cx = $sx - $a/2.0 - 5.0;
$cy = $sy;
} elsif ( $rot == "90" ) {
printf(OUT "%i/fs = +y\n", $p);
printf(OUT "%i/ss = -x\n", $p);
$cx = $sx;
$cy = $sy - $a/2.0;
} elsif ( $rot == "180" ) {
printf(OUT "%i/fs = +x\n", $p);
printf(OUT "%i/ss = +y\n", $p);
$cx = $sx - $a/2.0;
$cy = $sy - $b;
} elsif ( $rot == "270" ) {
printf(OUT "%i/fs = -y\n", $p);
printf(OUT "%i/ss = +x\n", $p);
$cx = $sx - $b;
$cy = $sy - $a/2.0 - 5.0;
}
printf(OUT "%i/corner_x = %5.2f\n", $p, $cx+$global_x);
printf(OUT "%i/corner_y = %5.2f\n", $p, $cy+$global_y);
printf(OUT "%i/no_index = 0\n", $p);
printf(OUT "\n");
#if ( ($y < 20) && ($x < 1) ) {
printf(STDERR "%f %f %f %f\n", $cx, $cy, $sc1, $sc2);
#}
$x++;
$p++;
if ( $x == 8 ) {
$x = 0;
$y++;
}
$quadrant++;
if ( $quadrant == 4 ) {
$quadrant = 0;
}
}
|
