Skip to content
Snippets Groups Projects
Commit e647e89b authored by Vardan Tadevosyan's avatar Vardan Tadevosyan
Browse files

Added/modified comments to the HMS calorimeter header files

parent 2ab14e51
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1521 additions and 0 deletions
examples/PARAM/52949/gbeam.param 0 → 100644
+ 96
− 0
View file @ e647e89b
; 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
examples/PARAM/52949/gdebug.param 0 → 100644
+ 3
− 0
View file @ e647e89b
gmisc_min_peds = 500
gelossdebug = 0 ;eloss debug (0=off,1=on)
examples/PARAM/52949/general.param 0 → 100644
+ 40
− 0
View file @ e647e89b
#real raddeg
raddeg=3.14159265/180
part_mass=0.000511
; hms/sosflags.param include spectrometer offsets and options.
#include "PARAM/52949/hcana.param"
#include "PARAM/52949/genflags.param"
#include "PARAM/52949/hmsflags.param"
#include "PARAM/52949/sosflags.param"
#include "PARAM/52949/gdebug.param"
#include "PARAM/52949/hdebug.param"
#include "PARAM/52949/sdebug.param"
#include "PARAM/52949/htracking.param"
#include "PARAM/52949/stracking.param"
#include "PARAM/52949/gscalers.param"
#include "PARAM/52949/gbeam.param"
#include "PARAM/52949/gtarget.param"
#include "PARAM/52949/hdc_offsets.param"
#include "PARAM/52949/hdc.pos"
#include "PARAM/52949/hhodo.pos"
#include "PARAM/52949/hcal.pos"
#include "PARAM/52949/sdc_offsets.param"
#include "PARAM/52949/sdc.pos"
#include "PARAM/52949/shodo.pos"
#include "PARAM/52949/scal.pos"
#include "PARAM/52949/hdeadwires.param"
#include "PARAM/52949/hdc.param.vt.52927"
#include "PARAM/52949/hms.driftmap.vt.52949"
#include "PARAM/52949/hhodo.param.vt"
#include "PARAM/52949/hcal.param.vt.52949"
#include "PARAM/52949/hcer.param.vt"
#include "PARAM/52949/haero.param"
#include "PARAM/52949/sdc.param.vt.53000"
#include "PARAM/52949/sos.driftmap.vt.52930"
#include "PARAM/52949/shodo.param.vt"
#include "PARAM/52949/scal.param.vt.52843"
#include "PARAM/52949/scer.param.vt"
examples/PARAM/52949/genflags.param 0 → 100644
+ 12
− 0
View file @ e647e89b
;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
examples/PARAM/52949/gscalers.param 0 → 100644
+ 30
− 0
View file @ e647e89b
; These BCM's (1 and 2) numbers are from Jan 05 experiments (VT)
;point to scaler indices for BCM's
gbcm1_index = 321
gbcm2_index = 322
gbcm3_index = 323
gunser_index = 324
gbcm1_gain = 0.000328449 ; microA/Hz (New Value)
gbcm2_gain = 0.000381451 ; microA/Hz (New Value)
gbcm3_gain = 0.00043343 ; microA/Hz (Old, Value)
gbcm1_offset = 250403. ; Hz (New Value)
gbcm2_offset = 250123. ; Hz (New Value)
gbcm3_offset = 245437. ; Hz (Old, Value)
gunser_offset = 514529. ; Hz (Old, Value)
gunser_gain = 0.00025001 ; microA/Hz (New Value)
;define the threshold current at which beam is "on". (in microamps)
g_beam_on_thresh_cur = 1.5, 1.5
bcm_for_threshold_cut = 1
;point to hardware clock (1MHz scaler, scaler channel 325)
gclock_index = 325
gclock_rate = 1.e+6
examples/PARAM/52949/gtarget.param 0 → 100644
+ 170
− 0
View file @ e647e89b
; 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
;
examples/PARAM/52949/haero.param 0 → 100644
+ 12
− 0
View file @ e647e89b
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.,
examples/PARAM/52949/hcal.param.vt.52949 0 → 100644
+ 57
− 0
View file @ e647e89b
; 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 52949, 435527 events processed
hcal_pos_gain_cor= 0.425, 0.355, 0.416, 0.385, 0.247, 0.512, 0.552, 0.378, 0.502, 0.345, 0.379, 0.390, 0.000,
0.375, 0.428, 0.340, 0.387, 0.449, 0.485, 0.225, 0.445, 0.303, 0.269, 0.332, 0.653, 0.558,
0.609, 0.630, 0.605, 0.857, 0.631, 0.658, 0.466, 0.547, 0.714, 0.733, 0.523, 0.636, 0.000,
0.802, 0.600, 0.573, 0.720, 0.724, 0.822, 0.851, 0.744, 0.832, 0.614, 0.699, 0.738, 0.000,
hcal_neg_gain_cor= 0.422, 0.269, 0.288, 0.367, 0.353, 0.050, 0.184, 0.449, 0.235, 0.334, 0.203, 0.170, 0.000,
0.401, 0.353, 0.389, 0.511, 0.332, 0.337, 0.454, 0.304, 0.470, 0.363, 0.338, 0.371, 1.156,
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,
examples/PARAM/52949/hcal.pos 0 → 100644
+ 81
− 0
View file @ e647e89b
; 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)
examples/PARAM/52949/hcana.param 0 → 100644
+ 17
− 0
View file @ e647e89b
;
; 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
hcal_num_layers = 4
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"
hcal_layer_names = "1pr 2ta 3ta 4ta"
hhodo_plane_names = "1x 1y 2x 2y"
examples/PARAM/52949/hcer.param.vt 0 → 100644
+ 24
− 0
View file @ e647e89b
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
examples/PARAM/52949/hdc.param.vt.52927 0 → 100644
+ 44
− 0
View file @ e647e89b
;---------------------------------------------------------------------
; 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
; hms dc tdc minimum tdc value array allowed for a good hit
hdc_tdc_min_win = 2750, 2750, 2750, 2750, 2750, 2750
2750, 2750, 2750, 2750, 2750, 2750
; hms dc tdc maximum tdc value array allowed for a good hit
hdc_tdc_max_win = 3500, 3500, 3500, 3500, 3500, 3500
3500, 3500, 3500, 3500, 3500, 3500
; hms drift chamber tdc's time per channel
hdc_tdc_time_per_channel = 0.50
; hms zero time for drift chambers !DECREASING this number moves the hdtime plots to LOWER time.
hdc_plane_time_zero = (1670+14+18-2-0)
(1670+12+20-2-0)
(1670+13+18-0-2)
(1670+13+20-2-0)
(1670+12+20-4-0)
(1670+14+16-0+2-2)
(1670+15+16+2-4)
(1670+11.5+16)
(1670+13+12+2-2)
(1670+13+12+2-0)
(1670+10.5+16-0)
(1670+13+18-0-2)
; 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
examples/PARAM/52949/hdc.pos 0 → 100644
+ 183
− 0
View file @ e647e89b
; 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
examples/PARAM/52949/hdc_offsets.param 0 → 100644
+ 125
− 0
View file @ e647e89b
; 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
examples/PARAM/52949/hdeadwires.param 0 → 100644
+ 10
− 0
View file @ e647e89b
; 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)?
examples/PARAM/52949/hdebug.param 0 → 100644
+ 117
− 0
View file @ e647e89b
; 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
examples/PARAM/52949/hhodo.param.vt 0 → 100644
+ 185
− 0
View file @ e647e89b
; 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
new variable for picking good hits for tof fitting
; this should not be set tight until you are ready to fit
; tof and you figured out good values
htof_tolerance = 100.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
; use htofusinginvadc=1 if want invadc_offset
; invadc_linear, and invadc_adc to be used
htofusinginvadc=0
examples/PARAM/52949/hhodo.pos 0 → 100644
+ 88
− 0
View file @ e647e89b
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/52949/hms.driftmap.vt.52949 0 → 100644
+ 175
− 0
View file @ e647e89b
; 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.0002,0.0003,0.0003,0.0004
0.0005,0.0006,0.0007,0.0008,0.0009,0.0011,0.0013,0.0015,0.0019,0.0027
0.0038,0.0056,0.0087,0.0138,0.0212,0.0309,0.0423,0.0554,0.0697,0.0853
0.1022,0.1201,0.1384,0.1581,0.1781,0.1987,0.2200,0.2417,0.2637,0.2863
0.3088,0.3318,0.3544,0.3769,0.3993,0.4219,0.4444,0.4669,0.4888,0.5107
0.5325,0.5543,0.5757,0.5972,0.6182,0.6393,0.6595,0.6796,0.6995,0.7191
0.7382,0.7569,0.7750,0.7929,0.8099,0.8266,0.8422,0.8571,0.8706,0.8831
0.8946,0.9055,0.9156,0.9247,0.9332,0.9410,0.9480,0.9546,0.9606,0.9662
0.9709,0.9754,0.9793,0.9826,0.9854,0.9879,0.9897,0.9914,0.9928,0.9939
0.9947,0.9953,0.9959,0.9965,0.9969,0.9972,0.9975,0.9977,0.9980,0.9982
0.9984,0.9985,0.9986,0.9988,0.9988,0.9990,0.9991,0.9992,0.9993,0.9994
0.9994,0.9995,0.9995,0.9996,0.9996,0.9996,0.9997,0.9997,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,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.0001,0.0002,0.0002,0.0002,0.0003
0.0003,0.0004,0.0004,0.0005,0.0006,0.0008,0.0010,0.0012,0.0017,0.0025
0.0038,0.0061,0.0104,0.0175,0.0276,0.0408,0.0571,0.0754,0.0952,0.1165
0.1384,0.1614,0.1845,0.2080,0.2315,0.2549,0.2786,0.3018,0.3247,0.3479
0.3706,0.3928,0.4148,0.4366,0.4583,0.4795,0.5006,0.5214,0.5419,0.5623
0.5825,0.6026,0.6225,0.6423,0.6619,0.6818,0.7011,0.7200,0.7388,0.7571
0.7751,0.7926,0.8096,0.8258,0.8414,0.8561,0.8701,0.8828,0.8946,0.9054
0.9150,0.9238,0.9320,0.9390,0.9455,0.9513,0.9564,0.9612,0.9654,0.9692
0.9728,0.9760,0.9789,0.9815,0.9838,0.9856,0.9874,0.9888,0.9901,0.9912
0.9923,0.9932,0.9939,0.9946,0.9951,0.9956,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.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
hwc1u1fract=0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.0001,0.0001
0.0001,0.0001,0.0001,0.0001,0.0002,0.0002,0.0003,0.0004,0.0007,0.0013
0.0026,0.0055,0.0108,0.0194,0.0312,0.0457,0.0623,0.0806,0.0999,0.1202
0.1415,0.1638,0.1860,0.2089,0.2316,0.2545,0.2774,0.3004,0.3229,0.3456
0.3680,0.3908,0.4128,0.4349,0.4565,0.4777,0.4988,0.5197,0.5403,0.5609
0.5812,0.6012,0.6212,0.6410,0.6605,0.6798,0.6987,0.7180,0.7364,0.7547
0.7726,0.7902,0.8071,0.8233,0.8388,0.8536,0.8672,0.8799,0.8915,0.9022
0.9115,0.9202,0.9278,0.9348,0.9411,0.9470,0.9526,0.9575,0.9621,0.9665
0.9705,0.9742,0.9775,0.9804,0.9829,0.9850,0.9868,0.9884,0.9898,0.9910
0.9920,0.9928,0.9936,0.9943,0.9949,0.9954,0.9959,0.9963,0.9966,0.9970
0.9973,0.9976,0.9978,0.9980,0.9982,0.9984,0.9985,0.9987,0.9988,0.9989
0.9990,0.9991,0.9992,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,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.0001,0.0002,0.0002,0.0002,0.0003,0.0003,0.0004,0.0006,0.0009,0.0016
0.0031,0.0058,0.0107,0.0192,0.0308,0.0454,0.0621,0.0808,0.1007,0.1217
0.1434,0.1654,0.1878,0.2105,0.2331,0.2558,0.2784,0.3009,0.3237,0.3462
0.3682,0.3905,0.4120,0.4341,0.4555,0.4764,0.4976,0.5184,0.5391,0.5596
0.5801,0.6000,0.6198,0.6397,0.6591,0.6785,0.6980,0.7173,0.7362,0.7548
0.7729,0.7907,0.8078,0.8243,0.8402,0.8552,0.8691,0.8820,0.8936,0.9040
0.9135,0.9220,0.9297,0.9367,0.9431,0.9489,0.9541,0.9591,0.9635,0.9677
0.9716,0.9750,0.9783,0.9810,0.9834,0.9855,0.9874,0.9890,0.9903,0.9915
0.9926,0.9934,0.9941,0.9947,0.9953,0.9957,0.9962,0.9965,0.9969,0.9972
0.9974,0.9977,0.9979,0.9981,0.9983,0.9984,0.9986,0.9987,0.9988,0.9990
0.9991,0.9992,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,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.0000,0.0001,0.0001,0.0001,0.0001
0.0001,0.0001,0.0002,0.0002,0.0002,0.0003,0.0004,0.0007,0.0011,0.0021
0.0040,0.0074,0.0132,0.0216,0.0329,0.0468,0.0632,0.0815,0.1011,0.1225
0.1446,0.1673,0.1904,0.2143,0.2375,0.2607,0.2840,0.3071,0.3298,0.3528
0.3752,0.3971,0.4190,0.4410,0.4628,0.4846,0.5053,0.5261,0.5466,0.5673
0.5872,0.6074,0.6272,0.6469,0.6662,0.6854,0.7043,0.7231,0.7413,0.7598
0.7773,0.7949,0.8118,0.8279,0.8432,0.8575,0.8711,0.8837,0.8951,0.9053
0.9147,0.9231,0.9307,0.9376,0.9438,0.9495,0.9547,0.9596,0.9639,0.9681
0.9717,0.9750,0.9779,0.9805,0.9829,0.9848,0.9865,0.9880,0.9893,0.9905
0.9916,0.9925,0.9933,0.9940,0.9947,0.9952,0.9957,0.9962,0.9965,0.9969
0.9973,0.9975,0.9978,0.9980,0.9982,0.9984,0.9985,0.9987,0.9989,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.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
hwc1x2fract=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.0008,0.0011,0.0017,0.0028
0.0047,0.0080,0.0133,0.0208,0.0308,0.0422,0.0554,0.0700,0.0858,0.1025
0.1206,0.1393,0.1592,0.1795,0.2005,0.2220,0.2439,0.2661,0.2886,0.3108
0.3336,0.3564,0.3790,0.4013,0.4233,0.4457,0.4682,0.4904,0.5123,0.5342
0.5555,0.5772,0.5985,0.6195,0.6404,0.6609,0.6814,0.7012,0.7212,0.7401
0.7589,0.7771,0.7948,0.8116,0.8279,0.8437,0.8581,0.8715,0.8840,0.8957
0.9064,0.9163,0.9252,0.9335,0.9413,0.9483,0.9549,0.9607,0.9661,0.9709
0.9752,0.9792,0.9826,0.9855,0.9878,0.9897,0.9913,0.9927,0.9937,0.9946
0.9954,0.9960,0.9965,0.9969,0.9972,0.9975,0.9977,0.9979,0.9981,0.9983
0.9985,0.9986,0.9987,0.9988,0.9989,0.9990,0.9991,0.9992,0.9993,0.9993
0.9994,0.9995,0.9995,0.9996,0.9996,0.9996,0.9997,0.9997,0.9997,0.9998
0.9998,0.9998,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
hwc2x1fract=0.0001,0.0001,0.0002,0.0002,0.0003,0.0004,0.0005,0.0006
0.0008,0.0010,0.0012,0.0014,0.0016,0.0021,0.0028,0.0041,0.0069,0.0119
0.0191,0.0288,0.0403,0.0533,0.0678,0.0835,0.0999,0.1171,0.1352,0.1539
0.1733,0.1935,0.2140,0.2350,0.2565,0.2785,0.3007,0.3228,0.3450,0.3673
0.3896,0.4125,0.4350,0.4576,0.4801,0.5030,0.5252,0.5475,0.5698,0.5914
0.6130,0.6344,0.6553,0.6764,0.6972,0.7179,0.7381,0.7575,0.7769,0.7958
0.8139,0.8310,0.8474,0.8625,0.8765,0.8892,0.9008,0.9113,0.9205,0.9292
0.9368,0.9440,0.9505,0.9565,0.9622,0.9674,0.9721,0.9763,0.9799,0.9832
0.9860,0.9884,0.9904,0.9918,0.9932,0.9942,0.9950,0.9956,0.9962,0.9966
0.9970,0.9974,0.9977,0.9979,0.9981,0.9983,0.9985,0.9986,0.9988,0.9989
0.9990,0.9991,0.9992,0.9993,0.9994,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.9998,0.9999
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
hwc2y1fract=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.0011
0.0021,0.0042,0.0082,0.0153,0.0258,0.0396,0.0562,0.0754,0.0967,0.1193
0.1428,0.1670,0.1909,0.2151,0.2392,0.2632,0.2869,0.3104,0.3344,0.3577
0.3805,0.4031,0.4250,0.4469,0.4682,0.4896,0.5107,0.5315,0.5522,0.5724
0.5926,0.6129,0.6331,0.6529,0.6726,0.6922,0.7115,0.7305,0.7490,0.7675
0.7853,0.8027,0.8195,0.8354,0.8507,0.8654,0.8792,0.8918,0.9034,0.9138
0.9231,0.9315,0.9388,0.9451,0.9510,0.9562,0.9609,0.9651,0.9689,0.9723
0.9754,0.9783,0.9807,0.9828,0.9847,0.9864,0.9878,0.9891,0.9903,0.9914
0.9923,0.9931,0.9938,0.9944,0.9949,0.9954,0.9959,0.9963,0.9966,0.9969
0.9973,0.9975,0.9977,0.9979,0.9981,0.9983,0.9985,0.9986,0.9988,0.9989
0.9990,0.9991,0.9992,0.9993,0.9994,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
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.0003,0.0004,0.0005,0.0009,0.0019
0.0042,0.0094,0.0188,0.0315,0.0476,0.0664,0.0872,0.1095,0.1325,0.1561
0.1801,0.2041,0.2281,0.2519,0.2757,0.2990,0.3222,0.3453,0.3678,0.3899
0.4119,0.4341,0.4556,0.4766,0.4975,0.5186,0.5393,0.5597,0.5796,0.5994
0.6193,0.6391,0.6585,0.6778,0.6971,0.7159,0.7345,0.7528,0.7709,0.7888
0.8061,0.8230,0.8390,0.8546,0.8690,0.8827,0.8951,0.9061,0.9161,0.9249
0.9328,0.9396,0.9458,0.9514,0.9562,0.9608,0.9649,0.9687,0.9721,0.9753
0.9782,0.9806,0.9829,0.9848,0.9866,0.9882,0.9894,0.9905,0.9916,0.9924
0.9931,0.9938,0.9943,0.9949,0.9954,0.9958,0.9961,0.9965,0.9968,0.9971
0.9973,0.9976,0.9978,0.9980,0.9982,0.9984,0.9986,0.9987,0.9989,0.9990
0.9991,0.9992,0.9993,0.9993,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,1.0000,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.0000,0.0001,0.0001,0.0001
0.0001,0.0002,0.0002,0.0002,0.0003,0.0004,0.0005,0.0007,0.0010,0.0017
0.0032,0.0065,0.0131,0.0235,0.0372,0.0542,0.0733,0.0947,0.1164,0.1394
0.1622,0.1857,0.2091,0.2328,0.2560,0.2793,0.3025,0.3253,0.3477,0.3703
0.3922,0.4141,0.4355,0.4566,0.4778,0.4986,0.5189,0.5396,0.5598,0.5801
0.6002,0.6199,0.6393,0.6589,0.6780,0.6969,0.7159,0.7343,0.7529,0.7711
0.7888,0.8059,0.8226,0.8388,0.8544,0.8687,0.8820,0.8941,0.9051,0.9147
0.9233,0.9310,0.9377,0.9438,0.9493,0.9543,0.9590,0.9631,0.9669,0.9704
0.9737,0.9767,0.9793,0.9818,0.9838,0.9856,0.9872,0.9886,0.9898,0.9909
0.9919,0.9927,0.9934,0.9940,0.9945,0.9951,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.9990,0.9991,0.9992,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
hwc2y2fract=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.0005,0.0007,0.0013
0.0023,0.0045,0.0088,0.0163,0.0271,0.0410,0.0578,0.0769,0.0977,0.1203
0.1436,0.1676,0.1916,0.2161,0.2403,0.2638,0.2876,0.3108,0.3342,0.3572
0.3804,0.4032,0.4254,0.4470,0.4686,0.4898,0.5111,0.5318,0.5525,0.5731
0.5936,0.6138,0.6339,0.6536,0.6735,0.6929,0.7122,0.7309,0.7495,0.7679
0.7860,0.8037,0.8205,0.8368,0.8524,0.8667,0.8803,0.8929,0.9044,0.9146
0.9240,0.9320,0.9394,0.9458,0.9516,0.9569,0.9614,0.9656,0.9694,0.9728
0.9759,0.9787,0.9812,0.9833,0.9850,0.9867,0.9882,0.9894,0.9906,0.9916
0.9925,0.9932,0.9939,0.9945,0.9950,0.9955,0.9959,0.9964,0.9967,0.9970
0.9973,0.9976,0.9978,0.9979,0.9981,0.9983,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.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
hwc2x2fract=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.0008,0.0010,0.0012,0.0017,0.0027
0.0042,0.0071,0.0121,0.0198,0.0299,0.0422,0.0564,0.0724,0.0892,0.1078
0.1270,0.1471,0.1677,0.1893,0.2112,0.2336,0.2564,0.2794,0.3021,0.3253
0.3481,0.3711,0.3937,0.4163,0.4389,0.4611,0.4829,0.5048,0.5266,0.5482
0.5697,0.5904,0.6112,0.6321,0.6523,0.6727,0.6927,0.7125,0.7317,0.7508
0.7695,0.7877,0.8050,0.8222,0.8385,0.8539,0.8682,0.8814,0.8936,0.9047
0.9147,0.9239,0.9322,0.9399,0.9468,0.9531,0.9588,0.9640,0.9690,0.9732
0.9770,0.9803,0.9832,0.9858,0.9879,0.9897,0.9910,0.9923,0.9934,0.9942
0.9949,0.9956,0.9961,0.9965,0.9969,0.9972,0.9975,0.9978,0.9980,0.9982
0.9984,0.9985,0.9987,0.9988,0.9989,0.9990,0.9991,0.9992,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.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
examples/PARAM/52949/hmsflags.param 0 → 100644
+ 52
− 0
View file @ e647e89b
; 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)
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment