Add parameters for DAQ07 data

examples/PARAM/daq07/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

examples/PARAM/daq07/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.13 ; 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.00 ; (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)
+
+

examples/PARAM/daq07/htracking.param

+; hms pattern recognition and tracking code parameters
+; hmax_pr_hits is the maximum number of decoded hits allowed in a chamber for
+;              pattern recognition to proceed. This protexts against unusual
+;              kami kaze events which hang the analyzer to a long time.
+;
+        hmax_pr_hits =    25, 25
+; maximum number of possible focal plane tracks, must be <= 20
+        hntracks_max_fp = 10
+; if h_remove_sppt_if_one_y_plane=1 then remove space point if space point
+;    contains only one y plane.
+; Recommend h_remove_sppt_if_one_y_plane=0.  
+        h_remove_sppt_if_one_y_plane=0
+; hminhit(i)   minimum number of hits in each space point (one for each chamber)
+	hmin_hit =     5, 5
+; hmin_combos(i) minimum number of pairs in each space point ( one per chamber)
+; should be 3/3 for 4/6 tracking, and 4/4 for 5/6 tracking (JRA)
+	hmin_combos = 4, 4 
+; hspace_point_criterion(i) minimum separation of distinct space points
+; recommended value for hspace_point_criterion ( based on meeting
+; Feb 2003) 
+       hspace_point_criterion = 1.0, 1.0
+
+        hstub_max_xpdiff = .05
+;
+; The following four criteria are applied to the stub fit results at the
+; focal plane position to link stubs to tracks
+; hxt_track_criterion  
+; hyt_track_criterion
+; hxpt_track_criterion
+; hypt_track_criterion
+; recommended value for hspace_point_criterion ( based on meeting
+; Feb 2003) 
+       hxt_track_criterion  = 10.
+       hyt_track_criterion  = 4.
+       hxpt_track_criterion = 0.2
+       hypt_track_criterion = 0.5
+;
+; hsingle_stub            If .ne. 0 make a track of each single stub
+  hsingle_stub = 0
+; turn on/off small angle approximation to left/right decision for y
+  hSmallAngleApprox = 1
+
+; Crude track selection criteria in h_physics
+; hsel_chi2_fpperdegmax   maximum chi2 per degree of freedom for track
+  hsel_chi2_fpperdegmax = 2000000.
+; hsel_dedx1min           minimum dedx to chose best track in scin 1
+  hsel_dedx1min = -100000
+; hsel_dedx1max           maximum dedx to chose best track in scin 1
+  hsel_dedx1max = 1000000.
+; hsel_betamin            minimum track beta to chose best track 
+  hsel_betamin  = -100000
+; hsel_betamax            maximum track beta to chose best track 
+  hsel_betamax  = 100000.
+; hsel_etmin              minimum track et to chose best track 
+  hsel_etmin = -100000.
+; hsel_etmax              maximum track et to chose best track 
+  hsel_etmax = 100000.
+; hsel_ndegreesmin        minimum number of degrees of freedom for track
+  hsel_ndegreesmin = 1
+
+
+; hstat_maxchisq          chisquared limit for tracks used to measure hodo. eff.
+  hstat_maxchisq = 10.
+; hsel_using_scin         uses scintillator for track selection
+  hsel_using_scin = 0
+; hstat_slop              distance from center of scin. to count as expecting hit
+  hstat_slop = 2.
+; hstat_mineff            warning level for scin. effic.
+  hstat_mineff = .95
+; hstat_cal_maxchisq      chisquared limit for tracks used to measure cal. eff.
+  hstat_cal_maxchisq = 10.
+; hstat_cal_slop          distance from center of block to count as expecting hit
+  hstat_cal_slop = 2.
+; hcer_min_eff            warning level for cerenkov effic.
+  hcer_min_eff = .98
+; hcer_threshold          threshold for good hit (for measureing effic.)
+  hcer_threshold = 0.5
+; hdc_min_eff             warning level for dc raw effic.(1 per plane)
+  hdc_min_eff = .96,.98,.96,.98,.98,.94
+                .96,.98,.98,.98,.98,.96
+; hdc_min_wire_eff        warning level for wire effic, measured using tracking
+  hdc_min_wire_eff = .85
+; hdc_min_plane_eff        warning level for plane effic, measured using tracking
+  hdc_min_plane_eff = .95,.95,.95,.95,.95,.95
+                      .95,.95,.95,.95,.95,.95
+;Derek added this for scintillator based fiducial cuts.  Everything from
+;loscin to hiscin inclusive will be included as "good" scintillator hits.
+;  hxloscin  = 5, 5
+;  hxhiscin  = 12, 12
+;  hyloscin  = 4, 4
+;  hyhiscin  = 7, 7
+  hxloscin  = 4, 4
+  hxhiscin  = 13, 13
+  hyloscin  = 4, 4
+  hyhiscin  = 7, 7
+;  htrack_eff_test_scin_planes is the number of planes nec needed to 
+;  set sweet spot to true. 4 is extra clean, 3 is good enough for e-'s.
+  htrack_eff_test_num_scin_planes = 4

examples/PARAM/daq07/online07.param

+#real raddeg
+raddeg=3.14159265/180
+
+;  This file contains the sos and hms parameters needed for all the analysis
+; 
+begin parm kinematics
+#include "PARAM/daq07/general.param"
+#include "PARAM/daq07/gscalers.param"
+#include "PARAM/daq07/gbeam.param"
+end parm kinematics
+
+begin parm hms_parameters
+#include "PARAM/daq07/hdeadwires.param"
+#include "PARAM/daq07/hdc.param"
+#include "PARAM/daq07/hdriftmap.param"
+#include "PARAM/daq07/hhodo.param"
+#include "PARAM/daq07/hcal.param"
+#include "PARAM/daq07/hcer.param"
+#include "PARAM/daq07/haero.param"
+end parm hms_parameters
+
+begin parm sos_parameters
+#include "PARAM/daq07/sdc.param"
+#include "PARAM/daq07/sdriftmap.param"
+#include "PARAM/daq07/shodo.param"
+#include "PARAM/daq07/scal.param"
+#include "PARAM/daq07/scer.param"
+end parm sos_parameters