1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Beam Emittance

  1. Nov 21, 2013 #1

    I want to write a very simple code to plot the beam emittance knowing x,y,z and vx,vy,vz of each particle.
    I have generated with a tracking code a very simple beam (cone distribution with 1° divergence) and I want to compare the emittance of this beam with the one calculated by a different code.
    I'm using matlab and this is the code I have written. It's not working properly, but I'm not able to spot the error. COuld you please help me?

    Code (Text):

    A= Vel ; %%%Vel is a matrix of n rows and columns: [x,y,z,vx,vy,vz]
    %%%%%   evaluate the momentum of the reference particle

    c=3e8; %[m/s] light speed
    e=1.602176e-19; %[C] electron charge

    E0=938.27; %[MeV] proton rest energy
    M0=E0/c^2; %[MeV/c^2] proton rest mass
    mkg=M0*10^6*e; %[kg] mass in kg e

    T=5; %[MeV] Kinetic energy (it is a monochromatic beam)

    Etot= E0+T; %Total Energy

    gamma = Etot/E0; %gamma relativistic

    beta= sqrt( 1 - (1/gamma^2)); %beta relativistic = v/c


    %%%%% momentum (not taking in account particle mass%%%%

    angX=sin(A(:,4)./V0); %%%angle in radians
    angX=angX.*10^3; %%% angle in mradians

    angY=sin(A(:,5)./V0); %%%angle in radians
    angY=angY.*10^3; %%% angle in mradians
  2. jcsd
  3. Nov 22, 2013 #2


    User Avatar
    2017 Award

    Staff: Mentor

    What is wrong?

    Did you check variables like v? Your protons are non-relativistic, you can check their speed with a non-relativistic calculation.
    What is A?
    The sin of something will never give you an angle. Do you mean the arcsin there?
  4. Nov 23, 2013 #3
    I check variables and their dimensions. I try with classic and relativistic formulas... what is wrong is the plot. Only the Vx-Vy plot seems to be ok.

    The sin was wrong, you are right, but actually nothing change as I hav very small angles...

    Attached Files:

  5. Nov 23, 2013 #4


    User Avatar
    2017 Award

    Staff: Mentor

    Looks like an odd distribution of your particles, unrelated to your code.
    Your code just scales everything in A.
  6. Nov 24, 2013 #5
    I don't know, the particles should be distributed around a circle.
    I have the right prifile in xy ad vx-vy, but in the x-vx I have a straight line.
    What I sent you is what I get in the focal point of a quadrupole
  7. Nov 24, 2013 #6


    User Avatar
    2017 Award

    Staff: Mentor

    I would not expect lines at all - you should have a 2D-distribution, otherwise your phase space volume is zero.
    There is nothing fundamentally wrong with a v_x,x-distribution that is very stretched. It just means you have to focus in that direction again.
  8. Nov 24, 2013 #7
    The v_x,x distribution is not stretched, it is a line in my plot... that's why I'm concerned...
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook