|
- Timestamp:
-
Dec 30, 2018, 1:06:54 PM (6 years ago)
- Author:
-
mosel
- Comment:
-
--
Legend:
- Unmodified
- Added
- Removed
- Modified
-
v13
|
v14
|
|
39 | 39 | |
40 | 40 | == How to compute cross sections from the perturbative weights == |
41 | | The output file `FinalEvents.dat` contains all the events generated. Per event all the four-momenta of final state particles are listed together with the incoming neutrino energy and various other useful properties (see documentation for `FinalEvents.dat`). As an example we consider here the calculation of the differential cross section dsigma/dE for the outgoing lepton. In `FinalEvents.dat` the lines with the particle number 902 contain all the muon kinematics as well as the perweight. In order to produce a spectrum one first has to bin the muon energies into energy bins. This binning process must preserve the connection between energy and perweight. Then all the perweights in a given energy bin are summed and divided by the bin width to obtain the differential cross section. If the GiBUU run used - for better statistics - a number of runs >1 at the same energies, then this number of runs has to be divided out to obtain the final differential cross section. |
| 41 | The output file `FinalEvents.dat` contains all the events generated. Per event all the four-momenta of final state particles are listed together with the incoming neutrino energy and various other useful properties (see documentation for `FinalEvents.dat`). |
| 42 | |
| 43 | As an example we consider here the calculation of the CC inclusive differential cross section dsigma/dE_mu for a neutrino-induced reaction on a nucleus; E_mu is the energy of the outgoing muon. In `FinalEvents.dat` the lines with the particle number 902 contain all the muon kinematics as well as the perweight. In order to produce a spectrum one first has to bin the muon energies into energy bins. This binning process must preserve the connection between energy and perweight. Then all the perweights in a given energy bin are summed and divided by the bin width to obtain the differential cross section. If the GiBUU run used - for better statistics - a number of runs >1 at the same energies, then this number of runs has to be divided out to obtain the final differential cross section. |
42 | 44 | |
43 | 45 | All cross sections in GiBUU, both the precomputed ones and the reconstructed ones, are given per nucleon. The units are 10^-38^ cm^2^ for neutrinos and 10^-33^ cm^2^ for electrons. |
|