QEDbase

This is QEDbase.jl, a julia package which provides the general data structures for calculations in relativistic particle physics.

This package is part of the QuantumElectrodynamics.jl library. For the description of the interoperability with other packages of QuantumElectrodynamics.jl see docs.

• Generic interface for Lorentz vectors
• concrete implementations of static and mutable Lorentz-vector/four-momentum types
• general Dirac bi-spinors, its adjoint counterpart as well as Dirac matrices
• particle spinors, i.e. solutions of Dirac's equation in momentum space
• Dirac's gamma matrices

To install the current stable version of QEDbase.jl you may use the standard julia package manager within the julia REPL

julia> using Pkg

julia> Pkg.add("QEDbase")

or you use the Pkg prompt by hitting ] within the Julia REPL and then type

(@v1.10) pkg> add QEDbase

To install the locally downloaded package on Windows, change to the parent directory and type within the Pkg prompt

(@v1.10) pkg> add ./QEDbase.jl

One can define a static four momentum component wise:

julia> using QEDbase; using QEDcore

juila> mass = rand()*10

julia> px,py,pz = rand(3)

julia> E = sqrt(px^2 + py^2 + pz^2 + mass^2) # on-shell condition

julia> P = SFourMomentum(E,px,py,pz)

Such SFourMomentum behaves like a four element static array (with all the standard arithmetics), but with the dot product exchanged with the Minkowski product

julia> @assert dot(P,P) == P*P == getMass2(P) == P[1]^2 - sum(P[1:].^2)

Furthermore, the Lorentz-vector interface provides a lot of properties for such a SFourMomentum, e.g.

julia> @assert isapprox(getRapidity(mom), 0.5*log((E+pz)/(E-pz)))
julia> @assert isapprox(getPlus(mom), 0.5*(E+pz))
julia> @assert isapprox(getPerp(mom), px^2 + py^2)

and a lot more (see here for a complete list). There is also a mutable version of a four vector in QEDcore.jl, where the Lorentz-vector interface provides setters to different properties as well (see here for details).

After installation it might be necessary to check if everything works properly. For that you can run the unittests by typing within the julia REPL

julia> using Pkg
julia> using QEDbase
julia> Pkg.test("QEDbase")

...

Testing Running tests...
Test Summary: | Pass  Total
QEDbase.jl    |  468    468
     Testing QEDbase tests passed

If you see the last line, you can assume that QEDbase.jl works properly for you.

If you want to contribute to QEDbase.jl feel free to do so by opening an issue or send me a pull request. In order to keep the packages within QuantumElectrodynamics.jl coherent, consider visiting the general contribution guide of QuantumElectrodynamics.jl.

This work was partly funded by the Center for Advanced Systems Understanding (CASUS) that is financed by Germany’s Federal Ministry of Education and Research (BMBF) and by the Saxon Ministry for Science, Culture and Tourism (SMWK) with tax funds on the basis of the budget approved by the Saxon State Parliament.

The core code of the package QEDbase.jl is developed by a small team at the Center for Advanced Systems Understanding (CASUS), namely

• Uwe Hernandez Acosta (u.hernandez@hzdr.de)

This package would not be possible without many contributions done from the community as well. For that we want send big thanks to:

• my Mate supplier
• the guy who will show me how to include the most recent contributers on github

MIT © Uwe Hernandez Acosta