Project

General

Profile

Kaon LT Meetings » mtg_24jan25.txt

Garth Huber, 01/30/2024 02:25 PM

 
1
                 Jan 25/24 PionLT/KaonLT Analysis Meeting Notes
2
                 ----------------------------------------------
3
                          (Notes by GH & Richard Trotta)
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 - Ali Usman, Muhammad Junaid, Vijay Kumar, Nacer Hamdi, Garth Huber
13
JLab - Dave Gaskell
14
CSULA - Konrad Aniol
15
FIU - Pete Markowitz
16
CUA - Tanja Horn, Richard Trotta
17

    
18
Junaid
19
------
20
Efficiency studies
21
- CSV files of eff values created for all data
22
- working on scripts for lumi, heep, and finalizing plotting scripts
23
- Plotting scripts for efficiencies read in csv of efficiencies and plot from
24
  there
25
- updating REPORT FILES but Heep efficiencies are calculating incorrectly.  The
26
  issue is the error bar of efficiencies, not the efficiencies themselves.
27
  - Ali: there might be a bug in the TEMPLATE file, such as the error formula
28
    somehow dividing by zero.  If efficiencies themselves are okay, then DEF
29
    files ok.
30

    
31
Richard
32
-------
33
Issues replaying Dummy runs
34
- dtype error
35
  - python cut script gave error for some runs and not processed
36
  - fixed code to catch error when there is an empty ROOT tree.  Code now says
37
    which trees are empty
38
- CT cut issue
39
  - met w/Stephen on Jan 24
40
  - CT offset was not applied to Dummy runs.  SK missed a few Dummy runs in
41
    standard.kinematics
42
  - Vijay found a similar issue for some runs
43

    
44
Weird features in some t-distributions
45
- some (e,e'K+) runs had a spiky t-distrubtion
46
  - getting events lower than the -t_min for K+ in that setting, which is a
47
    sudden jump in distribution from the allowed region
48
  - most likely due to insufficient PID cut, both pi+ and proton leakthrough,
49
    some evidence of both pi+n and p+omega peaks
50
  - after tightening HGC cut, the t-distribution looks better
51
  - seems to have evidence of Lambda(1405) and other resonances in Q2=2.115
52
    center, low epsilon data
53
    - this is pleasant news, in the online analysis only the high epsilon data
54
      (w/o diamond cut) showed strong evidence of the higher hyperon states
55

    
56
- GH: suggests to generate *clean* pi+ and proton spectra, normalize them to the
57
  pi+n and p+omega peaks and *subtract* them to get an even cleaner K+
58
  distribution
59
  - absolutely must be finalized before starting iterations!
60
  - also critical to determine the correct PID efficiencies to apply to data
61

    
62
Fixed binning issue
63
- long-standing binning issue, some bins were always empty when binning over
64
  Q2-bar, W-bar, t, epsilon
65
- rewrote both binning scripts for average kinematics and yields
66
- now looks much better
67

    
68
Separated cross-section script
69
- scripts now working!
70
- now can debug things properly with average kinematics and yields
71
  - some issues were encountered in not treating t-bins in a consistent
72
    sequence
73
- now can generate Iteration-0 plots, then can try to do the next iteration
74
  - process in general works, but need to resolve PID issues first
75

    
76
- DG: your high epsilon data show a weird left-right phi asymmetry that is
77
  unphysical
78
  - need to look more closely at what is causing this
79

    
80
SIMC MM weirdness
81
- plot shown for Q2=3.0, W=3.14, low-eps, center
82
  - K+L MM distribution has a bizarrely long tail that is much bigger than the
83
    data, exactly the opposite situation than was encountered by Alicia, where
84
    the pi+n MM tail is under-predicted by SIMC
85
- GH: maybe it's a radiative tail flag that is wrong??
86

    
87
SIMC negative counts for some kinematics at highest Q2
88
- for Q2=5.5, the Iteration-0 script gives negative weight events
89
  - corresponds to sigT<0 for some kinematics
90
  - the initial fit parameters in the SIMC model are completely inappropriate
91
    for this Q2, and need adjustment
92
  - will set to zero initially
93
  - GH: if this doesn't work, then find by hand some reasonable fit parameter
94
    values to try initially
95

    
96
Ali
97
---
98
- N-Delta Transition GPDs white paper contribution for ECT* draft was sent to
99
  collaborators and when feedback is returned he will send to everyone
100

    
101
- Optimizing BSA code, looking at left and right settings
102
  - Q2=2.115, not every bin has enough statistics
103
  - If Events<Threshold cut applied, sets these bins to zero
104
    - the issue is that the SIDIS+rho fits failed if not enough events
105
    - GH: Bill had to do something similar in the omega analysis, as the
106
      rho+PhaseSpace fits would fail if not enough statistics, or if radiation
107
      corrections were too high.  This is all explained in his thesis if you
108
      want
109

    
110
  - weighted average of everything should cover empty bins for off-center
111
    settings
112
  - Made changes to weighted average scripts, should have plots for all
113
    settings combined soon
114

    
115
Dave
116
----
117
hsxpfp vs delta correction for high HMS momenta
118
- looked at a subset of runs provided by Ali
119
- applied correction: delta=delta+C0*hsxpfp
120
  - correction factor C0<0
121
- plot of C0 vs P_HMS shows a monotonic decrease to about 5 GeV/c, then
122
  is reduced in magnitude from 5 to 6.5 GeV/c
123
  - likely the reversal in trend is due to partial HMS saturation correction
124
    applied to Matrix Elements
125
- 6.59 GeV/c correlation is not totally linear, did a best fit linear
126
  correction, but some slight curvature remains
127
- Will provide table of numerical C0 values
128

    
129
- some discussion on how to implement this in hcana
130
  - if recalculate hcana values on the fly, then just recalc delta with above
131
    correction, and not make any modifications to hcana
132
  - Richard will implement into ltsep python framework to recalculate delta
133
  - fortunately, this does not affect the PionLT data, as the effect
134
    disappeared when the power supply was replaced
135

    
136
- Final step will to be to optimize SIMC resolution for pi+n to get better
137
  agreement between data and MC
138

    
139
Vijay
140
-----
141
Busy with making talk for conference in India
142
- EIC workshop in Bangalore
143
- Will send around draft slides in next few days
144
- Leaving for India on Feb 2.  Talk is Feb 9
145

    
146
Nacer
147
-----
148
SHMS DC calibration
149
- resolved all bugs for KaonLT calib w/Junaid
150
- new t0 offsets needed
151
- Need to compare parameterization every few runs
152
  - systematically checking groups of runs.  6638-6691 shown
153

    
154
- Pre-calibration X-T and new T0 (red) is usually better than post-calibration X-T
155
  - this is because it over-calibrates (i.e. over-fits the wiggles in the data)
156
    in particular when there are low stats
157

    
158
- Basically only need to calibrate t0 and that's it. This is just tedious
159
  because need to check how it varies 
160
  - set up to check once per 8hr shift of data
161
  - ~115 runs still to go, should be done around mid to late February
162
  - batch jobs failing and such take a while
163

    
164

    
165
Next Meeting
166
------------
167
- Thur Feb 1 @ 15:00 Eastern/14:00 Regina
168
  - KaonLT will go first
169

    
(362-362/513)