1
|
Mar 23/23 PionLT/KaonLT Analysis Meeting Notes
|
2
|
----------------------------------------------
|
3
|
(Notes by GH)
|
4
|
|
5
|
Today: PionLT will be discussed first
|
6
|
|
7
|
Please remember to post your slides at:
|
8
|
https://redmine.jlab.org/projects/kltexp/wiki/Kaon_LT_Meetings
|
9
|
|
10
|
Present
|
11
|
-------
|
12
|
Regina - Garth Huber, Ali Usman, Alicia Postuma, Muhammad Junaid, Nathan
|
13
|
Heinrich, Vijay Kumar
|
14
|
ODU - Carlos Yero
|
15
|
JLab - Dave Gaskell
|
16
|
CUA - Richard Trotta, Tanja Horn
|
17
|
Ohio - Julie Roche
|
18
|
|
19
|
Carlos Discussion
|
20
|
-----------------
|
21
|
Carlos has agreed to attend and present some of his Lumi scaler analysis
|
22
|
results, so he goes first
|
23
|
|
24
|
Carlos discussed scaler analysis results with Richard, and also redid some
|
25
|
KaonLT lumi scaler analysis
|
26
|
- plots of T1(SHMS-3/4)/mC, Runs 5157,56,55,54 positive polarity runs
|
27
|
- sees steep rise with current to 55uA
|
28
|
- plots of T3(HMS-ELREAL)/mC, same runs (electron polarity)
|
29
|
- sees flat dependence with current
|
30
|
- also plots T1 vs scaler rate, T3 vs rate
|
31
|
- expects an increase with rate because of increasing randoms and scalers of
|
32
|
course can't eliminate them
|
33
|
|
34
|
Carlos shares his findings from CaFe analysis:
|
35
|
- motivation is to understand why there is a drop in yield (counts/mC) with
|
36
|
beam current observed in CaFe results. The yield has been corrected for
|
37
|
HMS/SHMS track eff, total LiveTime
|
38
|
|
39
|
His slides on relative yields rate dependence are available at:
|
40
|
https://hallcweb.jlab.org/DocDB/0012/001209/001/rate_dependence_slides_cyero_Mar09_2023.pdf
|
41
|
|
42
|
Slides on Event Type Definitions at:
|
43
|
https://hallcweb.jlab.org/DocDB/0012/001208/001/event_type_definitions_COIN_DAQ_cyero_Mar16_2023.pdf
|
44
|
|
45
|
The slides on event type definitions might be helpful when trying to
|
46
|
understand how to properly select which event types based on whether
|
47
|
singles or coincidence was taken, and how to extract the singles and
|
48
|
coincidence separately. This is sort of complementary to Jacob's EDTM
|
49
|
studies. It might help in selecting the right event types when calculating
|
50
|
EDTM live times for singles, in which case you would need to apply the
|
51
|
proper event type cuts.
|
52
|
|
53
|
- in CaFe trigger config, T5 (coin) is formed by overlap of SHMS (T2) and HMS
|
54
|
(T3)
|
55
|
- 1st hit in look-back window gets selected and is used by the analyzer to
|
56
|
form a coincidence time (i.e. T2-T3)
|
57
|
- it is needed to set a cut on the raw TDC spectrum for the relevant
|
58
|
channels (T2, T3), otherwise an accidental at arrives first in the
|
59
|
look-back window will "block" the real coincidence
|
60
|
- generic parameter file: PARAM/TRIG/tcoin.param
|
61
|
- coincidence time (corrected for path length) without this cut looks
|
62
|
drastically worse than one with properly set cuts
|
63
|
- also, the number of events with good CoinTime is higher with good cuts, as
|
64
|
fewer coincidences are blocked by accidentals
|
65
|
|
66
|
- the relative yield from Be9 target versus current is MUCH different after the
|
67
|
cut is optimized
|
68
|
- BEFORE: ~7% decrease in yield is observed to 60uA
|
69
|
- AFTER CUT OPTIMIZATION: the relative yield is nearly flat with current, at
|
70
|
most a ~1%(+/-1%) effect
|
71
|
|
72
|
- Conclusion: it is EXTREMELY IMPORTANT to verify that these cuts are set
|
73
|
correctly (see his example TDC spectra) as otherwise there will be a large
|
74
|
rate dependence and it will be hard to diagnose what is wrong!
|
75
|
|
76
|
Richard Continues this Discussion
|
77
|
---------------------------------
|
78
|
- Nscaler = sum(trig scaler) - EDTM_scaler
|
79
|
- Yield_scaler = Nscaler/Q
|
80
|
- runs 5154-58, 5298, all same setting
|
81
|
|
82
|
Compare Richard & Carlos results
|
83
|
- two main differences
|
84
|
1) when using same current cuts as Carlos, linear dep is restored, same as
|
85
|
Carlos
|
86
|
2) HMS: RT gets linear increase ~5%, while CY sees flat
|
87
|
- also RT gets rates up to ~100kHz, while CY has rates to 60kHz FOR SAME
|
88
|
RUNS
|
89
|
- CT data are middle 4 of RT's 6 runs
|
90
|
|
91
|
** RT, CY will meet to compare exact HMS scaler values on Friday, and see where
|
92
|
the rate discrepancy is coming from. Ali will join online
|
93
|
|
94
|
** Nathan will check PionLT data to see that trigger cuts that Carlos described
|
95
|
are set correctly
|
96
|
** it is important for RT, VK can confirm that these cuts are set correctly on
|
97
|
KaonLT and Summer-2019 data
|
98
|
- check both TDC and yield vs current for a solid target
|
99
|
- CY: also be sure that EDTM-LT cal uses EvtType >= 4, to get all COINS when
|
100
|
counting EDTM-accepted
|
101
|
|
102
|
Richard Updates
|
103
|
---------------
|
104
|
10.6GeV carbon lumi scan analysis
|
105
|
- Lumi #3, 5351-58 (HMS only)
|
106
|
- plots of Yscaler, Ynotrack, Ytrack
|
107
|
- get a consistent HMS current trend, ~2% increase to 65uA
|
108
|
- DG: most reliable is Carbon Track, which is cleanest
|
109
|
- worries about pi- contamination at high rate in scaler-no-track
|
110
|
- Carbon Track seems flat within errors (Runs 5154-5298)
|
111
|
- TH: this looks good, please move to CryoTarget study
|
112
|
- GH: it would be good to check 8 or 6 GeV (whichever is more reliable) so
|
113
|
that we know things are consistent
|
114
|
|
115
|
New script that lists the cuts applied to a given run in python anaysis
|
116
|
- UTIL_KAONLT/bin/cut-check
|
117
|
- used to check what cuts are used to trim the tree created by python
|
118
|
analysis
|
119
|
- hopefully this will be helpful in reducing errors, as otherwise one has to
|
120
|
go through a lot of files to list all of the cuts being applied, which
|
121
|
sounds like a recipe for error
|
122
|
|
123
|
RT, AP, AU meeting tomorrow to discuss the large right SHMS setting rate
|
124
|
discrepancy between the two analyses reported last week
|
125
|
|
126
|
Nathan Updates
|
127
|
--------------
|
128
|
Aerogel Calibs
|
129
|
- last report: PMT4+ calib had a slope vs run#
|
130
|
- broke up into 3 sections
|
131
|
- 2021-22 run now 2 parts, 1st part is most of 1st beam energy
|
132
|
- summer 2022 is 3rd part
|
133
|
- takes simple mean for each period, not an error-wt average, since errors seem
|
134
|
unreliable
|
135
|
- excluded some runs where NPE histo fits bad
|
136
|
- put new params in hcana DBASE
|
137
|
- replayed data w/new calib has better pion-separation from zeros in Aerogel
|
138
|
compared to online replay
|
139
|
|
140
|
Still waiting on Junaid for RefTime cut check
|
141
|
|
142
|
GHP presentation
|
143
|
- hoping Tanja has scaled cross section values from Scaling Paper at x=0.311
|
144
|
otherwise he will have to read them off the plot, which is less accurate
|
145
|
- will circulate draft slides next week
|
146
|
|
147
|
Junaid Updates
|
148
|
--------------
|
149
|
DC Calibs
|
150
|
- found some errors in PionLT.database file and kinematics file
|
151
|
- some missing run numbers and typos
|
152
|
- SHMS DC look ok
|
153
|
- HMS DC calibs off, still needs work
|
154
|
- getting high #events at zero bin entry
|
155
|
- DG: what kine of runs are these? Production
|
156
|
- part of the trick is to apply electron cuts to electron arm, otherwise
|
157
|
get a lot of junk that distort the DC spectra
|
158
|
- calib code is awkward, as you have to "hard wire" the PID cuts
|
159
|
- MJ still checking things, HMS plane residuals improving
|
160
|
|
161
|
Ali Updates
|
162
|
-----------
|
163
|
First look at pion-PID for high Q2 data
|
164
|
- shows plots where progressively apply cuts
|
165
|
- get a clean sample after HG cut, no obvious K/p leakage
|
166
|
- CoinTime cut helps a lot too
|
167
|
- DG: only apply the RF cut if you have to
|
168
|
- if get clean spectra w/o it, only HGC, then RF cut only gives an
|
169
|
inefficiency
|
170
|
- AU will need to look at # clean neutrons vs cut
|
171
|
|
172
|
Vijay Updates
|
173
|
-------------
|
174
|
Working on pion-LT package
|
175
|
- almost ready
|
176
|
- debugged Fortran scripts, made changes for Summer 2019 data
|
177
|
- working on cross section function for SIMC, getting some errors as GH gave
|
178
|
him an old version of the routine and some variable names have changed
|
179
|
- DG: yes, some variable names in the structures have changed, see include
|
180
|
files
|
181
|
|
182
|
Next Meeting
|
183
|
------------
|
184
|
- Thur Mar 31 @ 16:00 Eastern/14:00 Regina/13:00 Pacific
|
185
|
- KaonLT will go first
|
186
|
|
187
|
|
188
|
|
189
|
|
190
|
|
191
|
|
192
|
|
193
|
|
194
|
|
195
|
|
196
|
|
197
|
|
198
|
|
199
|
|
200
|
|