diff --git a/CHANGELOG.md b/CHANGELOG.md index ceabd8b5b..4332cf998 100644 --- a/CHANGELOG.md +++ b/CHANGELOG.md @@ -9,6 +9,7 @@ ### Interface changes: ### Features that are deprecated and will be removed after this release + * EBL model from Finke et al. 2022 ### Removed features * AMRMagneticField - underlying library (saga) is no longer supported. diff --git a/CMakeLists.txt b/CMakeLists.txt index 0b0082f32..1ea3fe40f 100644 --- a/CMakeLists.txt +++ b/CMakeLists.txt @@ -308,7 +308,7 @@ endif(APPLE) # Download data files (interaction data, masses, decay data ...) # ---------------------------------------------------------------------------- OPTION(DOWNLOAD_DATA "Download CRPropa data files" ON) -set(CRPROPA_DATAFILE_VER "2023-10-20") +set(CRPROPA_DATAFILE_VER "2024-03-11") if(DOWNLOAD_DATA) message("-- Downloading data files from sciebo ~ 73 MB") file(DOWNLOAD diff --git a/include/crpropa/PhotonBackground.h b/include/crpropa/PhotonBackground.h index 7f7d60f8f..1048cda06 100644 --- a/include/crpropa/PhotonBackground.h +++ b/include/crpropa/PhotonBackground.h @@ -233,6 +233,19 @@ class IRB_Saldana21_lower: public TabularPhotonField { IRB_Saldana21_lower() : TabularPhotonField("IRB_Saldana21_lower", true) {} }; +/** + @class IRB_Finke22 + @brief Extragalactic background light model from Finke et al. 2022 + + Source info: + DOI:10.3847/1538-4357/ac9843 + https://iopscience.iop.org/article/10.3847/1538-4357/ac9843/pdf + */ +class IRB_Finke22: public TabularPhotonField { +public: + IRB_Finke22() : TabularPhotonField("IRB_Finke22", true) {} +}; + /** @class URB @brief Extragalactic background light model from Protheroe & Biermann 1996 diff --git a/test/testInteraction.cpp b/test/testInteraction.cpp index 19287c63a..f2fe49dbf 100644 --- a/test/testInteraction.cpp +++ b/test/testInteraction.cpp @@ -40,6 +40,8 @@ TEST(ElectronPairProduction, allBackgrounds) { epp.setPhotonField(irb); irb = new IRB_Stecker16_lower(); epp.setPhotonField(irb); + irb = new IRB_Finke22(); + epp.setPhotonField(irb); } TEST(ElectronPairProduction, energyDecreasing) { @@ -342,6 +344,8 @@ TEST(PhotoDisintegration, allBackgrounds) { pd.setPhotonField(irb); irb = new IRB_Stecker16_lower(); pd.setPhotonField(irb); + irb = new IRB_Finke22(); + pd.setPhotonField(irb); urb = new URB_Nitu21(); pd.setPhotonField(urb); } @@ -558,6 +562,8 @@ TEST(PhotoPionProduction, allBackgrounds) { ppp.setPhotonField(irb); irb = new IRB_Stecker16_lower(); ppp.setPhotonField(irb); + irb = new IRB_Finke22(); + ppp.setPhotonField(irb); ref_ptr urb = new URB_Protheroe96(); ppp.setPhotonField(urb); urb = new URB_Nitu21(); @@ -718,6 +724,8 @@ TEST(EMPairProduction, allBackgrounds) { em.setPhotonField(ebl); ebl = new IRB_Stecker16_lower(); em.setPhotonField(ebl); + ebl = new IRB_Finke22(); + em.setPhotonField(ebl); urb = new URB_Fixsen11(); em.setPhotonField(urb); urb = new URB_Nitu21(); @@ -821,6 +829,8 @@ TEST(EMDoublePairProduction, allBackgrounds) { em.setPhotonField(ebl); ebl = new IRB_Stecker16_lower(); em.setPhotonField(ebl); + ebl = new IRB_Finke22(); + em.setPhotonField(ebl); urb = new URB_Fixsen11(); em.setPhotonField(urb); urb = new URB_Nitu21(); @@ -925,6 +935,8 @@ TEST(EMTripletPairProduction, allBackgrounds) { em.setPhotonField(ebl); ebl = new IRB_Stecker16_lower(); em.setPhotonField(ebl); + ebl = new IRB_Finke22(); + em.setPhotonField(ebl); urb = new URB_Fixsen11(); em.setPhotonField(urb); urb = new URB_Nitu21(); @@ -1030,6 +1042,8 @@ TEST(EMInverseComptonScattering, allBackgrounds) { em.setPhotonField(ebl); ebl = new IRB_Stecker16_lower(); em.setPhotonField(ebl); + ebl = new IRB_Finke22(); + em.setPhotonField(ebl); urb = new URB_Fixsen11(); em.setPhotonField(urb); urb = new URB_Nitu21();