Problem Set 2

Problem 1

Draw the Dalitz plot in the same manner as problem 2 in Problem Set 1:

Minimum energy = m1 (momentum is conserved by m2 and m3 having equal and opposite momenta).  Maximum energy occurs when m2 and m3 are moving opposite m1 and have the minimum effective mass between them (ie. they are moving with the same velocity).  Then they act as a single particle with effective mass of m2+m3

eq1 = m0 == (Sqrt[(p^2) + (m1^2)]) + (Sqrt[(p^2) + (((m2 + m3))^2)])

solve1 = Simplify[Solve[eq1, p]]

Get the maximum energy:

Em = 
      PowerExpand[
        Simplify[
          (Sqrt[(p^2) + (m1^2)]) /. 
            solve1[[1]]]]

Choose eta' decay:

particles = {m0 -> 1270. , m1 -> 139.6, m2 -> 139.6, 
        m3 -> 497.7}

Maximum energy:

Emax = Em /. particles

Emin = m1 /. particles

Given an E1 between Emin and Emax, determine E2 as a function of Cos[theta]:

Clear[E1]

eq1 = m0 == E1 + E2 + E3

eq2 = eq1 /. E3 -> (Sqrt[(p3^2) + (m3^2)])

eq3 = eq2 /. p3 -> (Sqrt[(p1^2) + (p2^2) + 2 p1 p2 costh])

eq4 = eq3 /. p1 -> (Sqrt[(E1^2) - (m1^2)])

eq5 = Simplify[eq4 /. p2 -> (Sqrt[(E2^2) - (m2^2)])]

Now solve for E2 in terms of E1 and Cos[theta]:

solve2 = Solve[eq5, E2]

The two solutions correspond to the uper and lower branches of the curve if Cos[theta]=+/- 1.

Plot lower branch:

plot1 = 
      Plot[Evaluate[
          E2 /. solve2[[1]] /. 
              particles /. costh -> 1], {E1, Emin, Emax}, 
        PlotRange -> {{100, 650}, {100, 650}}]

Plot upper branch:

plot2 = 
      Plot[Evaluate[
          E2 /. solve2[[2]] /. 
              particles /. costh -> 1], {E1, Emin, Emax}, 
        PlotRange -> {{100, 650}, {100, 650}}]

Plot both and energy conservation triangle:

dalitzplot = 
      Show[plot1, plot2, 
        Graphics[
          Line[{{m1, m2}, {m0 - m3 - m2, m2}, {m1, m0 - m3 - m1}, {m1, 
                  m2}} /. particles]], AspectRatio -> 1, 
        AxesLabel -> {"Pi- Energy", "Pi+ Energy"}]

Calculate pi- energy in two body decay with K892:

s892 = 
      Solve[m0 == (Sqrt[(p^2) + (892. ^2)]) + (Sqrt[(p^2) + (m1^2)]) /. particles, p]

E892 = 
      (Sqrt[(p^2) + (m1^2)]) /. 
          s892[[1]] /. particles

Calculate pi energy limits:

E2min = 
      E2 /. solve2[[1]] /. 
            particles /. costh -> 1 /. E1 -> E892

E2max = 
      E2 /. solve2[[2]] /. 
            particles /. costh -> 1 /. E1 -> E892

Show[dalitzplot, Graphics[Line[{{E892, E2min}, {E892, E2max}}]]]

Projections:
	Pi- Plot = Delta function at 329.4
	Pi+ Plot = Flat distribution from 239.8 to 431.
	K0 Plot = Flat distribution from 509.6 to 700.8