Fraction of particles getting through?

[unscientific]

Homework Statement

An unstable ion of mass ##M##, energy ##E## emits a massless particle of energy ##E_\nu## at angle ##\theta##. In the rest frame of the ion, find ##E_\nu^*## and ##cos \theta^*##.

Ions are now accelerated to ##\gamma=100## and a detector with radius ##r=20m## is placed ##D=200 km## away coaxially. Show that ##\cos \theta^* = \frac{1-\gamma^2 \theta^2}{1+\gamma^2 \theta^2 - \frac{\theta^2}{2}}## where ##\theta \simeq \frac{r}{D}##. Assuming in rest frame of ion, neutrinos are emitted equally in all directions, find the fraction of neutrinos that get through the detector.

Homework Equations

The Attempt at a Solution

Part (a)
By boosting to CM frame, we have
[tex]E_\nu^* = \gamma E_\nu \left( 1 - \beta \cos \theta\right) [/tex]
[tex]\cos \theta^* = \frac{\cos \theta - \beta}{1-\beta \cos \theta}[/tex]

Part (b)
Expanding small angles, we have
[tex]\cos \theta^* = \frac{1-\gamma^2 \theta^2}{1+\gamma^2 \theta^2 - \frac{\theta^2}{2}}[/tex]

How do I find the fraction of particles getting through? I did this in the rest frame:
[tex]f = \frac{4\int_0^{\theta^*} \sin \theta d\theta \int_0^{\phi^*} d\phi}{4\pi} [/tex]
[tex]f = \frac{1}{\pi} \left(1 - cos \theta^*\right) \theta^*[/tex]

Do I divide this by ##\gamma## to account for time dilation to find rate of particles getting through?

 

[mfb]

In (a), I think everything is supposed to be expressed with the given quantities.

How did you get the expansion for small angles?

How do I find the fraction of particles getting through?

Compare the solid angle (ion rest frame) where the neutrino hits the detector to the full solid angle.
There is no rate to calculate.

 

[unscientific]

In (a), I think everything is supposed to be expressed with the given quantities.

How did you get the expansion for small angles?

Compare the solid angle (ion rest frame) where the neutrino hits the detector to the full solid angle.
There is no rate to calculate.

I expanded ##\beta(\gamma)## and ##\theta \simeq \frac{r}{D}##.

So I can do everything in the ion's rest frame? Then the fraction is simply ##f = \frac{4\int_0^{\theta^*} \sin \theta d\theta \int_0^{\phi^*} d\phi}{4\pi}##.

 

[mfb]

The ϕ integral is trivial, but the prefactor of 4 looks odd.

So I can do everything in the ion's rest frame?

Sure, that's where you have the isotropic emission.

 

[unscientific]

The ϕ integral is trivial, but the prefactor of 4 looks odd.
Sure, that's where you have the isotropic emission.

I used ##\int_{-\theta^*}^{\theta^*} = 2 \int_0^{\theta^*}##.

 

[mfb]

θ does not get negative at all, or I misunderstood your coordinate system.
(Even then, where does the other 2 come from?)

 

[unscientific]

θ does not get negative at all, or I misunderstood your coordinate system.
(Even then, where does the other 2 come from?)

The cylindrical detector is on axis with the beam. ##\theta \simeq \frac{r}{D}##. So I think it should be ##\int_{-\phi^*}^{\phi^*} d\phi \int_0^{\theta^*} d\theta##

 

[mfb]

You need a really weird coordinate system to get negative θ. "The other side" should differ by pi in the angle ϕ, but not in θ.

 

[unscientific]

You need a really weird coordinate system to get negative θ. "The other side" should differ by pi in the angle ϕ, but not in θ.

Here's how I understand it: ##\theta## is the angle from the centre of the 'square' to the edge, so the area of the path is ##(r \sin \theta d\phi) \times (r d\theta)##. If I do ##\int_0^\theta d\theta \int_0^\phi d\phi##, it only covers 1/4 of the patch?

 

[mfb]

There is no square, the detector is a disk. With usual conventions for spherical coordinates, the angle goes from 0 to θmax or from -pi/2 to -pi/2+θmax, where the lower value is exactly in the center of the disk.
What do you mean with "area of the path"?

 

[unscientific]

There is no square, the detector is a disk. With usual conventions for spherical coordinates, the angle goes from 0 to θmax or from -pi/2 to -pi/2+θmax, where the lower value is exactly in the center of the disk.
What do you mean with "area of the path"?

I meant the area of the shaded patch. My worry is that the angle ##\theta## is defined as the angle from the centre of the patch to the edge of the patch..and using limits (0,theta) and (0,phi) only covers 1/4 of the patch.

 

[mfb]

My worry is that the angle θ\theta is defined as the angle from the centre of the patch to the edge of the patch.

That is the only reasonable definition.

and using limits (0,theta) and (0,phi) only covers 1/4 of the patch.

It covers 1/2, and if you integrate phi from 0 to 2 pi you cover everything.

 

[unscientific]

That is the only reasonable definition.
It covers 1/2, and if you integrate phi from 0 to 2 pi you cover everything.

Surely the limits of ##\phi## can't be from to 2pi. Ok, so the limits for the fraction are: (-phi,phi) and (0,theta).

 

[mfb]

-pi to pi is possible as well, it depends on the choice of your coordinate system which you still did not reveal. The result of the integral has to be 2 pi, the definition of the angle is not so important.

 

[unscientific]

-pi to pi is possible as well, it depends on the choice of your coordinate system which you still did not reveal. The result of the integral has to be 2 pi, the definition of the angle is not so important.

I suppose that the origin is where the particles fly out, and the detector is situated ##200km## away.

 

[mfb]

Yes that is obvious, but how do you define the angles phi and theta?

 

[unscientific]

Yes that is obvious, but how do you define the angles phi and theta?

Theta is the angle subtended from point of emission to edge of detector and phi is the same, but sideways.

 

[mfb]

That is a very unconventional and problematic coordinate system that leads to ugly integrals.
Spherical coordinates are so much easier.

 

[unscientific]

That is a very unconventional and problematic coordinate system that leads to ugly integrals.
Spherical coordinates are so much easier.

Ok, how would the integrals with limits look like in the spherical coordinate system?

 

[mfb]

Theta is the angle to the central axis (polar angle), phi is the azimuthal angle. Phi runs from 0 to 2 pi, theta from 0 to the maximal theta.