Automated calculation of matrix elements and physics motivated - - PowerPoint PPT Presentation

automated calculation of matrix elements and physics
SMART_READER_LITE
LIVE PREVIEW

Automated calculation of matrix elements and physics motivated - - PowerPoint PPT Presentation

1 Automated calculation of matrix elements and physics motivated observables Z. Was , Institute of Nuclear Physics, Polish Academy of Sciences Krakow (1) Once computers arrived, for me it was year 1980, approach to phenomenology of


slide-1
SLIDE 1

1

Automated calculation of matrix elements and physics motivated observables

  • Z. Was∗,

∗Institute of Nuclear Physics, Polish Academy of Sciences Krakow

  • (1) Once computers arrived, for me it was year 1980, approach to phenomenology of

theory/model based predictions could change a lot.

  • (2) Numerous benefits became available. Drawbacks appeared as well. For example,

methods of special functions expansions seem to be not as widespread as in the past.

  • (3) I will concentrate on examples of my personal experience. I do not have any intensions

to be systematic and balanced. Better picture will hopefully appear from other talks, e.g. examples of special functions expansions.

  • (4) I will not focus on successes of the field. These are well known.
  • (5) I will review traps which turned out to be rewarding to me once resolved; often in an

unexpected way.

  • Z. Was

Hayama, Octber, 2016

slide-2
SLIDE 2

2

  • Encouraged by Simizu-sensei conference, I choose to say what I

always wanted, but never did.

  • I thought the talk will be easy to prepare....
  • In contrary, I found work frustrating, but rewarding.
  • My plan is to show several simple examples of challenges

resulting from complexity and how automated calculations were of help, but also a source of difficulties.

  • Older examples originate from my work in Shimizu-san Minami

Tateya group I visited in 1995.

  • Each example in principle require substantial introduction,

impossible to cover in one talk.

  • My slides will show, outcome of my crippled attempts.
  • Z. Was

Hayama, Octber, 2016

slide-3
SLIDE 3

1983 Shoonship my first algebraic manipulation program 3

  • At that time Poland was an isolated place, but with enormous in-flow of students to
  • research. In reality a lot of contacts existed, but it was not to be seen by me.
  • Access to computing was limited and in fact quite awkward: hopeless loss of time it

seemed.

  • One of my first project was to evaluate spin density matrix for the process

e+e− → τ +τ −γ at Petra/PEP energies Monte Carlo Simulation of the Process e+ e- —> tau+ tau-

Including Radiative O(alpha**3) QED Corrections, Mass and Spin S. Jadach, Z. Was (Jagiellonian U.). Mar 1984. Comput.Phys.Commun. 36 (1985) 191.

  • Thhis work was performed under guidance of Prof. S. Jadach.
  • Fantastic experience in looking at spin amplitudes as (reducible) representations of

(Lorentz×gauge) groups.

  • It was great that we could spend all necessary time to understand details of what we

were doing.

  • In this particular case, how to represent moderatly complicated formulas of spin states

into compact forms, exploiting geometrical properties of formulae.

  • Z. Was

Hayama, Octber, 2016

slide-4
SLIDE 4

1983 Shoonship my first algebraic manipulation program 4

To simplify and to understand amplitudes:

Figure 2 RS( + ) B 3 (
  • )
  • !
QMS B 3 (
  • )
  • RS(
  • )
R 1 ( ) RS 1 (e + ) B 3 ( e )
  • !
CMS 1 B 3 ( e )
  • RS
1 (e
  • )
R 3 () R 3 () R 3 () RS(e + ) B 3 ( e )
  • !
CMS B 3 ( e )
  • RS(e
  • )
(2a) RS( + ) B 3 (
  • )
  • !
QMS B 3 (
  • )
  • RS(
  • )
R 1 ( 2 ) QMS
  • B
3 ( )R 3 () CMS
  • R
1 ( 1 ) RS 1 (e + ) B 3 ( e )
  • !
CMS 1 B 3 ( e )
  • RS
1 (e
  • )
R 3 ( 1 ) R 3 ( 1 ) R 3 ( 1 ) RS(e + ) B 3 ( e )
  • !
CMS B 3 ( e )
  • RS(e
  • )
(2b) 35
  • Z. Was

Hayama, Octber, 2016

slide-5
SLIDE 5

Compact and intitive representations of τ +τ − spin density matrix

5 forms of τ +τ − spin density matrix:

  • Z. Was

Hayama, Octber, 2016

slide-6
SLIDE 6

1994 Structure of spin amplitudes

6

  • General idea: to identify in amplitudes, with the help of gauge invariance structures

responsable later for phase-space enhancements: collinear-soft etc. This is fundamental, specially from the point of view of Monte Carlo algorithm construction.

  • Discussions with Shimizu-san were important.
  • Z. Was Gauge invariance, infrared / collinear singularities and tree level matrix element

for e+ e- —> nu(e) anti-nu(e) gamma gamma Eur.Phys.J. C44 (2005) 489,

  • A. van Hameren, Z. Was, Gauge invariant sub-structures of tree-level double-emission

exact QCD spin amplitudes, Eur.Phys.J. C61 (2009) 33

  • Also in this case algebraic manipulation mehods were providing the reference

calculations, necessary to cross check results.

  • I was not able to find patterns automatically, but algebraic progams were essential for

checks.

  • Only some of the patterns appear naturally. Feynman diagrams 1 and 2 combined (next

slide) are the complete amplitude for νµ¯

νµ production.

  • Z. Was

Hayama, Octber, 2016

slide-7
SLIDE 7

1994 Structure of spin amplitudes

7

Figure 1: The Feynman diagrams for e+e− → ¯

νeνeγ.

1 e
  • e
+ Z
  • i
  • i
  • 2
e
  • e
+ Z
  • i
  • i
3 e
  • e
W e +
  • e
4 e
  • e
W
  • e
e +
  • 5
e
  • e
W W
  • e
+
  • e
  • Z. Was

Hayama, Octber, 2016

slide-8
SLIDE 8

1994 Structure of spin amplitudes

8

  • The first two diagrams represent initial state QED bremsstrahlun amplitudes for νµ¯

νµ

pair production. It can be divided into parts, corresponding to β0, β1 of Yennie-Frautshi-Suura exponentiation.

  • Can separation be expanded to other cases, to higher orders, to terms of different

singularities/enhancements?

  • The answer seem to be always yes.
  • It is also important to observe that it extends to QCD, to scalar QED ...
  • I will sketch step for the calulation of single photon emission.
  • Slide 9 single photon emision in e+e− → νe¯

νe

  • Slide 10 double gluon emission in q¯

q → l+l−

  • Z. Was

Hayama, Octber, 2016

slide-9
SLIDE 9

1994 Structure of spin amplitudes

9

M1{I}

  • p

λ k1 σ1

  • = M0 + M1 + M2 + M3

M0 = eQe ¯ v(pb, λb) Mbd

{I}

pa + m − k1 −2k1pa ǫ⋆

σ1(k1) u(pa, λa)

+ eQe ¯ v(pb, λb) ǫ⋆

σ1 (k1) −pb + m+ k1

−2k1pb Mac

{I} u(pa, λa)

M1 = M1′ + M1′′ M1′ = +e ¯ v(pb, λb) Mbd,ac

{I}

u(pa, λa)ǫ⋆

σ1(k1) · (pc − pa)

1 ta − M2

W

1 tb − M2

W

, M1′′ = +e ¯ v(pb, λb) Mbd,ac

{I}

u(pa, λa)ǫ⋆

σ1(k1) · (pb − pd)

1 ta − M2

W

1 tb − M2

W

, M2 = +e ¯ v(pb, λb)gW eν

λb,λd ǫ⋆ σ1(k1) v(pd, λd)¯

u(pc, λc)gW eν

λc,λa k1 u(pa, λa)

1 ta − M2

W

1 tb − M2

W

M3 = −e ¯ v(pb, λb)gW eν

λb,λd k1 v(pd, λd)¯

u(pc, λc)gW eν

λc,λa ǫ⋆ σ1(k1) u(pa, λa)

1 ta − M2

W

1 tb − M2

W

,

(1)

  • Once manipulations completed, we separate the complete spin amplitude for the process e+e− → ¯

νeνeγ

into six individually QED gauge invariant parts. This conclusion is rather straightforward to check, replacing photon polarization vector with its four-momentum. Each of the obtained parts has well defined physical interpretation.

  • It is also easy to verify that the gauge invariance of each part can be preserved to the case of the extrapolation,

when because of additional photons, condition pa + pb = pc + pd + k1 is not valid.

  • Z. Was

Hayama, Octber, 2016

slide-10
SLIDE 10

QCD Eur.Phys.J. C61 (2009) 33

10

Ma,b = 1 2 ¯ v(p)

  • T aT bI(1,2) + T bT aI(2,1)

u(q) .

(2) For the T aT b-part, we find

I(1,2) = p·e1 p·k1 − k2·e1 k2·k1 − e /1k /1 2p·k1

  • J

/ k /2e /2 2q·k2 + k1·e2 k1·k2 − q·e2 q·k2

  • (3)

+ p·k2 p·k1 + p·k2 − k1·k2 p·e1 p·k1 − k2·e1 k2·k1 − e /1k /1 2p·k1 p·e2 p·k2 − k1·e2 k1·k2 − e /2k /2 2p·k2

  • J

/

(4)

+ J / q·k1 q·k1 + q·k2 − k1·k2 q·e1 q·k1 − k2·e1 k2·k1 − k /1e /1 2q·k1 q·e2 q·k2 − k1·e2 k1·k2 − k /2e /2 2q·k2

  • (5)

+ J /

  • 1 −

p·k2 p·k1 + p·k2 − k1·k2 − q·k1 q·k1 + q·k2 − k1·k2 k1·e2 k1·k2 k2·e1 k1·k2 − e1·e2 k1·k2

  • (6)

− 1 4 1 p·k1 + p·k2 − k1·k2 e /1k /1e /2k /2 − e /2k /2e /1k /1 k1·k2

  • J

/

(7)

− 1 4 J / 1 q·k1 + q·k2 − k1·k2 k /1e /1k /2e /2 − k /2e /2k /1e /1 k1·k2

  • .

(8) The part proportional to T bT a is obtained by a permutation of the momenta and polarization vectors of the gluons.

  • Z. Was

Hayama, Octber, 2016

slide-11
SLIDE 11

1998 Unexpected features,

11

  • The main purpose of my 1996 visit at KEK in MinamiTateya group, was to work on Grace

spin amplitudes (Comput.Phys.Commun. 153 (2003) 106).

  • Our KORALW Monte Carlo used Grace spin amplitudes for the e+e− → 4 fermion

processes.

  • Monte Carlo itegration of phase space regions where collinear configurations were

present, resulted in numerical difficulties. Abnormal features appeared. This required careful and painful work to avoid ‘trivial’ mistakes. Kind of faked ‘New Physics’, phenomnenon.

  • Let me show rather unexpected, at a time, example from the publication: Four quark final

state in W pair production: Case of signal and background, T. Ishikawa, Y. Kurihara, M. Skrzypek, Z. Was, Eur. Phys. J. C4 (1998) 75.

  • Interplay of theoretical effects and selection cuts can be confusing:
  • Z. Was

Hayama, Octber, 2016

slide-12
SLIDE 12

2→ 4 fermion processes at 195 GeV

12 W-pair production and decay, veto cut on 2 jets.

Figure 4: The

dσ2 dMs¯

s differential distribution of the “visible” s¯

s jets where c¯ c jets escape detection. The centre-of-mass energy is 195 GeV. Input parameters of type 2: CC-03 (thick line); and type 4: CC-43 (thin line). See Appendices A, B for a complete definition of all input parameters.

50 100 150 2.5 · 10−4 5.0 · 10−4 7.5 · 10−4 10.0 · 10−4 Ms¯

s [GeV]

dσ2 dMs¯

s[pb]

CC-43 CC-03

  • Z. Was

Hayama, Octber, 2016

slide-13
SLIDE 13

2→ 4 fermion processes at 195 GeV

13 effect due to spin: sometimes negigible sometimes not.

Figure 3: The

dσ2 dMs¯

s differential distribution of the “visible” s¯

s jets where c¯ c jets escape detection. The centre-of-mass energy is 195 GeV. Input parameters of type 1: CC-03 no spin correlation (thin line); and type 2: CC-03 spin correlations switched

  • n (thick line). See Appendices A, B for a complete definition of all input parameters.

50 100 150 .25 · 10−4 .50 · 10−4 .75 · 10−4 1.00 · 10−4 1.25 · 10−4 1.50 · 10−4 Ms¯

s [GeV]

dσ2 dMs¯

s[pb]

no-spin spin on

  • Z. Was

Hayama, Octber, 2016

slide-14
SLIDE 14

Higgs parity in H → ττ

14

  • I will show another example where complex observables need to be defined.
  • Important is slide nr 29: in case of H → ττ τ → 3πν we may want to mesaure

simultaneously 4 or 16 angles.

  • Each providing some of CP effect ...
  • ...but all of them correlated and under pressure from backgrounds.
  • These angles are extension of single acoplanarity angle which is used in case

H → τ +τ − τ ± → π±π0ν. Observable of mulidimensional nature cab be controlled

with ML techniques.

  • Risk of biasing.
  • I will skip some slides of introduction, we have no time to present that.
  • Z. Was

Hayama, Octber, 2016

slide-15
SLIDE 15

Higgs parity in H → ττ

15 The Higgs boson’s parity

  • H/A parity information can be extracted from the correlations between τ + and

τ − spin components which are further reflected in correlations between the τ

decay products in the plane transverse to the τ +τ − axes.

  • The decay probability

Γ(H/A → τ +τ −) ∼ 1 − sτ +

sτ −

  • ± sτ +

⊥ sτ − ⊥

is sensitive to the τ ± polarization vectors sτ − and sτ + (defined in their respective rest frames). The symbols ,⊥ denote components parallel/transverse to the Higgs boson momentum as seen from the respective

τ ± rest frames.

  • This spin case is technically easy, because ’Higgs spin’ is blind on Higgs origin.
  • Z. Was, M. Worek, Acta Phys. Polon. B33 (2002) 1875.
  • Z. Was

Hayama, Octber, 2016

slide-16
SLIDE 16

General formula for tau production and decay.

16 Formalism for τ +τ − : nothing changes

  • Because narrow τ width approximation can be obviously used for phase space,

cross-section for the process f ¯

f → τ +τ −Y ; τ + → X+¯ ν; τ − → νν reads: dσ =

  • spin

|M|2dΩ =

  • spin

|M|2dΩprod dΩτ + dΩτ −

  • This formalism is fine, but because of over 20 τ decay channels we have over

400 distinct processes. Also picture of production and decay are mixed.

  • Below only τ spin indices are explicitly written:

M =

2

  • λ1λ2=1

Mprod

λ1λ2 Mτ + λ1 Mτ − λ2

  • Cross section can be re-written into core formula of spin algorithms

dσ =

  • spin

|Mprod|2

spin

|Mτ +|2

spin

|Mτ −|2 wt dΩprod dΩτ + dΩτ −

  • Z. Was

Hayama, Octber, 2016

slide-17
SLIDE 17

αQED π

≃ 0.2% precision level

17 General formalism for semileptonic decays

  • Matrix element used in TAUOLA for semileptonic decay

τ(P, s) → ντ(N)X M =

G √ 2 ¯

u(N)γµ(v + aγ5)u(P)Jµ

  • Jµ the current depends on the momenta of all hadrons

|M|2 = G2 v2+a2

2

(ω + Hµsµ) ω = P µ(Πµ − γvaΠ5

µ)

Hµ =

1 M (M 2δν µ − PµP ν)(Π5 ν − γvaΠν)

Πµ = 2[(J∗ · N)Jµ + (J · N)J∗

µ − (J∗ · J)Nµ]

Π5µ = 2 Im ǫµνρσJ∗

ν JρNσ

γva = −

2va v2+a2

ˆ ω = 2 v2−a2

v2+a2 mνM(J∗ · J)

ˆ Hµ = −2 v2−a2

v2+a2 mν Im ǫµνρσJ∗ ν JρPσ

  • Z. Was

Hayama, Octber, 2016

slide-18
SLIDE 18

αQED π

≃ 0.2% precision level

18

H

iggs B

  • son Parity
  • Decay probability in formalism of Kramer et al.

Γ(H/A0 → τ +τ −) ∼ 1 − sτ +

sτ −

  • ± sτ +

⊥ sτ − ⊥

  • sτ is the τ polarization vectors.
  • / ⊥ denote components parallel / transverse to the Higgs boson momentum.
  • The spin weight is given by the following formula

wt = 1

4

  • 1 + 3

ij=1 Rijhihj

R33 = −1, R11 = ±1, R22 = ±1

  • Components for pure scalar and pseudoscalar Higgs boson respectively.
  • Z. Was

Hayama, Octber, 2016

slide-19
SLIDE 19

General idea

19 Density matrix

Only transverse spin correlations between τ + and τ − are different for scalar and pseudoscalar Higgs

  • The correlations can not be measured directly
  • One need to measure distributions of τ decay products
  • Precisely their transverse (to τ direction in Higgs boson rest frame) momenta
  • Most sensitive to spin is τ ± → π±ν
  • The largest branching ratio (25 % ) has τ ± → π±π0ν and we can look on

transverse spin correlations of ρ± → π±π0 decays.

  • Z. Was

Hayama, Octber, 2016

slide-20
SLIDE 20

Scalar or Pseudoscalar?

20

Pure Scalar And Pseudoscalar H

iggs B

  • son
  • Case of τ → ρντ decay, BR(τ → ρντ ) = 25%
  • The polarimeter vector is given by the formula where q for π± − π0, N for ντ .

hi = N

  • 2(q · N)qi − q2N i

q · N = (Eπ± − Eπ0)mτ

  • Acoplanarity of ρ+ and ρ− decay prod. (in ρ+ρ− r.f.) and events separation.

π π π π ρ ρ ϕ∗

− + − +

y1y2 > 0 ; y1y2 < 0 (in τ ± r.f.’s) y1 = Eπ+−Eπ0

Eπ++Eπ0 ;

y2 = Eπ−−Eπ0

Eπ−+Eπ0 .

  • Z. Was

Hayama, Octber, 2016

slide-21
SLIDE 21

Scalar or Pseudoscalar ?

21

R

esults Without Smearing

'
  • N
ev ts bin 0:5 1:0 1:5 2:0 2:5 3:0 0:0 '
  • N
ev ts bin 0:5 1:0 1:5 2:0 2:5 3:0 0:0
  • The ρ+ρ− decay products’ acoplanarity distribution without any smearing .
  • Selection y1y2 > 0 is used in the left plot, y1y2 < 0 is used for the right plot.
  • Thick line denote the case of the scalar Higgs and thin lines the pseudoscalar.
  • Complete spin correlations of h → τ +τ −, τ ± → ρ±ν, ρ± → π±π0 incl.
  • Z. Was

Hayama, Octber, 2016

slide-22
SLIDE 22

Scalar or Pseudoscalar ?

22

Phenomenology O

f General Case

  • Higgs boson Yukawa coupling expresed with the help of the

scalar–pseudoscalar mixing angle φ

¯ τN(cos φ + i sin φγ5)τ

  • Decay probability for the mixed scalar–pseudoscalar case

Γ(hmix → τ +τ −) ∼ 1 − sτ +

sτ −

  • + sτ +

⊥ R(2φ) sτ − ⊥

  • R(2φ) − operator for the rotation by angle 2φ around the direction.

R11 = R22 = cos 2φ R12 = −R21 = sin 2φ

  • Pure scalar case is reproduced for φ = 0.
  • For φ = π/2 we reproduce the pure pseudoscalar case.
  • Z. Was

Hayama, Octber, 2016

slide-23
SLIDE 23

Scalar or Pseudoscalar ?

23

O

bservable For M ixed Scalar−Pseudoscalar Case

  • For mixing angle φ, transverse component of τ + spin polarization vector is

correlated with the one of τ − rotated by angle 2φ.

  • Acoplanarity 0 < ϕ∗ < 2π is of physical interest, not just arc cos n− · n+.
  • Distinguish between the two cases 0 < ϕ∗ < π and 2π − ϕ∗
  • If no separation made the parity effect would wash itself out.
  • x

z y ρ+ π+

− −

∗ ρ π ϕ

Normal to planes: n± = pπ± × pπ0 Find the sign of

pπ− · n+

Negative

0 < ϕ∗ < π

Otherwise

2π − ϕ∗

  • Z. Was

Hayama, Octber, 2016

slide-24
SLIDE 24

Scalar or Pseudoscalar ?

24

R

esults For M ixed Scalar−Pseudoscalar Case

N ev ts bin '
  • 1:0
2:0 3:0 4:0 5:0 6:0 0:0

N ev ts bin '
  • 1:0
2:0 3:0 4:0 5:0 6:0 0:0

  • Only events where the signs of y1 and y2 are the same whether calculated

using the method without or with the help of the τ impact parameter.

  • Detector-like set-up is included (SIMDET).
  • The thick line corresponds to a scalar Higgs boson, the thin line to a mixed one.

Precision on φ ∼ 6 ◦, for 1ab−1 and 350 GeV CMS.

  • Z. Was

Hayama, Octber, 2016

slide-25
SLIDE 25

Hadronic currents: source of th. uncertainty

25

  • Improvements for ρ channel are technically straightforward: single real function to be

fitted: Jµ = (pπ± − pπ0)µFV (Q2) + (pπ± + pπ0)µFS(Q2) (FS ≃ 0).

  • For 3-scalar states: 4 complex function 3 variables each. Role of theoretical

assumptions is essential. Agreement on 1-dim distribution is a consistency check.

  • No go for model independent measurements? Not necessarily. Use of all dimensions for

data distributions: invariant masses Q2, s1, s2 as arguments of form-factors. Angular asymmetries help to separate currents: scalar Jµ

4 ∼ Qµ = (p1 + p2 + p3)µ, vector

1 ∼ (p1 − p3)µ|⊥Q and Jµ 2 ∼ (p2 − p3)µ|⊥Q and finally pseudovector

5 ∼ ǫ(µ, p1, p2, p3). Dependence on hadronic currents remain in calculation of

polarimetric vectors.

  • Model independent methods, template methods, neural networks, multidimensional
  • signatures. It was easier for Cleo. There, τ’s were produced nearly at rest, ντ

four-momentum was easy to reconstruct.

  • Fitting in complex situation is ... well complex ! Instead of acoplanarity angle in

a1 − a1 case we have 16 such angles.

  • Z. Was

Hayama, Octber, 2016

slide-26
SLIDE 26

ML techiques are needed.

26

from ρ± to a±

1 case.

  • 1. In case of τ → ρν there was one decay plane to define and sign of CP

sensitive sinusoid was dependent on sign of y+y−.

  • 2. In case of τ → a1ν four planes can be defined. Two for a1 → πρ0 and

another two for ρ0 → π+π− decays.

  • 3. We end up with 4 (or 16) angular distributions at number of yi like variables.
  • 4. That means meny sub categories to define sample ...
  • 5. All distributions are correlated.
  • 6. Methods of Maschine Learning necessary, to evaluate sensitivity of

mult-dimensional signatures → Brian.

  • Z. Was

Hayama, Octber, 2016

slide-27
SLIDE 27

Neural Network for CP parity of Higgs,arXiv:1608.02609 27

Acoplanarity angles of oriented half decay planes: ϕ∗

ρ0ρ0 (left), ϕ∗ a1ρ0 (middle) and ϕ∗ a1a1

(right), for events grouped by the sign of y+

ρ0y− ρ0, y+ a1y− ρ0 and y+ a1y− a1 respectively. Three

CP mixing angles φCP = 0.0 (scalar), 0.2 and 0.4. Note scale, effect on individual plot is so much smaller now. But up to 16 plots like that have to be measured, correlations

  • understood. But physics model depends on 1 parameter only and effect of φCP , the Higgs

mixing scalar pseudoscalar angle, is always a linear shift.

ρ ρ *

ϕ 1 2 3 4 5 6 Entries 120 125 130 135 140 145 150 155 160

3

10 ×

rest-frame ρ ρ

±

π 3 →

± 1

, a ν

± 1

a →

±

τ > 0

ρ

  • y

ρ +

y Scalar =0.2

CP

φ =0.4

CP

φ

1

a ρ *

ϕ 1 2 3 4 5 6 Entries 120 125 130 135 140 145 150 155 160

3

10 ×

rest-frame

1

a ρ

±

π 3 →

± 1

, a ν

± 1

a →

±

τ > 0

1

a

  • y

ρ +

y Scalar =0.2

CP

φ =0.4

CP

φ

  • 1

a

+ 1

a *

ϕ 1 2 3 4 5 6 Entries 120 125 130 135 140 145 150 155 160

3

10 ×

rest-frame

1

a

1

a

±

π 3 →

± 1

, a ν

± 1

a →

±

τ > 0

1

a

  • y

1

a +

y Scalar =0.2

CP

φ =0.4

CP

φ

  • Z. Was

Hayama, Octber, 2016

slide-28
SLIDE 28

Similar/supplementary projects.

28

Results relevant for fitting and for τ leptons.

  • 1. W production at LHC: lepton angular distributions and reference frames for probing hard

QCD, E. Richter-Was and Z. Was, arXiv:1609.02536

  • 2. Potential for optimizing Higgs boson CP measurement in H to tau tau decay at LHC and

ML techniques , R. Jozefowicz, E. Richter-Was and Z. Was,arXiv:1608.02609

  • 3. Separating electroweak and strong interactions in Drell−Yan processes at LHC: leptons

angular distributions and reference frames, E. Richter-Was and Z. Was, Eur.Phys.J. C76 (2016) 473

  • 4. “. Production of tau lepton pairs with high pT jets at the LHC and the TauSpinner

reweighting algorithm”, J. Kalinowski, W. Kotlarski, E. Richter-Was and Z. Was, arXiv:1604.00964

  • 5. “TauSpinner Program for Studies on Spin Effect in tau Production at the LHC”,
  • Z. Czyczula, T. Przedzinski and Z. Was, Eur. Phys. J. C 72, 1988 (2012)
  • Z. Was

Hayama, Octber, 2016

slide-29
SLIDE 29

Warning message

29

  • Biases in art, Giuseppe Arcimboldo (1572 - 1593).
  • Result depend on model assump-
  • tions. Models inspired with results ...

Fitting setup → biases.

  • Our algorithms are far less elaborate

than human eye/brain.

  • That may look worrisome.
  • Z. Was

Hayama, Octber, 2016

slide-30
SLIDE 30

We are not alone with the problem

30

Figure 2: Artificial Neural Networks have spurred remarkable recent progress in image classification and speech

recognition. But even though these are very useful tools based on well-known mathematical methods, we actually understand surprisingly little of why certain models work and others don’t. From http://googleresearch.blogspot.com/2015/06/inceptionism-going-deeper-into-neural.html Pattern recognition is an active field and deep concern and not only for us.

  • Z. Was

Hayama, Octber, 2016

slide-31
SLIDE 31

Summary

31

  • I have presented scattered results where use of computer algebraic methods or

pattern recognition techniques (Mashine Learning) was necessary.

  • My experience with such approaches started in 1996 in MinamiTateya group.
  • Working on my talk was inspiring to myself. Also, it was not easy to select slides

for a coherent presentation.

  • In fact, I am not sure if I was able to send the message: computer algebra

methods → correct huge expressions → loss of control on what should come

  • ut → how to understand/interpret/use → we are not alone with such

difficulties → how to avoid detection of non-existent...

  • Manpower/training is an essential issue for continuity of projects.
  • The challenges are more for newcomers, who may have missed long

years of rewarding failures.

  • Z. Was

Hayama, Octber, 2016