1 cdaq  1.1 *     hms_geometry.cmn
  2           *
  3           *     This include file has all the geometrical coefficients for the
  4           *     HMS wire chambers.
  5           *
  6           *     d.f. geesaman        1 September 1993
  7           *     modifed    dfg       14 Feb 1994
  8           *                change HPLANE_PARAM(2,)  to hdc_zpos
  9           *                change HPLANE_PARAM(3,)  to hdc_alpha_angle
 10           *                change HPLANE_PARAM(4,)  to hdc_beta_angle
 11           *                change HPLANE_PARAM(5,)  to hdc_gamma_angle
 12           *                change HPLANE_PARAM(6,)  to hdc_pitch
 13           *                change HPLANE_PARAM(7,)  to hdc_nrwire
 14           *                change HPLANE_PARAM(8,)  to hdc_central_wire
 15           *                change HPLANE_PARAM(9,)  to hdc_sigma
 16           *                change HPLANE_LABEL      to hdc_plane_name
 17           *                add     hdc_xcenter
 18           *                        hdc_ycenter

 19 cdaq  1.2 * $Log: hms_geometry.cmn,v $

 20 cdaq  1.8 * Revision 1.7  1994/09/13  19:19:37  cdaq
 21           * (JRA) Add global Wire chamber offsets
 22           *

 23 cdaq  1.7 * Revision 1.6  1994/08/05  17:54:15  cdaq
 24           * (SAW) Add makereg directive with required include files
 25           *       Add "CTPTYPE=parm" directive for auto generation of CTP reg calls
 26           *

 27 cdaq  1.6 * Revision 1.5  1994/07/27  19:15:18  cdaq
 28           * (DFG) Add hdrift_x0_break and hdrift_t0_break (probably obsoleted by DJM's
 29           *       Change to h_drift_dist_calc
 30           *

 31 cdaq  1.5 * Revision 1.4  1994/06/14  03:21:25  cdaq
 32           * (DFG) Added hdc_tdc_time_zero
 33           *

 34 cdaq  1.4 * Revision 1.3  1994/04/13  17:33:06  cdaq
 35           * (DFG) Added dummy and dmytst
 36           *

 37 cdaq  1.3 * Revision 1.2  1994/03/24  18:36:06  cdaq
 38           * (DFG) Additional parameters
 39           *

 40 cdaq  1.2 * Revision 1.1  1994/02/22  14:46:27  cdaq
 41           * Initial revision
 42           *                    

 43 cdaq  1.6 *     The following include statments must precede the inclusion of this
 44           *     file in each routine that uses it.  The *%% syntax is also a
 45           *     directive to makereg to tell it to include the code in the program
 46           *     that it generates.

 47 cdaq  1.1 *

 48 cdaq  1.6 *%%   include 'gen_data_structures.cmn'

 49 cdaq  1.1 *

 50 cdaq  1.6 *     CTPTYPE=parm             ! Probably don't really want the following registered

 51 cdaq  1.1 *
 52                 real*4 hzchi,hzpsi       ! geometrical coefficients defining z-z0
 53                 real*4 hxchi,hxpsi       ! x and y for each wire plane. 
 54                 real*4 hychi,hypsi
 55           *
 56                 real*4 hz0               ! z coordinate of intersection of chamber with
 57                                          ! z axis
 58                 real*4 hpsi0,hchi0,hphi0 ! psi, chi and phi  coordinates  where a 
 59                                          ! chamber normal passing through the origin
 60                                          ! intersects the chamber. Used in stub fits
 61                 
 62                 real*4 hstubcoef         ! coefficents used in stub fits
 63                                          ! note these contain one power of sigma
 64           
 65                 real*4 htanbeta,hsinbeta,hcosbeta
 66           *
 67                 real*4 hxsp              ! coefficents used in space point fits.
 68                 real*4 hysp
 69           *
 70                 real*8 hplane_coeff      ! coefficients used in final track fit
 71           *
 72 cdaq  1.1       integer*4 HNUM_PLANE_COEFF   ! number of plane track fit coefficients
 73                 parameter (HNUM_PLANE_COEFF=9)
 74                 integer*4 HNUM_RAY_PARAM     ! number of ray parameters
 75                 parameter (HNUM_RAY_PARAM=4)
 76           *
 77                 common/HMS_GEOMETRY/

 78 cdaq  1.2      &     hplane_coeff(HNUM_PLANE_COEFF,HMAX_NUM_DC_PLANES),
 79                &     hzpsi(HMAX_NUM_DC_PLANES),hzchi(HMAX_NUM_DC_PLANES),
 80                &     hxpsi(HMAX_NUM_DC_PLANES),hxchi(HMAX_NUM_DC_PLANES),
 81                &     hypsi(HMAX_NUM_DC_PLANES),hychi(HMAX_NUM_DC_PLANES),
 82                &     hz0(HMAX_NUM_DC_PLANES),hpsi0(HMAX_NUM_DC_PLANES),
 83                &     hchi0(HMAX_NUM_DC_PLANES),hphi0(HMAX_NUM_DC_PLANES),
 84                &     hstubcoef(HMAX_NUM_DC_PLANES,HNUM_RAY_PARAM),
 85                &     hxsp(HMAX_NUM_DC_PLANES),hysp(HMAX_NUM_DC_PLANES),
 86                &     htanbeta(HMAX_NUM_DC_PLANES),hsinbeta(HMAX_NUM_DC_PLANES),
 87                &     hcosbeta(HMAX_NUM_DC_PLANES)

 88 cdaq  1.1 *
 89                 real*4    hdc_zpos
 90                 real*4    hdc_alpha_angle
 91                 real*4    hdc_beta_angle
 92                 real*4    hdc_gamma_angle
 93                 real*4    hdc_pitch
 94                 real*4    hdc_nrwire
 95                 real*4    hdc_central_wire
 96                 real*4    hdc_sigma
 97                 real*4    hdc_xcenter
 98                 real*4    hdc_ycenter

 99 cdaq  1.7       real*4    hdc_center

100 cdaq  1.1       integer*4 hdc_chamber_planes
101                 character*16 hdc_plane_name
102                 common/HMS_PLANE_PARAMETERS/

103 cdaq  1.2      &     hdc_zpos(HMAX_NUM_DC_PLANES),
104                &     hdc_alpha_angle(HMAX_NUM_DC_PLANES),
105                &     hdc_beta_angle(HMAX_NUM_DC_PLANES),
106                &     hdc_gamma_angle(HMAX_NUM_DC_PLANES),
107                &     hdc_pitch(HMAX_NUM_DC_PLANES),
108                &     hdc_central_wire(HMAX_NUM_DC_PLANES),
109                &     hdc_nrwire(HMAX_NUM_DC_PLANES),
110                &     hdc_sigma(HMAX_NUM_DC_PLANES),

111 cdaq  1.8      &     hdc_xcenter(HMAX_NUM_CHAMBERS),
112                &     hdc_ycenter(HMAX_NUM_CHAMBERS),

113 cdaq  1.2      &     hdc_chamber_planes(HMAX_NUM_DC_PLANES),

114 cdaq  1.7      &     hdc_plane_name(HMAX_NUM_DC_PLANES),
115                &     hdc_center(HMAX_NUM_DC_PLANES)

116 cdaq  1.1 *
117           *
118                 real*4 hlocrayzt
119                 parameter (hlocrayzt=0.)
120           *

121 cdaq  1.6 *     CTPTYPE=parm

122 cdaq  1.1 *

123 cdaq  1.2 *     parameter file variables. Separate by type to make it easy to add
124           *     at the end
125           *     REAL*4
126           *

127 cdaq  1.1       real*4 hdrift_velocity        ! hms drift velocity in cm/ns

128 cdaq  1.2       real*4 hdc_tdc_time_per_channel ! hms drift chamber tdc calibration
129                 real*4 hdc_1_zpos
130                 real*4 hdc_2_zpos

131 cdaq  1.3       real*4 dummy

132 cdaq  1.4       real*4 hdc_plane_time_zero        ! zero time for drift chamber planes

133 cdaq  1.5       real*4 hdrift_t0_break            ! start of linear drift region in t
134                 real*4 hdrift_x0_break            ! start of linear drift region in x

135 cdaq  1.2       common/HMS_CHAMBER_READOUT_REAL/
136                &     hdrift_velocity,
137                &     hdc_tdc_time_per_channel,
138                &     hdc_1_zpos,

139 cdaq  1.3      &     hdc_2_zpos,

140 cdaq  1.4      &     dummy,

141 cdaq  1.5      &     hdc_plane_time_zero(HMAX_NUM_DC_PLANES),
142                &     hdrift_t0_break(HMAX_NUM_DC_PLANES),
143                &     hdrift_x0_break(HMAX_NUM_DC_PLANES)

144 cdaq  1.2 *

145 cdaq  1.3 *    INTEGER*4
146                 integer*4 hdc_wire_counting     ! readout numbering order
147                 integer*4 dmytst

148 cdaq  1.2 *
149                 common/HMS_CHAMBER_READOUT_INT/

150 cdaq  1.3      &     hdc_wire_counting(HMAX_NUM_DC_PLANES),
151                &     dmytst
152           
153