1 cdaq 1.1 *****************begin: hms_data_structures.cmn*************************
2 *
3 * include file hms_data_structures.cmn
4 *
|
5 jones 1.9 * Author: D. F. Geesaman 1 September 1993
|
6 cdaq 1.1 *
|
7 cdaq 1.2 * $Log: hms_data_structures.cmn,v $
|
8 cdaq 1.12.12.2 * Revision 1.12.12.1 2004/06/28 17:33:40 cdaq
9 * new variables like x_si, etc.
10 *
|
11 cdaq 1.12.12.1 * Revision 1.12 2003/09/05 20:39:18 jones
12 * Merge in online03 changes (mkj)
13 *
|
14 jones 1.12 * Revision 1.11 2003/04/01 13:55:08 jones
15 * Add variables hntracks_max_fp and h_remove_sppt_if_one_y_plane to
16 * hms_tracking.cmn
17 *
18 * Revision 1.10.2.2 2003/07/15 19:04:03 cdaq
19 * add hsinplane
20 *
21 * Revision 1.10.2.1 2003/04/10 12:36:56 cdaq
22 * comments added
23 *
|
24 jones 1.11 * Revision 1.10 2003/01/02 18:34:51 jones
25 * change HNTRACKS_MAX from 5 to 20
26 *
|
27 jones 1.10 * Revision 1.9 2002/12/20 21:52:34 jones
28 * Modified by Hamlet for new HMS aerogel
29 *
|
30 jones 1.9 *
31 * Revision 1.9 2002/09/26
32 * (Hamlet) Add structures for HMS Aerogel
33 *
34 * Revision 1.8 1999/02/23 19:09:29 csa
35 * Change some physics vars
36 *
|
37 csa 1.8 * Revision 1.7 1999/01/21 21:40:31 saw
38 * Extra shower counter tube modifications
39 *
|
40 saw 1.7 * Revision 1.6 1998/12/17 22:02:41 saw
41 * Support extra set of tubes on HMS shower counter
42 *
|
43 saw 1.6 * Revision 1.5 1996/09/04 15:49:31 saw
44 * (SAW) Double of all ' 's in comments
45 *
|
46 saw 1.5 * Revision 1.4 1996/01/24 16:18:59 saw
47 * (JRA) Cosmetic change
48 *
|
49 saw 1.4 * Revision 1.3 1996/01/17 16:00:02 cdaq
50 * (JRA) ADD HCER_RAW_ADC and HSCIN_FPTIME
51 *
|
52 cdaq 1.3 * Revision 1.2 1995/08/11 16:28:46 cdaq
53 * (CC) Add structure for # of photoelectrons in Cerenkov
54 *
|
55 cdaq 1.2 * Revision 1.1 1995/05/22 18:41:15 cdaq
56 * Initial revision
|
57 cdaq 1.1 *
|
58 cdaq 1.2 *
|
59 cdaq 1.1 * CTPTYPE=parm
60 *
61 * HTHETA_LAB HMS LAB ANGLE THETA (RADIANS)
62 * HPHI_LAB HMS LAB ANGLE PHI
63 * HPCENTRAL HMS CENTRAL MOMENUTM (GEV)
64 * HBFIELD HMS B FIELD INCLUDING SIGN
65 *
66 REAL*4 HTHETA_LAB ! HMS LAB ANGLE THETA (RADIANS)
67 REAL*4 HPHI_LAB ! HMS LAB ANGLE PHI
68 REAL*4 HPCENTRAL ! HMS CENTRAL MOMENUTM (GEV)
|
69 cdaq 1.12.12.2 REAL*4 HPCENTRAL_SET ! HMS CENTRAL MOMENUTM (GEV) (no offset)
|
70 cdaq 1.1 REAL*4 HBFIELD ! HMS B FIELD INCLUDING SIGN
71 REAL*4 HPARTMASS ! EXPECTED MASS OF DETECTED PARTICLE IN HMS
|
72 jones 1.9
|
73 cdaq 1.1 COMMON/HMS_SPECTROMETER/
74 & HTHETA_LAB,
75 & HPHI_LAB,
76 & HPCENTRAL,
|
77 cdaq 1.12.12.2 & HPCENTRAL_SET,
|
78 cdaq 1.1 & HBFIELD,
79 & HPARTMASS
80 *
81 * CTPTYPE=event
82 *
83 * HMS DECODED DATA
84 *
85 * FILLED BY G_decode_event_by_banks
86 * HMS DRIFT CHAMBER HITS
87 * EACH CHAMBER HIT (TDC VALUE) HAS A
88 * PLANE NUMBER
89 * WIRE NUMBER
90 * TDC VALUE
91 *
92 * THE TOTAL NUMBER OF HITS IS GIVEN IN HDC_RAW_TOT_HITS
93 *
94 INTEGER*4 HMAX_DC_HITS ! MAXIUM NUMBER OF DRIFT CHAMBER HITS
95 INTEGER*4 HMAX_NUM_DC_PLANES ! MAX NUMBER OF HMS DRIFT CHAMBER PLANES
96 INTEGER*4 HMAX_NUM_CHAMBERS ! NUMBER OF HMS DRIFT CHAMBERS
97 PARAMETER(HMAX_DC_HITS=3600)
98 PARAMETER(HMAX_NUM_DC_PLANES=12)
99 cdaq 1.1 PARAMETER(HMAX_NUM_CHAMBERS=2)
100 INTEGER*4 HDC_RAW_PLANE_NUM
101 INTEGER*4 HDC_RAW_WIRE_NUM
102 INTEGER*4 HDC_RAW_TDC
103 INTEGER*4 HDC_RAW_TOT_HITS
104 COMMON/HMS_RAW_DC/
105 1 HDC_RAW_PLANE_NUM(HMAX_DC_HITS),
106 2 HDC_RAW_WIRE_NUM(HMAX_DC_HITS),
107 3 HDC_RAW_TDC(HMAX_DC_HITS),
108 4 HDC_RAW_TOT_HITS
109 *
110 * HMS DECODED DC DATA
111 *
112 * FILLED BY H_TRANS_DC TRANSLATION ROUTINE
113 * HMS DRIFT CHAMBER HITS
114 * EACH CHAMBER HIT (TDC VALUE) HAS A
115 * PLANE NUMBER
116 * WIRE NUMBER
117 * TDC VALUE
118 * DRIFT TIME
119 * DRIFT DISTANCE
120 cdaq 1.1 * GENERALIZED COORDINATE OF HIT WIRE CENTER
121 * GENERALIZED COORDINATE OF HIT PERPENDICULAR TO WIRE DIRECTION
122 *
123 * THE TOTAL NUMBER OF HITS IN EACH PLANE IS GIVEN IN HDC_HITS_PER_PLANE(I)
124 * THE TOTAL NUMBER OF HITS IS GIVEN IN HDC_TOT_HITS
125 *
126 INTEGER*4 HDC_PLANE_NUM ! copied from HMS_RAW_DC
127 INTEGER*4 HDC_WIRE_NUM ! copied from HMS_RAW_DC
128 INTEGER*4 HDC_TDC ! copied from HMS_RAW_DC
129 INTEGER*4 HDC_TOT_HITS
130 REAL*4 HDC_DRIFT_TIME
131 REAL*4 HDC_DRIFT_DIS
132 REAL*4 HDC_WIRE_CENTER
133 REAL*4 HDC_WIRE_COORD
134 COMMON/HMS_DECODED_DC/
135 1 HDC_DRIFT_TIME(HMAX_DC_HITS),
136 2 HDC_DRIFT_DIS(HMAX_DC_HITS),
137 3 HDC_WIRE_CENTER(HMAX_DC_HITS),
138 4 HDC_WIRE_COORD(HMAX_DC_HITS),
139 5 HDC_PLANE_NUM(HMAX_DC_HITS),
140 6 HDC_WIRE_NUM(HMAX_DC_HITS),
141 cdaq 1.1 7 HDC_TDC(HMAX_DC_HITS),
142 9 HDC_TOT_HITS
143 *
144 * HMS RAW SCINTILLATOR HITS
145 * FILLED BY g_decode_evebt_by_banks
146 * EACH SCINTILLATOR HIT IS SPECIFIED BY A
147 * SCINTILLATOR PLANE NUMBER
148 * SCINTILLATOR COUNTER NUMBER
149 * SCINTILLATOR ADC AT POSITIVE GENERALIZED COORDINATE END
150 * SCINTILLATOR ADC AT NEGATIVE GENERALIZED COORDINATE END
151 * SCINTILLATOR TDC AT POSITIVE GENERALIZED COORDINATE END
152 * SCINTILLATOR TDC AT NEGATIVE GENERALIZED COORDINATE END
153 * THE TOTAL NUMBER OF HITS IS GIVEN BY HSCIN_ALL_TOT_HITS
154 INTEGER*4 HMAX_ALL_SCIN_HITS ! MAXIMUM TOTAL NUMBER OF SCIN HITS
155 PARAMETER (HMAX_ALL_SCIN_HITS=53) ! Number of elements + 1
156 INTEGER*4 HNUM_ALL_SCIN_PLANES ! TOTAL NUMBER OF SCIN PLANES
157 PARAMETER (HNUM_ALL_SCIN_PLANES=4)
158 INTEGER*4 HSCIN_ALL_PLANE_NUM
159 INTEGER*4 HSCIN_ALL_COUNTER_NUM
160 INTEGER*4 HSCIN_ALL_ADC_POS
161 INTEGER*4 HSCIN_ALL_ADC_NEG
162 cdaq 1.1 INTEGER*4 HSCIN_ALL_TDC_POS
163 INTEGER*4 HSCIN_ALL_TDC_NEG
164 INTEGER*4 HSCIN_ALL_TOT_HITS
165 *
166 COMMON/HMS_RAW_SCIN/
167 & HSCIN_ALL_PLANE_NUM(HMAX_ALL_SCIN_HITS),
168 & HSCIN_ALL_COUNTER_NUM(HMAX_ALL_SCIN_HITS),
169 & HSCIN_ALL_ADC_POS(HMAX_ALL_SCIN_HITS),
170 & HSCIN_ALL_ADC_NEG(HMAX_ALL_SCIN_HITS),
171 & HSCIN_ALL_TDC_POS(HMAX_ALL_SCIN_HITS),
172 & HSCIN_ALL_TDC_NEG(HMAX_ALL_SCIN_HITS),
173 & HSCIN_ALL_TOT_HITS
174 *
175 * HMS REAL SCINTILLATOR HITS (Hits with no TDC data stripped out)
176 * FILLED BY H_strip_scin (which is called by h_trans_scin)
177 * EACH SCINTILLATOR HIT IS SPECIFIED BY A
178 * SCINTILLATOR PLANE NUMBER
179 * SCINTILLATOR COUNTER NUMBER
180 * SCINTILLATOR ADC AT POSITIVE GENERALIZED COORDINATE END
181 * SCINTILLATOR ADC AT NEGATIVE GENERALIZED COORDINATE END
182 * SCINTILLATOR TDC AT POSITIVE GENERALIZED COORDINATE END
183 cdaq 1.1 * SCINTILLATOR TDC AT NEGATIVE GENERALIZED COORDINATE END
184 * THE TOTAL NUMBER OF HITS IS GIVEN BY HSCIN_TOT_HITS
185 INTEGER*4 HMAX_SCIN_HITS ! MAXIMUM TOTAL NUMBER OF SCIN HITS
186 PARAMETER (HMAX_SCIN_HITS=53) ! Number of elements + 1
187 INTEGER*4 HNUM_SCIN_PLANES ! TOTAL NUMBER OF SCIN PLANES
188 PARAMETER (HNUM_SCIN_PLANES=4)
189 INTEGER*4 HNUM_SCIN_ELEMENTS
190 PARAMETER (HNUM_SCIN_ELEMENTS=16)
191 INTEGER*4 HSCIN_PLANE_NUM
192 INTEGER*4 HSCIN_COUNTER_NUM
193 REAL*4 HSCIN_ADC_POS
194 REAL*4 HSCIN_ADC_NEG
195 INTEGER*4 HSCIN_TDC_POS
196 INTEGER*4 HSCIN_TDC_NEG
197 INTEGER*4 HSCIN_TOT_HITS
|
198 saw 1.5 INTEGER*4 HSCIN_SING_COUNTER ! DJM''s registered hit counter
|
199 cdaq 1.1 *
200 COMMON/HMS_REAL_SCIN/
201 & HSCIN_PLANE_NUM(HMAX_SCIN_HITS),
202 & HSCIN_COUNTER_NUM(HMAX_SCIN_HITS),
203 & HSCIN_ADC_POS(HMAX_SCIN_HITS),
204 & HSCIN_ADC_NEG(HMAX_SCIN_HITS),
205 & HSCIN_TDC_POS(HMAX_SCIN_HITS),
206 & HSCIN_TDC_NEG(HMAX_SCIN_HITS),
207 & HSCIN_TOT_HITS,
208 & HSCIN_SING_COUNTER(HNUM_SCIN_PLANES)
209 *
210 *
211 * DECODED SCIN HITS
212 * FILLED BY H_TRANS_SCIN
213 * THIS USES SCINTILLATOR INFORMATION ONLY TO CALCUATE
214 * HSCIN_APPROX_HIT_COORD GENERALIZED COORDINATE ALONG HIT
215 * SCINTILLATOR DETERMINED FROM TDC INFO.
216 * HSCIN_COR_ADC CORRECTED PULSE HEIGHT AT HIT
217 * HSCIN_COR_TIME CORRECTED TIME AT HIT
218 * THE TOTAL NUMBER OF HITS IN EACH PLANE IS GIVEN BY HSCIN_HITS_PER_PLANE
219 * THE HSTART_TIME IS DEFINED AS THE TIME IF THE TRACK HAD GONE THROUGH
220 cdaq 1.1 * THE CENTER OF S1X. IT IS USED BY H_TRANS_DC FOR THE DRIFT TIME
221 * CALCULATION.
222 * HTWO_GOOD_TIMES IS A LOGICAL VARIABLE WHICH IS TRUE IF BOTH ENDS
223 * OF THE SCINTILLATOR ARE HIT AND FALSE IF ONLY ONE END IS HIT.
224 * HGOOD_START_TIME IS TRUE IF A START TIME WAS FOUND
225 *
226 REAL*4 HSCIN_ZPOS(HMAX_SCIN_HITS)
227 REAL*4 HSCIN_CENTER_COORD(HMAX_SCIN_HITS)
228 REAL*4 HSCIN_DEC_HIT_COORD(HMAX_SCIN_HITS)
229 REAL*4 HSCIN_WIDTH(HMAX_SCIN_HITS)
230 REAL*4 HSCIN_SLOP(HMAX_SCIN_HITS)
231 REAL*4 HSCIN_COR_ADC(HMAX_SCIN_HITS)
232 REAL*4 HSCIN_COR_TIME(HMAX_SCIN_HITS)
233 REAL*4 HSTART_TIME
234 INTEGER*4 HSTART_HITNUM
235 INTEGER*4 HSTART_HITSIDE
236 INTEGER*4 HSCIN_HITS_PER_PLANE(HNUM_SCIN_PLANES)
237 LOGICAL*4 HTWO_GOOD_TIMES(HMAX_SCIN_HITS)
238 LOGICAL*4 HGOOD_START_TIME
239
240 COMMON/HMS_DECODED_SCIN/
241 cdaq 1.1 & HSCIN_ZPOS,
242 & HSCIN_CENTER_COORD,
243 & HSCIN_WIDTH,
244 & HSCIN_SLOP,
245 & HSCIN_DEC_HIT_COORD,
246 & HSCIN_COR_ADC,
247 & HSCIN_COR_TIME,
248 & HSTART_TIME,
249 & HSCIN_HITS_PER_PLANE,
250 & HSTART_HITNUM,
251 & HSTART_HITSIDE,
252 & HTWO_GOOD_TIMES,
253 & HGOOD_START_TIME
254 *
255 *
256 * HMS CALORIMETER HITS
257 *
258 * ALLOW FOR NO SPARCIFICATION OF SHOWER COUNTER ELEMENTS
259 *
260 * EACH COUNTER HAS A
261 * COLUMN NUMBER
262 cdaq 1.1 * ROW NUMBER
|
263 saw 1.6 * ADC VALUE AT Positive End (This end if only one end used)
264 * ADC VALUE AT Negative End
|
265 cdaq 1.1 *
266 * THE TOTAL NUMBER OF SHOWER HITS IS HCAL_TOT_HITS
267 *
268 INTEGER*4 HMAX_CAL_BLOCKS ! TOTAL NUMBER OF SHOWER BLOCKS
269 PARAMETER (HMAX_CAL_BLOCKS=52)
270 INTEGER*4 HMAX_CAL_ROWS !Number of calorimeter rows
271 PARAMETER (HMAX_CAL_ROWS=13)
272 INTEGER*4 HMAX_CAL_COLUMNS !Number of calorimeter columns
273 PARAMETER (HMAX_CAL_COLUMNS=4)
274 INTEGER*4 HCAL_TOT_HITS
|
275 saw 1.7 INTEGER*4 HCAL_POS_HITS
276 INTEGER*4 HCAL_NEG_HITS
|
277 cdaq 1.1 INTEGER*4 HCAL_COLUMN(HMAX_CAL_BLOCKS)
278 INTEGER*4 HCAL_ROW(HMAX_CAL_BLOCKS)
|
279 saw 1.6 INTEGER*4 HCAL_ADC_POS(HMAX_CAL_BLOCKS)
280 INTEGER*4 HCAL_ADC_NEG(HMAX_CAL_BLOCKS)
|
281 cdaq 1.1 INTEGER*4 HCAL_ADC(HMAX_CAL_BLOCKS)
|
282 saw 1.6 equivalence (hcal_adc, hcal_adc_pos) ! For old code
283 REAL*4 HCAL_REALADC_POS(HMAX_CAL_BLOCKS)
284 REAL*4 HCAL_REALADC_NEG(HMAX_CAL_BLOCKS)
|
285 cdaq 1.1 REAL*4 HCAL_REALADC(HMAX_CAL_BLOCKS)
|
286 saw 1.6 equivalence (hcal_realadc, hcal_realadc_pos)
|
287 cdaq 1.1
288 COMMON/HMS_RAW_CAL/
289 & HCAL_COLUMN,
290 & HCAL_ROW,
|
291 saw 1.6 & HCAL_ADC_POS,
292 & HCAL_ADC_NEG,
293 & HCAL_REALADC_POS,
294 & HCAL_REALADC_NEG,
|
295 saw 1.7 & HCAL_TOT_HITS,
296 & HCAL_POS_HITS,
297 & HCAL_NEG_HITS
|
298 cdaq 1.1 *
299 * DECODED CALORIMETER QUANTITIES
300 * FILLED BY H_TRANS_CAL USING ONLY THE CALORIMETER INFORMATION
301 *
302 * X COORDINATE OF BLOCK CENTER
303 * Z COORDINATE OF BLOCK CENTER
304 * ENERGY DEPOSITED IN THE BLOCK
305 * ENERGY DEPOSITED IN COLUMN #1
306 * #2
307 * #3
308 * #4
309 * TOTAL ENERGY IN THE CALORIMETER
310 INTEGER*4 HNHITS_CAL !NUMBER OF CALORIMETER HITS ABOVE THRESHOLD
311 REAL*4 HBLOCK_XC !X COORDINATE OF BLOCK CENTER
312 REAL*4 HBLOCK_ZC !Z COORDINATE OF BLOCK CENTER
313 REAL*4 HBLOCK_DE !ENERGY DEPOSITION IN THE BLOCKS
314 REAL*4 HCAL_E1 !ENERGY DEPOSITION IN COLUMN #1
315 REAL*4 HCAL_E2 ! #2
316 REAL*4 HCAL_E3 ! #3
317 REAL*4 HCAL_E4 ! #4
318 REAL*4 HCAL_ET !TOTAL ENERGY IN THE CALORIMETER
|
319 saw 1.7 REAL*4 HCAL_E1_POS
320 REAL*4 HCAL_E1_NEG
321 REAL*4 HCAL_E2_POS
322 REAL*4 HCAL_E2_NEG
323 REAL*4 HBLOCK_DE_POS
324 REAL*4 HBLOCK_DE_NEG
|
325 cdaq 1.1 COMMON/HMS_DECODED_CAL/
326 & HBLOCK_XC(HMAX_CAL_BLOCKS),
327 & HBLOCK_ZC(HMAX_CAL_BLOCKS),
328 & HBLOCK_DE(HMAX_CAL_BLOCKS),
329 & HCAL_E1,
330 & HCAL_E2,
331 & HCAL_E3,
332 & HCAL_E4,
333 & HCAL_ET,
|
334 saw 1.7 & HCAL_E1_POS,
335 & HCAL_E1_NEG,
336 & HCAL_E2_POS,
337 & HCAL_E2_NEG,
338 & HBLOCK_DE_POS(HMAX_CAL_BLOCKS),
339 & HBLOCK_DE_NEG(HMAX_CAL_BLOCKS),
|
340 cdaq 1.1 & HNHITS_CAL
341 *
342 * HMS CERENKOV HITS
343 * FILLED BY G_decode_event_by_banks
344 *
345 * THERE ARE TWO CERENKOV PHOTOTUBES. EACH HIT HAS
346 * TUBE NUMBER
347 * ADC VALUE
348 * THE TOTAL NUMBER OF PHOTOTUBE HITS IS HCER_TOT_HITS
349 * HOWEVER TO MAKE THE DECODERS SYMMETRICAL TO A SCINT WE MUST ADD
350 * A DUMMY PLANE NUMBER.
351 INTEGER*4 HMAX_CER_HITS
352 PARAMETER(HMAX_CER_HITS=2)
353 INTEGER*4 HCER_TOT_HITS
354 INTEGER*4 HCER_TUBE_NUM
|
355 cdaq 1.3 INTEGER*4 HCER_RAW_ADC
|
356 cdaq 1.1 INTEGER*4 HCER_PLANE
357 COMMON/HMS_RAW_CER/
358 & HCER_TUBE_NUM(HMAX_CER_HITS),
|
359 cdaq 1.3 & HCER_RAW_ADC(HMAX_CER_HITS),
|
360 cdaq 1.1 & HCER_PLANE(HMAX_CER_HITS),
361 & HCER_TOT_HITS
362 *
|
363 cdaq 1.2 * DECODED CERENKOV QUANTITIES
364 * FILLED BY H_TRANS_CER
365 *
366 INTEGER*4 HCER_NUM_HITS ! NUMBER OF CERENKOV HITS ABOVE THRESHOLD
367 REAL*4 HCER_NPE ! ADC CONVERTED TO NUMBER OF PHOTOELECTRONS
|
368 saw 1.5 REAL*4 HCER_NPE_SUM ! SUM OVER TUBES OF NPE''S
|
369 cdaq 1.3 REAL*4 HCER_ADC ! PED SUBTRACTED ADC FOR EACH *TUBE*
|
370 cdaq 1.2 COMMON/HCER_DECODED_CER/
371 & HCER_NUM_HITS,
372 & HCER_NPE(HMAX_CER_HITS),
|
373 cdaq 1.3 & HCER_NPE_SUM,
374 & HCER_ADC(HMAX_CER_HITS)
|
375 cdaq 1.2 *
|
376 jones 1.9 *..........................................................................
377 * (Last revison 01 Dec'02; Hamlet)
378 *
379 * HMS AEROGEL HITS
380 * FILLED BY G_decode_event_by_banks
381 *
382 * THERE ARE 16 AEROGEL PHOTOTUBES. We will pair tubes on the left and
383 * right side so that are 8 "counters".
384 * POS pmt's on the positive side of Y axis (as for Hodoscope)
385 *
386 * TUBE Row
387 * ADC Left VALUE
388 * ADC Right VALUE
389 * THE TOTAL NUMBER OF PHOTOTUBE HITS IS HAERO_TOT_HITS
390 * HOWEVER TO MAKE THE DECODERS SYMMETRICAL TO A SCINT WE MUST ADD
391 * A DUMMY PLANE NUMBER.
392 INTEGER*4 HMAX_AERO_HITS
393 PARAMETER(HMAX_AERO_HITS=8)
394 INTEGER*4 HNUM_AERO_BLOCKS
395 PARAMETER(HNUM_AERO_BLOCKS=8)
396 INTEGER*4 HAERO_TOT_HITS
397 jones 1.9 INTEGER*4 HAERO_PLANE
398 INTEGER*4 HAERO_PAIR_NUM
399 INTEGER*4 HAERO_ADC_POS
400 INTEGER*4 HAERO_ADC_NEG
401 INTEGER*4 HAERO_TDC_POS
402 INTEGER*4 HAERO_TDC_NEG
403 COMMON/HMS_RAW_AERO/
404 & HAERO_PLANE(HMAX_AERO_HITS),
405 & HAERO_PAIR_NUM(HMAX_AERO_HITS),
406 & HAERO_ADC_POS(HMAX_AERO_HITS),
407 & HAERO_ADC_NEG(HMAX_AERO_HITS),
408 & HAERO_TDC_POS(HMAX_AERO_HITS),
409 & HAERO_TDC_NEG(HMAX_AERO_HITS),
410 & HAERO_TOT_HITS
411
412 *
413 *..........................................................................
|
414 cdaq 1.1 *
415 *
416 * HMS DETECTOR TRACK QUANTITIES
417 * FILLED BY H_TRACK SUBROUTINE
418 *
|
419 jones 1.11 INTEGER*4 HNTRACKS_MAX ! size of array
420 PARAMETER (HNTRACKS_MAX=20) ! SET MAXIMUM TO 20
|
421 cdaq 1.1 INTEGER*4 HNTRACKHITS_MAX ! MAXIMUM NUMBER OF HITS IN EACH TRACK
422 PARAMETER (HNTRACKHITS_MAX=16) ! SET MAXIMUM TO 15
423 INTEGER*4 HNTRACKS_FP ! NUMBER OF FOCAL PLANE TRACKS FOUND
424 * ALL THE FOLLOWING VARIABLES ARE ARRAYS
425 *
426 REAL*4 HX_FP ! X POSITION OF TRACK IN FOCAL PLANE
427 REAL*4 HY_FP ! Y POSITION OF TRACK IN FOCAL PLANE
428 REAL*4 HZ_FP ! Z POSITION OF TRACK IN FOCAL PLANE
429 REAL*4 HXP_FP ! X SLOPE OF TRACK IN FOCAL PLANE
430 REAL*4 HYP_FP ! Y SLOPE OF TRACK IN FOCAL PLANE
431 REAL*4 HCHI2_FP ! FIT QUALITY IN FOCAL PLANE
432 REAL*4 HCHI2PERDOF_FP ! fp chi2 per degree of freedom
433 INTEGER*4 HNFREE_FP ! NUMBER OF DEGREES OF FREEDOM IN FIT
434 INTEGER*4 HNTRACK_HITS ! LIST OF HITS ON EACH TRACK
435 REAL*4 HDEL_FP ! FOCAL PLANE ERROR MATRIX
436 REAL*4 HX_FP_rot ! x in rotated focal plane
437 REAL*4 HY_FP_rot ! y in rotated f plane
438 REAL*4 HXP_FP_rot ! x slope in rotated f plane
439 REAL*4 HYP_FP_rot ! y slope in rot f plane
440
441 * THE FOCAL PLANE ERROR MATRIX IS A 4 BY 4 BY HNTRACK_MAX ARRAY
442 cdaq 1.1 * THE FOUR INDECIES FOR EACH TRACK ARE X, XP, Y, YP
443 * THE ERROR MATRIX FOR EACH TRACK IS SYMMETRIC ABOUT THE DIAGONAL
444 * DELXX DELXXP DELXY DELXYP
445 * DELXPX DELXPXP DELXPY DELXPYP
446 * DELYPX DELYXP DELYY DELYYP
447 * DELYPX DELYPXP DELYPY DELYPYP
448 COMMON/HMS_FOCAL_PLANE/
449 1 HX_FP(HNTRACKS_MAX),
450 2 HY_FP(HNTRACKS_MAX),
451 3 HZ_FP(HNTRACKS_MAX),
452 4 HXP_FP(HNTRACKS_MAX),
453 5 HYP_FP(HNTRACKS_MAX),
454 6 HCHI2_FP(HNTRACKS_MAX),
455 7 HDEL_FP(4,4,HNTRACKS_MAX),
456 8 HNTRACK_HITS(HNTRACKS_MAX,HNTRACKHITS_MAX+1),
457 9 HNFREE_FP(HNTRACKS_MAX),
458 A HNTRACKS_FP,HCHI2PERDOF_FP(HNTRACKS_MAX),
459 > HX_FP_rot(HNTRACKS_MAX),
460 > HY_FP_rot(HNTRACKS_MAX),
461 > HXP_FP_rot(HNTRACKS_MAX),
462 > HYP_FP_rot(HNTRACKS_MAX)
|
463 jones 1.11
|
464 cdaq 1.1
465
466 * HMS TARGET QUANTITIES
467 * ASSUME THE SAME NUMBER OF MAXIMUM TRACKS AS IN THE FOCAL PLANE
468 * HNTRACKS_MAX
469 *
470 * FILLED BY H_TARG_TRANS SUBROUTINE
471 INTEGER*4 HNTRACKS_TAR ! NUMBER OF TARGET TRACKS FOUND
472 * ALL THE FOLLOWING VARIABLES ARE ARRAYS
473 *
474 REAL*4 HX_TAR ! X POSITION OF TRACK AT TARGET
475 REAL*4 HY_TAR ! Y POSITION OF TRACK AT TARGET
476 REAL*4 HZ_TAR ! Z POSITION OF TRACK AT TARGET
477 REAL*4 HXP_TAR ! X SLOPE OF TRACK AT TARGET
478 REAL*4 HYP_TAR ! Y SLOPE OF TRACK AT TARGET
479 REAL*4 HDELTA_TAR ! FRACTION TRACK MOMENTUM
480 REAL*4 HP_TAR ! MOMENTUM OF TRACK AT TARGET
481 REAL*4 HCHI2_TAR ! FIT QUALITY AT TARGET
482 INTEGER*4 HNFREE_TAR ! NUMBER OF DEGREES OF FREEDOM IN FIT
483 REAL*4 HDEL_TAR ! TARGET ERROR MATRIX
484 * THE FOCAL PLANE ERROR MATRIX IS A 5 BY 5 BY HNTRACK_MAX ARRAY
485 cdaq 1.1 * THE FIVE INDECIES FOR EACH TRACK ARE X, Y , XP, YP and P
486 * THE ERROR MATRIX FOR EACH TRACK IS SYMMETRIC ABOUT THE DIAGONAL
487 * DELXX DELXXP DELXY DELXYP DELXP
488 * DELXPX DELXPXP DELXPY DELXPYP DELXPP
489 * DELYX DELYXPY DELYY DELYYP DELYP
490 * DELYPX DELYPXP DELYPY DELYPYP DELYPP
491 * DELPX DELPXP DELPY DELPXP DELPP
492 INTEGER*4 HLINK_TAR_FP ! LINK POINTER TO FOCAL PLANE TRACK NUMBER
493 COMMON/HMS_TARGET/
494 1 HX_TAR(HNTRACKS_MAX),
495 2 HY_TAR(HNTRACKS_MAX),
496 3 HZ_TAR(HNTRACKS_MAX),
497 4 HXP_TAR(HNTRACKS_MAX),
498 5 HYP_TAR(HNTRACKS_MAX),
499 6 HDELTA_TAR(HNTRACKS_MAX),
500 7 HP_TAR(HNTRACKS_MAX),
501 8 HCHI2_TAR(HNTRACKS_MAX),
502 9 HDEL_TAR(5,5,HNTRACKS_MAX),
503 A HNFREE_TAR(HNTRACKS_MAX),
504 B HLINK_TAR_FP(HNTRACKS_MAX),
505 C HNTRACKS_TAR
506 cdaq 1.1 *
507 *
508 * HMS_TRACK_TESTS
509 *
510 * PARTICLE ID INFORMATION FILLED BY H_TOF and H_CAL
511 * THIS STORES THE RESULTS OF SHOWER AND SCINTILLATOR CALCULATIONS
512 * FOR EACH OF THE TRACKS GENERATED BY H_TRACK
513 *
514 INTEGER*4 HNBLOCKS_CAL(HNTRACKS_MAX)
515 REAL*4 HTRACK_E1(HNTRACKS_MAX)
516 REAL*4 HTRACK_E2(HNTRACKS_MAX)
517 REAL*4 HTRACK_E3(HNTRACKS_MAX)
518 REAL*4 HTRACK_E4(HNTRACKS_MAX)
519 REAL*4 HTRACK_ET(HNTRACKS_MAX)
520 REAL*4 HTRACK_PRESHOWER_E(HNTRACKS_MAX)
|
521 saw 1.7 REAL*4 HTRACK_E1_POS(HNTRACKS_MAX)
522 REAL*4 HTRACK_E1_NEG(HNTRACKS_MAX)
523 REAL*4 HTRACK_E2_POS(HNTRACKS_MAX)
524 REAL*4 HTRACK_E2_NEG(HNTRACKS_MAX)
|
525 cdaq 1.1 *
526 INTEGER*4 HSCIN_HIT(HNTRACKS_MAX,HMAX_SCIN_HITS)
527 INTEGER*4 HNUM_SCIN_HIT(HNTRACKS_MAX)
528 INTEGER*4 HNUM_PMT_HIT(HNTRACKS_MAX)
529 REAL*4 HDEDX(HNTRACKS_MAX,HMAX_SCIN_HITS)
530 REAL*4 HBETA(HNTRACKS_MAX)
531 REAL*4 HBETA_CHISQ(HNTRACKS_MAX)
532 REAL*4 HTIME_AT_FP(HNTRACKS_MAX)
|
533 cdaq 1.3 REAL*4 HSCIN_FPTIME(HNTRACKS_MAX,HMAX_SCIN_HITS)
|
534 cdaq 1.1 *
535 COMMON/HMS_TRACK_TESTS/
|
536 saw 1.7 1 HTRACK_E1, ! ENERGY IN FIRST COLUMN NEAR TRACK
|
537 cdaq 1.1 2 HTRACK_E2, ! SECOND COLUMN
538 3 HTRACK_E3, ! THIRD COLUMN
539 4 HTRACK_E4, ! FOURTH COLUMN
540 5 HTRACK_ET, ! TOTAL SHOWER ENERGY ALONG TRACK
541 6 HTRACK_PRESHOWER_E, ! PRESHOWER ENERGY-Note definition is exp. dep.
542 7 HDEDX, !
543 8 HBETA, ! VELOCITY OF TRACK
544 9 HBETA_CHISQ, ! CHISQ OF FIT TO BETA OF TRACK
545 A HTIME_AT_FP, !
546 B HNBLOCKS_CAL, ! NUMBER OF SHOWER BLOCKS ON EACH TRACK
547 C HSCIN_HIT, ! ARRAY OF SCIN HITS ASSOCIATED WITH
548 D ! EACH TRACK
549 E HNUM_SCIN_HIT, ! NUMBER OF HITS FOR EACH TRACK
|
550 cdaq 1.3 F HNUM_PMT_HIT, ! NUMBER OF PMT HITS FOR EACH TRACK
|
551 saw 1.7 G HSCIN_FPTIME,
552 H HTRACK_E1_POS,
553 I HTRACK_E1_NEG,
554 J HTRACK_E2_POS,
555 K HTRACK_E2_NEG
|
556 cdaq 1.1 *
557 * HMS SINGLES PHYSICS COMMON BLOCKS
558 * THESE ARE FILLED BY H_PHYSICS
559 *
560 *
561 REAL*4 HSP ! Lab momentum of chosen track in GeV/c
562 REAL*4 HSENERGY ! Lab total energy of chosen track in GeV
|
563 csa 1.8 REAL*4 HSCORRP ! electron momentum corrected for eloss
564 REAL*4 HSCORRE ! electron energy corrected for eloss
|
565 cdaq 1.1 REAL*4 HSDELTA ! Spectrometer delta of chosen track
566 REAL*4 HSTHETA ! Lab Scattering angle in radians
567 REAL*4 HSPHI ! Lab Azymuthal angle in radians
|
568 jones 1.12 REAL*4 HSINPLANE ! In plane angle
|
569 cdaq 1.1 REAL*4 HSMINV ! Invariant Mass of remaing hadronic system
570 REAL*4 HSZBEAM ! Lab Z coordinate of intersection of beam
571 ! track with spectrometer ray
572 REAL*4 HSDEDX(4) ! DEDX of chosen track in each scin plane
573 REAL*4 HSBETA ! BETA of chosen track
574 REAL*4 HSTRACK_ET ! Total shower energy of chosen track
575 REAL*4 HSTRACK_PRESHOWER_E ! preshower of chosen track
|
576 saw 1.7 REAL*4 HSTRACK_E1_POS ! Shower energy E1_POS of chosen track
577 REAL*4 HSTRACK_E1_NEG ! Shower energy E1_NEG of chosen track
578 REAL*4 HSTRACK_E2_POS ! Shower energy E2_POS of chosen track
579 REAL*4 HSTRACK_E2_NEG ! Shower energy E2_NEG of chosen track
|
580 cdaq 1.1 REAL*4 HSTIME_AT_FP
581 REAL*4 HSX_FP ! X focal plane position
582 REAL*4 HSY_FP
583 REAL*4 HSXP_FP
584 REAL*4 HSYP_FP
585 REAL*4 HSCHI2PERDEG ! CHI2 per degree of freedom of chosen track.
586 REAL*4 HSX_TAR
587 REAL*4 HSY_TAR
588 REAL*4 HSXP_TAR
589 REAL*4 HSYP_TAR
590 REAL*4 HSBETA_CHISQ
591 *
592 real*4 hsmass2 ! Mass squared
593 real*4 hst ! invariant t
594 real*4 hsu ! invariant u
595 real*4 hseloss
596 real*4 hsq3 ! Lab frame momentum transfer
|
597 csa 1.8 * real*4 hsthetaq, hsphiq ! Direction of q3
|
598 cdaq 1.1 real*4 hsbigq2 ! Q**2
599 real*4 hsx ! fraction of nucleon p carried by quark
|
600 saw 1.5 real*4 hsy ! fraction of lepton''s E lost in lab
|
601 cdaq 1.1 real*4 hsw ! Invariant mass of recoil system
602 real*4 hsw2 ! Invariant mass**2 of recoil system
|
603 cdaq 1.12.12.1 real*4 hsnu ! Ebeam-HSenergy
604 real*4 hsx_bj ! Bjorken x
605 real*4 y_scale ! scaling variable y
606 real*4 yscale ! real function that calculates Y
607 real*4 xsi_scale ! nachtmann variable (=x at Q2=1)
608 real*4 yscale_eps ! separation energy
|
609 cdaq 1.1 *
610 INTEGER*4 HSNUM_FPTRACK ! Index of focal plane track chosen
611 INTEGER*4 HSNUM_TARTRACK ! Index of target track chosen
612 INTEGER*4 HSID_LUND ! LUND particle ID code -- not yet filled
613 INTEGER*4 HSNFREE_FP
614 INTEGER*4 HSNUM_SCIN_HIT ! # OF SCINTILLATORS ON TRACK
615 INTEGER*4 HSNUM_PMT_HIT ! # OF HODOSCOPE PMTS ON TRACK
616 *
617 COMMON/HMS_PHYSICS_R4/
618 & HSP,
619 & HSENERGY,
620 & HSDELTA,
621 & HSTHETA,
622 & HSPHI,
|
623 jones 1.12 & HSINPLANE,
|
624 cdaq 1.1 & HSMINV,
625 & HSZBEAM,
626 & HSDEDX,
627 & HSBETA,
628 & HSTRACK_ET,
629 & HSTRACK_PRESHOWER_E,
|
630 saw 1.7 & HSTRACK_E1_POS,
631 & HSTRACK_E1_NEG,
632 & HSTRACK_E2_POS,
633 & HSTRACK_E2_NEG,
|
634 cdaq 1.1 & HSTIME_AT_FP,
635 & HSX_FP ,
636 & HSY_FP ,
637 & HSXP_FP ,
638 & HSYP_FP ,
639 & HSCHI2PERDEG ,
640 & HSX_TAR ,
641 & HSY_TAR ,
642 & HSXP_TAR ,
643 & HSYP_TAR ,
644 & HSBETA_CHISQ,
|
645 csa 1.8 & hsmass2,
646 & hst,
647 & hsu,
648 & hseloss,
649 & hsq3,
650 & hsbigq2,
651 & hsx,
652 & hsy,
653 & hsw,
654 & hsw2,
655 & hscorrp,
|
656 cdaq 1.12.12.1 & hscorre,
657 & hsnu,
658 & hsx_bj,
659 & y_scale,
660 & xsi_scale
|
661 csa 1.8 * & hsthetaq,hsphiq
|
662 cdaq 1.1 *
663 COMMON/HMS_PHYSICS_I4/
664 & HSNUM_FPTRACK,
665 & HSNUM_TARTRACK,
666 & HSID_LUND,
667 & HSNFREE_FP,
668 & HSNUM_SCIN_HIT,
669 & HSNUM_PMT_HIT
670 *
671 *
672 * Non-Hits data with HMS gates or starts.
673 * (Energy SUMS, logic timing,...)
674 *
675 INTEGER HMAX_MISC_HITS
|
676 cdaq 1.3 PARAMETER(HMAX_MISC_HITS=100)
|
677 cdaq 1.1 INTEGER*4 HMISC_TOT_HITS
|
678 saw 1.4 INTEGER*4 HMISC_RAW_ADDR1 ! "Plane" (1=TDC,2=ADC)
|
679 cdaq 1.1 INTEGER*4 HMISC_RAW_ADDR2 ! "Counter"
680 INTEGER*4 HMISC_RAW_DATA
681 COMMON/H_RAW_MISC/
682 & HMISC_RAW_ADDR1(HMAX_MISC_HITS),
683 & HMISC_RAW_ADDR2(HMAX_MISC_HITS),
684 & HMISC_RAW_DATA(HMAX_MISC_HITS),
685 & HMISC_TOT_HITS
686 *
687 *
688 *
689 *******************end: hms_data_structures.cmn*************************
690 * Local Variables:
691 * mode: fortran
692 * fortran-continuation-string: "&"
693 * comment-column: 35
694 * End:
|
695 jones 1.9
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716 jones 1.9
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737 jones 1.9
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758 jones 1.9
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779 jones 1.9
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800 jones 1.9
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821 jones 1.9
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842 jones 1.9
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863 jones 1.9
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884 jones 1.9
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905 jones 1.9
906
907
908
909
910
911
|
912 cdaq 1.1
|