|
1
|
Jan 15-16/26 PionLT/KaonLT Analysis Meeting Notes
|
|
2
|
-------------------------------------------------
|
|
3
|
(Notes by GH)
|
|
4
|
|
|
5
|
Today: KaonLT will be discussed first
|
|
6
|
|
|
7
|
Thursday: Present
|
|
8
|
-----------------
|
|
9
|
Regina - Garth Huber, Muhammad Junaid, Alicia Postuma, Nermin Sadoun,
|
|
10
|
Nacer Hamdi
|
|
11
|
CUA - Sameer Jain, Tanja Horn, Chi Kin Tam
|
|
12
|
Virginia - Richard Trotta
|
|
13
|
Ohio - Julie Roche
|
|
14
|
FIU - Pete Markowitz
|
|
15
|
|
|
16
|
Chi Kin
|
|
17
|
-------
|
|
18
|
KaonLT High Q2 LT-sep
|
|
19
|
- diamond cut adjustments for:
|
|
20
|
Q2=2.1 W=2.95; Q2=3.0 W=3.14; Q2=4.4 W=2.74; Q2=5.5 W=3.02
|
|
21
|
- compared data and SIMC for high, low epsilon, take the minimum overlap
|
|
22
|
region
|
|
23
|
- previous studies only looked at center SHMS settings, now overlaying
|
|
24
|
left-right-center SHMS for both epsilons, and made some adjustments
|
|
25
|
|
|
26
|
Q2=3.0, W=2.32
|
|
27
|
- found an issue with high epsilon RIGHT setting
|
|
28
|
- diamond is smaller than center & left SHMS, missing low W part of diamond
|
|
29
|
- will reduce the size of the diamond for all settings at this Q2-W to
|
|
30
|
compensate
|
|
31
|
- Richard: this was the very first KaonLT setting taken, something may have
|
|
32
|
been off
|
|
33
|
- *NB* will take another look at the run kinematics to see if we can
|
|
34
|
identify what caused this
|
|
35
|
|
|
36
|
- Q2=4.4 new iterations
|
|
37
|
- compared with Richard's results from 11/06/25
|
|
38
|
- Richard had sigT ~8, while CKT has ~60
|
|
39
|
- *NB* Richard and Chi Kin will cross-check
|
|
40
|
|
|
41
|
Sameer
|
|
42
|
------
|
|
43
|
KaonLT CoinTime Blocking
|
|
44
|
- corrected all rates for BeamOnTime
|
|
45
|
- had a question on whether the beam current cuts have been set correctly
|
|
46
|
- looked at the run sheet target current, in some cases this was different
|
|
47
|
than BCM2 value, which was then used
|
|
48
|
- showed some plots of CoinBlock vs corrected CoinRate
|
|
49
|
- trends with rate generally looked a lot better now
|
|
50
|
- Q2=3.0 W=3.14: SHMS Right CoinBlock are low compared to Left and Center
|
|
51
|
- Garth: it could either be that the rate is calculated too low and the
|
|
52
|
point is really further to the right, or that the CoinBlock is wrong and
|
|
53
|
should really be higher
|
|
54
|
- *NB* Tanja: SJ should check if the rate is too low, or if the cuts needed
|
|
55
|
for the correction are off for this setting
|
|
56
|
|
|
57
|
Nacer
|
|
58
|
-----
|
|
59
|
KaonLT Q2=0.5 LT-sep
|
|
60
|
- tried a variety of L/T parameterizations
|
|
61
|
- the best Data/MC ratios are found with
|
|
62
|
sigT = p1/Q2 + |t|**p2/(Q2+p3)**2 and L=LT=TT=0
|
|
63
|
- tried to add back a function for sigL, such as
|
|
64
|
sigL = p4/Q2 + p5/|t|**p6
|
|
65
|
which results in the ratios getting worse
|
|
66
|
- also, sigT has a drop in cross section at low -t and he hasn't been able to
|
|
67
|
find a funtion which describes this while getting a good fit at low -t
|
|
68
|
- Junaid: what do the kinematic plot comparisons look like?
|
|
69
|
- Garth: more important than the stability of the Data/MC ratios is the
|
|
70
|
stability of the resulting cross sections
|
|
71
|
- *NB* please make a spreadsheet comparing sigL,T for a model w/sigL and a
|
|
72
|
model w/o sigL and let's see to what degree the cross sections are stable
|
|
73
|
- if they are stable within errors, and the kinematic plots look good, then
|
|
74
|
that's good enough. If the cross sections are not stable, then we need
|
|
75
|
to figure out next steps.
|
|
76
|
|
|
77
|
- K+ absorption correction
|
|
78
|
- Junaid helped getting the Geant4 code to work
|
|
79
|
- removed NGC, the other cuts in the script seem applicable
|
|
80
|
- get a corretion of 4.2% at 2.5 GeV/c
|
|
81
|
- *NB* it would be good to cross check with Vijay's correction
|
|
82
|
|
|
83
|
- using Nathan's LH2 boiling factor
|
|
84
|
- *NB* should cross check it with what Chi Kin/Richard use
|
|
85
|
|
|
86
|
- Sameer has KaonLT Q2=0.5 runs, will look at CoinBlock for these setting
|
|
87
|
- Nacer will communicate on cuts and other issues
|
|
88
|
|
|
89
|
Alicia
|
|
90
|
------
|
|
91
|
KaonLT u-channel background subtraction
|
|
92
|
- found a work around for the weird Pythia issues at Q2=3.0 W=2.32
|
|
93
|
- double bump in Python MM distribution
|
|
94
|
- this setting also has worse MM resolution than the others
|
|
95
|
- previously, found that adjustments to Pythia input parameters made no
|
|
96
|
difference
|
|
97
|
- however, applying acceptance cuts makes the MM distribution smoother
|
|
98
|
- plots of ssxpfp,ssypfp for MM<0.9(rho-omega) and MM>0.9(phasespace) regions
|
|
99
|
show some interesting correlations
|
|
100
|
- the idea would be to adjust the ssxpfp cut to get better agreement with
|
|
101
|
MM distribution shape
|
|
102
|
- change to ssxpfp<0.02 instead of <0.04
|
|
103
|
- resulting Pythia MM distribution is smoother
|
|
104
|
- Note: the previous MM upper edge issue (Data and MC had different limist) was
|
|
105
|
traced to the diamond cut and fixed
|
|
106
|
- *NB* Garth: to confirm that applying a different cut to Pythia MC than to
|
|
107
|
data and other MC doesn't create new issues, we should look at the Data vs
|
|
108
|
MC comparison of variables other than MM
|
|
109
|
- from the MM fit of MC to data, find the MC scale factor
|
|
110
|
- apply this scale factor to plots of other variables
|
|
111
|
- make plots of subtracted data vs omega MC, and plots of data vs sum of
|
|
112
|
normalized MCs
|
|
113
|
|
|
114
|
|
|
115
|
Firday: Present
|
|
116
|
---------------
|
|
117
|
Regina - Garth Huber, Alicia Postuma, Nathan Heinrich, Nacer Hamdi,
|
|
118
|
Nermin Sadoun, Vijay Kumar
|
|
119
|
JMU - Gabriel Niculesu, Ionana Niculescu
|
|
120
|
Glasgow - Kathleen Ramage, Rachel Montgomery
|
|
121
|
CUA - Sameer Jain, Tanja Horn, Chi Kin Tam
|
|
122
|
Virginia - Richard Trotta
|
|
123
|
FIU - Pete Markowitz
|
|
124
|
|
|
125
|
|
|
126
|
Richard
|
|
127
|
-------
|
|
128
|
KaonLT High Q2 LT-sep
|
|
129
|
- Data/MC plots for Q2=3.0 W=2.32 and Q2=5.5, using exactly the same parameters
|
|
130
|
and functions as Q2=4.4
|
|
131
|
- ratios are surprisingly close to 1
|
|
132
|
- uncertainties now include background subtraction uncertainties, so are much
|
|
133
|
larger in some cases
|
|
134
|
- errors look quite large for larger Q2 settings, need to reduce t-binning and
|
|
135
|
change upper and lower t-ranges
|
|
136
|
- Q2=4.4: need to go from 4 to 2 or 3 t-bins, also 15 to 9 phi bins
|
|
137
|
- Q2=5.5: need to go from 2 to 1 t-bin, reduce phi bins too
|
|
138
|
|
|
139
|
- kinematic plots: Q2=4.4 HMS-delta has some disagreement in shape, the others
|
|
140
|
are fairly reasonable
|
|
141
|
|
|
142
|
To Do:
|
|
143
|
- need to deal with the diamond cut identified by Chi Kin
|
|
144
|
- need to address Q2=4.4 cross section discrepancy between CKT and RT
|
|
145
|
|
|
146
|
Vijay
|
|
147
|
-----
|
|
148
|
PionLT Low Q2 systematic studies
|
|
149
|
- changing CoinTime cut by +/-10%, last week found that low epsilon sigUNS was
|
|
150
|
much more sensitive than mid and high epsilon
|
|
151
|
- given this is sensitivity, Garth asked VK to look at this vs t
|
|
152
|
- obtained separated cross sections for the +/-10% cut variation
|
|
153
|
- Garth had asked to try with no iterations to model
|
|
154
|
- Q2=0.375: sigL changed by 18%
|
|
155
|
- *NB* Tanja: suggests to do 1-2 iterations to see if this sensitivity
|
|
156
|
remains
|
|
157
|
|
|
158
|
- shows some fits of PKT model to Q2=0.375,0.425 data
|
|
159
|
- model shows a turn over in sigL at lowest -t that is not exhibited by data
|
|
160
|
- we have some discusion about this
|
|
161
|
- more studies are clearly needed to assess the reliability of the model to
|
|
162
|
these kinematics
|
|
163
|
|
|
164
|
Nermin
|
|
165
|
------
|
|
166
|
PionLT LD2 setting up for PID studies
|
|
167
|
- created a run list of LD2 data for all 4 beam energies
|
|
168
|
- had some issues submitting iFarm runs, which has been fixed
|
|
169
|
- produced ROOT files were corruped, need to look again at production scripts
|
|
170
|
|
|
171
|
Kathleen
|
|
172
|
--------
|
|
173
|
PionLT setting up for data analysis
|
|
174
|
- looking at previous Heep data replay raw trees
|
|
175
|
- using Junaid's script: singyield_heep.py for Run 10600
|
|
176
|
- looking into why H_gtr_xptar,yptar are missing
|
|
177
|
- will compare W vs yield for SIMC and Coin Data
|
|
178
|
|
|
179
|
|
|
180
|
Next Meetings
|
|
181
|
------------------
|
|
182
|
- Thurs: Jan 22 @ 16:00 Eastern/15:00 Regina
|
|
183
|
- PionLT will go first
|
|
184
|
|
|
185
|
- Fri: Jan 23 @ 11:00 Eastern/10:00 Regina
|
|
186
|
- we will continue where we left off
|
|
187
|
|
|
188
|
|
|
189
|
|
|
190
|
|
|
191
|
|
|
192
|
|
|
193
|
|
|
194
|
|
|
195
|
|
|
196
|
|