Turn ON radiative splitting for annihilation

- Split annihilation in flight and at rest events for primary and fat positrons.
nrc-cnrc · Aug 3, 2023 · 4214601 · 4214601
1 parent fe362fd
commit 4214601
Showing 1 changed file with 9 additions and 3 deletions.
diff --git a/HEN_HOUSE/egs++/ausgab_objects/egs_radiative_splitting/egs_radiative_splitting.h b/HEN_HOUSE/egs++/ausgab_objects/egs_radiative_splitting/egs_radiative_splitting.h
@@ -142,11 +142,17 @@ class EGS_RADIATIVE_SPLITTING_EXPORT EGS_RadiativeSplitting : public EGS_AusgabO
         bool is_phat = (app->top_p.wt - wthin) > epsilon;
         bool is_primary = app->top_p.latch == 0 ? true : false;
 
-        /* Split primary and fat electrons */
-        if ( iarg == EGS_Application::BeforeBrems && (is_primary || is_phat) ){
+        /* Split radiative events ONLY for primary and fat electrons */
+        if ( iarg == EGS_Application::BeforeBrems       ||
+             iarg == EGS_Application::BeforeAnnihFlight ||
+             iarg == EGS_Application::BeforeAnnihRest   && 
+             (is_primary || is_phat)                    ){
            app->setRadiativeSplitting(nsplit);
         }
-        else if ( iarg == EGS_Application::AfterBrems ){
+        /* Avoids higher order splitting of radiative events */        
+        else if ( iarg == EGS_Application::AfterBrems       ||
+                  iarg == EGS_Application::AfterAnnihFlight ||
+                  iarg == EGS_Application::AfterAnnihRest   ){
            app->setRadiativeSplitting(1);
            app->setLatch(app->getNpOld()+1,1);
         }