(file) Return to hms_data_structures.cmn CVS log (file) (dir) Up to [HallC] / Analyzer / INCLUDE

  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           

Analyzer/Replay: Mark Jones, Documents: Stephen Wood
Powered by
ViewCVS 0.9.2-cvsgraph-1.4.0