Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Get constant neutrino loss fraction #145

Draft
wants to merge 2 commits into
base: develop
Choose a base branch
from
Draft

Get constant neutrino loss fraction #145

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
01780d3
Get constant neutrino loss fraction
Nov 5, 2024
e684915
update
Nov 7, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 5 additions & 1 deletion artisoptions_kilonova_lte.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@

#include "constants.h"

constexpr int MPKTS = 15000;
constexpr int MPKTS = 100000;

constexpr auto GRID_TYPE = GridType::CARTESIAN3D;
constexpr int CUBOID_NCOORDGRID_X = 50;
Expand Down Expand Up @@ -151,5 +151,9 @@ constexpr auto PARTICLE_THERMALISATION_SCHEME = ThermalisationScheme::DETAILED;

constexpr auto GAMMA_THERMALISATION_SCHEME = ThermalisationScheme::DETAILED;

constexpr float BETAMINUS_ENERGY_GAMMA_SPLITUP = 0.45;

constexpr float BETAMINUS_ENERGY_ELECTRON_SPLITUP = 0.2;

// NOLINTEND(modernize*,misc-unused-parameters)
#endif // ARTISOPTIONS_H
6 changes: 6 additions & 0 deletions decay.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -999,6 +999,12 @@ void init_nuclides(const std::vector<int> &custom_zlist, const std::vector<int>
// Read in data for gamma ray lines and make a list of them in energy order.
gammapkt::init_gamma_data();

// manipulate betaminus decay splitup ratios here
for (int i = 0; i < static_cast<int>(nuclides.size()); i++) {
nuclides[i].endecay_gamma = BETAMINUS_ENERGY_GAMMA_SPLITUP * nuclides[i].endecay_q[DECAYTYPE_BETAMINUS];
nuclides[i].endecay_electron = BETAMINUS_ENERGY_ELECTRON_SPLITUP * nuclides[i].endecay_q[DECAYTYPE_BETAMINUS];
}

// TODO: generalise this to all included nuclides
printout("decayenergy(NI56), decayenergy(CO56), decayenergy_gamma(CO56): %g, %g, %g\n",
nucdecayenergytotal(28, 56) / MEV, nucdecayenergytotal(27, 56) / MEV, nucdecayenergygamma(27, 56) / MEV);
Expand Down
5 changes: 4 additions & 1 deletion gammapkt.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -236,7 +236,10 @@ void init_xcom_photoion_data() {
__host__ __device__ auto choose_gamma_ray(const int nucindex) -> double {
// Routine to choose which gamma ray line it'll be.

const double E_gamma = decay::nucdecayenergygamma(nucindex); // Average energy per gamma line of a decay
double E_gamma = 0.; // Average energy per gamma line of a decay
for (ptrdiff_t n = 0; n < std::ssize(gamma_spectra[nucindex]); n++) {
E_gamma += gamma_spectra[nucindex][n].probability * gamma_spectra[nucindex][n].energy;
}

const double zrand = rng_uniform();
double runtot = 0.;
Expand Down
Loading