1 of 36

CLAS12 Double-Pion (N ππ) �Electroproduction Analysis

Krishna Neupane

USC, Columbia, SC

2 of 36

Physics Analysis: Double-Pion Channel Cross-Section

2

  • Particle Identification:
  • - Electron pid cuts ✅
  • - Hadron pid cuts ✅
  • Event Selection
  • Detector Inefficiency Cuts
  • Cross-Section Calculations
  • Corrections
  • Simulations ☑️
  • Holes Filling

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 of 36

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 of 36

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:

  • 50 MeV/c W bin size
  • 9 Q2 bins
  • 7 bins for invariant masses
  • 10 bins for θ
  • 6 bins for φ
  • 10 bins for α

5 of 36

Simulations

  • Twopeg event generator available on Jlab’s osg portal (Dr. Iuliia skorodumina)
  • 1.35 GeV < W< 2.15 GeV & 1.95 GeV2 < Q2 < 9.0 GeV2
  • 45nA background merging files
  • (gemc 5.10) (pass2) + gemc 5.4

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 of 36

6

Acceptance Correction Factor A

7 of 36

7

8 of 36

8

Background based on: MMSQ peaks for ep->e’p’pi+pi- topology

Simulations

Experiment

9 of 36

9

 

 

10 of 36

10

11 of 36

11

 

 

12 of 36

12

13 of 36

13

Background Subtraction

14 of 36

14

Exclusive topology

Missing pim topology

Experiment

Simulations

15 of 36

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 of 36

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 of 36

17

Bin-Centering Corrections:1060-W-1.65 GeV, 4.2-Q2-5.0 GeV2

18 of 36

18

CLAS6 and Twopeg cross-sections

19 of 36

19

 

 

 

a and b will be calculated from the fit

Twopeg Cross-sections

20 of 36

20

The fit gives us:

a = 0.31660

b = 1.18085

 

21 of 36

21

Extending the fit Up to Q2 = 30.0 GeV2

22 of 36

Hole Fractions

23 of 36

23

Efficiency cut

1.90-W-1.95 GeV 3.5-Q2-4.2 GeV2

24 of 36

24

Preliminary Results

25 of 36

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 of 36

26

Preliminary Results

27 of 36

27

Preliminary Results

28 of 36

28

Preliminary Results

29 of 36

29

Preliminary Results

30 of 36

30

Preliminary Results

31 of 36

31

Preliminary Results

32 of 36

32

Preliminary Results

33 of 36

Tasks to be done

  • Systematic and statistical Error analysis
  • 45, 50, 55 nA beam current adjustments

33

34 of 36

Conclusions

  • Particle identification cuts have been implemented
  • Based on missing π- topology, two-pion events are selected
  • Using TWOPEG event generator, acceptance correction factor are applied and remaining holes are filled
  • Energy/momentum corrections, radiative effect corrections, background subtractions are applied
  • Smearing MC data, Detector efficiency studies are in done
  • Double-pion channel total and nine 1D preliminary Cross-sections are extracted
  • Cross-sections are compared with CLAS6 analysis for overlapping regions and they are matching within the error bar limits
  • Systematic studies are needed to quantify uncertainties and extract precise Cross-Sections

34

35 of 36

35

36 of 36

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:

  • 50 MeV/c W bin size
  • 8 Q2 bins
  • 6 bins for invariant masses
  • 10 bins for θ
  • 6 bins for φ
  • 8 bins for α

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