CLAS12 Double-Pion (N ππ) �Electroproduction Analysis
Krishna Neupane
USC, Columbia, SC
Physics Analysis: Double-Pion Channel Cross-Section
2
Exp. Data: - Using Run Group A, Fall2018 data with inbending configuration (pass2), 126 runs same as Inclusive analysis
- The beam energy is 10.6041 GeV
- 1.4 GeV < W< 2.15 GeV & 2.0 GeV2 < Q2 < 9.0 GeV2
Simulations: - TWOPEG event generator is used (gemc 5.4/5.10) (pass2)
Topology: missing pim topology (ep->e’p’pi+)
Done: Energy loss, Momentum, ΔP ✅
Acceptance Correction
Detector Efficiency from Experimental Data ✅
Bin centering Corrections ✅
Radiative Effects ✅
Background subtraction ✅
In Progress: Empty target subtraction
Next: Error Analysis (Systematic Uncertainties)
ep-->e’p’ π+π-
3
cut_mid_0.95 : {38.12, 38.12, 38.12, 38.12, 38.12, 38.12, };
cut_tight_0.9 : {41.88, 40.62, 41.88, 41.88, 41.88, 41.88, };
cut_loose_0.99 : {36.88, 33.12, 35.62, 35.62, 35.62, 35.62, };
Electron PID
4
Double-Pion Channel Cross-Section
Electron Scattering Cross-Section:
Hadronic Cross-Section:
Single Differential Cross-Section:
The goal of the analysis is to extract the virtual photoproduction cross section σv of the reaction γvp → p′π+π−. This virtual photoproduction cross section σv is five-fold differential and in the single-photon exchange approximation connected with the seven-fold differential electron scattering cross section σe via :
Nine 1D experimental yields for 1.65-W-1.70 GeV 3.5-Q2-4.2 GeV2
1) Mπ-π+, Mπ+p’, θπ-, φπ- and α (pπ-)(p’π+) (ie. απ-) [π−, π+, p′]
2) Mp’π+, Mπ+π-, θp’ , φp’ and α(pp’)(π+π-) (ie. αp’ ) [p′, π−, π+]
3) Mπ+π-, Mπ-p’ , θπ+, φπ+ and α (pπ+)(p’π-) (ie. απ+) [π−, π+, p′]
Binning:
Simulations
5
Generated no of events in each 2D, W-Q2 bins (reference 108 events per bin)
gemc 5.10
2-2.4: 2.6 x107
2.4-3.0: 3.8 x107
3-3.5: 3.2 x107
3.5-4.2: 4.3 x107
4.2:5.0: 4.9 x107
5.0-6.0: 6.1 x107
6.0-7.0: 6.1 x107
7.0-8.0: 6.0 x107
8.0-9.0: 5.5 x107
gemc 5.4 Q2 < 3.3 to 8.5 GeV2
3.5-4.2: 2.9 x107
4.2:5.0: 3.4 x107
5.0-6.0: 4.2 x107
6.0-7.0: 4.2 x107
7.0-8.0: 3.9 x107
6
Acceptance Correction Factor A
7
8
Background based on: MMSQ peaks for ep->e’p’pi+pi- topology
Simulations
Experiment
9
10
11
12
13
Background Subtraction
14
Exclusive topology
Missing pim topology
Experiment
Simulations
15
Bin Centering Corrections
w_mid (GeV) | q2_mid(GeV) | hadronic_mid (GeV) | Total (Sum) | average value | middle bin value | Had corr | Had. bc error | w_corr | w_bc_error | q2_corr | q2_bc_error | Total_corr | Total_err |
1.425 | 2.2 | 1.0925 | 215.51146 | 19.59195 | 18.95441 | 0.96746 | 0.00361 | 0.83775 | 0.00348 | 0.97047 | 0.00401 | 0.78655586 | 0.00546374 |
1.425 | 2.2 | 1.122 | 658.24889 | 59.84081 | 58.41135 | 0.97611 | 0.00283 | 0.91781 | 0.00264 | 0.97249 | 0.00278 | 0.87123776 | 0.00434321 |
1.425 | 2.2 | 1.152 | 1167.53088 | 106.13917 | 105.57159 | 0.99465 | 0.00262 | 1.00427 | 0.00262 | 0.97431 | 0.00254 | 0.97323549 | 0.0044109 |
1.425 | 2.2 | 1.1815 | 1761.56491 | 160.14226 | 161.98165 | 1.01149 | 0.00267 | 1.01436 | 0.00263 | 0.97337 | 0.00253 | 0.99869222 | 0.00451584 |
1.425 | 2.2 | 1.211 | 1937.80174 | 176.16379 | 179.93126 | 1.02139 | 0.00276 | 1.01781 | 0.00267 | 0.97192 | 0.00258 | 1.01038952 | 0.00465548 |
1.425 | 2.2 | 1.241 | 1276.81974 | 116.07452 | 116.64342 | 1.0049 | 0.00289 | 0.99656 | 0.00278 | 0.9714 | 0.00278 | 0.97280187 | 0.00478979 |
1.425 | 2.2 | 1.2705 | 457.43673 | 41.58516 | 36.83623 | 0.8858 | 0.00366 | 0.88691 | 0.00326 | 0.96832 | 0.00388 | 0.76073628 | 0.00519525 |
16
1.675 | 3.85 | 1.1105 | 295.8804 | 26.89822 | 26.10299 | 0.97044 | 0.00266 | 0.99269 | 0.00287 | 0.96723 |
1.675 | 3.85 | 1.176 | 1120.22391 | 101.83854 | 99.77551 | 0.97974 | 0.00202 | 0.99092 | 0.00208 | 0.96747 |
1.675 | 3.85 | 1.2415 | 1566.16702 | 142.37882 | 147.51986 | 1.03611 | 0.002 | 0.99367 | 0.00193 | 0.96476 |
1.675 | 3.85 | 1.3065 | 863.56068 | 78.50552 | 76.01157 | 0.96823 | 0.00186 | 0.98888 | 0.00191 | 0.9671 |
1.675 | 3.85 | 1.372 | 605.14949 | 55.01359 | 54.72749 | 0.9948 | 0.00203 | 0.99092 | 0.00202 | 0.9667 |
1.675 | 3.85 | 1.4375 | 459.05753 | 41.7325 | 42.96744 | 1.02959 | 0.00242 | 1.00453 | 0.00234 | 0.96562 |
1.675 | 3.85 | 1.5025 | 171.0792 | 15.55265 | 14.76409 | 0.9493 | 0.00326 | 0.94789 | 0.00325 | 0.96558 |
1.525 | 4.6 | 1.1 | 31.88747 | 2.89886 | 2.791 | 0.96279 | 0.00268 | 0.98622 | 0.00303 | 0.96895 |
1.525 | 4.6 | 1.144 | 137.22181 | 12.47471 | 11.76909 | 0.94344 | 0.00189 | 0.99714 | 0.00205 | 0.96633 |
1.525 | 4.6 | 1.188 | 350.19645 | 31.83604 | 32.49528 | 1.02071 | 0.00186 | 1.00617 | 0.00186 | 0.96731 |
1.525 | 4.6 | 1.232 | 419.38118 | 38.12556 | 38.9836 | 1.02251 | 0.00183 | 1.00791 | 0.0018 | 0.96468 |
1.525 | 4.6 | 1.276 | 281.50358 | 25.59123 | 25.35616 | 0.99081 | 0.00186 | 1.00532 | 0.00189 | 0.96642 |
1.525 | 4.6 | 1.3195 | 153.30083 | 13.93644 | 13.32792 | 0.95634 | 0.00209 | 0.99783 | 0.00217 | 0.96491 |
1.525 | 4.6 | 1.363 | 56.13451 | 5.10314 | 4.74391 | 0.92961 | 0.00293 | 0.9734 | 0.00289 | 0.96485 |
2.025 | 5.5 | 1.1355 | 269.44525 | 24.49502 | 23.08156 | 0.9423 | 0.00239 | 0.99929 | 0.00292 | 0.96427 |
2.025 | 5.5 | 1.251 | 470.19436 | 42.74494 | 43.35199 | 1.0142 | 0.00195 | 0.98636 | 0.00194 | 0.96808 |
2.025 | 5.5 | 1.3665 | 308.47078 | 28.0428 | 27.62633 | 0.98515 | 0.00174 | 1.00207 | 0.00177 | 0.96627 |
2.025 | 5.5 | 1.4815 | 278.22608 | 25.29328 | 25.39429 | 1.00399 | 0.00181 | 1.00419 | 0.00182 | 0.96828 |
2.025 | 5.5 | 1.597 | 307.06639 | 27.91513 | 28.13875 | 1.00801 | 0.00206 | 1.00374 | 0.00206 | 0.96774 |
2.025 | 5.5 | 1.7125 | 283.70257 | 25.79114 | 27.04574 | 1.04864 | 0.00253 | 1.00934 | 0.00243 | 0.9655 |
2.025 | 5.5 | 1.8275 | 78.02528 | 7.09321 | 5.59966 | 0.78944 | 0.00222 | 0.97962 | 0.00323 | 0.96534 |
w_mid (GeV) | q2_mid(GeV) | hadronic_mid (GeV) | Total (Sum) | average value | middle bin value | Had corr | Had. bc error | w_corr | w_bc_error | q2_corr | q2_bc_error | Total_corr | |
17
Bin-Centering Corrections:1060-W-1.65 GeV, 4.2-Q2-5.0 GeV2
18
CLAS6 and Twopeg cross-sections
19
a and b will be calculated from the fit
Twopeg Cross-sections
20
The fit gives us:
a = 0.31660
b = 1.18085
21
Extending the fit Up to Q2 = 30.0 GeV2
Hole Fractions
23
Efficiency cut
1.90-W-1.95 GeV 3.5-Q2-4.2 GeV2
24
Preliminary Results
25
Preliminary Results: Best and worse comparison
Nine 1D cross-sections yields for 1.85-W-1.90 GeV 2.0-Q2-2.4 GeV2
Nine 1D cross-sections yields for 1.60-W-1.55 GeV 2.0-Q2-2.4 GeV2
26
Preliminary Results
27
Preliminary Results
28
Preliminary Results
29
Preliminary Results
30
Preliminary Results
31
Preliminary Results
32
Preliminary Results
Tasks to be done
33
Conclusions
34
35
36
Two-Pion Channel Cross-Section
Electron scattering Cross-Section:
Total cross-section
Δ N = L σ
From
dN/dt = (dN/dt)beamNtargetσ / A
Theta pim is angle betn
Three mom of initial
Proton and three mom
Of final pim.
1) Mπ-π+, Mπ+p’, θπ-, φπ- and α (pπ-)(p’π+) (ie. απ-) [π−, π+, p′]
2) Mp’π+, Mπ+π-, θp’ , φp’ and α(pp’)(π+π-) (ie. αp’ ) [p′, π−, π+]
3) Mπ+π-, Mπ-p’ , θπ+, φπ+ and α (pπ+)(p’π-) (ie. απ+) [π−, π+, p′]
Binning:
Where:
ΔN, Q are no of two-pion events inside 7-differential bin and charge on faraday cup with full an empty target, A, E, R are correction factors, ΔW, ΔQ2 are kinematical bins, L is luminosity,
∆τ = ∆Mpπ+∆Mπ+π-∆(−cos(θπ-))∆φπ-∆απ- is an element of the hadronic 5-dimensional phase space
36