Skip to content
Snippets Groups Projects
Commit a0a8b193 authored by Eric Pooser's avatar Eric Pooser
Browse files

Commit the revised directory structure

parent 60588185
No related branches found
No related tags found
No related merge requests found
Showing
with 1514 additions and 0 deletions
;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
; 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
;
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.,
; 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,
; 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)
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
;---------------------------------------------------------------------
; HMS_TRACKING
; CTP parameter file containing all tracking parameters for the HMS
;----------------------------------------------------------------------
; sigma of wire chamber resolution for each plane
hdc_sigma = 0.020
0.020
0.020
0.020
0.020
0.020
0.020
0.020
0.020
0.020
0.020
0.020
hdc_tdc_min_win = 2000,2000,2000,2000,2000,2000
2000,2000,2000,2000,2000,2000
hdc_tdc_max_win = 5500,5000,5500,5500,5000,5500
5500,5000,5500,5500,5000,5500
; hms drift chamber tdc's time per channel
hdc_tdc_time_per_channel = -0.10
; hms zero time for drift chambers !DECREASING this number moves the hdtime plots to LOWER time.
hdc_plane_time_zero = -213.0,-213.0,-213.0,-213.0,-213.0,-213.0
-213.0,-213.0,-213.0,-213.0,-213.0,-213.0
; Dave Abbott's wire velocity correction
hdc_wire_velocity = 12.0
hdc_central_time = 7,9,3,4,6,5
7,5,3,4,6,6
; 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
; 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
; Custom parameter file which should loaded when aiming to
; analyze HMS drift chamber data with no tracking.
; Increase HMS TDC window size
hdc_tdc_min_win = -50000,-50000,-50000,-50000,-50000,-50000
-50000,-50000,-50000,-50000,-50000,-50000
hdc_tdc_max_win = 50000,50000,50000,50000,50000,50000
50000,50000,50000,50000,50000,50000
; Remove hodoscope dependence on tracking
hsel_using_scin = 0
; hsel_using_prune = 0
; Adjust HMS DC plane t = 0
; hdc_plane_time_zero = 1500.0,1500.0,1500.0,1500.0,1500.0,1500.0
; 1500.0,1500.0,1500.0,1500.0,1500.0,1500.0
; 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)?
; Lookup table
;number of bins in Meek's time to distance lookup table
hdriftbins=138
;number of 1st bin in Meek's table in ns
hdrift1stbin=-24
;bin size in ns of Meek's table
hdriftbinsz=2
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
;
; Parameters that were built into Fortran analyzer that we want
; to pass as parameters so that the resulting code can be more generic.
;
; Scintillator parameters
href_npaddles = 0
href_nsperchan = 0.1
href_offset = 0.0
;
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
# 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
; 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
; 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)
; 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 = 35, 35
; 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 = 0.05
; 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 = 1
; 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 = 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 = 3
This diff is collapsed.
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
;---------------------------------------------------------------------
; SHMS_TRACKING
; CTP parameter file containing all tracking parameters for the HMS
;----------------------------------------------------------------------
; sigma of wire chabmer resolution for each plane
pdc_sigma = 0.02, 0.02, 0.02, 0.02, 0.02, 0.02
0.02, 0.02, 0.02, 0.02, 0.02, 0.02
; tdc window, wide limits for now
pdc_tdc_min_win = -55000,-55000,-55000,-55000,-55000,-55000
-55000,-55000,-55000,-55000,-55000,-55000
pdc_tdc_max_win = 55000,55000,55000,55000,55000,55000
55000,55000,55000,55000,55000,55000
; conversion from tdc lciks to ns
pdc_tdc_time_per_channel = -0.10 ; CAEN1190 TDC step size = 100ps
; zero time, decreasing moves sdtime plots to lower time
pdc_plane_time_zero = 120.0,120.0,120.0,120.0,120.0,120.0
120.0,120.0,120.0,120.0,120.0,120.0
; wire velocity
pdc_wire_velocity = 13.0
pdc_central_time = 1,1,1,1,1,1,1,1,1,1,1,1
ppartmass = 0.139
phodo_num_planes = 1
; Scintillator parameters
pref_npaddles = 6
pref_nsperchan = 0.1
pref_offset = 0.0
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment