@@ -10,12 +10,12 @@ A detailed description can be found [here](https://git.kern.phys.au.dk/oliskir/b
- Program that calculates ξ and Θ: [spincoeff.f](https://git.kern.phys.au.dk/oliskir/bedepe/blob/master/spincoeff.f)
# Compile
First you must install CERNLIB. In Ubuntu this is easily accomplished with the command sudo apt-get install cernlib.
Bedepe may then be compiled with the command gfortran bedepe.out bedepe.f -static -lpacklib -lkernlib -lmathlib.
This creates the executable file bedepe.out.
First you must install CERNLIB. In Ubuntu this is easily accomplished with the command ```sudo apt-get install cernlib```.
Bedepe may then be compiled with the command ```gfortran bedepe.out bedepe.f -static -lpacklib -lkernlib -lmathlib```.
This creates the executable file ```bedepe.out```.
# Output and input
The program simulates beta-delayed particle emission (neutron, proton, alpha). The output file (bedepe.dat) has seven columns with the following information:
The program simulates beta-delayed particle emission (neutron, proton, alpha). The output file (```bedepe.dat```) has seven columns with the following information:
1. Event number
2. Kinetic energy of beta
3. Kinetic energy of delayed particle
...
...
@@ -26,11 +26,11 @@ The program simulates beta-delayed particle emission (neutron, proton, alpha). T
Quantities are measured in the rest frame of the decaying nucleus. Angles are in degrees and energies are in keV. The maximum weight assigned to any event is given at the end of the output file.
The program also produced a histogram file (bedepe.his), which can be viewed in PAW, with the following histograms:
The program also produced a histogram file (```bedepe.his```), which can be viewed in PAW, with the following histograms:
- Histogram #10: Recoil energy
- Histograms #11-15: Column 2-6 as defined above
The following information must be provided in the input file (input0):
The following information must be provided in the input file (```input0``):
1. Number of events
2. Q-value
3. Particle-decay threshold in daughter nucleus
...
...
@@ -41,7 +41,8 @@ The following information must be provided in the input file (input0):
8. ξ
9. Θ
The example input file is for simulating beta-delayed alpha emission from 8B assuming 3.0 MeV excitation energy in 8Be. a3 is the so-called triple-correlation coefficient which equals -1 for pure Gamow-Teller decays, +1 for pure Fermi decays and takes on values between -1 and +1 for mixed decays. ξ and Θ are spin-dependent coefficients given by the these formulas.
The example input file is for simulating beta-delayed alpha emission from 8B assuming 3.0 MeV excitation energy in 8Be. a3 is the so-called triple-correlation coefficient which equals -1 for pure Gamow-Teller decays, +1 for pure Fermi decays and takes on values between -1 and +1 for mixed decays.
ξ and Θ are spin-dependent coefficients given by the [these formulas](https://git.kern.phys.au.dk/oliskir/bedepe/blob/master/spincoeff.pdf).
# Comments
Feel free to contact me with questions or comments. Already now, I can point to one way in which the simulation program may be improved, namely, by including Coulomb effects in the beta phase-space function.
