Initial files

.gitignore

+*~
+ROOTFILES/*

DBASE/general.param

+#real raddeg
+raddeg=3.14159265/180
+;
+#include "PARAM/hcana.param"
+; hmsflags.param include spectrometer offsets and options.
+#include "PARAM/genflags.param"
+#include "PARAM/hmsflags.param"
+
+#include "PARAM/gdebug.param"
+#include "PARAM/hdebug.param"
+#include "PARAM/htracking.param"
+
+#include "PARAM/gtarget.param"
+#include "PARAM/gbeam.param"
+;  
+#include "PARAM/hdc_offsets.param"
+#include "PARAM/hdc.pos"
+#include "PARAM/hhodo.pos"
+#include "PARAM/hcal.pos"
+#include "PARAM/hhodo.param"
+#include "PARAM/haero.param"
+#include "PARAM/hdc.param"
+#include "PARAM/hdriftmap.param"
+#include "PARAM/hcal.param"
+#include "PARAM/hcer.param"

DBASE/standard.database

+0-99999
+g_ctp_parm_filename ="PARAM/general.param"
+g_decode_map_filename ="MAPS/hms.map"
+g_ctp_kinematics_filename ="DBASE/standard.kinematics"

PARAM/gbeam.param

+;                    BPM calibration constants
+;                    =========================
+  guse_bpm_in_recon  = 0   ; if 1 use bpm information for reconstruction
+  guse_bpmc = 1  ; if 1 use all 3 bpm info, if 0 use bpm A and B 
+  gbpm_sample  = 500 ; number of events, which are used for average 
+                     ; beam position ( <5000 ). Optimal value dep. on rate
+
+; the names are based on Paul Gueye's 'Status of the actual Beam Position
+;                   Monitors in the Hall C Beamline', December 1, 1995.
+; (mkj) apr-4-03 kappa,alpha and off for 3rd bpm are guesses,
+; need to find right values.
+;
+;  gbpm_kappa    = 1.85 ,1.85,1.85     ; sensitivity in cm
+gbpm_kappa    = 2.00 ,1.95,1.85  ; modified mkj 4/9/03
+  gbpm_alpha_x = 1.8192,0.7330 ,0.7   ; calibration gain: we may get them from the
+  gbpm_alpha_y = 1.0063,0.8935 ,0.7   ;  EPICS events, right now that's according to Paul's note
+
+; the following offsets get added to the calculated positions
+;  gbpm_x_off = -0.011,-0.004,0.      ; in cm: survey according to Paul's note
+;  gbpm_y_off = +0.052,+0.056,0.
+; From Dahlberg survey DT_C853 Apr 7, 2003
+gbpm_x_off = 0.017+.412,0.061-.353,0.037
+gbpm_y_off = 0.037+.381,-0.04-.211,0.028
+
+; average beam positions: only used if 'guse_bpm_in_recon=0' 
+;  gbeam_xoff  = +0.18
+  gbeam_xoff  = +0.00
+  gbeam_xpoff = +0.00
+  gbeam_yoff  = +0.00
+  gbeam_ypoff = +0.00
+
+; spectrometers would like to see this positions (from optics studies) 
+; we treat SOS and HMS the same, however, we may change this if necessary
+;  gspec_xoff  = +0.18
+  gspec_xoff  = +0.00
+  gspec_xpoff = +0.00
+  gspec_yoff  = +0.00
+  gspec_ypoff = +0.00
+
+; Pedestals for BPM ADCs: from cosmic run #10933, 10/3/96, book XI-137
+  gbpm_xp_ped = 431.9,331.0,0.  ; these pedestals have to be determined
+  gbpm_xm_ped = 514.7,350.6,0.  ; from runs with no beam in the cavities,
+  gbpm_yp_ped = 406.9,358.8,0.  ; e.g. cosmic runs.
+  gbpm_ym_ped = 499.7,292.7,0. ;
+
+; positions of BPMs relative to target (from Paul's note)
+  gbpm_zpos = 345.5,166.3,0.   ; cm
+
+
+;             Fast Raster calibration constants
+;             =================================
+
+; Various fast raster quantities: gUse* are flags
+
+  gusefr               = 1  ; if 1 correct for FRY in reconstruction
+
+  guse_frdefault       = 1  ; if 1 do no phase correction (default)
+                                ;if 0 apply phase correction 
+;      gfr_cal_mom      = 3.245  ; = beam momnetum during calibration run 
+;      gfrx_adcpercm    = 11328. ; = FR channels per cm deflection on target
+;      gfry_adcpercm    = 10708. ; from run 9981, August 25, 1996, book X-22
+
+      gfr_cal_mom      = 2.038  ; = beam momnetum during calibration run 
+      gfrx_adcpercm    = 4364.7 ; = FR channels per cm deflection on target
+      gfry_adcpercm    = 5471.9 ; from harp scan (H00A), June 19, 2004
+
+
+      gfrx_dphase      = 0.0    ; phase shift (never measured for x so far) 
+      gFrx_synccut     = 0.     ; should be zero if applied on centered distr.
+      gfrx_adcmax      = 1000   ; ADC amplitude in channels.
+      gfrx_maxsize     = 0.1    ; fast raster amplitude in centimeter.
+
+      gfry_dphase      = 5.8    ; phase shift 
+      gFry_synccut     = 0.     ; should be zero if applied on centered distr.
+      gfry_adcmax      = 1000   ; ADC amplitude in channels.
+      gfry_maxsize     = 0.1    ; fast raster amplitude in centimeter.
+
+; The latest FR phase analysis from spring '96 showed, that there is no
+; measurable phase shift. During early running (E91-13, E89-12) the FRY-phase
+; was determined to be 5.8 degree.
+
+; positions of FR magnets relative to target
+  gfrx_dist  = 2111   ; cm
+  gfry_dist  = 2071   ; cm
+
+
+
+
+
+
+
+
+
+
+
+

PARAM/gdebug.param

+  gmisc_min_peds = 500
+
+  gelossdebug = 0	;eloss debug (0=off,1=on)

PARAM/genflags.param

+;default settings for general option flags.
+  gen_run_starting_event = 1
+  gen_run_stopping_event = 0
+  gen_run_hist_dump_interval = 40000
+
+; hack_enable             run hack_anal if .ne. 0
+  hack_enable = 0
+
+; gen_eloss_enable         (0=disabled) total_eloss.f
+  gen_eloss_enable = 1
+
+

PARAM/gtarget.param

+; MEC - updated for E02-109/E04-001  1/8/05
+;
+; This is the target parameter file.  Numbers 1-10 are solid targets,
+; 11-16 are cryotargets, and 17-20 are dummy/optics.
+;
+;The target order is:
+;
+;       Au_1%         C_2%        Fe_1%
+;       C_1%          C_5%        C_6%
+;       Cu_6%         Be_3%       Au_6%
+;       ---
+;       H (4cm)       ---
+;       ---           ---
+;       D (4cm)       ---
+;       Dummy (4cm)
+;       optics (8cm)
+;       optics (4 foil)
+;       NO target
+
+
+; non-existant targets have A=Z=M=...=0
+; dummy targets have A=Z=M=...=0
+; empty/dummy cells have A=Z=M=...= AS ALUMINUM (MAY NEED TO CHANGE FOR DUMMY DATA!!!)
+;
+;
+; 1=tuna can, 2=beer can, 21=solid
+ gtarg_type   =  21.     ,    21.     ,   21.
+	         21.     ,    21.     ,   21.
+                 21.     ,    21.     ,   21.
+                 21.
+                 1.      ,     1.
+                 1.      ,     1.
+                 1.      ,     1.
+                 21.     ,    21.
+                 21.     ,    21.
+
+  gtarg_z     = 79.0     ,     6.0    ,   26.0
+                 6.0     ,     6.0    ,    6.0
+                29.0     ,     0.0    ,   79.0
+                 0.0
+                 1.0     ,     0.0
+                 0.0     ,     0.0
+                 1.0     ,     0.0
+                 13.0    ,     0.0
+		 0.0     ,     0.0
+
+  gtarg_a     = 197.0    ,    12.0    ,    57.0
+                12.0     ,    12.0    ,    12.0
+                64.0     ,     0.0    ,   197.0
+                 0.0
+                 1.0     ,     0.0
+                 0.0     ,     0.0
+                 2.0     ,     0.0
+                 27.0    ,     0.0
+		 0.0     ,     0.0
+
+  gtarg_mass  = 196.9237    ,   12.0107   ,    56.85
+                12.0107     ,   12.0107   ,    12.0107 
+                63.546      ,   0.0  ,        196.9237    
+                 0.0
+                 1.00727647 ,    0.0 
+                 0.0   ,    0.0
+                 2.01355322,     0.0
+                 26.98       ,   0.0
+                 0.0        ,    0.0
+
+
+; Next is the radiation length in percent...  SOLID=APPROXIMATE, CRYO=GUESSES!!!!!
+  gtarg_lrad  =  1.0      ,    2.0     ,     1.0
+                 1.0      ,    5.2     ,     6.4
+                 6.0      ,    3.0     ,     6.0
+                 0.0 
+                 0.46     ,    0.0
+                 0.0      ,    0.0
+                 0.55     ,    0.0
+                 2.2      ,    0.0
+                 0.0      ,    0.0
+
+
+; Next is the thickness of the target in g/cm^2...
+; NOTE THE 3HE,4HE USE THICKNESS FROM KAON EXPERIMENT, WHICH HAD DIFFERENT CELLS!!!
+
+  gtarg_thick =  0.0417 ,    0.6706    ,     0.119
+                 0.346  ,    1.79      ,     2.735
+                 0.7986 ,    0.5       ,     0.3795
+                 0.0
+                 0.289  ,    0.0
+                 0.0    ,    0.0
+                 0.648  ,    0.0
+                 0.526  ,    0.0
+		 0.0    ,    0.0
+
+
+; Finally, the density of the target in g/cm^3.
+
+; Densities taken from http://pdg.lbl.gov/AtomicNuclearProperties
+;  used "Carbon (compact)" for Carbon - could be way off
+;  3He and 4He come from taking the total thickness (above) / 4cm - thickness could be way off
+
+  gtarg_dens  =  19.320    ,  2.265   ,     7.87
+                 2.265     ,  2.265   ,     2.265
+                 8.960     ,  0.0     ,     19.320
+                 0.0
+                 0.0723    ,  0.0
+                 0.0       ,  0.0
+                 0.167     ,  0.0
+		 2.265     ,  0.0
+                 0.0       ,  0.0
+
+; Target angle:  =90 deg. is perpendicluar to beam
+;                <90 deg. is face towards HMS
+;                >90 deg. is face towards SOS 
+; for E99118 target was rotated 20.3deg. to face HMS =>69.7deg. = 1.2164945rad
+  gtarg_theta =  1.570796327 
+
+
+;general target energy loss parameters. Vansyoc Mar.98
+
+;target cell****************
+; JRA: 09/19/04: Updated based on Meekins' report, June 30, 2004
+; 		gcell_radius was 1.0 - don't know why.  Changed to 2.0
+;  DID NOT UPDATED WALL/END/FRONT THICKNESSES - Don't know the units :(
+
+      gcell_radius = 2.008
+      gz_cell 	 =  13.0
+      ga_cell 	 =  27.0
+      gcell_den 	=   2.7
+      gwall_thk	 =  0.03429
+      gend_thk  	=   0.03429
+      gfront_thk  =  0.020574   
+
+;Air gap between the chamber and the entrance window
+      gair_dens 	=   0.00121
+      gair_thk 	=   0.018
+      gair_a 	=   14.68
+      gair _z 	=   7.32
+
+
+;HMS********
+;    HMS scattering chamber window specs.********** 
+      hscat_win_thk  = 0.109728
+      hscat_win_den =  2.70
+      hscat_win_z   =  13.0
+      hscat_win_a   =  27.0
+;    HMS entrance window specs.********************
+      hdet_ent_thk  =  0.049098
+      hdet_ent_den  =  0.878636
+      hdet_ent_z    =  2.67
+      hdet_ent_a    =  4.67
+;SOS********
+;    SOS scattering chamber window specs.**********
+      sscat_win_thk =  0.054864
+      sscat_win_den =  2.70
+      sscat_win_z   =  13.0
+      sscat_win_a   =  27.0
+;    SOS entrance window specs.********************
+      sdet_ent_thk  =  0.016421
+      sdet_ent_den  =  0.862
+      sdet_ent_z    =  2.67
+      sdet_ent_a    =  4.67
+;eloss debug   0 = off,1 = on
+;      gelossdebug = 0
+;use old cryo geometry? (1)
+;      old_tgeom = 1
+
+;
+
+
+
+

PARAM/haero.param

+     haero_pos_gain   = 1./168.3, 1./233.9, 1./181.3, 1./169.2, 1./177.6, 1./155.0, 1./169.5, 1./149.7,
+
+     haero_neg_gain   = 1./248.8, 1./186.4, 1./236.4, 1./172.0,
+                        1./205.3, 1./198.1, 1./227.9, 1./229.2,
+
+     haero_pos_ped_limit =  1000., 1000., 1000., 1000.,
+                            1000., 1000., 1000., 1000.,
+
+     haero_neg_ped_limit =  1000., 1000., 1000., 1000.,
+                            1000., 1000., 1000., 1000.,
+

PARAM/hcal.param

+; Slop in x position.  Allowed distance between track and edge of block (in cm)
+; too tight djm hcal_slop = 7.5
+hcal_slop = 7.5
+
+;Turn on HMS cal. fiducial volume cut. 0="no cut"
+;Default hcal_fv_test=0
+hcal_fv_test = 0
+
+hcal_pos_cal_const =0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+
+hcal_neg_cal_const =0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+                    0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001,0.001
+
+hcal_pos_gain_ini  =1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+
+hcal_neg_gain_ini  =1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+
+hcal_neg_gain_cur  =1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+                    1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000
+
+;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+hcal_pos_ped_limit =1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+
+hcal_neg_ped_limit =1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+                    1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000,1000
+
+;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+; Calibration constants for run   53049,  71823 events processed
+
+hcal_pos_gain_cor= 0.289, 0.233, 0.347, 0.388, 0.294, 0.284, 0.346, 0.392, 0.356, 0.361, 0.320, 0.312, 0.473,
+                   0.158, 0.353, 0.311, 0.368, 0.340, 0.349, 0.325, 0.366, 0.317, 0.297, 0.332, 0.557, 0.400,
+                   0.630, 0.641, 0.611, 0.893, 0.613, 0.692, 0.468, 0.566, 0.744, 0.731, 0.499, 0.590, 0.751,
+                   0.415, 0.585, 0.624, 0.721, 0.690, 0.782, 0.857, 0.777, 0.834, 0.602, 0.741, 0.739, 0.686,
+hcal_neg_gain_cor= 0.577, 0.385, 0.347, 0.400, 0.294, 0.182, 0.440, 0.371, 0.329, 0.312, 0.271, 0.296, 0.262,
+                   0.604, 0.409, 0.427, 0.559, 0.447, 0.423, 0.363, 0.403, 0.457, 0.334, 0.321, 0.416, 0.531,
+                   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,
+                   0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000,

PARAM/hcal.pos

+; neglect extra tubes for now.
+hcal_num_neg_columns = 2
+
+; Z position of front of shower counter layers
+; guesses, based on the fact that the back scin plane is supposedly at 322.0
+;  (add 25 to each from original value) - JRA
+      hcal_1pr_zpos = (350.0-11.31)
+      hcal_2ta_zpos = (361.0-11.31)
+      hcal_3ta_zpos = (372.0-11.31)
+      hcal_4ta_zpos = (383.0-11.31)
+; Thickness of shower counter blocks, blocks are 10 * 10 * 70 cm^3
+      hcal_1pr_thick = 10.0
+      hcal_2ta_thick = 10.0
+      hcal_3ta_thick = 10.0
+      hcal_4ta_thick = 10.0
+; Number of shower counter blocks per layer
+      hcal_1pr_nr = 13
+      hcal_2ta_nr = 13
+      hcal_3ta_nr = 13
+      hcal_4ta_nr = 13
+; X,Y positions of shower counter blocks
+      hcal_1pr_left  =  35.0
+      hcal_1pr_right = -35.0
+      hcal_1pr_top   = (-69.0-1.4)
+                       (-59.0-1.4)
+                       (-49.0-1.4)
+                       (-39.0-1.4)
+                       (-29.0-1.4)
+                       (-19.0-1.4)
+                       ( -9.0-1.4)
+                       (  1.0-1.4)
+                       ( 11.0-1.4)
+                       ( 21.0-1.4)
+                       ( 31.0-1.4)
+                       ( 41.0-1.4)
+                       ( 51.0-1.4)
+      hcal_2ta_left  =  35.0
+      hcal_2ta_right = -35.0
+      hcal_2ta_top   = (-69.0-1.4)
+                       (-59.0-1.4)
+                       (-49.0-1.4)
+                       (-39.0-1.4)
+                       (-29.0-1.4)
+                       (-19.0-1.4)
+                       ( -9.0-1.4)
+                       (  1.0-1.4)
+                       ( 11.0-1.4)
+                       ( 21.0-1.4)
+                       ( 31.0-1.4)
+                       ( 41.0-1.4)
+                       ( 51.0-1.4)
+      hcal_3ta_left  =  35.0
+      hcal_3ta_right = -35.0
+      hcal_3ta_top   = (-69.0-1.4)
+                       (-59.0-1.4)
+                       (-49.0-1.4)
+                       (-39.0-1.4)
+                       (-29.0-1.4)
+                       (-19.0-1.4)
+                       ( -9.0-1.4)
+                       (  1.0-1.4)
+                       ( 11.0-1.4)
+                       ( 21.0-1.4)
+                       ( 31.0-1.4)
+                       ( 41.0-1.4)
+                       ( 51.0-1.4)
+      hcal_4ta_left  =  35.0
+      hcal_4ta_right = -35.0
+      hcal_4ta_top   = (-69.0-1.4)
+                       (-59.0-1.4)
+                       (-49.0-1.4)
+                       (-39.0-1.4)
+                       (-29.0-1.4)
+                       (-19.0-1.4)
+                       ( -9.0-1.4)
+                       (  1.0-1.4)
+                       ( 11.0-1.4)
+                       ( 21.0-1.4)
+                       ( 31.0-1.4)
+                       ( 41.0-1.4)
+                       ( 51.0-1.4)

PARAM/hcana.param

+;
+; Parameters that were built into Fortran analyzer that we want
+; to pass as parameters so that the resulting code can be more generic.
+;
+
+hhodo_num_planes = 4
+hhodo_plane_names = "1x 1y 2x 2y"
+
+hcal_num_layers = 4
+
+rraster_num_signals = 4
+htofusinginvadc=0
+ hcer_tot_pmts = 2
+# Exclusion band width for the calorimeter's fiducial volume.
+hcal_fv_delta = 5.
+
+# Constants for the coordiante correction of the calorimeter energy depositions
+hcal_a_cor = 200.
+hcal_b_cor = 8000.
+hcal_c_cor = 64.36, 64.36	# for positive and negative sides
+hcal_d_cor =  1.66,  1.66
+
+hcal_layer_names = "1pr 2ta 3ta 4ta"
+
+haero_num_pairs = 8
+
+# Names of planes so that parameter names can be constructed
+hdc_plane_names = "1x1 1y1 1u1 1v1 1y2 1x2 2x1 2y1 2u1 2v1 2y2 2x2"
+
+# The following were defined in REPLAY.PARAM
+h_recon_coeff_filename =    'DATFILES/hms_recon_coeff.dat'  ;hms optics matrix
+s_recon_coeff_filename =    'DATFILES/sos_recon_coeff.dat'  ;SOs optics matrix
+
+# The following are set to zero to replicate historical ENGINE behavior
+# For new analyses they should be set to 1.  If not defined here,
+# hcana will default 1, the new and correct behaviour.
+
+# If 1, Let a hit have different L/R assignment for different space points
+# instead of L/R assignment from first sp it appears in.
+hdc_fix_lr = 0
+# If 1, don't do the the propagation along the wire each time the hit
+# appears in a space point.  (Which means the correction accumulates)
+hdc_fix_propcorr = 0
+
+# SOS parameters
+shodo_num_planes = 4
+shodo_plane_names = "1x 1y 2x 2y"
+
+scal_num_layers = 4
+sdbg_init_cal=0
+# Exclusion band width for the calorimeter's fiducial volume.
+# (saw) Don't know what this should be.  Copied it from HMS.
+scal_fv_delta = 5.
+
+# Constants for the coordiante correction of the calorimeter energy depositions
+scal_a_cor = 400.
+scal_b_cor = 12000.
+scal_c_cor = -87.1628, -100.	# The positive side constants reproduce
+scal_d_cor =  1.65054,    3.	# correction in Engine to accuracy better 0.005.
+
+scal_layer_names = "1pr 2ta 3ta 4ta"
+
+# Names of planes so that parameter names can be constructed
+sdc_plane_names = "1u1 1u2 1x1 1x2 1v1 1v2 2u1 2u2 2x1 2x2 2v1 2v2"
+
+# The following were defined in REPLAY.PARAM
+
+# Fortran ENGINE only had this as a parameter for HMS.  Need it here
+# because same code used for both spectrometers
+sntracks_max_fp = 10
+
+# The following are set to zero to replicate historical ENGINE behavior
+# For new analyses they should be set to 1.  If not defined here,
+# hcana will default 1, the new and correct behaviour.
+
+# If 1, Let a hit have different L/R assignment for different space points
+# instead of L/R assignment from first sp it appears in.
+sdc_fix_lr = 0
+# If 1, don't do the the propagation along the wire each time the hit
+# appears in a space point.  (Which means the correction accumulates)
+sdc_fix_propcorr = 0
+
+# Total number of PMTs in Gas Cherenkov detector.
+scer_tot_pmts = 4
+

PARAM/hcer.param

+hcer_width = 50, 50
+
+hcer_adc_to_npe = 1/145.0, 1/122.0
+
+
+hcer_chi2max = 50.
+hcer_beta_min = 0.8
+hcer_beta_max = 1.2
+hcer_et_min = 0.95
+hcer_et_max = 5.0
+hcer_mirror_zpos = 230
+; hcer_region: 8 values for each region (1 per mirror + sum)
+; central x,y,dx,dy values and x,y,dx,dy half widths.
+hcer_region =   30,  -30,   0,
+                 0,    0,   0,
+                 0,    0,   0,
+                 0,    0,   0,
+                20,   20,  60,
+                30,   30,  30,
+                .1,   .1,  .1,
+                .1,   .1,  .1
+
+hcer_ped_limit =  1000,1000
+	       

PARAM/hdc.pos

+; Number of planes installed in HMS detector setup
+      hdc_num_planes = 12
+; Number of chambers installed in HMS detector setup
+      hdc_num_chambers = 2
+; Names of each wire plane
+;
+; Z positions of various planes in HMS chambers
+; hdc_n_zpos is the surveyed Z position of the center of chamber n.
+      hdc_1_zpos = (-40.611-11.312+0.003)
+      hdc_2_zpos = (+40.611-11.312-0.008)
+      
+      hdc_zpos   = hdc_1_zpos - 3.6000
+                   hdc_1_zpos - 1.8000
+                   hdc_1_zpos - 0.0000
+                   hdc_1_zpos + 1.8000
+                   hdc_1_zpos + 3.6000
+                   hdc_1_zpos + 5.4000
+                   hdc_2_zpos - 3.6000
+                   hdc_2_zpos - 1.8000
+                   hdc_2_zpos - 0.0000
+                   hdc_2_zpos + 1.8000
+                   hdc_2_zpos + 3.6000
+                   hdc_2_zpos + 5.4000
+; Angle alpha of wires in wire chamber planes
+        hdc_alpha_angle =  (90. - 0.071)*raddeg
+                           (0.0 - 0.071)*raddeg
+                           (74.925 - 0.071)*raddeg
+                           (105.075 - 0.071)*raddeg
+                           (0.0 - 0.071)*raddeg
+                           (90.0 - 0.071)*raddeg
+                           (89.90814 - 0.153)*raddeg
+                           (0.01611 - 0.153)*raddeg
+                           (74.85 - 0.153)*raddeg
+                           (105.05 - 0.153)*raddeg
+                           (0.01611 - 0.153)*raddeg
+                           (89.90814 - 0.153)*raddeg
+;        hdc_alpha_angle = ( 90.1541 - 0.071)*raddeg
+;                          ( -0.0408 - 0.071)*raddeg
+;                          ( 75.1267 - 0.071)*raddeg
+;                          (105.1267 - 0.071)*raddeg
+;                          ( -0.0408 - 0.071)*raddeg
+;                          ( 90.1541 - 0.071)*raddeg
+;                          ( 90.1095 - 0.153)*raddeg
+;                          ( -0.0370 - 0.153)*raddeg
+;                          ( 75.0904 - 0.153)*raddeg
+;                          (105.0904 - 0.153)*raddeg
+;                          ( -0.0370 - 0.153)*raddeg
+;                          ( 90.1095 - 0.153)*raddeg
+;
+; Angle beta of wires in wire chamber planes
+      hdc_beta_angle =    -0.041*raddeg
+                          -0.041*raddeg
+                          -0.041*raddeg
+                          -0.041*raddeg
+                          -0.041*raddeg
+                          -0.041*raddeg
+                          +0.054*raddeg
+                          +0.054*raddeg
+                          +0.054*raddeg
+                          +0.054*raddeg
+                          +0.054*raddeg
+                          +0.054*raddeg
+; Angle gamma of wires in wire chamber planes
+      hdc_gamma_angle =   +0.0376*raddeg
+                          +0.0376*raddeg
+                          +0.0376*raddeg
+                          +0.0376*raddeg
+                          +0.0376*raddeg
+                          +0.0376*raddeg
+                          -0.145*raddeg
+                          -0.145*raddeg
+                          -0.145*raddeg
+                          -0.145*raddeg
+                          -0.145*raddeg
+                          -0.145*raddeg
+; Pitch
+      hdc_pitch = 1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+                  1.000252
+; Number of wires per plane
+      hdc_nrwire = 113
+                    52
+                   107
+                   107
+                    52
+                   113
+                   113
+                    52
+                   107
+                   107
+                    52
+                   113
+; X,Y position of center of wire chamber
+; increasing dc2 coordinates decreases dpos histograms.
+; new values according to henk and rolf
+      hdc_xcenter = (1.670 - 0.0155 - 0.020)
+                    (2.758 - 0.0155 + 0.040)
+;      hdc_xcenter = (1.670 - 0.0155)
+;                    (2.758 - 0.0155)
+;      hdc_xcenter = 1.670 + 0.15 + 0.0155
+;                    2.758 + 0.16 + 0.0155
+;  Y CENTER CHANGED TO MATCH OPTICS DATA(offset of 1.1 cm in y)
+; new values included according to henk and rolf
+      hdc_ycenter = (1.443 - 0.1699 - 0.0279 - 0.04)
+                    (2.753 - 0.1699 + 0.0279 + 0.04)
+;      hdc_ycenter = (1.443 - 0.1699 - 0.0279 - 0.035)
+;                    (2.753 - 0.1699 + 0.0279 + 0.035)
+;      hdc_ycenter = (1.443 - 0.1699 - 0.0279)
+;                    (2.753 - 0.1699 + 0.0279)
+
+;      hdc_ycenter = 0.343 + 1.1 - 0.12 + 0.1699 + 0.0279
+;                    1.653 + 1.1 - 0.22 + 0.1699 - 0.0279
+;      hdc_ycenter = 0.343  THESE ARE THE ORIGINAL Y CENTERS FROM SURVEY DATA
+;                    1.653
+;      hdc_xcenter = 0.0
+;                    1.22
+;      hdc_ycenter = 0.0
+;                    1.24
+
+; Wire number of center of wire chamber
+; assume 1st wire of 1st dc of a set is closest to edge where counting starts
+       hdc_central_wire = 57.257, 26.240, 54.001, 53.999, 26.760 , 56.743
+                          57.244, 26.242, 53.998, 54.002, 26.758 , 56.756
+
+;       hdc_central_wire = 57.25, 26.245, 54.0, 54.0, 26.755 , 56.75
+;                          57.24, 26.245, 54.0, 54.0, 26.755 , 56.76
+;       hdc_central_wire = 57.25, 26.25, 54.0, 54.0, 26.75 , 56.75
+;                          57.24, 26.25, 54.0, 54.0, 26.75 , 56.76
+; hdc_chamber_planes    array giving the chamber number for each plane
+      hdc_chamber_planes = 1
+                           1
+                           1
+                           1
+                           1
+                           1
+                           2
+                           2
+                           2
+                           2
+                           2
+                           2
+; The hms readout numbers some planes in reverse order.
+; The following array is a flag on the order number.
+; If hdc_wire_counting(plane) = 0
+;  the wire center is at (wire - hdc_central_wire) * pitch
+; If hdc_wire_counting(plane) = 1
+;  the wire center is at ( hdc_nrwire + 1 - wire - hdc_central_wire) * pitch
+        hdc_wire_counting = 1
+                            1
+                            0
+                            1
+                            0
+                            0
+                            1
+                            1
+                            0
+                            1
+                            0
+                            0
+; The velocity correction is the distance from the center of the wire divided
+; by the velocity of propagation times hdc_drifttime_sign(pln).  +/-1
+; for disc. card at +/- coord. (i.e. top = -x direction, so top readout is +1)
+       hdc_drifttime_sign = -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1
+                            -1

PARAM/hdc_offsets.param

+; cable id's for HMS Drift Chambers
+; array is 113,12 (wires,planes) but most planes have less than 113 wires)
+
+hdc_card_no=
+; Plane 1, wires 1-113, 16 per card (+one extra wire)
+              1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
+              2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
+              3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3
+              4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4
+              5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5
+              6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6
+              7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7
+             12
+; plane 2, wires 1-52, 13 per card
+              8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8
+              9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9
+             10,10,10,10,10,10,10,10,10,10,10,10,10
+             11,11,11,11,11,11,11,11,11,11,11,11,11
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+; plane 3, wires 1-107, 15 or 16 per card (12 has one x wire)
+             12,12,12,12,12,12,12,12,12,12,12,12,12,12
+             13,13,13,13,13,13,13,13,13,13,13,13,13,13,13
+             14,14,14,14,14,14,14,14,14,14,14,14,14,14,14,14
+             15,15,15,15,15,15,15,15,15,15,15,15,15,15,15
+             16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16
+             17,17,17,17,17,17,17,17,17,17,17,17,17,17,17
+             18,18,18,18,18,18,18,18,18,18,18,18,18,18,18,18
+             -1,-1,-1,-1,-1,-1
+; plane 4, wires 1-107, 15 or 16 per card (19 has one x' wire)
+             19,19,19,19,19,19,19,19,19,19,19,19,19,19
+             20,20,20,20,20,20,20,20,20,20,20,20,20,20,20
+             21,21,21,21,21,21,21,21,21,21,21,21,21,21,21,21
+             22,22,22,22,22,22,22,22,22,22,22,22,22,22,22
+             23,23,23,23,23,23,23,23,23,23,23,23,23,23,23,23
+             24,24,24,24,24,24,24,24,24,24,24,24,24,24,24
+             25,25,25,25,25,25,25,25,25,25,25,25,25,25,25,25
+             -1,-1,-1,-1,-1,-1
+; plane 5, wires 1-52, 13 per card
+             26,26,26,26,26,26,26,26,26,26,26,26,26
+             27,27,27,27,27,27,27,27,27,27,27,27,27
+             28,28,28,28,28,28,28,28,28,28,28,28,28
+             29,29,29,29,29,29,29,29,29,29,29,29,29
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+; Plane 6, wires 1-113, 16 per card (+one extra wire)
+             30,30,30,30,30,30,30,30,30,30,30,30,30,30,30,30
+             31,31,31,31,31,31,31,31,31,31,31,31,31,31,31,31
+             32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32
+             33,33,33,33,33,33,33,33,33,33,33,33,33,33,33,33
+             34,34,34,34,34,34,34,34,34,34,34,34,34,34,34,34
+             35,35,35,35,35,35,35,35,35,35,35,35,35,35,35,35
+             36,36,36,36,36,36,36,36,36,36,36,36,36,36,36,36
+             19
+; Chamber 2
+; Plane 7, wires 1-113, 16 per card (+one extra wire)
+             37,37,37,37,37,37,37,37,37,37,37,37,37,37,37,37
+             38,38,38,38,38,38,38,38,38,38,38,38,38,38,38,38
+             39,39,39,39,39,39,39,39,39,39,39,39,39,39,39,39
+             40,40,40,40,40,40,40,40,40,40,40,40,40,40,40,40
+             41,41,41,41,41,41,41,41,41,41,41,41,41,41,41,41
+             42,42,42,42,42,42,42,42,42,42,42,42,42,42,42,42
+             43,43,43,43,43,43,43,43,43,43,43,43,43,43,43,43
+             48
+; plane 8, wires 1-52, 13 per card
+             44,44,44,44,44,44,44,44,44,44,44,44,44
+             45,45,45,45,45,45,45,45,45,45,45,45,45
+             46,46,46,46,46,46,46,46,46,46,46,46,46
+             47,47,47,47,47,47,47,47,47,47,47,47,47
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+; plane 9, wires 1-107, 15 or 16 per card (48 has one x wire)
+             48,48,48,48,48,48,48,48,48,48,48,48,48,48
+             49,49,49,49,49,49,49,49,49,49,49,49,49,49,49
+             50,50,50,50,50,50,50,50,50,50,50,50,50,50,50,50
+             51,51,51,51,51,51,51,51,51,51,51,51,51,51,51
+             52,52,52,52,52,52,52,52,52,52,52,52,52,52,52,52
+             53,53,53,53,53,53,53,53,53,53,53,53,53,53,53
+             54,54,54,54,54,54,54,54,54,54,54,54,54,54,54,54
+             -1,-1,-1,-1,-1,-1
+; plane 10, wires 1-107, 15 or 16 per card (55 has one x' wire)
+             55,55,55,55,55,55,55,55,55,55,55,55,55,55
+             56,56,56,56,56,56,56,56,56,56,56,56,56,56,56
+             57,57,57,57,57,57,57,57,57,57,57,57,57,57,57,57
+             58,58,58,58,58,58,58,58,58,58,58,58,58,58,58
+             59,59,59,59,59,59,59,59,59,59,59,59,59,59,59,59
+             60,60,60,60,60,60,60,60,60,60,60,60,60,60,60
+             61,61,61,61,61,61,61,61,61,61,61,61,61,61,61,61
+             -1,-1,-1,-1,-1,-1
+; plane 11, wires 1-52, 13 per card
+             62,62,62,62,62,62,62,62,62,62,62,62,62
+             63,63,63,63,63,63,63,63,63,63,63,63,63
+             64,64,64,64,64,64,64,64,64,64,64,64,64
+             65,65,65,65,65,65,65,65,65,65,65,65,65
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+             -1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1
+; Plane 12, wires 1-113, 16 per card (+one extra wire)
+             66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66
+             67,67,67,67,67,67,67,67,67,67,67,67,67,67,67,67
+             68,68,68,68,68,68,68,68,68,68,68,68,68,68,68,68
+             69,69,69,69,69,69,69,69,69,69,69,69,69,69,69,69
+             70,70,70,70,70,70,70,70,70,70,70,70,70,70,70,70
+             71,71,71,71,71,71,71,71,71,71,71,71,71,71,71,71
+             72,72,72,72,72,72,72,72,72,72,72,72,72,72,72,72
+             55
+;
+;       Cable Delays on a per/card basys
+;
+hdc_card_delay =  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0, 0, 0, 0, 0, 0, 0, 0, 0
+                  0, 0
+

PARAM/hdeadwires.param

+; list of dead wires in the HMS (plane,wire).  Dead means <20% efficiency.
+hdc_num_deadwires = 6
+
+hdc_deadwire_plane =  2,   6,   6,   7,   9,  10
+
+hdc_deadwire_num =   43,  18,  60,  97,  30,   8
+
+; rough efficiency  22%,  3%, .4%,  0%,  0%,  7%
+
+; pln9,wire99 is HOT (used to be that 9/#100 was DEAD)?

PARAM/hdebug.param

+; hms tracking code debug flags. 
+; when these flags are non-zero, dump routines are called to 
+; print out the data banks for each event.
+;
+;  hluno     fortran lun for output
+      hluno = 66
+; hdebugdumptof           Dump timing information for fitting offsets
+      hdebugdumptof = 0
+; hdebugdumpcal           Dump calorimeter information for fitting gains
+      hdebugdumpcal = 0
+; number of pedestal events required to overwrite orig. peds.
+  hhodo_min_peds = 500
+  hcal_min_peds  = 500
+  hcer_min_peds  = 500
+  haero_min_peds  = 500
+  hmisc_min_peds  = 500
+
+; DEBUGGING OUTPUT FLAGS
+; ----------------------
+; hdebugprintrawdc        Dump HMS_RAW_DC Bank
+      hdebugprintrawdc = 0
+; hdebugprintdecodeddc    Dump HMS_DECODED_DC BAnk
+      hdebugprintdecodeddc = 0
+; hdebugflagpsi               Dump calculated coordinate positions
+      hdebugflagpsi =  0
+; hdebugflaggeoemtry          Dump calculated geometrical parameters
+      hdebugflaggeometry = 0
+; hdebugflagpr                Dump intermediate pattern recognition results
+      hdebugflagpr = 0
+; hdebugflagstubchisq         Dump chi2 of left-right fits
+      hdebugstubchisq = 0
+; hdebugflagstubs             Dump stub fit       
+      hdebugflagstubs = 0
+; hdebuglinkstubs             Dump linked hits
+      hdebuglinkstubs = 0
+; hdebugtrackprint            Dump focal plane track fit results
+      hdebugtrackprint = 0
+; hdebugtartrackprint         Dump tracks at target
+      hdebugtartrackprint = 0
+; tof debuging 
+;
+;  hdebugprintscinraw         Dump HMS_RAW_SCIN
+    hdebugprintscinraw = 0
+;  hdebugprintscindec         Dump HMS_DECODED_SCIN
+    hdebugprintscindec = 0
+;  hdebugprinttoftracks       Dump  hms_scin_tof
+    hdebugprinttoftracks = 0
+;  hdebugprinttracktests      Dump HMS_TRACK_TESTS
+    hdebugprinttracktests = 0
+;
+;  hms calorimeter debug flags
+;  hlun_dbg_cal                lun
+   hlun_dbg_cal    = hluno
+;  hdbg_raw_cal        call h_prt_cal_raw
+   hdbg_raw_cal = 0
+;  hdbg_sparsified_cal call h_prt_cal_sparsified
+   hdbg_sparsified_cal = 0
+;  hdbg_decoded_cal    call h_prt_cal_decoded
+   hdbg_decoded_cal = 0
+;  hdbg_clusters_cal   call h_prt_cal_clusters
+   hdbg_clusters_cal = 0
+;  hdbg_tracks_cal     call h_prt_cal_tracks
+   hdbg_tracks_cal = 0
+;  hdbg_tests_cal      call h_prt_cal_tests    HMS_TRACK_TESTS
+   hdbg_tests_cal = 0
+;
+;  hard wired histograms. If these flags .ne. 0 then the histogram
+;     blocks are filled
+;
+; hturnon_scin_raw_hist       Histogram HMS raw hodoscope data (each block)
+  hturnon_scin_raw_hist = 1
+; hturnon_decoded_dc_hist     Histogram HMS_DECODED_DC
+  hturnon_decoded_dc_hist = 1
+; hturnon_focal_plane_hist    Histogram HMS_FOCAL_PLANE 
+  hturnon_focal_plane_hist = 1
+; hturnon_target_hist         HISTOGRAM HMS_TARGET
+  hturnon_target_hist = 1
+;
+; Bypass paramters to bypass code elementes in h_reconstruction
+;      Stored in hmsbypass_switches.cmn
+;      Code element is bypassed if switch .ne. 0
+; hbypass_trans_dc         bypass h_trans_dc
+  hbypass_trans_dc   = 0 
+; hbypass_track            bypass h_track
+  hbypass_track      = 0
+; hbypass_targ_trans       bypass h_targ_trans
+  hbypass_targ_trans = 0
+; hbypass_dc_eff           bypass h_dc_eff and h_dc_eff_shutdown
+  hbypass_dc_eff     = 0
+; hbypass_track_eff        bypass h_tracking efficiency code
+  hbypass_track_eff  = 0
+; hbypass_track_eff_files  bypass h_tracking efficiency output files
+  hbypass_track_eff_files  = 1
+
+; hbypass_trans_scin       bypass h_trans_scin
+  hbypass_trans_scin = 0
+; hbypass_tof              bypass h_tof
+  hbypass_tof        = 0
+; hbypass_scin_eff         bypass h_scin_eff and h_scin_eff_shutdown
+  hbypass_scin_eff   = 0
+
+; hbypass_trans_cer        bypass h_trans_cer
+  hbypass_trans_cer  = 0
+; hbypass_cer              bypass h_cer
+  hbypass_cer        = 0
+; hbypass_cer_eff          bypass h_cer_eff and h_cer_eff_shutdown
+  hbypass_cer_eff    = 0
+
+; hbypass_trans_cal        bypass h_trans_cal
+  hbypass_trans_cal  = 0
+; hbypass_cal              bypass h_cal
+  hbypass_cal        = 0
+; hbypass_cal_eff          bypass h_cal_eff and h_cal_eff_shutdown
+  hbypass_cal_eff    = 0
+
+; hbypass_physics          bypass h_physics
+  hbypass_physics    = 0

PARAM/hhodo.param

+; hstart_time_center  center of allowed time window                             
+   hstart_time_center = 32.                                                     
+; hstart_time_slop    1/2 width of time window                                  
+   hstart_time_slop = 25.                                                       
+; hscin_tdc_min       minimum tdc value in hms scin                             
+    hscin_tdc_min = 0                                                           
+; hscin_tdc_max       maximum allowed tdc value                                 
+   hscin_tdc_max = 4000                                                         
+; hscin_tdc_to_time   scin tdc time per channel                                 
+   hscin_tdc_to_time = 0.0259                                                   
+; tof and you figured out good values
+   htof_tolerance = 3.0
+;                                                                               
+; hms_tof_params                                                                
+; hnum_scin_counters, hhodo_zpos, hhodo_center_coord, hhodo_width               
+; are all calculated within h_init_scin                                         
+;    hhodo_pos_coord and hhodo_neg_coord are not yet used                       
+;                                                                               
+     hhodo_slop     =    2.,  2.,  4.,  4.                                      
+;                                                                               
+     hhodo_vel_light =  14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+                        14.0,14.8,15.5,15.8                                     
+;                                                                               
+      hhodo_pos_sigma = .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+;                                                                               
+      hhodo_neg_sigma = .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                        .3,.3,.3,.3,.3,.3,.3,.3                                 
+                                                                                
+  hhodo_pos_minph = 125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                                                                                
+  hhodo_neg_minph = 125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+                    125,135,125,120                                             
+  hhodo_pos_phc_coeff = -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                        -1.45,-1.1,-1.,-1.2                                     
+                                                                                
+  hhodo_neg_phc_coeff = -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                        -1.3,-1.1,-1.2,-1.4                                     
+                                                                                
+                                                                                
+; csa 9/8/98 -- I had to hand-twaddle a few of the values                       
+; based on (relative) offsets of older hhodo.param                              
+                                                                                
+                                                                                
+hhodo_pos_time_offset =   6.7992,   0.0000,  -7.1126,   4.7475
+                         -1.1436,   0.4026,  -0.1443,   8.8108
+                          5.1911,   1.7935,  -7.0357,   7.7590
+                          5.6115,   0.1231,  -4.3477,  12.8552
+                          1.2354,  -6.4896,   8.6617,   9.3857
+                          1.8675,  -0.3651,  -4.7731,   7.9034
+                          6.4738,  -1.0463,   1.1786,   8.9740
+                          4.5395,   2.3368,  -4.4623,   4.5968
+                          2.6241,  -2.4480,  -6.4964,   6.2390
+                          4.2056,   0.0000,  -0.8924,  -1.2482
+                          6.1303,   0.0000,  -7.9503,   0.0000
+                          5.5326,   0.0000,  -0.3564,   0.0000
+                          4.3238,   0.0000,   0.8133,   0.0000
+                          1.7467,   0.0000,   1.9021,   0.0000
+                          4.1333,   0.0000,  -5.4523,   0.0000
+                          0.9416,   0.0000,   1.0158,   0.0000
+hhodo_neg_time_offset =  -2.1515,   0.0000,   2.5771,  11.2323
+                         -1.6547,  -2.2888,  10.1063,  13.3595
+                         -3.8050,  -4.1933,   3.4133,  15.0118
+                          0.7925,  -1.9592,  10.9311,  15.1057
+                         -1.2621,  -4.6538,   2.3474,  10.2752
+                         -2.9459,   1.2711,   9.6479,   7.8567
+                          1.2439,  -3.4946,   3.7283,  10.7897
+                         -4.8205,   1.1743,   8.5908,  11.1753
+                         -6.3324,  -5.2290,   0.2444,  14.7184
+                          1.8168,   0.0000,  11.2387,  15.5869
+                          0.4240,   0.0000,   8.2657,   0.0000
+                         -0.9945,   0.0000,   9.2381,   0.0000
+                         -2.3217,   0.0000,   5.1855,   0.0000
+                         -5.1166,   0.0000,  10.0781,   0.0000
+                          0.2537,   0.0000,   6.2455,   0.0000
+                          7.0554,   0.0000,  12.8714,   0.0000
+                                                                                
+;                                                                               
+hhodo_pos_ped_limit = 1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+;                                                                               
+hhodo_neg_ped_limit = 1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                      1000,1000,1000,1000,1000,1000,1000,1000                   
+                                                                                

PARAM/hhodo.pos

+      hpathlength_central = 2500
+; Z positions of hodoscopes
+      hscin_1x_zpos =  (89.14-11.31)
+      hscin_1y_zpos =  (108.83-11.31)
+      hscin_2x_zpos =  (310.13-11.31)
+      hscin_2y_zpos =  (329.82-11.31)
+      hscin_1x_dzpos = 2.12
+      hscin_1y_dzpos = 2.12
+      hscin_2x_dzpos = 2.12
+      hscin_2y_dzpos = 2.12
+      hscin_1x_size = 8.0
+      hscin_1y_size = 8.0
+      hscin_2x_size = 8.0
+      hscin_2y_size = 8.0
+      hscin_1x_spacing = 7.5
+      hscin_1y_spacing = 7.5
+      hscin_2x_spacing = 7.5
+      hscin_2y_spacing = 7.5
+; Number of hodoscope paddles per layer
+      hscin_1x_nr = 16
+      hscin_1y_nr = 10
+      hscin_2x_nr = 16
+      hscin_2y_nr = 10
+; X,Y positions of hodoscope paddles
+      hscin_1x_left  =  37.750
+      hscin_1x_right = -37.750
+      hscin_1x_offset=  -1.3
+      hscin_1x_center= -56.250
+                       -48.750
+                       -41.250
+                       -33.750
+                       -26.250
+                       -18.750
+                       -11.250
+                        -3.750
+                         3.750
+                        11.250
+                        18.750
+                        26.250
+                        33.750
+                        41.250
+                        48.750
+                        56.250
+      hscin_1y_top   = -60.250
+      hscin_1y_bot   =  60.250
+      hscin_1y_offset=  -1.3
+      hscin_1y_center=  33.750
+                        26.250
+                        18.750
+                        11.250
+                         3.750
+                        -3.750
+                       -11.250
+                       -18.750
+                       -26.250
+                       -33.750
+      hscin_2x_left  =  37.750
+      hscin_2x_right = -37.750
+      hscin_2x_offset=  -0.6
+      hscin_2x_center= -56.250
+                       -48.750
+                       -41.250
+                       -33.750
+                       -26.250
+                       -18.750
+                       -11.250
+                        -3.750
+                         3.750
+                        11.250
+                        18.750
+                        26.250
+                        33.750
+                        41.250
+                        48.750
+                        56.250
+      hscin_2y_top   = -60.250
+      hscin_2y_bot   =  60.250
+      hscin_2y_offset=  -2.4
+      hscin_2y_center=  33.750
+                        26.250
+                        18.750
+                        11.250
+                         3.750
+                        -3.750
+                       -11.250
+                       -18.750
+                       -26.250
+                       -33.750

PARAM/hms.driftmap

+hwc1x1fract=0.0000,0.0001,0.0001,0.0002,0.0003,0.0004,0.0004,0.0005
+0.0007,0.0008,0.0010,0.0012,0.0015,0.0020,0.0028,0.0041,0.0060,0.0086
+0.0121,0.0167,0.0230,0.0311,0.0412,0.0529,0.0658,0.0800,0.0951,0.1114
+0.1284,0.1462,0.1645,0.1838,0.2036,0.2241,0.2451,0.2667,0.2884,0.3099
+0.3319,0.3541,0.3764,0.3989,0.4213,0.4438,0.4662,0.4888,0.5113,0.5334
+0.5555,0.5777,0.5996,0.6208,0.6419,0.6628,0.6838,0.7045,0.7245,0.7446
+0.7639,0.7826,0.8007,0.8181,0.8346,0.8500,0.8646,0.8784,0.8910,0.9029
+0.9134,0.9234,0.9324,0.9405,0.9479,0.9547,0.9608,0.9663,0.9712,0.9756
+0.9795,0.9827,0.9854,0.9878,0.9897,0.9912,0.9925,0.9936,0.9945,0.9952
+0.9957,0.9963,0.9967,0.9971,0.9974,0.9977,0.9979,0.9981,0.9983,0.9984
+0.9986,0.9987,0.9989,0.9990,0.9991,0.9992,0.9992,0.9993,0.9994,0.9995
+0.9995,0.9995,0.9996,0.9996,0.9997,0.9997,0.9997,0.9998,0.9998,0.9998
+0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc1y1fract=0.0000,0.0001,0.0001,0.0002,0.0002,0.0003,0.0004,0.0005
+0.0005,0.0006,0.0008,0.0009,0.0011,0.0014,0.0019,0.0025,0.0034,0.0048
+0.0071,0.0109,0.0172,0.0263,0.0385,0.0530,0.0697,0.0879,0.1078,0.1289
+0.1506,0.1727,0.1951,0.2180,0.2406,0.2634,0.2863,0.3092,0.3320,0.3544
+0.3769,0.3986,0.4206,0.4421,0.4634,0.4851,0.5059,0.5269,0.5479,0.5685
+0.5889,0.6096,0.6299,0.6503,0.6704,0.6900,0.7097,0.7289,0.7478,0.7663
+0.7846,0.8022,0.8192,0.8357,0.8517,0.8663,0.8803,0.8931,0.9046,0.9149
+0.9246,0.9328,0.9403,0.9471,0.9529,0.9582,0.9630,0.9672,0.9710,0.9745
+0.9777,0.9802,0.9825,0.9845,0.9864,0.9880,0.9894,0.9906,0.9916,0.9924
+0.9933,0.9940,0.9946,0.9951,0.9957,0.9961,0.9965,0.9969,0.9972,0.9975
+0.9978,0.9980,0.9982,0.9984,0.9986,0.9987,0.9989,0.9990,0.9991,0.9992
+0.9993,0.9994,0.9995,0.9995,0.9996,0.9997,0.9997,0.9998,0.9998,0.9998
+0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000,1.0000,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc1u1fract=0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0001,0.0001
+0.0001,0.0001,0.0002,0.0002,0.0003,0.0006,0.0013,0.0025,0.0043,0.0071
+0.0108,0.0160,0.0239,0.0345,0.0472,0.0618,0.0783,0.0958,0.1143,0.1339
+0.1539,0.1743,0.1957,0.2175,0.2394,0.2614,0.2839,0.3061,0.3280,0.3503
+0.3724,0.3944,0.4163,0.4382,0.4595,0.4808,0.5023,0.5232,0.5442,0.5652
+0.5861,0.6064,0.6270,0.6474,0.6676,0.6880,0.7077,0.7272,0.7464,0.7654
+0.7838,0.8016,0.8187,0.8350,0.8502,0.8644,0.8778,0.8901,0.9014,0.9119
+0.9214,0.9297,0.9374,0.9441,0.9504,0.9559,0.9608,0.9653,0.9693,0.9730
+0.9764,0.9793,0.9818,0.9841,0.9861,0.9877,0.9891,0.9903,0.9914,0.9923
+0.9930,0.9936,0.9942,0.9948,0.9952,0.9957,0.9961,0.9965,0.9968,0.9971
+0.9975,0.9977,0.9979,0.9981,0.9983,0.9985,0.9987,0.9988,0.9990,0.9991
+0.9992,0.9993,0.9994,0.9994,0.9995,0.9996,0.9996,0.9997,0.9997,0.9997
+0.9998,0.9998,0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999
+0.9999,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc1v1fract=0.0000,0.0000,0.0000,0.0000,0.0001,0.0001,0.0001,0.0001
+0.0002,0.0002,0.0003,0.0003,0.0004,0.0005,0.0006,0.0009,0.0014,0.0025
+0.0046,0.0087,0.0152,0.0250,0.0373,0.0519,0.0688,0.0872,0.1066,0.1271
+0.1481,0.1695,0.1911,0.2126,0.2344,0.2566,0.2792,0.3017,0.3238,0.3460
+0.3680,0.3900,0.4117,0.4334,0.4545,0.4761,0.4973,0.5184,0.5390,0.5597
+0.5801,0.6006,0.6215,0.6418,0.6621,0.6822,0.7018,0.7213,0.7403,0.7590
+0.7775,0.7958,0.8132,0.8303,0.8465,0.8616,0.8758,0.8885,0.9002,0.9109
+0.9205,0.9291,0.9368,0.9437,0.9498,0.9554,0.9602,0.9650,0.9691,0.9728
+0.9761,0.9791,0.9817,0.9839,0.9859,0.9875,0.9890,0.9901,0.9912,0.9921
+0.9929,0.9936,0.9942,0.9947,0.9952,0.9956,0.9960,0.9964,0.9967,0.9971
+0.9973,0.9976,0.9978,0.9981,0.9982,0.9984,0.9986,0.9988,0.9989,0.9990
+0.9991,0.9992,0.9993,0.9994,0.9995,0.9995,0.9996,0.9997,0.9997,0.9997
+0.9998,0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc1y2fract=0.0000,0.0000,0.0000,0.0001,0.0001,0.0001,0.0002,0.0002
+0.0003,0.0003,0.0004,0.0005,0.0007,0.0009,0.0012,0.0016,0.0026,0.0043
+0.0073,0.0121,0.0198,0.0299,0.0427,0.0578,0.0749,0.0937,0.1139,0.1347
+0.1561,0.1781,0.2003,0.2230,0.2455,0.2679,0.2901,0.3122,0.3344,0.3568
+0.3789,0.4010,0.4232,0.4450,0.4663,0.4873,0.5085,0.5298,0.5506,0.5714
+0.5921,0.6128,0.6332,0.6529,0.6729,0.6925,0.7119,0.7314,0.7501,0.7685
+0.7863,0.8038,0.8209,0.8371,0.8523,0.8667,0.8800,0.8925,0.9039,0.9140
+0.9230,0.9312,0.9388,0.9456,0.9516,0.9569,0.9617,0.9660,0.9697,0.9732
+0.9763,0.9789,0.9814,0.9835,0.9853,0.9870,0.9884,0.9895,0.9906,0.9915
+0.9923,0.9931,0.9938,0.9944,0.9949,0.9954,0.9959,0.9963,0.9966,0.9970
+0.9973,0.9975,0.9977,0.9980,0.9982,0.9984,0.9986,0.9987,0.9989,0.9990
+0.9992,0.9993,0.9994,0.9995,0.9995,0.9996,0.9997,0.9997,0.9998,0.9998
+0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc1x2fract=0.0000,0.0000,0.0001,0.0001,0.0001,0.0002,0.0002,0.0003
+0.0003,0.0004,0.0005,0.0006,0.0007,0.0009,0.0012,0.0018,0.0029,0.0048
+0.0076,0.0120,0.0183,0.0265,0.0369,0.0489,0.0625,0.0771,0.0930,0.1097
+0.1272,0.1458,0.1650,0.1846,0.2051,0.2256,0.2469,0.2688,0.2909,0.3128
+0.3351,0.3573,0.3798,0.4024,0.4247,0.4470,0.4691,0.4913,0.5136,0.5353
+0.5571,0.5792,0.6009,0.6223,0.6434,0.6643,0.6853,0.7060,0.7262,0.7460
+0.7651,0.7841,0.8026,0.8201,0.8363,0.8521,0.8665,0.8803,0.8930,0.9044
+0.9151,0.9247,0.9337,0.9419,0.9492,0.9558,0.9617,0.9671,0.9718,0.9761
+0.9797,0.9828,0.9854,0.9876,0.9895,0.9910,0.9923,0.9933,0.9941,0.9949
+0.9955,0.9959,0.9964,0.9967,0.9971,0.9973,0.9976,0.9978,0.9980,0.9983
+0.9984,0.9986,0.9987,0.9988,0.9990,0.9991,0.9992,0.9993,0.9993,0.9994
+0.9994,0.9995,0.9995,0.9996,0.9996,0.9997,0.9997,0.9997,0.9998,0.9998
+0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2x1fract=0.0000,0.0001,0.0002,0.0002,0.0003,0.0003,0.0004,0.0005
+0.0006,0.0007,0.0008,0.0010,0.0012,0.0014,0.0017,0.0022,0.0030,0.0044
+0.0069,0.0107,0.0160,0.0233,0.0321,0.0431,0.0554,0.0692,0.0839,0.1000
+0.1171,0.1348,0.1536,0.1728,0.1927,0.2131,0.2343,0.2559,0.2775,0.2991
+0.3213,0.3438,0.3664,0.3889,0.4113,0.4338,0.4566,0.4789,0.5017,0.5242
+0.5464,0.5681,0.5901,0.6114,0.6331,0.6551,0.6761,0.6973,0.7178,0.7382
+0.7584,0.7776,0.7963,0.8145,0.8318,0.8481,0.8634,0.8775,0.8905,0.9025
+0.9135,0.9234,0.9326,0.9407,0.9478,0.9544,0.9606,0.9659,0.9708,0.9750
+0.9787,0.9819,0.9846,0.9870,0.9890,0.9906,0.9919,0.9930,0.9939,0.9947
+0.9953,0.9959,0.9964,0.9968,0.9971,0.9975,0.9977,0.9979,0.9982,0.9983
+0.9985,0.9986,0.9987,0.9989,0.9990,0.9991,0.9992,0.9993,0.9993,0.9994
+0.9995,0.9995,0.9996,0.9996,0.9997,0.9997,0.9998,0.9998,0.9998,0.9998
+0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2y1fract=0.0000,0.0000,0.0001,0.0001,0.0001,0.0001,0.0001,0.0002
+0.0002,0.0003,0.0003,0.0004,0.0005,0.0006,0.0008,0.0011,0.0017,0.0027
+0.0047,0.0083,0.0142,0.0233,0.0356,0.0507,0.0682,0.0876,0.1083,0.1304
+0.1533,0.1769,0.2002,0.2235,0.2469,0.2702,0.2935,0.3167,0.3394,0.3618
+0.3847,0.4074,0.4295,0.4515,0.4733,0.4946,0.5158,0.5368,0.5576,0.5783
+0.5991,0.6201,0.6404,0.6610,0.6811,0.7011,0.7206,0.7397,0.7587,0.7775
+0.7958,0.8131,0.8303,0.8467,0.8624,0.8770,0.8904,0.9028,0.9140,0.9238
+0.9326,0.9404,0.9473,0.9535,0.9591,0.9639,0.9681,0.9719,0.9752,0.9782
+0.9808,0.9831,0.9850,0.9867,0.9881,0.9894,0.9905,0.9916,0.9924,0.9932
+0.9938,0.9944,0.9950,0.9955,0.9959,0.9962,0.9965,0.9969,0.9971,0.9974
+0.9976,0.9978,0.9980,0.9983,0.9984,0.9986,0.9987,0.9988,0.9990,0.9991
+0.9992,0.9993,0.9994,0.9995,0.9995,0.9996,0.9997,0.9997,0.9998,0.9998
+0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2u1fract=0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0001,0.0001
+0.0001,0.0001,0.0001,0.0002,0.0002,0.0002,0.0003,0.0004,0.0006,0.0012
+0.0024,0.0052,0.0107,0.0196,0.0320,0.0474,0.0652,0.0852,0.1066,0.1288
+0.1509,0.1737,0.1965,0.2199,0.2431,0.2663,0.2895,0.3123,0.3350,0.3577
+0.3800,0.4019,0.4240,0.4457,0.4668,0.4882,0.5093,0.5302,0.5510,0.5716
+0.5920,0.6121,0.6323,0.6525,0.6723,0.6919,0.7113,0.7306,0.7497,0.7685
+0.7869,0.8046,0.8222,0.8393,0.8555,0.8704,0.8844,0.8973,0.9087,0.9193
+0.9283,0.9363,0.9435,0.9498,0.9553,0.9602,0.9646,0.9686,0.9721,0.9752
+0.9779,0.9804,0.9827,0.9846,0.9863,0.9878,0.9891,0.9902,0.9912,0.9920
+0.9927,0.9933,0.9939,0.9944,0.9949,0.9952,0.9956,0.9959,0.9963,0.9966
+0.9969,0.9972,0.9974,0.9977,0.9979,0.9981,0.9983,0.9984,0.9986,0.9988
+0.9989,0.9990,0.9991,0.9993,0.9994,0.9995,0.9995,0.9996,0.9996,0.9997
+0.9997,0.9998,0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999
+0.9999,0.9999,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2v1fract=0.0000,0.0000,0.0000,0.0000,0.0001,0.0001,0.0001,0.0001
+0.0001,0.0002,0.0002,0.0003,0.0003,0.0004,0.0006,0.0008,0.0012,0.0020
+0.0038,0.0077,0.0148,0.0252,0.0386,0.0549,0.0738,0.0941,0.1152,0.1374
+0.1601,0.1827,0.2055,0.2288,0.2519,0.2750,0.2981,0.3207,0.3434,0.3662
+0.3881,0.4102,0.4315,0.4530,0.4746,0.4960,0.5171,0.5381,0.5585,0.5791
+0.5993,0.6196,0.6395,0.6594,0.6790,0.6989,0.7177,0.7370,0.7558,0.7744
+0.7926,0.8108,0.8278,0.8440,0.8595,0.8745,0.8879,0.9005,0.9115,0.9218
+0.9307,0.9385,0.9454,0.9514,0.9566,0.9614,0.9657,0.9695,0.9730,0.9760
+0.9789,0.9813,0.9834,0.9853,0.9869,0.9883,0.9895,0.9906,0.9915,0.9922
+0.9929,0.9936,0.9941,0.9947,0.9952,0.9955,0.9959,0.9963,0.9967,0.9970
+0.9973,0.9976,0.9978,0.9980,0.9982,0.9984,0.9986,0.9987,0.9989,0.9991
+0.9992,0.9993,0.9994,0.9995,0.9996,0.9996,0.9997,0.9997,0.9998,0.9998
+0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000,1.0000,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2y2fract=0.0000,0.0000,0.0000,0.0001,0.0001,0.0001,0.0002,0.0002
+0.0002,0.0003,0.0004,0.0005,0.0006,0.0007,0.0010,0.0014,0.0022,0.0035
+0.0058,0.0097,0.0163,0.0258,0.0382,0.0539,0.0715,0.0912,0.1122,0.1341
+0.1573,0.1803,0.2038,0.2275,0.2514,0.2749,0.2984,0.3215,0.3442,0.3667
+0.3893,0.4112,0.4332,0.4549,0.4767,0.4985,0.5191,0.5403,0.5612,0.5817
+0.6025,0.6231,0.6436,0.6642,0.6839,0.7036,0.7235,0.7427,0.7617,0.7803
+0.7984,0.8159,0.8328,0.8491,0.8645,0.8788,0.8919,0.9043,0.9152,0.9250
+0.9337,0.9414,0.9482,0.9541,0.9595,0.9641,0.9683,0.9721,0.9753,0.9782
+0.9807,0.9829,0.9849,0.9866,0.9881,0.9894,0.9905,0.9914,0.9922,0.9929
+0.9936,0.9942,0.9948,0.9952,0.9956,0.9960,0.9964,0.9967,0.9970,0.9972
+0.9975,0.9977,0.9979,0.9981,0.9983,0.9985,0.9987,0.9988,0.9989,0.9990
+0.9991,0.9993,0.9994,0.9994,0.9995,0.9996,0.9996,0.9997,0.9997,0.9998
+0.9998,0.9998,0.9998,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000
+hwc2x2fract=0.0000,0.0001,0.0001,0.0001,0.0002,0.0002,0.0003,0.0003
+0.0004,0.0005,0.0006,0.0007,0.0009,0.0012,0.0016,0.0023,0.0039,0.0064
+0.0100,0.0151,0.0223,0.0316,0.0430,0.0561,0.0707,0.0864,0.1032,0.1212
+0.1397,0.1592,0.1791,0.1996,0.2206,0.2422,0.2638,0.2856,0.3081,0.3308
+0.3534,0.3759,0.3990,0.4214,0.4438,0.4662,0.4885,0.5107,0.5324,0.5543
+0.5762,0.5979,0.6196,0.6411,0.6621,0.6828,0.7033,0.7238,0.7438,0.7632
+0.7824,0.8007,0.8182,0.8352,0.8514,0.8666,0.8806,0.8934,0.9052,0.9160
+0.9256,0.9346,0.9426,0.9498,0.9561,0.9617,0.9670,0.9715,0.9755,0.9790
+0.9820,0.9846,0.9869,0.9888,0.9904,0.9917,0.9929,0.9938,0.9946,0.9952
+0.9958,0.9962,0.9966,0.9970,0.9973,0.9976,0.9978,0.9981,0.9982,0.9984
+0.9986,0.9987,0.9989,0.9990,0.9991,0.9992,0.9993,0.9994,0.9994,0.9995
+0.9996,0.9996,0.9996,0.9997,0.9997,0.9997,0.9998,0.9998,0.9998,0.9998
+0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,1.0000
+1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000

PARAM/hmsflags.param

+; Offset in the spectrometer momentum.
+  hmomentum_factor = 0.000 ; leave 0 so it will have no effect
+                           ; use hpcentral_offset if one needs to offset central momentum
+;
+; The hdelta_offset,htheta_offset,hphi_offset effect the
+; reconstructed target quantities. Used in h_targ_trans.f 
+; For transport x is in the dispersive direction with + down
+;               y is in the horizontal + towards small angles.
+;               z is along the central optics axis.
+; In transport coordinates phi = hyptar = dy/dz and theta = hxptar = dx/dz 
+;    but for unknown reasons the yp offset is named  htheta_offset
+;    and  the xp offset is named  hphi_offset
+; Do not to change these values, since these are the zero order
+;  CMOP matrix elements. If you do change then your hms sieve
+;  plots will be screwed up.
+  hdelta_offset = 0.       ; (%)   hdelta_tar = hdelta_tar + hdelta_offset
+  htheta_offset = 4.83e-5  ; (rad) hyp_tar = hyp_tar + htheta_offset
+  hphi_offset   = -4.73e-3 ; (rad) hxp_tar = hxp_tar + hphi_offset
+;
+
+;saturation correction flag
+  genable_hms_satcorr = 2000 ;(0=disabled) - h_satcorr.f
+                          ; a correction to hsdelta event by event
+                          ; for a problem in setting Q3 current.
+                          ; There was an unknown zero offset in the Q3 current.
+                          ; The magnet setting code field00.f partially 
+                          ; fixes this problem. T. Horn in 2003 determined
+                          ; the corrections to delta.
+                          ; Data taken with fields set by field99.f or earlier should set to 1999.
+                          ; Data taken with fields set by field00.f or later should set to 2000.
+; central field  correction
+  genable_hms_fieldcorr = 1 ; (1=disabled) - h_fieldcorr.f 
+                            ; Need to enable for experiments before Jan 2002 .
+			    ; experiments using field02 and field03 should disable. 
+;
+; The following offsets are applied to the central kinematic variables
+;  in h_apply_offsets.f  . These might be modified by an experiment
+;  after doing calibration with elastic ep.
+; The values below are from T. Horn 2003 analysis
+  hpcentral_offset = -0.25   ; sets hpcentral = hpcentral * ( 1. + hpcentral_offset / 100. )
+		    ; experiments earlier than April 2003 
+                    ; should use  about -0.3 which is the best
+                    ; estimate based on several previous experiments.
+  hthetacentral_offset = 0.0002 ; (rad) 
+                         ;htheta_lab=htheta_lab + hthetacentral_offset/degree 
+; This offset is determined from elastic ep data.
+; Must be added to ssxptar when used in calculating lab angles.
+; Example is in h_physics.f
+;
+  h_oopcentral_offset = 0.0011  ; (rad)
+
+