Analyzing p K data OverComplete Measurements and partialwave

  • Slides: 37
Download presentation
Analyzing p(γ, K+)Λ data: (Over)Complete Measurements and partial-wave analyses

Analyzing p(γ, K+)Λ data: (Over)Complete Measurements and partial-wave analyses

4 topics •

4 topics •

The data

The data

Data: visualized NEW CLAS

Data: visualized NEW CLAS

New data: visualized

New data: visualized

CLAS overview

CLAS overview

Experiments

Experiments

Partial-wave analysis and amplitude extraction

Partial-wave analysis and amplitude extraction

Differential cross section

Differential cross section

Multipoles

Multipoles

L_max and ambiguities

L_max and ambiguities

Chiang and Tabakin Single-(s, t) amplitudes

Chiang and Tabakin Single-(s, t) amplitudes

Case for real data: entropy

Case for real data: entropy

Real data AND models

Real data AND models

Comparing data sets

Comparing data sets

More interesting: the future • Projects that are on my to-do list: • Tackle

More interesting: the future • Projects that are on my to-do list: • Tackle model-independent PWA in a Bayesian framework. (finished proof of principle) • Finite-energy Sum Rules in eta photoproduction. (main research) • Further down the road • Improving the Regge-plus-resonance model for electroproduction. • Extending the Vrancx-Ryckebusch models (above the resonance region).

Bayesian approach • Chi-square: point estimate of the parameters • What is the uncertainty

Bayesian approach • Chi-square: point estimate of the parameters • What is the uncertainty of your results? • NOT (!!!) Minuit error bars (see Sandorfi plots) • You can do bootstrapping • Bayesian framework: works with distributions. • Has been applied to similar problems and is well-accepted. • You get a very (very!) good idea of what the errors are. • You can properly determine the L_max from your data (evidence).

Probability theory •

Probability theory •

Observable level

Observable level

No overfitting or unnecessary model dependence + NO STARTING VALUES!

No overfitting or unnecessary model dependence + NO STARTING VALUES!

VERY first results for P. O. C. 1. Identification of outliers 2. Identification of

VERY first results for P. O. C. 1. Identification of outliers 2. Identification of ambiguities 3. Allows for smooth connections (continuous distributions!) Violin plot for Lmax = 1, Re E 0+ These distributions can be projected to observables Results with Multi. Nest: Gives you ALL maxima

About yesterday Make cuts here and see which multipoles would be affected

About yesterday Make cuts here and see which multipoles would be affected

Evidence chart

Evidence chart

Observables (only maxima)

Observables (only maxima)

Pros and cons • Pros • Posterior samples from which one can calculate •

Pros and cons • Pros • Posterior samples from which one can calculate • Marginals • Probabilities with multipoles with evidence • Distributions for the observables • Not just ‘visual’ inspection of proper Lmax, but Occam’s razer naturally built in. • You have a better idea of how significant your results are than with a single Minuit minimization. • Cons • Computational time (limits to Lmax = 5 if no phases are fixed). • Maybe I did not consider the proper reaction for this.

Collaboration with V. Mathieu Finite-Energy Sum Rules Unphysical Regge Low-energy Models Im s Re

Collaboration with V. Mathieu Finite-Energy Sum Rules Unphysical Regge Low-energy Models Im s Re s

Dispersive integral Nucleon pole Regge Unphysical domain (2) Physical domain

Dispersive integral Nucleon pole Regge Unphysical domain (2) Physical domain

Continuation of the multipoles

Continuation of the multipoles

Legendre polynomials

Legendre polynomials

Legendre polynomials L=0 J = 1/2

Legendre polynomials L=0 J = 1/2

Legendre polynomials

Legendre polynomials

Legendre polynomials

Legendre polynomials

2 CGLN amplitudes (t = -0, 5 Ge. V )

2 CGLN amplitudes (t = -0, 5 Ge. V )

Invariant amplitudes (t = -0, 5 Ge. V 2)

Invariant amplitudes (t = -0, 5 Ge. V 2)

Fin

Fin

Different than penalty terms Regularization or ‘ridge regression’

Different than penalty terms Regularization or ‘ridge regression’