1
|
Feb 22/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 - Nathan Heinrich, Garth Huber, Ali Usman, Alicia Postuma,
|
13
|
Muhammad Junaid, Nacer Hamdi
|
14
|
CSULA - Konrad Aniol
|
15
|
CUA - Richard Trotta, Tanja Horn
|
16
|
York - Stephen Kay
|
17
|
FIU - Pete Markowitz
|
18
|
|
19
|
Nathan
|
20
|
------
|
21
|
Writing Lumi analysis script for PionLT
|
22
|
- modifying script from Richard, which always assumed an active HMS and SHMS
|
23
|
trigger for each run, as some PionLT Lumi scans were single arm only
|
24
|
- had to make similar edits earlier for CoinLumi script, so confident of the
|
25
|
needed changes
|
26
|
- edits made, still testing
|
27
|
|
28
|
Script for plots for SHMS NIM paper
|
29
|
- started to modify old plotting script
|
30
|
- looked briefly so far at the draft NIM paper text
|
31
|
- once the plotting code has been run, will look up needed details for the
|
32
|
paper
|
33
|
|
34
|
Junaid
|
35
|
------
|
36
|
PionLT efficiency studies
|
37
|
- investigating outlier runs identified last week, about 30 runs in total
|
38
|
- will make some decisions soon on which runs to discard
|
39
|
- then will post efficiency slides on RedMine
|
40
|
|
41
|
Ali
|
42
|
---
|
43
|
SIMC resolution study
|
44
|
- on hold until Alicia finishes HMS delta-correction study
|
45
|
|
46
|
Setting up for pi+n LT-analysis
|
47
|
- replayed 3 of 5 settings so far (both epsilons)
|
48
|
- of these, analysis completed for 2 of 3 settings
|
49
|
- will give Richard clean pi+n data for subtraction
|
50
|
|
51
|
- updated SIMC w/ physics_iterate.f model
|
52
|
- compared to Vijay's code, almost identical except for different Q2-settings
|
53
|
|
54
|
- looking at Richard's LT-sep scripts, hoping to have some initial ratios soon
|
55
|
|
56
|
Info re HGC cut issue
|
57
|
- no HGC cut intended to be used for pi+n analysis, get a small K+L leakage,
|
58
|
but in MM region well above pi+n
|
59
|
- this means pi+n analysis will also not need an HGC efficiency to be applied
|
60
|
- on the other hand, pi+Delta0 region is affected by K+L leakage, and will
|
61
|
need an HGC cut as well as an HGC efficiency to be applied
|
62
|
- the pi+ sample for Richard will not have an HGC cut applied
|
63
|
- will give skimmed ROOT files w/o HGC cut to RT, he will have to apply the
|
64
|
HGC cut before subtracting pi+ leakage sample from K+ data
|
65
|
|
66
|
Richard
|
67
|
-------
|
68
|
Comparing SIMC model to preliminary data
|
69
|
- looked at 4 of 6 settings (including Q2=0.5)
|
70
|
- Q2=0.5 has lots of issues that he will flag for Nacer to eventually look at
|
71
|
- standard.kinematics issues:
|
72
|
- need to confirm that CoinTime offset is set correctly for EVERY run
|
73
|
- SHMS_RF_Offset was also not set for some low Q2 runs
|
74
|
|
75
|
Initial look at Q2=2.115, W=2.95
|
76
|
- all LT, TT terms set to zero for now in model
|
77
|
- sigL model: adjusted exp(|t|-0.2) factor to exp(|t|+0.2), seems to give more
|
78
|
stable ratios vs t
|
79
|
- sigT model: compared p5+p6*log(Q2) with alternative p5*log(Q2)+p6/Q4
|
80
|
|
81
|
- also compared to Tanja's global K+ analysis in PRC 85, 018202 (2012)
|
82
|
- checked W-factor and sigT models in this paper
|
83
|
- W dependence: tried changing W-factor=1/(W2-Mp2)^2 to 1/(W2-Mp2)^3
|
84
|
- Tanja's paper has the power as 2.5
|
85
|
- when changing the W-factor, found that the sig_UNS expt error bars changed
|
86
|
a lot, which should have not happened
|
87
|
- this indicates a coding error that needs fixing
|
88
|
|
89
|
- then tried sigT=p5/(1+p6*Q2) and old W-factor
|
90
|
- manually changing p5 from 120 to 25 gave better ratios
|
91
|
- doing things by hand right now, instead of using MINUIT to get new params
|
92
|
|
93
|
Initial look at Q2=3.0, W=3.4
|
94
|
- same procedure, changing sigT and sigL factors
|
95
|
- lowering sigT p5 param improves ratios, just as at Q2=2.115
|
96
|
|
97
|
- *WARNING* data have NO MM-cut applied at all, not only do the exp yields
|
98
|
contain pi+ and proton leakthrough, they also include all K+ contributions
|
99
|
from inclusive spectrum above Lambda peak
|
100
|
- right, now sees a weird phi Left-Right asymmetry in the data
|
101
|
- the experimental yields will CHANGE SIGNIFICANTLY once Richard has a clean
|
102
|
sample
|
103
|
|
104
|
- Ali offers to generate for Richard also a clean proton sample for subtraction
|
105
|
|
106
|
Initial look at Q2=3.0, W=2.32
|
107
|
- looked at using W-factor with ^3 exponent
|
108
|
- getting negative sigL for lowest -t bin, CAVEATS in data as above
|
109
|
|
110
|
Compared initial results vs x_B
|
111
|
- early indications that a systematic change in the model vs x_B will be
|
112
|
needed, will be interesting if this holds up in more thorough analysis
|
113
|
|
114
|
GH: as a reminder, we will eventually have to look again at the residual
|
115
|
Q2-dependence of the Heep-Coin ratios, and resolve this issue before we can
|
116
|
finalize the data and publish
|
117
|
- this should happen after Richard graduates
|
118
|
|
119
|
Alicia
|
120
|
------
|
121
|
Implementing HMS delta-correction, hsdelta_NEW=hsdelta-p0*hsxpfp
|
122
|
- took MM calc from Richard's reconstruction script and implemented calc for
|
123
|
data
|
124
|
- as a check, compared new MM calc w/p0=0 to old MM calc
|
125
|
- overall very similar, but a slight difference between them
|
126
|
- plot difference in #events vs Mx, most pi+n events within +/-2 MeV compared
|
127
|
to old MM calc
|
128
|
- best guess is that there is some rounding error difference in the HMS
|
129
|
momentum calc compared to hcana when applying offsets
|
130
|
- GH: as a coding check, can try comparing old and new MM w/o offset
|
131
|
|
132
|
- then applied Dave's p0 values and compared the old-new data
|
133
|
- MM gets wider instead of narrower, something still seems wrong
|
134
|
- GH: suggest to make a plot of hsdelta_NEW vs hsxpfp to see what kind of
|
135
|
correlation you get, should be a straight vertical band if applied
|
136
|
correctly
|
137
|
- Alicia will send her script to Richard for review
|
138
|
|
139
|
Nacer
|
140
|
-----
|
141
|
Low Q2 KaonLT calibs
|
142
|
- finished DC and Hodo timewalk calcs
|
143
|
- pull request to GitHub was accepted
|
144
|
- did a full replay, including all of Ali and Richard's improvements
|
145
|
|
146
|
- plot of EDTM vs Coin_trig_rate
|
147
|
- sees sharp drop off above 2.7kHz trig rate, as expected, due to buffered
|
148
|
mode data acquisition (good news)
|
149
|
- plot of COIN_CPULT vs COIN_trig_rate
|
150
|
- similar behavior, looks good
|
151
|
- plots of HMS(e), SHMS(pi) SING_TRK_Eff vs S1X hodo rate
|
152
|
- no HMS outliers
|
153
|
- a few SHMS outliers that need investigation
|
154
|
|
155
|
To Do:
|
156
|
- check other plots
|
157
|
- move to K+ PID
|
158
|
- also should get Lumi scripts from Richard and look at Low-Q2 Lumi scans
|
159
|
|
160
|
Next Meeting
|
161
|
------------
|
162
|
- Thur Feb 29 @ 15:00 Eastern/14:00 Regina
|
163
|
- KaonLT will go first
|
164
|
|