1 cdaq 1.1 *****************begin: hms_data_structures.cmn*************************
2 *
3 * include file hms_data_structures.cmn
4 *
5 * Author: D. F. Geesaman 1 September 1993
6 *
|
7 cdaq 1.2 * $Log: hms_data_structures.cmn,v $
|
8 saw 1.5 * Revision 1.4 1996/01/24 16:18:59 saw
9 * (JRA) Cosmetic change
10 *
|
11 saw 1.4 * Revision 1.3 1996/01/17 16:00:02 cdaq
12 * (JRA) ADD HCER_RAW_ADC and HSCIN_FPTIME
13 *
|
14 cdaq 1.3 * Revision 1.2 1995/08/11 16:28:46 cdaq
15 * (CC) Add structure for # of photoelectrons in Cerenkov
16 *
|
17 cdaq 1.2 * Revision 1.1 1995/05/22 18:41:15 cdaq
18 * Initial revision
|
19 cdaq 1.1 *
|
20 cdaq 1.2 *
|
21 cdaq 1.1 * CTPTYPE=parm
22 *
23 * HTHETA_LAB HMS LAB ANGLE THETA (RADIANS)
24 * HPHI_LAB HMS LAB ANGLE PHI
25 * HPCENTRAL HMS CENTRAL MOMENUTM (GEV)
26 * HBFIELD HMS B FIELD INCLUDING SIGN
27 *
28 REAL*4 HTHETA_LAB ! HMS LAB ANGLE THETA (RADIANS)
29 REAL*4 HPHI_LAB ! HMS LAB ANGLE PHI
30 REAL*4 HPCENTRAL ! HMS CENTRAL MOMENUTM (GEV)
31 REAL*4 HBFIELD ! HMS B FIELD INCLUDING SIGN
32 REAL*4 HPARTMASS ! EXPECTED MASS OF DETECTED PARTICLE IN HMS
33 COMMON/HMS_SPECTROMETER/
34 & HTHETA_LAB,
35 & HPHI_LAB,
36 & HPCENTRAL,
37 & HBFIELD,
38 & HPARTMASS
39 *
40 * CTPTYPE=event
41 *
42 cdaq 1.1 * HMS DECODED DATA
43 *
44 * FILLED BY G_decode_event_by_banks
45 * HMS DRIFT CHAMBER HITS
46 * EACH CHAMBER HIT (TDC VALUE) HAS A
47 * PLANE NUMBER
48 * WIRE NUMBER
49 * TDC VALUE
50 *
51 * THE TOTAL NUMBER OF HITS IS GIVEN IN HDC_RAW_TOT_HITS
52 *
53 INTEGER*4 HMAX_DC_HITS ! MAXIUM NUMBER OF DRIFT CHAMBER HITS
54 INTEGER*4 HMAX_NUM_DC_PLANES ! MAX NUMBER OF HMS DRIFT CHAMBER PLANES
55 INTEGER*4 HMAX_NUM_CHAMBERS ! NUMBER OF HMS DRIFT CHAMBERS
56 PARAMETER(HMAX_DC_HITS=3600)
57 PARAMETER(HMAX_NUM_DC_PLANES=12)
58 PARAMETER(HMAX_NUM_CHAMBERS=2)
59 INTEGER*4 HDC_RAW_PLANE_NUM
60 INTEGER*4 HDC_RAW_WIRE_NUM
61 INTEGER*4 HDC_RAW_TDC
62 INTEGER*4 HDC_RAW_TOT_HITS
63 cdaq 1.1 COMMON/HMS_RAW_DC/
64 1 HDC_RAW_PLANE_NUM(HMAX_DC_HITS),
65 2 HDC_RAW_WIRE_NUM(HMAX_DC_HITS),
66 3 HDC_RAW_TDC(HMAX_DC_HITS),
67 4 HDC_RAW_TOT_HITS
68 *
69 * HMS DECODED DC DATA
70 *
71 * FILLED BY H_TRANS_DC TRANSLATION ROUTINE
72 * HMS DRIFT CHAMBER HITS
73 * EACH CHAMBER HIT (TDC VALUE) HAS A
74 * PLANE NUMBER
75 * WIRE NUMBER
76 * TDC VALUE
77 * DRIFT TIME
78 * DRIFT DISTANCE
79 * GENERALIZED COORDINATE OF HIT WIRE CENTER
80 * GENERALIZED COORDINATE OF HIT PERPENDICULAR TO WIRE DIRECTION
81 *
82 * THE TOTAL NUMBER OF HITS IN EACH PLANE IS GIVEN IN HDC_HITS_PER_PLANE(I)
83 * THE TOTAL NUMBER OF HITS IS GIVEN IN HDC_TOT_HITS
84 cdaq 1.1 *
85 INTEGER*4 HDC_PLANE_NUM ! copied from HMS_RAW_DC
86 INTEGER*4 HDC_WIRE_NUM ! copied from HMS_RAW_DC
87 INTEGER*4 HDC_TDC ! copied from HMS_RAW_DC
88 INTEGER*4 HDC_TOT_HITS
89 REAL*4 HDC_DRIFT_TIME
90 REAL*4 HDC_DRIFT_DIS
91 REAL*4 HDC_WIRE_CENTER
92 REAL*4 HDC_WIRE_COORD
93 COMMON/HMS_DECODED_DC/
94 1 HDC_DRIFT_TIME(HMAX_DC_HITS),
95 2 HDC_DRIFT_DIS(HMAX_DC_HITS),
96 3 HDC_WIRE_CENTER(HMAX_DC_HITS),
97 4 HDC_WIRE_COORD(HMAX_DC_HITS),
98 5 HDC_PLANE_NUM(HMAX_DC_HITS),
99 6 HDC_WIRE_NUM(HMAX_DC_HITS),
100 7 HDC_TDC(HMAX_DC_HITS),
101 9 HDC_TOT_HITS
102 *
103 * HMS RAW SCINTILLATOR HITS
104 * FILLED BY g_decode_evebt_by_banks
105 cdaq 1.1 * EACH SCINTILLATOR HIT IS SPECIFIED BY A
106 * SCINTILLATOR PLANE NUMBER
107 * SCINTILLATOR COUNTER NUMBER
108 * SCINTILLATOR ADC AT POSITIVE GENERALIZED COORDINATE END
109 * SCINTILLATOR ADC AT NEGATIVE GENERALIZED COORDINATE END
110 * SCINTILLATOR TDC AT POSITIVE GENERALIZED COORDINATE END
111 * SCINTILLATOR TDC AT NEGATIVE GENERALIZED COORDINATE END
112 * THE TOTAL NUMBER OF HITS IS GIVEN BY HSCIN_ALL_TOT_HITS
113 INTEGER*4 HMAX_ALL_SCIN_HITS ! MAXIMUM TOTAL NUMBER OF SCIN HITS
114 PARAMETER (HMAX_ALL_SCIN_HITS=53) ! Number of elements + 1
115 INTEGER*4 HNUM_ALL_SCIN_PLANES ! TOTAL NUMBER OF SCIN PLANES
116 PARAMETER (HNUM_ALL_SCIN_PLANES=4)
117 INTEGER*4 HSCIN_ALL_PLANE_NUM
118 INTEGER*4 HSCIN_ALL_COUNTER_NUM
119 INTEGER*4 HSCIN_ALL_ADC_POS
120 INTEGER*4 HSCIN_ALL_ADC_NEG
121 INTEGER*4 HSCIN_ALL_TDC_POS
122 INTEGER*4 HSCIN_ALL_TDC_NEG
123 INTEGER*4 HSCIN_ALL_TOT_HITS
124 *
125 COMMON/HMS_RAW_SCIN/
126 cdaq 1.1 & HSCIN_ALL_PLANE_NUM(HMAX_ALL_SCIN_HITS),
127 & HSCIN_ALL_COUNTER_NUM(HMAX_ALL_SCIN_HITS),
128 & HSCIN_ALL_ADC_POS(HMAX_ALL_SCIN_HITS),
129 & HSCIN_ALL_ADC_NEG(HMAX_ALL_SCIN_HITS),
130 & HSCIN_ALL_TDC_POS(HMAX_ALL_SCIN_HITS),
131 & HSCIN_ALL_TDC_NEG(HMAX_ALL_SCIN_HITS),
132 & HSCIN_ALL_TOT_HITS
133 *
134 * HMS REAL SCINTILLATOR HITS (Hits with no TDC data stripped out)
135 * FILLED BY H_strip_scin (which is called by h_trans_scin)
136 * EACH SCINTILLATOR HIT IS SPECIFIED BY A
137 * SCINTILLATOR PLANE NUMBER
138 * SCINTILLATOR COUNTER NUMBER
139 * SCINTILLATOR ADC AT POSITIVE GENERALIZED COORDINATE END
140 * SCINTILLATOR ADC AT NEGATIVE GENERALIZED COORDINATE END
141 * SCINTILLATOR TDC AT POSITIVE GENERALIZED COORDINATE END
142 * SCINTILLATOR TDC AT NEGATIVE GENERALIZED COORDINATE END
143 * THE TOTAL NUMBER OF HITS IS GIVEN BY HSCIN_TOT_HITS
144 INTEGER*4 HMAX_SCIN_HITS ! MAXIMUM TOTAL NUMBER OF SCIN HITS
145 PARAMETER (HMAX_SCIN_HITS=53) ! Number of elements + 1
146 INTEGER*4 HNUM_SCIN_PLANES ! TOTAL NUMBER OF SCIN PLANES
147 cdaq 1.1 PARAMETER (HNUM_SCIN_PLANES=4)
148 INTEGER*4 HNUM_SCIN_ELEMENTS
149 PARAMETER (HNUM_SCIN_ELEMENTS=16)
150 INTEGER*4 HSCIN_PLANE_NUM
151 INTEGER*4 HSCIN_COUNTER_NUM
152 REAL*4 HSCIN_ADC_POS
153 REAL*4 HSCIN_ADC_NEG
154 INTEGER*4 HSCIN_TDC_POS
155 INTEGER*4 HSCIN_TDC_NEG
156 INTEGER*4 HSCIN_TOT_HITS
|
157 saw 1.5 INTEGER*4 HSCIN_SING_COUNTER ! DJM''s registered hit counter
|
158 cdaq 1.1 *
159 COMMON/HMS_REAL_SCIN/
160 & HSCIN_PLANE_NUM(HMAX_SCIN_HITS),
161 & HSCIN_COUNTER_NUM(HMAX_SCIN_HITS),
162 & HSCIN_ADC_POS(HMAX_SCIN_HITS),
163 & HSCIN_ADC_NEG(HMAX_SCIN_HITS),
164 & HSCIN_TDC_POS(HMAX_SCIN_HITS),
165 & HSCIN_TDC_NEG(HMAX_SCIN_HITS),
166 & HSCIN_TOT_HITS,
167 & HSCIN_SING_COUNTER(HNUM_SCIN_PLANES)
168 *
169 *
170 * DECODED SCIN HITS
171 * FILLED BY H_TRANS_SCIN
172 * THIS USES SCINTILLATOR INFORMATION ONLY TO CALCUATE
173 * HSCIN_APPROX_HIT_COORD GENERALIZED COORDINATE ALONG HIT
174 * SCINTILLATOR DETERMINED FROM TDC INFO.
175 * HSCIN_COR_ADC CORRECTED PULSE HEIGHT AT HIT
176 * HSCIN_COR_TIME CORRECTED TIME AT HIT
177 * THE TOTAL NUMBER OF HITS IN EACH PLANE IS GIVEN BY HSCIN_HITS_PER_PLANE
178 * THE HSTART_TIME IS DEFINED AS THE TIME IF THE TRACK HAD GONE THROUGH
179 cdaq 1.1 * THE CENTER OF S1X. IT IS USED BY H_TRANS_DC FOR THE DRIFT TIME
180 * CALCULATION.
181 * HTWO_GOOD_TIMES IS A LOGICAL VARIABLE WHICH IS TRUE IF BOTH ENDS
182 * OF THE SCINTILLATOR ARE HIT AND FALSE IF ONLY ONE END IS HIT.
183 * HGOOD_START_TIME IS TRUE IF A START TIME WAS FOUND
184 *
185 REAL*4 HSCIN_ZPOS(HMAX_SCIN_HITS)
186 REAL*4 HSCIN_CENTER_COORD(HMAX_SCIN_HITS)
187 REAL*4 HSCIN_DEC_HIT_COORD(HMAX_SCIN_HITS)
188 REAL*4 HSCIN_WIDTH(HMAX_SCIN_HITS)
189 REAL*4 HSCIN_SLOP(HMAX_SCIN_HITS)
190 REAL*4 HSCIN_COR_ADC(HMAX_SCIN_HITS)
191 REAL*4 HSCIN_COR_TIME(HMAX_SCIN_HITS)
192 REAL*4 HSTART_TIME
193 INTEGER*4 HSTART_HITNUM
194 INTEGER*4 HSTART_HITSIDE
195 INTEGER*4 HSCIN_HITS_PER_PLANE(HNUM_SCIN_PLANES)
196 LOGICAL*4 HTWO_GOOD_TIMES(HMAX_SCIN_HITS)
197 LOGICAL*4 HGOOD_START_TIME
198
199 COMMON/HMS_DECODED_SCIN/
200 cdaq 1.1 & HSCIN_ZPOS,
201 & HSCIN_CENTER_COORD,
202 & HSCIN_WIDTH,
203 & HSCIN_SLOP,
204 & HSCIN_DEC_HIT_COORD,
205 & HSCIN_COR_ADC,
206 & HSCIN_COR_TIME,
207 & HSTART_TIME,
208 & HSCIN_HITS_PER_PLANE,
209 & HSTART_HITNUM,
210 & HSTART_HITSIDE,
211 & HTWO_GOOD_TIMES,
212 & HGOOD_START_TIME
213 *
214 *
215 * HMS CALORIMETER HITS
216 *
217 * ALLOW FOR NO SPARCIFICATION OF SHOWER COUNTER ELEMENTS
218 *
219 * EACH COUNTER HAS A
220 * COLUMN NUMBER
221 cdaq 1.1 * ROW NUMBER
222 * ADC VALUE
223 *
224 * THE TOTAL NUMBER OF SHOWER HITS IS HCAL_TOT_HITS
225 *
226 INTEGER*4 HMAX_CAL_BLOCKS ! TOTAL NUMBER OF SHOWER BLOCKS
227 PARAMETER (HMAX_CAL_BLOCKS=52)
228 INTEGER*4 HMAX_CAL_ROWS !Number of calorimeter rows
229 PARAMETER (HMAX_CAL_ROWS=13)
230 INTEGER*4 HMAX_CAL_COLUMNS !Number of calorimeter columns
231 PARAMETER (HMAX_CAL_COLUMNS=4)
232 INTEGER*4 HCAL_TOT_HITS
233 INTEGER*4 HCAL_COLUMN(HMAX_CAL_BLOCKS)
234 INTEGER*4 HCAL_ROW(HMAX_CAL_BLOCKS)
235 INTEGER*4 HCAL_ADC(HMAX_CAL_BLOCKS)
236 REAL*4 HCAL_REALADC(HMAX_CAL_BLOCKS)
237
238 COMMON/HMS_RAW_CAL/
239 & HCAL_COLUMN,
240 & HCAL_ROW,
241 & HCAL_ADC,
242 cdaq 1.1 & HCAL_TOT_HITS,
243 & HCAL_REALADC
244 *
245 * DECODED CALORIMETER QUANTITIES
246 * FILLED BY H_TRANS_CAL USING ONLY THE CALORIMETER INFORMATION
247 *
248 * X COORDINATE OF BLOCK CENTER
249 * Z COORDINATE OF BLOCK CENTER
250 * ENERGY DEPOSITED IN THE BLOCK
251 * ENERGY DEPOSITED IN COLUMN #1
252 * #2
253 * #3
254 * #4
255 * TOTAL ENERGY IN THE CALORIMETER
256 INTEGER*4 HNHITS_CAL !NUMBER OF CALORIMETER HITS ABOVE THRESHOLD
257 REAL*4 HBLOCK_XC !X COORDINATE OF BLOCK CENTER
258 REAL*4 HBLOCK_ZC !Z COORDINATE OF BLOCK CENTER
259 REAL*4 HBLOCK_DE !ENERGY DEPOSITION IN THE BLOCKS
260 REAL*4 HCAL_E1 !ENERGY DEPOSITION IN COLUMN #1
261 REAL*4 HCAL_E2 ! #2
262 REAL*4 HCAL_E3 ! #3
263 cdaq 1.1 REAL*4 HCAL_E4 ! #4
264 REAL*4 HCAL_ET !TOTAL ENERGY IN THE CALORIMETER
265 COMMON/HMS_DECODED_CAL/
266 & HBLOCK_XC(HMAX_CAL_BLOCKS),
267 & HBLOCK_ZC(HMAX_CAL_BLOCKS),
268 & HBLOCK_DE(HMAX_CAL_BLOCKS),
269 & HCAL_E1,
270 & HCAL_E2,
271 & HCAL_E3,
272 & HCAL_E4,
273 & HCAL_ET,
274 & HNHITS_CAL
275 *
276 * HMS CERENKOV HITS
277 * FILLED BY G_decode_event_by_banks
278 *
279 * THERE ARE TWO CERENKOV PHOTOTUBES. EACH HIT HAS
280 * TUBE NUMBER
281 * ADC VALUE
282 * THE TOTAL NUMBER OF PHOTOTUBE HITS IS HCER_TOT_HITS
283 * HOWEVER TO MAKE THE DECODERS SYMMETRICAL TO A SCINT WE MUST ADD
284 cdaq 1.1 * A DUMMY PLANE NUMBER.
285 INTEGER*4 HMAX_CER_HITS
286 PARAMETER(HMAX_CER_HITS=2)
287 INTEGER*4 HCER_TOT_HITS
288 INTEGER*4 HCER_TUBE_NUM
|
289 cdaq 1.3 INTEGER*4 HCER_RAW_ADC
|
290 cdaq 1.1 INTEGER*4 HCER_PLANE
291 COMMON/HMS_RAW_CER/
292 & HCER_TUBE_NUM(HMAX_CER_HITS),
|
293 cdaq 1.3 & HCER_RAW_ADC(HMAX_CER_HITS),
|
294 cdaq 1.1 & HCER_PLANE(HMAX_CER_HITS),
295 & HCER_TOT_HITS
296 *
|
297 cdaq 1.2 * DECODED CERENKOV QUANTITIES
298 * FILLED BY H_TRANS_CER
299 *
300 INTEGER*4 HCER_NUM_HITS ! NUMBER OF CERENKOV HITS ABOVE THRESHOLD
301 REAL*4 HCER_NPE ! ADC CONVERTED TO NUMBER OF PHOTOELECTRONS
|
302 saw 1.5 REAL*4 HCER_NPE_SUM ! SUM OVER TUBES OF NPE''S
|
303 cdaq 1.3 REAL*4 HCER_ADC ! PED SUBTRACTED ADC FOR EACH *TUBE*
|
304 cdaq 1.2 COMMON/HCER_DECODED_CER/
305 & HCER_NUM_HITS,
306 & HCER_NPE(HMAX_CER_HITS),
|
307 cdaq 1.3 & HCER_NPE_SUM,
308 & HCER_ADC(HMAX_CER_HITS)
|
309 cdaq 1.2 *
|
310 cdaq 1.1 *
311 *
312 * HMS DETECTOR TRACK QUANTITIES
313 * FILLED BY H_TRACK SUBROUTINE
314 *
315 INTEGER*4 HNTRACKS_MAX ! NUMBER OF TRACKS ALLOWED IN FOCAL PLANE
316 PARAMETER (HNTRACKS_MAX=5) ! SET MAXIMUM TO 5
317 INTEGER*4 HNTRACKHITS_MAX ! MAXIMUM NUMBER OF HITS IN EACH TRACK
318 PARAMETER (HNTRACKHITS_MAX=16) ! SET MAXIMUM TO 15
319 INTEGER*4 HNTRACKS_FP ! NUMBER OF FOCAL PLANE TRACKS FOUND
320 * ALL THE FOLLOWING VARIABLES ARE ARRAYS
321 *
322 REAL*4 HX_FP ! X POSITION OF TRACK IN FOCAL PLANE
323 REAL*4 HY_FP ! Y POSITION OF TRACK IN FOCAL PLANE
324 REAL*4 HZ_FP ! Z POSITION OF TRACK IN FOCAL PLANE
325 REAL*4 HXP_FP ! X SLOPE OF TRACK IN FOCAL PLANE
326 REAL*4 HYP_FP ! Y SLOPE OF TRACK IN FOCAL PLANE
327 REAL*4 HCHI2_FP ! FIT QUALITY IN FOCAL PLANE
328 REAL*4 HCHI2PERDOF_FP ! fp chi2 per degree of freedom
329 INTEGER*4 HNFREE_FP ! NUMBER OF DEGREES OF FREEDOM IN FIT
330 INTEGER*4 HNTRACK_HITS ! LIST OF HITS ON EACH TRACK
331 cdaq 1.1 REAL*4 HDEL_FP ! FOCAL PLANE ERROR MATRIX
332 REAL*4 HX_FP_rot ! x in rotated focal plane
333 REAL*4 HY_FP_rot ! y in rotated f plane
334 REAL*4 HXP_FP_rot ! x slope in rotated f plane
335 REAL*4 HYP_FP_rot ! y slope in rot f plane
336
337 * THE FOCAL PLANE ERROR MATRIX IS A 4 BY 4 BY HNTRACK_MAX ARRAY
338 * THE FOUR INDECIES FOR EACH TRACK ARE X, XP, Y, YP
339 * THE ERROR MATRIX FOR EACH TRACK IS SYMMETRIC ABOUT THE DIAGONAL
340 * DELXX DELXXP DELXY DELXYP
341 * DELXPX DELXPXP DELXPY DELXPYP
342 * DELYPX DELYXP DELYY DELYYP
343 * DELYPX DELYPXP DELYPY DELYPYP
344 COMMON/HMS_FOCAL_PLANE/
345 1 HX_FP(HNTRACKS_MAX),
346 2 HY_FP(HNTRACKS_MAX),
347 3 HZ_FP(HNTRACKS_MAX),
348 4 HXP_FP(HNTRACKS_MAX),
349 5 HYP_FP(HNTRACKS_MAX),
350 6 HCHI2_FP(HNTRACKS_MAX),
351 7 HDEL_FP(4,4,HNTRACKS_MAX),
352 cdaq 1.1 8 HNTRACK_HITS(HNTRACKS_MAX,HNTRACKHITS_MAX+1),
353 9 HNFREE_FP(HNTRACKS_MAX),
354 A HNTRACKS_FP,HCHI2PERDOF_FP(HNTRACKS_MAX),
355 > HX_FP_rot(HNTRACKS_MAX),
356 > HY_FP_rot(HNTRACKS_MAX),
357 > HXP_FP_rot(HNTRACKS_MAX),
358 > HYP_FP_rot(HNTRACKS_MAX)
359
360
361 * HMS TARGET QUANTITIES
362 * ASSUME THE SAME NUMBER OF MAXIMUM TRACKS AS IN THE FOCAL PLANE
363 * HNTRACKS_MAX
364 *
365 * FILLED BY H_TARG_TRANS SUBROUTINE
366 INTEGER*4 HNTRACKS_TAR ! NUMBER OF TARGET TRACKS FOUND
367 * ALL THE FOLLOWING VARIABLES ARE ARRAYS
368 *
369 REAL*4 HX_TAR ! X POSITION OF TRACK AT TARGET
370 REAL*4 HY_TAR ! Y POSITION OF TRACK AT TARGET
371 REAL*4 HZ_TAR ! Z POSITION OF TRACK AT TARGET
372 REAL*4 HXP_TAR ! X SLOPE OF TRACK AT TARGET
373 cdaq 1.1 REAL*4 HYP_TAR ! Y SLOPE OF TRACK AT TARGET
374 REAL*4 HDELTA_TAR ! FRACTION TRACK MOMENTUM
375 REAL*4 HP_TAR ! MOMENTUM OF TRACK AT TARGET
376 REAL*4 HCHI2_TAR ! FIT QUALITY AT TARGET
377 INTEGER*4 HNFREE_TAR ! NUMBER OF DEGREES OF FREEDOM IN FIT
378 REAL*4 HDEL_TAR ! TARGET ERROR MATRIX
379 * THE FOCAL PLANE ERROR MATRIX IS A 5 BY 5 BY HNTRACK_MAX ARRAY
380 * THE FIVE INDECIES FOR EACH TRACK ARE X, Y , XP, YP and P
381 * THE ERROR MATRIX FOR EACH TRACK IS SYMMETRIC ABOUT THE DIAGONAL
382 * DELXX DELXXP DELXY DELXYP DELXP
383 * DELXPX DELXPXP DELXPY DELXPYP DELXPP
384 * DELYX DELYXPY DELYY DELYYP DELYP
385 * DELYPX DELYPXP DELYPY DELYPYP DELYPP
386 * DELPX DELPXP DELPY DELPXP DELPP
387 INTEGER*4 HLINK_TAR_FP ! LINK POINTER TO FOCAL PLANE TRACK NUMBER
388 COMMON/HMS_TARGET/
389 1 HX_TAR(HNTRACKS_MAX),
390 2 HY_TAR(HNTRACKS_MAX),
391 3 HZ_TAR(HNTRACKS_MAX),
392 4 HXP_TAR(HNTRACKS_MAX),
393 5 HYP_TAR(HNTRACKS_MAX),
394 cdaq 1.1 6 HDELTA_TAR(HNTRACKS_MAX),
395 7 HP_TAR(HNTRACKS_MAX),
396 8 HCHI2_TAR(HNTRACKS_MAX),
397 9 HDEL_TAR(5,5,HNTRACKS_MAX),
398 A HNFREE_TAR(HNTRACKS_MAX),
399 B HLINK_TAR_FP(HNTRACKS_MAX),
400 C HNTRACKS_TAR
401 *
402 *
403 * HMS_TRACK_TESTS
404 *
405 * PARTICLE ID INFORMATION FILLED BY H_TOF and H_CAL
406 * THIS STORES THE RESULTS OF SHOWER AND SCINTILLATOR CALCULATIONS
407 * FOR EACH OF THE TRACKS GENERATED BY H_TRACK
408 *
409 INTEGER*4 HNBLOCKS_CAL(HNTRACKS_MAX)
410 REAL*4 HTRACK_E1(HNTRACKS_MAX)
411 REAL*4 HTRACK_E2(HNTRACKS_MAX)
412 REAL*4 HTRACK_E3(HNTRACKS_MAX)
413 REAL*4 HTRACK_E4(HNTRACKS_MAX)
414 REAL*4 HTRACK_ET(HNTRACKS_MAX)
415 cdaq 1.1 REAL*4 HTRACK_PRESHOWER_E(HNTRACKS_MAX)
416 *
417 INTEGER*4 HSCIN_HIT(HNTRACKS_MAX,HMAX_SCIN_HITS)
418 INTEGER*4 HNUM_SCIN_HIT(HNTRACKS_MAX)
419 INTEGER*4 HNUM_PMT_HIT(HNTRACKS_MAX)
420 REAL*4 HDEDX(HNTRACKS_MAX,HMAX_SCIN_HITS)
421 REAL*4 HBETA(HNTRACKS_MAX)
422 REAL*4 HBETA_CHISQ(HNTRACKS_MAX)
423 REAL*4 HTIME_AT_FP(HNTRACKS_MAX)
|
424 cdaq 1.3 REAL*4 HSCIN_FPTIME(HNTRACKS_MAX,HMAX_SCIN_HITS)
|
425 cdaq 1.1 *
426 COMMON/HMS_TRACK_TESTS/
427 1 HTRACK_E1, ! ENERGY IN FIRST CALORIMETER COLUMN NEAR TRACK
428 2 HTRACK_E2, ! SECOND COLUMN
429 3 HTRACK_E3, ! THIRD COLUMN
430 4 HTRACK_E4, ! FOURTH COLUMN
431 5 HTRACK_ET, ! TOTAL SHOWER ENERGY ALONG TRACK
432 6 HTRACK_PRESHOWER_E, ! PRESHOWER ENERGY-Note definition is exp. dep.
433 7 HDEDX, !
434 8 HBETA, ! VELOCITY OF TRACK
435 9 HBETA_CHISQ, ! CHISQ OF FIT TO BETA OF TRACK
436 A HTIME_AT_FP, !
437 B HNBLOCKS_CAL, ! NUMBER OF SHOWER BLOCKS ON EACH TRACK
438 C HSCIN_HIT, ! ARRAY OF SCIN HITS ASSOCIATED WITH
439 D ! EACH TRACK
440 E HNUM_SCIN_HIT, ! NUMBER OF HITS FOR EACH TRACK
|
441 cdaq 1.3 F HNUM_PMT_HIT, ! NUMBER OF PMT HITS FOR EACH TRACK
442 G HSCIN_FPTIME
|
443 cdaq 1.1 *
444 * HMS SINGLES PHYSICS COMMON BLOCKS
445 * THESE ARE FILLED BY H_PHYSICS
446 *
447 *
448 REAL*4 HSP ! Lab momentum of chosen track in GeV/c
449 REAL*4 HSENERGY ! Lab total energy of chosen track in GeV
450 REAL*4 HSDELTA ! Spectrometer delta of chosen track
451 REAL*4 HSTHETA ! Lab Scattering angle in radians
452 REAL*4 HSPHI ! Lab Azymuthal angle in radians
453 REAL*4 HSMINV ! Invariant Mass of remaing hadronic system
454 REAL*4 HSZBEAM ! Lab Z coordinate of intersection of beam
455 ! track with spectrometer ray
456 REAL*4 HSDEDX(4) ! DEDX of chosen track in each scin plane
457 REAL*4 HSBETA ! BETA of chosen track
458 REAL*4 HSTRACK_ET ! Total shower energy of chosen track
459 REAL*4 HSTRACK_PRESHOWER_E ! preshower of chosen track
460 REAL*4 HSTIME_AT_FP
461 REAL*4 HSX_FP ! X focal plane position
462 REAL*4 HSY_FP
463 REAL*4 HSXP_FP
464 cdaq 1.1 REAL*4 HSYP_FP
465 REAL*4 HSCHI2PERDEG ! CHI2 per degree of freedom of chosen track.
466 REAL*4 HSX_TAR
467 REAL*4 HSY_TAR
468 REAL*4 HSXP_TAR
469 REAL*4 HSYP_TAR
470 REAL*4 HSBETA_CHISQ
471 *
472 real*4 hsmass2 ! Mass squared
473 real*4 hst ! invariant t
474 real*4 hsu ! invariant u
475 real*4 hseloss
476 real*4 hsq3 ! Lab frame momentum transfer
477 real*4 hsthetaq, hsphiq ! Direction of q3
478 real*4 hsbigq2 ! Q**2
479 real*4 hsx ! fraction of nucleon p carried by quark
|
480 saw 1.5 real*4 hsy ! fraction of lepton''s E lost in lab
|
481 cdaq 1.1 real*4 hsw ! Invariant mass of recoil system
482 real*4 hsw2 ! Invariant mass**2 of recoil system
483 *
484 INTEGER*4 HSNUM_FPTRACK ! Index of focal plane track chosen
485 INTEGER*4 HSNUM_TARTRACK ! Index of target track chosen
486 INTEGER*4 HSID_LUND ! LUND particle ID code -- not yet filled
487 INTEGER*4 HSNFREE_FP
488 INTEGER*4 HSNUM_SCIN_HIT ! # OF SCINTILLATORS ON TRACK
489 INTEGER*4 HSNUM_PMT_HIT ! # OF HODOSCOPE PMTS ON TRACK
490 *
491 COMMON/HMS_PHYSICS_R4/
492 & HSP,
493 & HSENERGY,
494 & HSDELTA,
495 & HSTHETA,
496 & HSPHI,
497 & HSMINV,
498 & HSZBEAM,
499 & HSDEDX,
500 & HSBETA,
501 & HSTRACK_ET,
502 cdaq 1.1 & HSTRACK_PRESHOWER_E,
503 & HSTIME_AT_FP,
504 & HSX_FP ,
505 & HSY_FP ,
506 & HSXP_FP ,
507 & HSYP_FP ,
508 & HSCHI2PERDEG ,
509 & HSX_TAR ,
510 & HSY_TAR ,
511 & HSXP_TAR ,
512 & HSYP_TAR ,
513 & HSBETA_CHISQ,
514 & hsmass2,
515 & hst,
516 & hsu,
517 & hseloss,
518 & hsq3,
519 & hsthetaq, hsphiq,
520 & hsbigq2,
521 & hsx,
522 & hsy,
523 cdaq 1.1 & hsw,
524 & hsw2
525 *
526 COMMON/HMS_PHYSICS_I4/
527 & HSNUM_FPTRACK,
528 & HSNUM_TARTRACK,
529 & HSID_LUND,
530 & HSNFREE_FP,
531 & HSNUM_SCIN_HIT,
532 & HSNUM_PMT_HIT
533 *
534 *
535 * Non-Hits data with HMS gates or starts.
536 * (Energy SUMS, logic timing,...)
537 *
538 INTEGER HMAX_MISC_HITS
|
539 cdaq 1.3 PARAMETER(HMAX_MISC_HITS=100)
|
540 cdaq 1.1 INTEGER*4 HMISC_TOT_HITS
|
541 saw 1.4 INTEGER*4 HMISC_RAW_ADDR1 ! "Plane" (1=TDC,2=ADC)
|
542 cdaq 1.1 INTEGER*4 HMISC_RAW_ADDR2 ! "Counter"
543 INTEGER*4 HMISC_RAW_DATA
544 COMMON/H_RAW_MISC/
545 & HMISC_RAW_ADDR1(HMAX_MISC_HITS),
546 & HMISC_RAW_ADDR2(HMAX_MISC_HITS),
547 & HMISC_RAW_DATA(HMAX_MISC_HITS),
548 & HMISC_TOT_HITS
549 *
550 *
551 *
552 *******************end: hms_data_structures.cmn*************************
553 * Local Variables:
554 * mode: fortran
555 * fortran-continuation-string: "&"
556 * comment-column: 35
557 * End:
558
|