1 |
* DATA ZTOF/55.,25.,-23/ |
2 |
* |
3 |
* TOF data |
4 |
* first index : 1 = A, 2 = B |
5 |
* second index : 1.. number of paddle |
6 |
* third index : 1 = TDC, 2 = ADC |
7 |
|
8 |
c INTEGER tof11(2,8,2),tof12(2,6,2) |
9 |
c INTEGER tof21(2,2,2),tof22(2,2,2) |
10 |
c INTEGER tof31(2,3,2),tof32(2,3,2) |
11 |
REAL tof11(2,8,2),tof12(2,6,2) |
12 |
REAL tof21(2,2,2),tof22(2,2,2) |
13 |
REAL tof31(2,3,2),tof32(2,3,2) |
14 |
|
15 |
c Normalized ADC value |
16 |
REAL adc_c(4,12) |
17 |
|
18 |
c k1 constans for calculation of beta |
19 |
REAL k1_s11s31(24),k1_s12s32(18),k1_S21S31(6),k1_S22S32(6) |
20 |
|
21 |
c middle y (or x) position of the upper and middle ToF-Paddle |
22 |
c to calculate theta-angle without track information |
23 |
REAL tof11_x(8), tof21_y(2), tof31_x(3) |
24 |
REAL tof12_y(6), tof22_x(2), tof32_y(3) |
25 |
|
26 |
c relation of time differences in a tof paddle and the location of the |
27 |
c passage of the particle |
28 |
c first index : 1... number of paddle |
29 |
c second index : 1 = offset, 2 = slope |
30 |
REAL y_coor_lin11(8,2),x_coor_lin12(6,2) |
31 |
REAL x_coor_lin21(2,2),y_coor_lin22(2,2) |
32 |
REAL y_coor_lin31(3,2),x_coor_lin32(3,2) |
33 |
|
34 |
c ADC(x) - correction |
35 |
c first index : 1 = left, 2 = right |
36 |
c second index : 1..8 number of paddle |
37 |
c third index: parameter as follows |
38 |
c left: xkorr=par(1)*exp(-x/par(2)) |
39 |
c right:xkorr=par(1)*exp(x/par(2)) |
40 |
c then:adc(x,right/left)=adc(x,right/left)/xkorr |
41 |
c normalize to x=0 (middle of paddle) |
42 |
REAL adcx11(2,8,2),adcx12(2,6,2) |
43 |
REAL adcx21(2,2,2),adcx22(2,2,2) |
44 |
REAL adcx31(2,3,2),adcx32(2,3,2) |
45 |
|
46 |
|
47 |
c Time Walk - correction |
48 |
c first index : 1 = left, 2 = right |
49 |
c second index : 1..8 number of paddle |
50 |
REAL tw11(2,8),tw12(2,6) |
51 |
REAL tw21(2,2),tw22(2,2) |
52 |
REAL tw31(2,3),tw32(2,3) |
53 |
|
54 |
|
55 |
c value for status of each PM-data |
56 |
c first index : 1 = left, 2 = right |
57 |
c second index : 1... number of paddle |
58 |
INTEGER tof11_event(2,8),tof12_event(2,6) |
59 |
INTEGER tof21_event(2,2),tof22_event(2,2) |
60 |
INTEGER tof31_event(2,3),tof32_event(2,3) |
61 |
|
62 |
C---- |
63 |
|
64 |
INTEGER left,right,top,bot,mid,itdc,iadc |
65 |
INTEGER tdc_ev, adc_ev, none_ev, none_find |
66 |
INTEGER offset, slope |
67 |
INTEGER tof11_i,tof21_i,tof31_i,tof11_j,tof21_j,tof31_j |
68 |
INTEGER tof12_i,tof22_i,tof32_i,tof12_j,tof22_j,tof32_j |
69 |
* REAL adcoff(4,2,52),a,b |
70 |
c REAL delta_x12,delta_y12,delta_x23,delta_y23,delta_x13,delta_y13 |
71 |
REAL xtofpos(3),ytofpos(3) |
72 |
REAL xkorr, secure |
73 |
INTEGER i, j, ihelp |
74 |
LOGICAL check |
75 |
REAL theta12,theta13,theta23 |
76 |
|
77 |
REAL beta_a(5) |
78 |
|
79 |
INTEGER tof_i_flag(6),tof_j_flag(6) |
80 |
|
81 |
|
82 |
DATA TOFfst / 0/ |
83 |
|
84 |
PARAMETER (tofarm23 = 45.) |
85 |
PARAMETER (tofarm13 = 81.) |
86 |
PARAMETER (c2_13=108.) ! = 2.*tofarm13*0.01/(3.E08*50.E-12) |
87 |
PARAMETER (c2_23=60.) ! = 2.*tofarm23*0.01/(3.E08*50.E-12) |
88 |
|
89 |
|
90 |
* S11 8 paddles 33.0 x 5.1 cm |
91 |
* S12 6 paddles 40.8 x 5.5 cm |
92 |
* S21 2 paddles 18.0 x 7.5 cm |
93 |
* S22 2 paddles 15.0 x 9.0 cm |
94 |
* S31 3 paddles 15.0 x 6.0 cm |
95 |
* S32 3 paddles 18.0 x 5.0 cm |
96 |
|
97 |
DATA tof11_x/ -17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85/ |
98 |
DATA tof12_y/ -13.75,-8.25,-2.75,2.75,8.25,13.75/ |
99 |
DATA tof21_y/ -3.75,3.75/ |
100 |
DATA tof22_x/ -4.5,4.5/ |
101 |
DATA tof31_x/ -6.0,0.,6.0/ |
102 |
DATA tof32_y/ -5.0,0.0,5.0/ |
103 |
|
104 |
integer ch11a(8),hb11a(8),ch11b(8),hb11b(8) |
105 |
integer ch12a(6),hb12a(6),ch12b(6),hb12b(6) |
106 |
integer ch21a(2),hb21a(2),ch21b(2),hb21b(2) |
107 |
integer ch22a(2),hb22a(2),ch22b(2),hb22b(2) |
108 |
integer ch31a(3),hb31a(3),ch31b(3),hb31b(3) |
109 |
integer ch32a(3),hb32a(3),ch32b(3),hb32b(3) |
110 |
|
111 |
|
112 |
data ch11a / 4, 4, 4, 4, 1, 1, 2, 2/ |
113 |
data hb11a / 1, 3, 5, 7, 10, 12, 2, 4/ |
114 |
data ch11b / 2, 2, 2, 2, 1, 1, 1, 1/ |
115 |
data hb11b / 6, 8, 12, 10, 8, 6, 4, 2/ |
116 |
|
117 |
data ch12a / 3, 3, 3, 3, 3, 3/ |
118 |
data hb12a / 2, 4, 6, 8, 10, 12/ |
119 |
data ch12b / 4, 4, 4, 4, 4, 4/ |
120 |
data hb12b / 12, 10, 8, 6, 4, 2/ |
121 |
|
122 |
data ch21a / 1, 1/ |
123 |
data hb21a / 1, 5/ |
124 |
data ch21b / 2, 1/ |
125 |
data hb21b / 9, 7/ |
126 |
|
127 |
data ch22a / 1, 1/ |
128 |
data hb22a / 3, 9/ |
129 |
data ch22b / 2, 1/ |
130 |
data hb22b / 11, 11/ |
131 |
|
132 |
data ch31a / 2, 3, 3/ |
133 |
data hb31a / 7, 9, 11/ |
134 |
data ch31b / 2, 2, 2/ |
135 |
data hb31b / 5, 3, 1/ |
136 |
|
137 |
data ch32a / 3, 3, 4/ |
138 |
data hb32a / 1, 5, 9/ |
139 |
data ch32b / 3, 3, 4/ |
140 |
data hb32b / 3, 7, 11/ |
141 |
|
142 |
common / tofcalc /k1_s11s31,k1_s12s32,k1_S21S31, |
143 |
& k1_S22S32,beta_a, |
144 |
& adcx11,adcx12,adcx21,adcx22,adcx31,adcx32, |
145 |
& tw11,tw12,tw21,tw22,tw31,tw32, |
146 |
& y_coor_lin11,x_coor_lin12,x_coor_lin21, |
147 |
& y_coor_lin22,y_coor_lin31,x_coor_lin32, |
148 |
& xtofpos,ytofpos,adc_c,tof_i_flag,tof_j_flag |