I put in the following values for the masses into the code in an attempt to recreate the parameters of MLB:
Particle 4: Mass = 170 GeV
Particle 3: Mass = 80 GeV
Particle 2: Mass = 0.0001 GeV
Particle A: Mass = 0.1 GeV
Particle B: Mass = 0.1GeV
To test the code I calculated the invariant mass of 2+A+B to make sure that it came out to 170 GeV as it should:
Pretty much spherically symmetric. Currently continuing the investigating on why the code isn't giving my triangles.
Particle 4: Mass = 170 GeV
Particle 3: Mass = 80 GeV
Particle 2: Mass = 0.0001 GeV
Particle A: Mass = 0.1 GeV
Particle B: Mass = 0.1GeV
To test the code I calculated the invariant mass of 2+A+B to make sure that it came out to 170 GeV as it should:
So I know the mass formulae that I derived for the 4-vectors are correct. Next I tried to make some triangles using the invariant mass of A+B:
Not quite a triangle unfortunately. The end point is in the correct place and I ran the code again using only 180 degree decay angles and it gave a single peak at 150 GeV. So at least I can be 100% confident that the Lorentz boost function is working properly.
I thought maybe there is a problem with the random number generator so I plotted the Theta/Phi distribution of 2 and B in the rest frame of 3:
I thought maybe there is a problem with the random number generator so I plotted the Theta/Phi distribution of 2 and B in the rest frame of 3:
No comments:
Post a Comment