|
1
|
Feb 19-20/26 PionLT/KaonLT Analysis Meeting Notes
|
|
2
|
-------------------------------------------------
|
|
3
|
(Notes by GH)
|
|
4
|
|
|
5
|
Today: PionLT will be discussed first
|
|
6
|
|
|
7
|
Thursday: Present
|
|
8
|
-----------------
|
|
9
|
Regina - Garth Huber, Nacer Hamdi, Muhammad Junaid, Vijay Kumar,
|
|
10
|
Alicia Postuma, Nermin Sadoun
|
|
11
|
JLab - Dave Gaskell
|
|
12
|
CSULA - Konrad Aniol
|
|
13
|
CUA - Tanja Horn, Chi Kin Tam, Sameer Jain
|
|
14
|
FIU - Pete Markowitz
|
|
15
|
|
|
16
|
|
|
17
|
Junaid
|
|
18
|
------
|
|
19
|
PionLT Q2=3.85 analysis
|
|
20
|
- started working on PhD thesis and on systematic studies in parallel
|
|
21
|
- will meet w/Nathan on systematic studies already completed
|
|
22
|
|
|
23
|
Vijay
|
|
24
|
-----
|
|
25
|
Low Q2 PionLT analysis
|
|
26
|
- presented new model systematics studies
|
|
27
|
- followed instructions in Blok et al. paper:
|
|
28
|
- varied sigL, sigT +/-10% then 2 iterations
|
|
29
|
- effect is very small overall, with 1st t-bin is more sensitive:
|
|
30
|
sigL only 0.004 ub/GeV2 variation
|
|
31
|
- apparently this is not what Tanja did, not sure why the text in the paper
|
|
32
|
is written this way. A change in the parameters should result in
|
|
33
|
negligible change after 2 iterations as the parameters return to their
|
|
34
|
optimal values
|
|
35
|
- Tanja suggests instead to change the functional form in some way and see
|
|
36
|
what that change is
|
|
37
|
- Garth: Vijay already has a slightly different parameterization (from
|
|
38
|
Nov/25), so all he needs to do is compare the latest paramterization
|
|
39
|
result with that one
|
|
40
|
- the difference is that the old one has no Q2-dep in the functional
|
|
41
|
form, while the new one does
|
|
42
|
|
|
43
|
Alicia
|
|
44
|
------
|
|
45
|
u-channel analysis
|
|
46
|
- comparing data and MC kinematics plots for Q2=3.0, W=3.14, high eps, center
|
|
47
|
- apply MC normalization params from MM shape study to other kinematic
|
|
48
|
variables
|
|
49
|
- two types of plots presented: sum of MC compared to data, subtract bkd MCs
|
|
50
|
from data and compare difference with omega MC
|
|
51
|
- encouraging level of agreement, but more checks needed
|
|
52
|
- ssypfp, ssxpfp distributions have a lot of zero events, need to check if
|
|
53
|
1-2 runs for the setting were replayed improperly
|
|
54
|
- removed the zero events for now, and compared shapes with MC
|
|
55
|
- need to recheck MC normalizations to data
|
|
56
|
|
|
57
|
- working in parallel on Phi-007 MM reconstruction resolution for Expt
|
|
58
|
Readiness Review w/ Henry Klest
|
|
59
|
|
|
60
|
Chi Kin
|
|
61
|
-------
|
|
62
|
KaonLT high Q2 analysis
|
|
63
|
- looking into discrepancy between his results and Richard's
|
|
64
|
- narrowed down that difference is not coming from cuts
|
|
65
|
- next possibility is that it might be the SIMC weights
|
|
66
|
- meeting tomorrow w/ Richard and will report again tomorrow
|
|
67
|
|
|
68
|
- started drafting paper on Overleaf
|
|
69
|
|
|
70
|
Sameer
|
|
71
|
------
|
|
72
|
PionLT luminosity analysis
|
|
73
|
- working with Nathan's lumi scripts to learn how to use them and apply to his
|
|
74
|
own analysis
|
|
75
|
|
|
76
|
Nacer
|
|
77
|
-----
|
|
78
|
Low Q2 KaonLT analysis
|
|
79
|
- further work on MC model optimization
|
|
80
|
- prior best model:
|
|
81
|
- functional form only for sigT, sigL,LT,TT=0 functional form
|
|
82
|
- ratios were fairly good, but some wiggles
|
|
83
|
- fit of function to sigL,T results vs t not very good
|
|
84
|
|
|
85
|
- new changes:
|
|
86
|
- added another t-bin, since sigL shows a large jump between bins 1,2
|
|
87
|
- New limits (8 bins):
|
|
88
|
0.070, 0.084, 0.090, 0.096, 0.102, 0.108, 0.116, 0.128, 0152
|
|
89
|
Old limits (7 bins):
|
|
90
|
0.070, 0.086, 0.094, 0.100, 0.106, 0.114, 0.126, 0.152
|
|
91
|
- has ~10k events/bin
|
|
92
|
- then used only sigT model for new 8-bins, Data/MC ratios changed more
|
|
93
|
than expected but still near 1
|
|
94
|
|
|
95
|
- change to simple functional forms for L,LT,TT same as Chi Kin
|
|
96
|
sigL = p2*exp(-|p3*t|)
|
|
97
|
sigLT= p4*exp(-|p5*t|)*sin(theta*)
|
|
98
|
sigTT= p6*exp(-|p7*t|)*sin^2(theta*)
|
|
99
|
- Data/MC ratios look very good, close to 1 and much flatter than before,
|
|
100
|
even the highest -t bins look good
|
|
101
|
- fit of function to sigT,LT,TT vs t also looks pretty good
|
|
102
|
- new 2nd t-bin has sigL slightly higher than zero by ~1.5 sigma, bins 3-8
|
|
103
|
remain statistically consistent with zero
|
|
104
|
- kinematic comparison plots (Data vs MC):
|
|
105
|
- HMS xptar looks better after adjusting HMS ME offset
|
|
106
|
- Dave: a few nagging discrepancies but not bad agreement overall
|
|
107
|
|
|
108
|
- Dave: still concerned that rise in sigL at low -t is anti-correlated with
|
|
109
|
drop in sigT for same bin, need to demonstrate that this is real, and not an
|
|
110
|
analysis artifact
|
|
111
|
- are the MM offsets different at low and high epsilon?
|
|
112
|
- the difference in offsets is ~2.5 MeV
|
|
113
|
- The MM offset is applied only to MM, and is not propagated to the other
|
|
114
|
kinematic variables. The concern is that MM is correlated with t, so
|
|
115
|
significantly different MM offsets at high and low epsilon could imply the
|
|
116
|
t-distributions are shifted by different amounts at high and low epsilon
|
|
117
|
and this could give rise to anti-correlated behavior between sigT,L
|
|
118
|
|
|
119
|
*NB* Garth will work with Nacer to estimate the implied shift in t by a 2.5
|
|
120
|
MeV MM shift
|
|
121
|
|
|
122
|
|
|
123
|
Friday: Present
|
|
124
|
---------------
|
|
125
|
Regina - Garth Huber, Ailcia Postuma, Nacer Hamdi, Nermin Sadoun,
|
|
126
|
Muhammad Junaid, Vijay Kumar
|
|
127
|
York - Stephen Kay
|
|
128
|
JMU - Gabriel Niculescu
|
|
129
|
CUA - Tanja Horn, Chi Kin Tam
|
|
130
|
Virginia - Richard Trotta
|
|
131
|
FIU - Pete Markowitz
|
|
132
|
Glasgow - Rachel Montgomery
|
|
133
|
|
|
134
|
|
|
135
|
Nermin
|
|
136
|
------
|
|
137
|
PionLT LD2 analysis
|
|
138
|
- met w/ Kathleen Ramage yesterday on analysis issues
|
|
139
|
|
|
140
|
- RF times offset correction for LD- Q2=1.6, W=3.08, high eps=9.16GeV
|
|
141
|
- before: RF time distribution had 2 peaks separated by 2ns, both consistent
|
|
142
|
with pi-p missing mass
|
|
143
|
- determined that some runs needed the RF time to be offset by 2ns compared
|
|
144
|
to the others
|
|
145
|
- the issue, as we learned during the experiment, is that sometimes Hall C
|
|
146
|
receives the even beam bunches, and sometimes the odd beam bunches, and
|
|
147
|
we have to correct the odd ones by 2ns to bring them into time with the
|
|
148
|
even ones
|
|
149
|
- the shifted runs were: 12043,44 which were in the middle of the setting
|
|
150
|
- after: RF time has 1 nice peak
|
|
151
|
- shows RF time plot w/o and w/ CT, Hcer, Hcal, Acceptance Cuts:
|
|
152
|
background is much reduced with these cuts applied
|
|
153
|
|
|
154
|
- NGC plots have most events near zero
|
|
155
|
Mean=0.134NPE #Entries=646774
|
|
156
|
- apply Pcal>0.85 cut to select electrons in SHMS
|
|
157
|
Mean=0.153NPE #Entries=305096
|
|
158
|
- ~50% of SHMS events are electrons
|
|
159
|
- Nermin looks up LD analysis report by Argha Das
|
|
160
|
- NGC mean is ~4 NPE for a different setting, much higher than here
|
|
161
|
- Garth: suggests to look at low epsilon data for same setting to see if
|
|
162
|
the issue is just these runs, or overall
|
|
163
|
- Nermin will look into the NGC calibration scripts
|
|
164
|
|
|
165
|
Rachel
|
|
166
|
------
|
|
167
|
PionLT LD2 analysis
|
|
168
|
- Kathleen looking at a different LD- setting than Nermin, sees a very large
|
|
169
|
#randoms, hard to tell where the prompt peak is
|
|
170
|
- Nacer: suggests to look at CoinTime with MM cut applied, it should make the
|
|
171
|
prompt peak stand out more
|
|
172
|
- Garth: the LD- settings have high random rates due to electrons in the
|
|
173
|
SHMS, it is easier to set up the analysis first on LD+ data for same
|
|
174
|
kinematics
|
|
175
|
|
|
176
|
Richard
|
|
177
|
-------
|
|
178
|
KaonLT high Q2 analysis
|
|
179
|
- looking into discrepancy between his results and Chi Kin's
|
|
180
|
- compared results bin-by-bin
|
|
181
|
- Kin's yields ~50% higher
|
|
182
|
- found 2 differences in cuts:
|
|
183
|
- CKT has optimized the diamond cut
|
|
184
|
- HGC hole cut also slightly different, Richard needs to confirm that x,y
|
|
185
|
axes are applied correctly (so the hole cut is where it is intended)
|
|
186
|
- SIMC weights were consistent, so that's not the issue
|
|
187
|
- believes the issue is caused by the order of cuts being applied
|
|
188
|
- in principle, cuts are supposed to be commutative (i.e. the order doesn't
|
|
189
|
matter), but possibly there is a NumPy array issue which screws this up
|
|
190
|
- Gabriel: is there a milestone plan on who checks what?
|
|
191
|
- Tanja: it's important to try to keep some things constant (such as
|
|
192
|
identical cuts between RT and CKT) so the cause of the difference can be
|
|
193
|
identified more readily
|
|
194
|
|
|
195
|
|
|
196
|
Next Meetings
|
|
197
|
------------------
|
|
198
|
- Thurs: Feb 26 @ 16:00 Eastern/15:00 Regina
|
|
199
|
- KaonLT will go first
|
|
200
|
|
|
201
|
- Fri: Feb 27 @ 11:00 Eastern/10:00 Regina
|
|
202
|
- we will continue where we left off
|
|
203
|
|
|
204
|
|
|
205
|
|
|
206
|
|
|
207
|
|
|
208
|
|
|
209
|
|
|
210
|
|
|
211
|
|