Solving Particles in a Tube: Neutrino Energy from Pion Mass, E & Angle

[stunner5000pt]

Homework Statement

One way to produce a (muon-) neutrino beam is to produce a beam of pions and direct them down an evacuated pipe letting the pions decay to a neutrino (i.e., Muon + neutrino) Assume the neutrino is massless, what is the expression for the energy of the neutrino as a function of the pion’s mass m, energy E, and of the angle of emission of the neutrino in the lab frame (i.e., the angle between the neutrino and pion directions in the lab)?
Assume that the particles are moving at relativistic speeds

2. The attempt at a solution
The attached diagram is not part of hte problem's given. The diagram is my own work. Is it wrong to assume that the pion is going to split in a way analogous to Compton scattering?
Until the diagram is approved here is the description. assume the right and upward to be positive. The pion moves to the right and splits into the muon which goes off at angle phi and the neutrino which goes off at angle theta.

the momentum and the energy of the pion is
p_{\pi} = \gamma_{\pi} m_{\pi} v_{\pi} (1)
E_{\pi} = \gamma m_{\pi} c^2 (2)

the momentum in the X and Y directions of the muon
p_{x} = \gamma_{\mu} m_{\mu} v_{\mu} \cos\phi (3)
p_{y} = \gamma_{\mu} m_{\mu} v_{\mu} \sin\phi (4)
Energy is
E_{\mu} = \gamma_{\mu} m_{\mu}c^2 (5)

The energy of the neutrino is E_{\nu}
the momentum in the X and Y directions of the neutrino
p_{x} = \frac{E_{\nu}}{c} \cos\theta (6)
p_{y} = -\frac{E_{\nu}}{c} \sin\theta (7)

so are my equations correct??

For conservation of momentum in the X direction
\gamma_{\pi} m_{\pi} v_{\pi} = \gamma_{\mu} m_{\mu} v_{\mu} \cos\phi + \frac{E_{\nu}}{c} \cos\theta

For conservation of momentum in the Y direction
\gamma_{\mu} m_{\mu} v_{\mu} \sin\phi = \frac{E_{\nu}}{c} \sin\theta
E_{\nu} = \frac{\gamma_{\mu} m_{\mu} v_{\mu}c}{sin\theta} \sin\phi

For conservation of energy
\gamma_{\pi} m_{\pi} c^2 = \gamma_{\mu} m_{\mu}c^2 + E_{\nu}
E_{\nu} = E_{\pi} - \gamma_{\mu} m_{\mu}c^2

But energy of the muon is E_{\mu}^2 = p_{\mu}^2 c^2 - m_{\mu}^2 c^4

E_{\nu} = E_{\pi}- (p_{\mu}^2 c^2+ m_{\mu}^2 c^4

using conservation of momentum we can get p_{\mu} = \frac{E_{\nu}\sin\theta}{c\sin\phi}

So i have the Energy of the neturino is terms of the eneryg of hte pion, the emission angle of the neutrino theta and the mass of the muon... How do i get rid of the mass of the muon and the emission angle of the muon... phi? I ve been trying for quite some time please help!

Thanks for any and all of your help!

 

[stunner5000pt]

I posted this separate because the equations above look really complicating when really they arent

For conservation of momentum in the X direction
\gamma_{\pi} m_{\pi} v_{\pi} = \gamma_{\mu} m_{\mu} v_{\mu} \cos\phi + \frac{E_{\nu}}{c} \cos\theta

For conservation of momentum in the Y direction
\gamma_{\mu} m_{\mu} v_{\mu} \sin\phi = \frac{E_{\nu}}{c} \sin\theta
E_{\nu} = \frac{\gamma_{\mu} m_{\mu} v_{\mu}c}{sin\theta} \sin\phi

For conservation of energy
\gamma_{\pi} m_{\pi} c^2 = \gamma_{\mu} m_{\mu}c^2 + E_{\nu}
E_{\nu} = E_{\pi} - \gamma_{\mu} m_{\mu}c^2

But energy of the muon is E_{\mu}^2 = p_{\mu}^2 c^2 - m_{\mu}^2 c^4

E_{\nu} = E_{\pi}- (p_{\mu}^2 c^2+ m_{\mu}^2 c^4

using conservation of momentum we can get p_{\mu} = \frac{E_{\nu}\sin\theta}{c\sin\phi}

So i have the Energy of the neturino is terms of the eneryg of hte pion, the emission angle of the neutrino theta and the mass of the muon... How do i get rid of the mass of the muon and the emission angle of the muon... phi? I ve been trying for quite some time please help!

 

[stunner5000pt]

OK i worked on it a bit more
here's what i have

E_{\nu} = E_{\pi} - E_{\mu}
E_{\nu} = E_{\pi} - p_{\mu}^2 c^2 + m_{\mu}^2 c^4

Conservation of energy in the X direction

p_{\pi} = p_{\nu}\cos\theta + p_{\mu}\cos\phi
\frac{E_{\pi}^2}{c^2} - m_{\pi}^2 c^2 = \frac{E_{nu}}{c}\cos\theta +\left(\frac{E_{\mu}^2}{c^2} - m_{\mu}^2 c^2\right) \cos\phi

Conservation of momentum in the Y direction
\frac{E_{\nu}}{c}\sin\theta = \left(\frac{E_{\mu}^2}{c^2}-m_{\mu}^2 c^2\right)\sin\phidividing the two above equations we get \tan\phi and we can solve for phi
One we know Phi we can solve for the momentum of the muon. Once we have that we can sub this into the energy conservation eqatuions above and we are done

i get this after some solving (also i make c=1, 'natural units')
E_{\nu} = E_{\pi} - E_{\nu}^2 \sin^2 \theta \left(\sin \left[\arctan\left(\frac{E_{\nu}\sin\theta}{E_{\pi}^2 - m_{\pi}^2 - E_{\nu}\cos\theta}}\right)\right]\right)^2 + m_{\mu}^2

But now i have the neutrino energy everywhere! Is there a simpler way of doing this?

 

[kdv]

OK i worked on it a bit more
here's what i have

E_{\nu} = E_{\pi} - E_{\mu}
E_{\nu} = E_{\pi} - p_{\mu}^2 c^2 + m_{\mu}^2 c^4

Conservation of energy in the X direction

p_{\pi} = p_{\nu}\cos\theta + p_{\mu}\cos\phi
\frac{E_{\pi}^2}{c^2} - m_{\pi}^2 c^2 = \frac{E_{nu}}{c}\cos\theta +\left(\frac{E_{\mu}^2}{c^2} - m_{\mu}^2 c^2\right) \cos\phi

Conservation of momentum in the Y direction
\frac{E_{\nu}}{c}\sin\theta = \left(\frac{E_{\mu}^2}{c^2}-m_{\mu}^2 c^2\right)\sin\phi


dividing the two above equations we get \tan\phi and we can solve for phi
One we know Phi we can solve for the momentum of the muon. Once we have that we can sub this into the energy conservation eqatuions above and we are done

i get this after some solving (also i make c=1, 'natural units')
E_{\nu} = E_{\pi} - E_{\nu}^2 \sin^2 \theta \left(\sin \left[\arctan\left(\frac{E_{\nu}\sin\theta}{E_{\pi}^2 - m_{\pi}^2 - E_{\nu}\cos\theta}}\right)\right]\right)^2 + m_{\mu}^2

But now i have the neutrino energy everywhere! Is there a simpler way of doing this?

Looks complicated.

I haven't checked your steps but have you thought about using four-momentum? soundsl ike it woudl be easy.

writing P_{muon} = P_{pion} - P_{\nu} and squaring both sides seems like it would directly give the result

 

[stunner5000pt]

Looks complicated.

I haven't checked your steps but have you thought about using four-momentum? soundsl ike it woudl be easy.

writing P_{muon} = P_{pion} - P_{\nu} and squaring both sides seems like it would directly give the result

unsure how to do this...
Is the 4 momentum like (E,px,py,pz)

The mometum of the pion is
P_{\pi}= (E_{\pi},p_{\pi},0,0)<br /> P_{\mu}= (E_{\mu},p_{\mu}\cosphi,p_{\mu}\sin\phi,0)&lt;br /&gt; P_{\mu}= (E_{\nuu},-p_{\nu}\cos\theta,-p_{\nu}\sin\phi,0)&amp;lt;br /&amp;gt; P_{muon} = P_{pion} - P_{\nu}&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; But wouldn&amp;amp;#039;t squaring the above give us&amp;lt;br /&amp;gt; P_{\mu}^2 = (P_{\pi} - P_{\nu})^2&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; the other temrs are fine except for the 2 P_{\pi}\cdot P_{\nu}&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; plus how will i get rid of the momentum (not 4 momentum) of the muon?

 

[kdv]

unsure how to do this...
Is the 4 momentum like (E,px,py,pz)

The mometum of the pion is
P_{\pi}= (E_{\pi},p_{\pi},0,0)<br /> P_{\mu}= (E_{\mu},p_{\mu}\cosphi,p_{\mu}\sin\phi,0)&lt;br /&gt; P_{\mu}= (E_{\nuu},-p_{\nu}\cos\theta,-p_{\nu}\sin\phi,0)&amp;lt;br /&amp;gt; P_{muon} = P_{pion} - P_{\nu}&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; But wouldn&amp;amp;#039;t squaring the above give us&amp;lt;br /&amp;gt; P_{\mu}^2 = (P_{\pi} - P_{\nu})^2&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; the other temrs are fine except for the 2 P_{\pi}\cdot P_{\nu}&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; plus how will i get rid of the momentum (not 4 momentum) of the muon?

&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; The square of the four-momentum of any particle gives P^2 = m^2 (that&amp;amp;#039;s not true for so-called virtual particles but that&amp;amp;#039; irrelevant here). So P_{\mu}^2 = m_{\mu}^2&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; Also, you should write the components of the four-momentum of the neutrino explicitly. this is where the direction of the neutrino relative to the pion will enter. Note that you don&amp;amp;#039;t have to introduce a phi and a theta angle. You only need to introduce one angle by working in the plane containing the three-momenta of all three particles.

 

[stunner5000pt]

The square of the four-momentum of any particle gives P^2 = m^2 (that's not true for so-called virtual particles but that' irrelevant here). So P_{\mu}^2 = m_{\mu}^2

If this is true then doesn't P_{\pi}^2 = m_{pi}^2? So then how do we get the energy of the pion into our solution

Also, you should write the components of the four-momentum of the neutrino explicitly. this is where the direction of the neutrino relative to the pion will enter. Note that you don't have to introduce a phi and a theta angle. You only need to introduce one angle by working in the plane containing the three-momenta of all three particles.

P_{\pi}= (E_{\pi},p_{\pi},0,0)
P_{\mu}= (E_{\mu},p_{\mu}\cosphi,p_{\mu}\sin\phi,0)

neutrino's four momentum
P_{\nu}= (E_{\nu},-p_{\nu}\cos\theta,-p_{\nu}\sin\phi,0)

consevaton of 4 momentum
P_{\nu} = P_{pion} - P_{\mu}

do u mean like that?

 

[kdv]

If this is true then doesn't P_{\pi}^2 = m_{pi}^2? So then how do we get the energy of the pion into our solution



P_{\pi}= (E_{\pi},p_{\pi},0,0)
P_{\mu}= (E_{\mu},p_{\mu}\cosphi,p_{\mu}\sin\phi,0)

neutrino's four momentum
P_{\nu}= (E_{\nu},-p_{\nu}\cos\theta,-p_{\nu}\sin\phi,0)

consevaton of 4 momentum
P_{\nu} = P_{pion} - P_{\mu}

do u mean like that?

I am not sure what your theta and your phi are. But you may choose the axis such that the components of the neutrino four-momentum are simply P_{\nu}= (E_{\nu},p_{\nu}\sin\theta,p_{\nu}\cos\theta,0) with theta being the angle between the three momentum of the decaying pion and the three-momentum of the neutrino.

You will get the energy of the pion when you calculate the product P_{\pi} \cdot P_{\nu}

 

[kdv]

consevaton of 4 momentum
P_{\nu} = P_{pion} - P_{\mu}

do u mean like that?

Like I said in one of my previous posts, start instead from
P_{\mu} = P_{\pi} - P_{\nu} and square both sides. This way you won't have to deal with the energy of the muon which you don't want in your equation (right? You want the energy of the neutrino interms of the energy of the pion, correct?)

 

[kdv]

I am not sure what your theta and your phi are. But you may choose the axis such that the components of the neutrino four-momentum are simply P_{\nu}= (E_{\nu},p_{\nu}\sin\theta,p_{\nu}\cos\theta,0) with theta being the angle between the three momentum of the decaying pion and the three-momentum of the neutrino.

You will get the energy of the pion when you calculate the product P_{\pi} \cdot P_{\nu}

Sorry, I had the cos and sin theta switched in the neutrino four momentum It should have been

P_{\nu}= (E_{\nu},p_{\nu}\cos\theta,p_{\nu}\sin\theta,0)

 

[stunner5000pt]

I am not sure what your theta and your phi are. But you may choose the axis such that the components of the neutrino four-momentum are simply P_{\nu}= (E_{\nu},p_{\nu}\sin\theta,p_{\nu}\cos\theta,0) with theta being the angle between the three momentum of the decaying pion and the three-momentum of the neutrino.

You will get the energy of the pion when you calculate the product P_{\pi} \cdot P_{\nu}

i see your point now

Using the diagram i have posted in the first post i had it set up like this
P_{\nu}= (E_{\nu},p_{\nu}\cos\theta,-p_{\nu}\sin\theta,0)
P_{\pi} = (E_{\pi},p_{\pi},0,0}

does that make sense using the diagram up there?

OK so using
P_{\mu} = P_{\pi} - P_{\nu}

and squaring both sides

m_{\mu}^2 = m_{\pi}^2 - 2 P_{\pi}\cdot P_{\nu} + E_{\nu}^2

and using the 4 momentum of the muon and the neutrino
m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

now what is the square of the 4 momentum of the neutrino? If, for a particle iwht mass, the square of the 4 momentum is the mass of the particle, then for a massless particle it should be equal to the E/c?

and after this should i rewrite all the momenta of the muon and neutrino in terms of energy?

And

Thanks for your help so far !

 

[stunner5000pt]

Oh and by the way when you said

P^2 = m^2

this is using 'natural' units right? where h=c=1?

otherwise it would be
P^2 = m^2 c^4?

 

[kdv]

Oh and by the way when you said

P^2 = m^2

this is using 'natural' units right? where h=c=1?

otherwise it would be
P^2 = m^2 c^4?

Yes...Being a particle physicist, I always unconsciously set c=1

 

[kdv]

i see your point now

Using the diagram i have posted in the first post i had it set up like this
P_{\nu}= (E_{\nu},p_{\nu}\cos\theta,-p_{\nu}\sin\theta,0)
P_{\pi} = (E_{\pi},p_{\pi},0,0}

does that make sense using the diagram up there?

Yes. Sorry, I had not even noticed the diagram!

OK so using
P_{\mu} = P_{\pi} - P_{\nu}

and squaring both sides

m_{\mu}^2 = m_{\pi}^2 - 2 P_{\pi}\cdot P_{\nu} + E_{\nu}^2

and using the 4 momentum of the muon and the neutrino
m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

now what is the square of the 4 momentum of the neutrino? If, for a particle iwht mass, the square of the 4 momentum is the mass of the particle, then for a massless particle it should be equal to the E/c?

No, the square of the four-momentum of a massless particle is simply 0!

and after this should i rewrite all the momenta of the muon and neutrino in terms of energy?

Right, you have to reexpress the magnitude of the three-momenta (which is what you have in your final expression) in terms of the energies. For the neutrino it's simply p_{\nu} = E_{\nu}/c since it's massless. For the pion, it will have the square root of E^2 minus m^2 but this is where you could use the fact that the pion is relativistic to neglect its mass compared to its total energy so that, in the relativistic limit, p_{\pi} \approx E_{\pi}/c also.

Thanks for your help so far !

You are welcome.

 

[stunner5000pt]

alright then

m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

rewriting p_{\nu} = E_{\nu} where c=1 for 'natural' units

and p_{\pi} = \sqrt{E_{\pi}^2 - m_{\pi}^2}

m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+E_{\nu}\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)
m_{\mu}^2 = m_{\pi}^2 - 2E_{\nu} \left(E_{\pi}+\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)
and then rearranging we get

E_{\nu} = \frac{m_{\pi}^2 - m_{\mu}^2}{2\left(E_{\pi}+\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)}

good?

 

[kdv]

alright then

m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

rewriting p_{\nu} = E_{\nu} where c=1 for 'natural' units

and p_{\pi} = \sqrt{E_{\pi}^2 - m_{\pi}^2}

m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+E_{\nu}\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)
m_{\mu}^2 = m_{\pi}^2 - 2E_{\nu} \left(E_{\pi}+\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)
and then rearranging we get

E_{\nu} = \frac{m_{\pi}^2 - m_{\mu}^2}{2\left(E_{\pi}+\sqrt{E_{\pi}^2 - m_{\pi}^2}\cos\theta\right)}

good?

Looks good!

 

[stunner5000pt]

Looks good!

Thats awesome

thank you so much

 

[kdv]

Thats awesome

thank you so much

You are welcome.

Isn't that a nice formula? very compact and sweet. (Again, if you use the fact that teh pion is ultrarelativistic, it simplifies even more since you may neglect the mass compared to the total energy of the pion)

and you see how using four-momentum conservation is a powerful approach. You can recover the result using separately energy and three-momentum conservation but it's a very long and painful calculation.

I am glad I could help.

Take care.

 

[kdv]

alright then

m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}+p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

Oops. sorry, I had read too quickly. There is a sign mistake. It should read


m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}-p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

recall that the dot product between two four vectors involves a relative minus sign.

sorry if I missed it earlier, I was hurrying up because I was in the middel of helping students.

 

[stunner5000pt]

Oops. sorry, I had read too quickly. There is a sign mistake. It should read


m_{\mu}^2 = m_{\pi}^2 - 2 \left(E_{\pi}E_{\nu}-p_{\pi}p_{\nu}\cos\theta\right) + P_{\nu}^2

recall that the dot product between two four vectors involves a relative minus sign.

sorry if I missed it earlier, I was hurrying up because I was in the middel of helping students.

alright thanks for that
but why is there negative sign?

 

[kdv]

alright thanks for that
but why is there negative sign?

That's because of the metric. In general (in the flat spacetime of special relativity),

A_\mu B^\mu = A_0 B_0 - \vec{A} \cdot \vec{B}