Added a maximum threshold of 1000 km/s for Disberg & Mandel (2025) LOGNORMAL kicks, matching paper

Author: Ilya Mandel

src/BaseStar.cpp

@@ -3904,10 +3904,16 @@ double BaseStar::DrawSNKickMagnitude(const double p_Sigma,
 
         case KICK_MAGNITUDE_DISTRIBUTION::LOGNORMAL: {                                          // LOGNORMAL
             // Only draw Disberg & Mandel (2025) kicks for CCSN or PPISN (if they receive a kick)
-	    // use Maxwellians with sigma set in CalculateSNKickMagnitude() for other SN types
+            // use Maxwellians with sigma set in CalculateSNKickMagnitude() for other SN types
             SN_EVENT thisSNevent = utils::SNEventType(m_SupernovaDetails.events.current);       // current SN event
-            if (thisSNevent == SN_EVENT::CCSN || (thisSNevent == SN_EVENT::PPISN && OPTIONS->NatalKickForPPISN()))
-                kickMagnitude = gsl_cdf_lognormal_Pinv(p_Rand, DISBERG_MANDEL_MU, DISBERG_MANDEL_SIGMA);
+            if (thisSNevent == SN_EVENT::CCSN || (thisSNevent == SN_EVENT::PPISN && OPTIONS->NatalKickForPPISN())) {
+                kickMagnitude = 1001.0;
+                double rand   = p_Rand;                                                         // makes it possible to adjust if p_Rand is too low, to avoid getting stuck
+                while( utils::Compare(kickMagnitude, 1000.0) > 0 ) {                            // maximum kick of 1000 km/s, following Disberg & Mandel (2025)
+                    kickMagnitude = gsl_cdf_lognormal_Pinv(p_Rand, DISBERG_MANDEL_MU, DISBERG_MANDEL_SIGMA);
+                    rand          = min(rand + p_Rand + 0.0001, 1.0);
+                }
+            }
             else
                 kickMagnitude = DrawKickMagnitudeDistributionMaxwell(p_Sigma, p_Rand);
             } break;

src/changelog.h

@@ -1591,6 +1591,8 @@
 //  03.20.06   IM - June 25, 2025       - Enhancements:
 //                                          - The MAXWELLIAN NS CCSN kick changed from the Hobbs value of 265 km/s to 217 km/s based on 48 younger than 10 Myr pulsars with proper motions from Disberg & Mandel (2025) sample; corrects Hobbs+ 2005 missing Jacobian
 //                                          - Implemented a LOGNORMAL NS CCSN kick magnitude distribution based on Disberg & Mandel, 2025
+//  03.20.07   IM - June 25, 2025       - Enhancement:
+//                                          - Added a maximum threshold of 1000 km/s for Disberg & Mandel (2025) LOGNORMAL kicks, matching paper
 
 const std::string VERSION_STRING = "03.20.06";