1 cdaq 1.1 *
2 * HMS calorimeter. Parameters from hms_positions.parm
3 *
|
4 cdaq 1.2 * $Log: hms_calorimeter.cmn,v $
|
5 saw 1.9 * Revision 1.8 1996/01/17 15:59:26 cdaq
6 * (JRA) Add normalized event quantities
7 *
|
8 cdaq 1.8 * Revision 1.7 1995/08/11 16:27:42 cdaq
9 * (JRA) Add accumulators for calorimeter
10 *
|
11 cdaq 1.7 * Revision 1.6 1995/05/22 19:06:29 cdaq
12 * (SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts
13 *
|
14 cdaq 1.6 * Revision 1.5 1995/03/13 19:01:46 cdaq
15 * (JRA) Move hmax_cal_rows and hmax_cal_columns to gen_data_structures
16 *
|
17 cdaq 1.5 * Revision 1.4 1995/01/04 17:09:32 cdaq
18 * (HGM) Increase hmax_cal_rows to 13
19 *
|
20 cdaq 1.4 * Revision 1.3 1994/10/11 19:44:44 cdaq
21 * (SAW) Add hcal_slop
22 *
|
23 cdaq 1.3 * Revision 1.2 1994/08/05 15:35:42 cdaq
24 * (SAW) Add makereg directive with required include files
25 * Add "CTPTYPE=event" directive for auto generation of CTP reg calls
26 *
|
27 cdaq 1.2 * Revision 1.1 1994/04/12 21:32:21 cdaq
28 * Initial revision
29 *
30 * The following include statments must precede the inclusion of this
31 * file in each routine that uses it. The *%% syntax is also a
32 * directive to makereg to tell it to include the code in the program
33 * that it generates.
|
34 cdaq 1.1 *
|
35 cdaq 1.6 *%% include 'hms_data_structures.cmn'
|
36 cdaq 1.2 *
37 * CTPTYPE=parm
38 *
|
39 cdaq 1.1 real*4 hcal_1pr_zpos ! Z positions of front of shower counter layers
40 real*4 hcal_2ta_zpos !
41 real*4 hcal_3ta_zpos !
42 real*4 hcal_4ta_zpos !
43 real*4 hcal_1pr_thick ! Thickness of shower counter blocks
44 real*4 hcal_2ta_thick !
45 real*4 hcal_3ta_thick !
46 real*4 hcal_4ta_thick !
47 integer*4 hcal_1pr_nr !Number of shower blocks per column
48 integer*4 hcal_2ta_nr !
49 integer*4 hcal_3ta_nr !
50 integer*4 hcal_4ta_nr !
51 real*4 hcal_1pr_left !Y position, column #1
52 real*4 hcal_1pr_right !
53 real*4 hcal_1pr_top ! X positions, column #1
54 real*4 hcal_2ta_left ! Y position, column #2
55 real*4 hcal_2ta_right !
56 real*4 hcal_2ta_top ! X positions, column#2
57 real*4 hcal_3ta_left ! Y position, column #3
58 real*4 hcal_3ta_right !
59 real*4 hcal_3ta_top ! X positions, column #3
60 cdaq 1.1 real*4 hcal_4ta_left ! Y position, column #4
61 real*4 hcal_4ta_right !
62 real*4 hcal_4ta_top ! X positions, column #4
|
63 cdaq 1.3 REAL*4 HCAL_SLOP !SLOP IN DISTANCE BETWEEN TRACK AND BLOCK
|
64 cdaq 1.1 common/hms_cal_parms/
65 & hcal_1pr_zpos,
66 & hcal_2ta_zpos,
67 & hcal_3ta_zpos,
68 & hcal_4ta_zpos,
69 & hcal_1pr_thick,
70 & hcal_2ta_thick,
71 & hcal_3ta_thick,
72 & hcal_4ta_thick,
73 & hcal_1pr_nr,
74 & hcal_2ta_nr,
75 & hcal_3ta_nr,
76 & hcal_4ta_nr,
77 & hcal_1pr_left,
78 & hcal_1pr_right,
|
79 saw 1.9 & hcal_1pr_top(HMAX_CAL_ROWS),
|
80 cdaq 1.1 & hcal_2ta_left,
81 & hcal_2ta_right,
|
82 saw 1.9 & hcal_2ta_top(HMAX_CAL_ROWS),
|
83 cdaq 1.1 & hcal_3ta_left,
84 & hcal_3ta_right,
|
85 saw 1.9 & hcal_3ta_top(HMAX_CAL_ROWS),
|
86 cdaq 1.1 & hcal_4ta_left,
87 & hcal_4ta_right,
|
88 saw 1.9 & hcal_4ta_top(HMAX_CAL_ROWS),
|
89 cdaq 1.3 & hcal_slop
|
90 cdaq 1.1 *
91 * HMS calorimeter. Geometrical constants filled by h_init_cal
92 *
93 real*4 hcal_block_xsize !
94 real*4 hcal_block_ysize !Block dimensions - 10*70*10 cm^3
95 real*4 hcal_block_zsize !
96 real*4 hcal_block_xc !
97 real*4 hcal_block_yc !X,Y,Z coordinates of block centers
98 real*4 hcal_block_zc !
99 real*4 hcal_xmin !
100 real*4 hcal_xmax !
101 real*4 hcal_ymin !Boundaries of the HMS
102 real*4 hcal_ymax !calorimeter stack
103 real*4 hcal_zmin !
104 real*4 hcal_zmax !
105 real*4 hcal_fv_xmin !
106 real*4 hcal_fv_xmax !
107 real*4 hcal_fv_ymin !Boundaries of the
108 real*4 hcal_fv_ymax !fiducial volume
109 real*4 hcal_fv_zmin !
110 real*4 hcal_fv_zmax !
111 cdaq 1.1 common/hms_geometry_cal/
112 & hcal_block_xsize,
113 & hcal_block_ysize,
114 & hcal_block_zsize,
|
115 saw 1.9 & hcal_block_xc(HMAX_CAL_BLOCKS),
116 & hcal_block_yc(HMAX_CAL_BLOCKS),
117 & hcal_block_zc(HMAX_CAL_BLOCKS),
|
118 cdaq 1.1 & hcal_xmin,hcal_xmax,
119 & hcal_ymin,hcal_ymax,
120 & hcal_zmin,hcal_zmax,
121 & hcal_fv_xmin,hcal_fv_xmax,
122 & hcal_fv_ymin,hcal_fv_ymax,
123 & hcal_fv_zmin,hcal_fv_zmax
124 *
125 * HMS calorimeter. Sparsified data filled by h_sparsify_cal
126 *
|
127 cdaq 1.2 * CTPTYPE=event
128 *
|
129 cdaq 1.1 integer*4 hcal_rows !Row number. Copied from hms_raw_cal
130 integer*4 hcal_cols !Column number. Copied from hms_raw_cal
131 real*4 hcal_adcs !Pulse height in channels(ADC-PED).
132 integer*4 hcal_num_hits !Total number of hits above threshold
133 common/hms_sparsified_cal/
|
134 saw 1.9 & hcal_rows(HMAX_CAL_BLOCKS),
135 & hcal_cols(HMAX_CAL_BLOCKS),
136 & hcal_adcs(HMAX_CAL_BLOCKS),
|
137 cdaq 1.1 & hcal_num_hits
138 *
139 * HMS calorimeter. Cluster data, filled by h_clusters_cal
140 *
141 integer*4 hnclusters_max !Number of clusters allowed in the calorimeter
142 parameter (hnclusters_max=5) !Set the maximum to 5
143 integer*4 hnclusters_cal !Number of clusters found
144 integer*4 hcluster_hit !Link pointer to cluster number for each hit
145 integer*4 hcluster_size !Number of hits in a cluster
146 real*4 hcluster_xc !X-coordinate of a cluster
147 real*4 hcluster_e1 !Energy deposition in column #1
148 real*4 hcluster_e2 ! #2
149 real*4 hcluster_e3 ! #3
150 real*4 hcluster_e4 ! #4
151 real*4 hcluster_et !Total energy deposition
152 common/hms_clusters_cal/
|
153 saw 1.9 & hcluster_hit(HMAX_CAL_BLOCKS),
|
154 cdaq 1.1 & hcluster_size(hnclusters_max),
155 & hcluster_xc(hnclusters_max),
156 & hcluster_e1(hnclusters_max),
157 & hcluster_e2(hnclusters_max),
158 & hcluster_e3(hnclusters_max),
159 & hcluster_e4(hnclusters_max),
160 & hcluster_et(hnclusters_max),
161 & hnclusters_cal
162 *
163 * HMS calorimeter. Calorimeter track quantities,filled by h_tracks_cal
164 *
165 real*4 htrack_xc !X,Y position of track on
166 real*4 htrack_yc ! calorimeter front surface
167 integer*4 hcluster_track !Link pointer to calorimeter cluster number
168 integer*4 hntracks_cal !Number of tracks for which a cluster was found
169 common/hms_tracks_cal/
|
170 saw 1.9 & htrack_xc(HNTRACKS_MAX),
171 & htrack_yc(HNTRACKS_MAX),
172 & hcluster_track(HNTRACKS_MAX),
|
173 cdaq 1.1 & hntracks_cal
|
174 cdaq 1.2 *
175 * CTPTYPE=parm
176 *
|
177 cdaq 1.1 *
178 * HMS calorimeter. ADC pedestals and thresholds
179 *
180 real*4 hcal_ped_mean !Mean pedestal value
181 real*4 hcal_ped_rms !Pedestal rms deviation
182 real*4 hcal_threshold !Threshold=3.*hcal_ped_rms
183 common/hms_cal_pedestals/
|
184 saw 1.9 & hcal_ped_mean(HMAX_CAL_BLOCKS),
185 & hcal_ped_rms(HMAX_CAL_BLOCKS),
186 & hcal_threshold(HMAX_CAL_BLOCKS)
|
187 cdaq 1.1 *
188 * HMS calorimeter. Calibration constants
189 *
190 real*4 hcal_cal_const !Calibration constants
191 common/hms_cal_const/
|
192 saw 1.9 & hcal_cal_const(HMAX_CAL_BLOCKS)
|
193 cdaq 1.1 *
194 * HMS calorimeter. Relative gains & correction factors
195 *
196 real*4 hcal_gain_ini !Relative gains during the last calibration
197 real*4 hcal_gain_cur !Current relative gains
198 real*4 hcal_gain_cor !Correction factor: cor=ini/cur
199 common/hms_cal_monitor/
|
200 saw 1.9 & hcal_gain_ini(HMAX_CAL_BLOCKS),
201 & hcal_gain_cur(HMAX_CAL_BLOCKS),
202 & hcal_gain_cor(HMAX_CAL_BLOCKS)
|
203 cdaq 1.1 *
204 * HMS calorimeter. Debuging LUN and flags
205 *
206 integer*4 hlun_dbg_cal
207 integer*4 hdbg_raw_cal
208 integer*4 hdbg_sparsified_cal
209 integer*4 hdbg_decoded_cal
210 integer*4 hdbg_clusters_cal
211 integer*4 hdbg_tracks_cal
212 integer*4 hdbg_tests_cal
213 common/hms_cal_flags/
214 & hlun_dbg_cal,
215 & hdbg_raw_cal,
216 & hdbg_sparsified_cal,
217 & hdbg_decoded_cal,
218 & hdbg_clusters_cal,
219 & hdbg_tracks_cal,
220 & hdbg_tests_cal
|
221 cdaq 1.7
222 *
223 * CTPTYPE=event
224 *
|
225 saw 1.9 integer*4 hcal_zero_sum(HMAX_CAL_BLOCKS)
226 integer*4 hcal_zero_sum2(HMAX_CAL_BLOCKS)
227 integer*4 hcal_zero_num(HMAX_CAL_BLOCKS)
228 real*4 hcal_zero_ave(HMAX_CAL_BLOCKS)
229 real*4 hcal_zero_sig(HMAX_CAL_BLOCKS)
230 real*4 hcal_zero_thresh(HMAX_CAL_BLOCKS)
|
231 cdaq 1.7
232 common /hms_cal_zero/
233 & hcal_zero_ave,
234 & hcal_zero_sig,
235 & hcal_zero_thresh,
236 & hcal_zero_num,
237 & hcal_zero_sum,
238 & hcal_zero_sum2
|
239 cdaq 1.8
240 *
241 * CTPTYPE=event
242 *
243 real*4 hscal_suma !normalized sum of layer A.
244 real*4 hscal_sumb !normalized sum of layer B.
245 real*4 hscal_sumc !normalized sum of layer C.
246 real*4 hscal_sumd !normalized sum of layer D.
247 real*4 hsprsum !normalized PR sum.
248 real*4 hsshsum !normalized total sum.
249 real*4 hsprtrk !normalized PR sum on track.
250 real*4 hsshtrk !normalized total sum on track.
251
252 common /hms_cal_normalized/
253 & hscal_suma,
254 & hscal_sumb,
255 & hscal_sumc,
256 & hscal_sumd,
257 & hsprsum,
258 & hsshsum,
259 & hsprtrk,
260 cdaq 1.8 & hsshtrk
|