-
Notifications
You must be signed in to change notification settings - Fork 4
/
prex.dat
72 lines (50 loc) · 1.92 KB
/
prex.dat
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
# (optional) Setup parameters for hamc for PREX
# Setup choices for HRS. Using prex2 collimator normally
hrs_setup warmseptum thrstrans useempiricalangle
#useprex2collim
#usepaulcollim
#usetrkcorr
#
#
# HRSTrans parameters. Here is tune P used in the 2019 PREX run
quad_fields q1 0.11505 q2 -0.136543 q3 -0.159619
# For references, here is standard HRS tune, this does not belong with septum
# quad_fields q1 0.1746 q2 -0.1385 q3 -0.1281
# turn off A_T hole for PREX-2
spreader_collim a 9 z 4 t 0.0001
# HRS momentum and angle (4.74 for PREX-II)
HRS_P0 0.95
HRS_angle 4.74
# If Angle_cut is turned on, there is a cut on scattering angle
# Ideally this should be commented out and let the collimator do the job.
# Angle_cut 3.0
HRS_arm left
# HRS_arm right
PREX_model horpb
# The HRS_acceptoff is normally commented out (or deleted)
# When uncommented it's a way of testing the acceptance model.
# Turn on one thing at a time to check the cuts at each point.
# 0 means ignore the acceptance limit, 1 means use it.
# 1 2 3 4 5 6 7 8 9 10
# target, collim septin, septout, q1, dipin, dipexit, q3in, q3exit, focal
# cut at collim
HRS_acceptoff 0 1 0 0 1 1 0 1 1 0
# HRS_acceptoff 0 0 0 0 0 0 0 0 0 0
# Beam parameters.
# Raster is in meters, typ 0.006 (full width, so +/- half this)
# normally raster 0.006 x 0.004
# current in uA.
beam E0 0.95 E0sigma 0.00011 rastered xraster 0.006 yraster 0.004 current 70
# range of angles (scatt, azimuth) in degrees
# phi=0 is up. phi=90 is into HRS, phi=180 is down.
out_angles thetamin 0.001 thetamax 20 phimin 0 phimax 180
radcor numer
# distance of detector to focal plane (meters)
detector dist 0
# uncomment the following line if you want to iterate
# iterate
# Here is the variable to iterate
# kick:track P0 0.003
# to turn off ntup (good for high stats).
# comment this out to receive an ntuple output.
# hamc_inout nontup