1
|
Sept 5/24 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, Alicia Postuma, Muhammad Junaid, Nathan Heinrich,
|
13
|
Ali Usman, Nacer Hamdi
|
14
|
Virginia - Richard Trotta
|
15
|
JLab - Dave Gaskell
|
16
|
CUA - Tanja Horn
|
17
|
Ohio - Julie Roche
|
18
|
CSULA - Konrad Aniol
|
19
|
|
20
|
Junaid
|
21
|
------
|
22
|
PionLT Heep Coin Hodo-3/4 effs
|
23
|
- new plot of 3/4 effs vs SHMS momentum
|
24
|
- planes 1,2,3 show a linear drop of eff with increasing momentum
|
25
|
- plane 4 shows a low eff, but very large variation (50-75%) at P~2.4 GeV/c
|
26
|
- GH: the drop off at low momentum is expected, but the linear drop at
|
27
|
high momentum is puzzling
|
28
|
- Dave: *NB* confident that there is an issue in the PID selection cuts used
|
29
|
to calculate these efficiencies, that needs to be sorted out first, and
|
30
|
then we can come back to this
|
31
|
|
32
|
- finds the 3/4 effs are sensitive to the ptof_tolerance value
|
33
|
- PionLT value=10, apparently set by C. Yero at start of PionLT run
|
34
|
- changing to 15 results in slightly LARGER efficiencies
|
35
|
- Ali: please recall the beta vs. cointime issue KaonLT encountered last
|
36
|
year, this was because the tolerance variable was set wrong (was 100, it
|
37
|
was supposed to be 2)
|
38
|
- a value of 2 should be applied to ALL data
|
39
|
- Junaid had tried this as well, results in even LOWER effs
|
40
|
|
41
|
More discussion about that to do about the 3/4 effs:
|
42
|
|
43
|
- Dave: NCER cut needs to be removed entirely, makes sense for negative
|
44
|
polarity SHMS, but no sense at all for positive polarity
|
45
|
- Junaid tried varying the CERCUT parameter, all the way down to zero
|
46
|
(default is 2pe) and it made NO difference to the efficiencies
|
47
|
- note Nathan did calibrate the NGC, so the pedestals should have a
|
48
|
reasonable value when applying npe=0 cut that Junaid tried
|
49
|
|
50
|
- Nathan thinks the CERCUT is legacy code that is not called, since for
|
51
|
example no NGC was installed for KaonLT. What did the code do then?
|
52
|
- we spend some time looking at the code, it seems that the HGC is being
|
53
|
called here instead of the NGC
|
54
|
- however, protons also never Cherenkov in the HGC either, so this cut also
|
55
|
should not be applied to Heep data effs
|
56
|
- Dave suggests to remove this code completely from hcana and see what
|
57
|
happens
|
58
|
- Dave looks a bit more at the code, and sees that in THcHodoEff.cxx both
|
59
|
Cherenkov and Calorimeter cuts are used
|
60
|
- *NB* Junaid will investigate the code more thoroughly
|
61
|
|
62
|
Nathan
|
63
|
------
|
64
|
PionLT Lumi studies
|
65
|
- new plots using Zach's code applied to 2021 data
|
66
|
- SHMS Carbon 12143-12154
|
67
|
- Scaler analysis shows a small linear rate dep
|
68
|
- ElLT not yet applied
|
69
|
- will look at fixing slopes on Carbon Scaler/NoTrack/Track analyses
|
70
|
|
71
|
- Mild antiboiling seen for LH2 12192-12198
|
72
|
- Dave and Tanja expect Coin data are not as reliable for determining
|
73
|
boiling corr than the 2022 dedicated singles runs
|
74
|
|
75
|
- One issue is what to use for TLT
|
76
|
- right now, Nathan picks the higher TLT from either EDTM or CPU*ElLT
|
77
|
methods
|
78
|
- GH: *NB* this is "cherry picking", much better to decide which method is
|
79
|
more reliable for a given run period, and then use it consistently for
|
80
|
both Heep and physics data
|
81
|
|
82
|
Richard
|
83
|
-------
|
84
|
KaonLT Ratio debugging for Q2=4.4, W=2.74
|
85
|
- everything is looking pretty good, except that Ratios converge around 1/3
|
86
|
instead of 1
|
87
|
- ratios much flatter than before
|
88
|
- comparisons of data and MC distribution shapes are reasonable
|
89
|
- Richard reports that this factor is consistent across all settings, not
|
90
|
just the Q2=4.4 shown today
|
91
|
- GH and Tanja agree there seems to be a missing factor of ~pi somewhere in
|
92
|
one section of the code that is consistent with the others. When this is
|
93
|
identified, things should look fairly good
|
94
|
|
95
|
- various studies to try to find the missing factor
|
96
|
- reran SIMC from scratch instead of recalculating weights -- ratio unchanged
|
97
|
- looking at histo binning in Exp Yield calc
|
98
|
- also will check SIMC w/o recon_hcana
|
99
|
- should change the distribution shapes, but not the normalization
|
100
|
|
101
|
- Tanja suggested to run for only a small t range (~0.1) to check if it's a
|
102
|
functional issue w/t-parameteriation
|
103
|
|
104
|
- Nacer: try applying your code to Vijay's data and hunt for the difference
|
105
|
between your plots and his
|
106
|
- Vijay should be back and available for a meeting with Richard next week
|
107
|
|
108
|
Nacer
|
109
|
-----
|
110
|
KaonLT Heep studies
|
111
|
- proton absorption correction
|
112
|
- reviewed John Matter's thesis and spreadsheet
|
113
|
- they compared calculated and experimentally determined proton absorption
|
114
|
results:
|
115
|
A_pred=8.56% spreadsheet
|
116
|
A_exp=9.03+/-0.07% by looking for missing protons other spectrometer in
|
117
|
Heep singles data
|
118
|
|
119
|
- modified their table, removed NGC, need to adjust aerogel density for actual
|
120
|
values used during KaonLT
|
121
|
- 0.2 density should correspond to n=1.03 aerogel
|
122
|
- Tanja: 1.015 density=0.08, 1.011 density=0.04
|
123
|
all trays should have the same thickness aerogel. The area will be
|
124
|
different though
|
125
|
- the correction factor ~7.5%, applying this factor changes the ratio from
|
126
|
near 1.0 to ~1.08, except for 10.6GeV data where R=1.00 after factor
|
127
|
|
128
|
- tried different MM cuts to exclude pi0
|
129
|
- ratios get a bit smaller for MM<0.10, 0.09, 0.08 cuts, but all above 1.0
|
130
|
except 10.6GeV
|
131
|
|
132
|
- HMS matrix element check, hcana default vs old 6 GeV matrix
|
133
|
- data vs MC shape difference gets a bit worse (8.2 GeV data)
|
134
|
- ratios also get a few % worse, which is a bit surprising, only expected
|
135
|
events to move around but not get lost
|
136
|
- clearly the default matrix elements are better, this issue closed
|
137
|
|
138
|
More discussion about that to do about the ratios with absorption corr included:
|
139
|
|
140
|
- Ali: asks about the momentum dependence to the interaction length in
|
141
|
spreadsheet
|
142
|
- Matter's thesis: cross section indep of momentum >1 GeV/c
|
143
|
- confirmed with their data for 4-10 GeV/c
|
144
|
- our data are 2-8 GeV/c should be fine
|
145
|
- Dave: the proton absorption correction is applied by other analyses, it is
|
146
|
not optional as there are no proton absorption effects in SIMC
|
147
|
|
148
|
- Ali: there is a cut applied on aero tray in the replay, irrespective of
|
149
|
whether an aero NPE cut is applied to the data
|
150
|
- *NB* need to be 100% sure the same cut is applied also to SIMC
|
151
|
|
152
|
- GH: *NB* the MM distribution clearly has a much larger tail in data compared
|
153
|
to MC, independent of the pi0 contamination. This tail is ~5%, so seems to
|
154
|
be the cause of the Ratio problem
|
155
|
- Dave: *NB* please try looking at Emiss, W for Data & MC, not so used to
|
156
|
looking at MM for Heep data, maybe it will show something
|
157
|
|
158
|
- Dave: *NB* it is important to look at some Single Arm Heep to see if there
|
159
|
are similar issues
|
160
|
- GH agrees this is an excellent idea. Since no protons are detected, no
|
161
|
absorption corr is needed to be applied, so it will give important
|
162
|
information
|
163
|
|
164
|
- *NB* GH's suspicion of what's wrong (added after meeting):
|
165
|
- the spreadsheet assumes that interacting protons are completely lost and
|
166
|
no trigger is made
|
167
|
- if the protons scatter but still give a trigger, they end up in the
|
168
|
Heep tail, and the spreadsheet correction is an over-estimate
|
169
|
- if so, we would need to apply a much tigher MM cut to remove the
|
170
|
interacting protons not included in the SIMC tail
|
171
|
- hopefully the W and EM plots will shed light on what cut to apply
|
172
|
- the Heep singles Data/MC ratios will also be critical to this
|
173
|
- Henk came up with a modified method to evaluate the absorption
|
174
|
correction, looking at the proportions of events in peak and tail, to
|
175
|
accont for this
|
176
|
- *NB* GH will review this info and see how to apply it here
|
177
|
|
178
|
Alicia
|
179
|
------
|
180
|
Update from PRL:
|
181
|
- referee reports received Aug 12,14,31
|
182
|
- queries sent to 2 referees after receiving 3rd referee report, 1 response
|
183
|
received so far
|
184
|
- we infer that one referee gives a report disagreeing with the other 2, and
|
185
|
they have been asked to comment on it
|
186
|
|
187
|
Garth
|
188
|
-----
|
189
|
Hall C Winter Workshop Jan 13-14/25
|
190
|
- at last meeting, discussed a session on "nearly final results" and what
|
191
|
KaonLT might have ready by then
|
192
|
- 3 talks seem likely:
|
193
|
- Vijay, Richard, Ali
|
194
|
- Richard and Vijay should coordinate e.g. on how best to show the various
|
195
|
systematic studies
|
196
|
- Tanja: what about general experiment updates?
|
197
|
- GH: this was not explicitly discussed. Since there are starting to be lots
|
198
|
of results nearing publication (unlike a few year ago), priority will be
|
199
|
given first them.
|
200
|
- once that part of the program is fleshed out, we will likely go back and
|
201
|
ask for more general updates
|
202
|
|
203
|
Next Meeting
|
204
|
-------------
|
205
|
- Thur Sept 12 @ 16:30 Eastern/14:30 Regina
|
206
|
- KaonLT will go first
|
207
|
- *NB* Please note the new time and new Zoom ID!
|
208
|
|
209
|
|
210
|
|
211
|
|
212
|
|
213
|
|
214
|
|
215
|
|
216
|
|
217
|
|