forked from oslocyclotronlab/usersort
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Zn68_sort.txt
145 lines (125 loc) · 5.45 KB
/
Zn68_sort.txt
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
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
# where to load the gainshifts from; one can also use:
# parameter gain_e = ...
# and so on
#gain file gainshifts_plain.dat
gain file gainshifts_si28_oscarcal.dat
# (constant fraction) cactus energy correction parameters for cactus time
parameter tnai_corr_enai = 0 1 0 0
parameter tnai_corr_enai = -0.838857 2192.67 230.828 -8.44176e-05
# (constant fraction) siri energy correction parameters for cactus time
parameter tnai_corr_esi = 0 1 0 0
# parameter tnai_corr_esi = 50.4471 -475514 4506.25 -0.00185164
################ Only when PPACs are used ##################################
# (constant fraction) engery/time correction parameters for PPAC
# The correction for the PPAC coutners are necessary due to SiRi, so
# usually the same corrections parmeters as in tnai_corr_esi should be used
# ### comment out this line(s) if no PPACs are used
#parameter tppac_corr_esi = 0 1 0 0
# parameter tppac_corr_esi = 50.4471 -475514 4506.25 -0.00185164
############################################################################
# Parameters for the range curve, 3He or alpha
#rangefile zrange_he3.dat
#rangefile zrange_a.dat
#rangefile zrange_d.dat
rangefile zrange_p.dat
# Cut of low-energy events by making a rectangle which is excluded
# in the down, left corner of the banana.
# Contains E-minimum 1, DE-minimum 1, E-minimum 2, DE-minimum 2.
# Including Z=1 particles:
parameter ede_rect = 500 250 30 500
# Excluding Z=1 particles:
#parameter ede_rect = 500 2500 2000 2500
# Thickness gate for proton banana
# parameter thick_range = 82 12 0
parameter thick_range = 130 20 0
# parameter thick_range = 140 20 0
# parameter thick_range = 276 75 0
# parameter thick_range = 230 25 0
# Thickness gate for 3He banana
#parameter thick_range = 48 8 0
# Thickness gate for alpha banana
#parameter thick_range = 165 10 0
#parameter thick_range = 48 5 0
# Fit of rkinz/qkinz column 1 as function of qkinz/rkinz column 6
# 240Pu(d,p) -- REALLY CHANGE THESE TO YOUR REACTION!!! WRITE IN KeV
parameter ex_from_ede = 15.169624e3 -0.992924 -0.001135e-3 \
1.5521e+4 -9.6197e-1 -1.4187e-6 \
1.5514e+4 -9.6096e-1 -1.4466e-6 \
1.5505e+4 -9.5981e-1 -1.4785e-6 \
1.5495e+4 -9.5852e-1 -1.5148e-6 \
1.5484e+4 -9.5707e-1 -1.5558e-6 \
1.5471e+4 -9.5545e-1 -1.6023e-6 \
1.5455e+4 -9.5363e-1 -1.6550e-6
# empirical excitation energy correction for the above, e.g. from known peaks
parameter ex_corr_exp = 0 1 \
0 1 \
0 1 \
0 1 \
0 1 \
0 1 \
0 1 \
0 1 \
# NaI time gates, for making the ALFNA matrices
# found by looking at the m_nai_t plot and projecting it
# down to the x-axis/m_nai_e_t_c x-axis plot
# remember to subtract the same amount of channels
# as the ones that are filled, and make sure you "hit" a top!
# parameter = lower cut of prompt , higher gate of prompt peak, \
# lower cut on backgr., higher cut backgr.
parameter nai_time_cuts = 193 205 \
144 156
############## Attention when no PPACs are used ###########################
# define id of NaI detectors and PPAC detectors (need to change usersort for more then four PPACs!)
# All other channel are considered as NaI detectors. In case you don't use fission detectors,
# define this parameters as "-1" - all id's are positive numbers, thus none will match -1.
# ### Corresponding line /without/ fission detectors:
# parameter channel_PPAC = -1 -1 -1 -1
# ### Corresponding line /with/ fission detectors: [adopt channel numbers]
parameter channel_PPAC = 4 12 30 31
################ Only when PPACs are used ##################################
# PPAC time gates
# time gates of the ppacs are found from the fission blob
# in the m_siri_e_t_ppac plot (fission blob in tPPAC vs E_SiRi gate)
# - time gates from the y-axis, and energy from the x-axis
# parameter = lower cut of prompt , higher gate of prompt peak, \
# lower cut on backgr., higher cut backgr.
# ### comment out this line(s) if no PPACs are used
# parameter ppac_time_cuts = 190 200 \
# 290 300
# Energy gate on fission
# the energy gate is found by looking at the
# m_siri_e_t_ppac plot (fission blob in tPPAC vs E_SiRi gate)
# - it should be roughly the fission barrier
# (Fission starts at some point, only events with
# a higher energy than excitation_energy_min can be true fission events)
# ### comment out this line(s) if no PPACs are used
# parameter fission_excitation_energy_min = 5000 # in keV
# Total efficiency of the PPACs in 4Pi
# parameter ppac_efficiency = 0.48
############################################################################
#data directory Datafiles
# data directory /path/to/data_dir
data directory /Users/frida/usersort
#maximum number of buffers to read for each file; for testing
#max_buffers 1
# Zn68 is 4,5,6,7
# list all data files
#data file sirius-20180123-102129.data
#data file sirius-20180123-160159.data
#data file sirius-20180123-115357.data
data file sirius-20180123-115357-big-000.data
data file sirius-20180123-112829.data
data file sirius-20180123-173618.data
data file sirius-20180123-173618-big-000.data
#data file sirius-20180124-164441-big-000.data
#data file sirius-20180124-164441.data
#data file sirius-20180124-224822.data
#data file sirius-20180125-091200.data #Si
# Write to MAMA matrix
#export mama h_ex_f0 h_ex_f0_MAMA
#export mama m_alfna_bg alfna_bg
export mama m_alfna alfna
# at the end, spectra are dumped to the root file
# root output file
export root offline_Zn68_pp_calsiri.root
export mama m_alfna alfna.m