diff --git a/History.md b/History.md
--- a/History.md
+++ b/History.md
@@ -11,6 +11,9 @@
 ===
 ## R02-02-0X
 
+16 Dec 2022 John Back
+* D89: Added probabilities for B_c -> V pi+ and V pi+ pi0 EvtBcVHad modes.
+
 16 Nov 2022 Tom Latham
 * T123: Provide documention of how to contribute bug reports, feature requests, new/modified code
 
diff --git a/src/EvtGenModels/EvtBcVHad.cpp b/src/EvtGenModels/EvtBcVHad.cpp
--- a/src/EvtGenModels/EvtBcVHad.cpp
+++ b/src/EvtGenModels/EvtBcVHad.cpp
@@ -126,7 +126,33 @@
 
 void EvtBcVHad::initProbMax()
 {
-    if ( out_code == 3 ) {
+    if ( out_code == 1 ) {
+        if ( idVector == EvtPDL::getId( "J/psi" ).getId() && whichfit == 1 &&
+             getNDaug() == 2 )
+            setProbMax( 500. );
+        else if ( idVector == EvtPDL::getId( "J/psi" ).getId() &&
+                  whichfit == 2 && getNDaug() == 2 )
+            setProbMax( 300. );
+        else if ( idVector == EvtPDL::getId( "psi(2S)" ).getId() &&
+                  whichfit == 1 && getNDaug() == 2 )
+            setProbMax( 17. );
+        else if ( idVector == EvtPDL::getId( "psi(2S)" ).getId() &&
+                  whichfit == 2 && getNDaug() == 2 )
+            setProbMax( 40. );
+    } else if ( out_code == 2 ) {
+        if ( idVector == EvtPDL::getId( "J/psi" ).getId() && whichfit == 1 &&
+             getNDaug() == 3 )
+            setProbMax( 10950. );
+        else if ( idVector == EvtPDL::getId( "J/psi" ).getId() &&
+                  whichfit == 2 && getNDaug() == 3 )
+            setProbMax( 4200. );
+        else if ( idVector == EvtPDL::getId( "psi(2S)" ).getId() &&
+                  whichfit == 1 && getNDaug() == 3 )
+            setProbMax( 500. );
+        else if ( idVector == EvtPDL::getId( "psi(2S)" ).getId() &&
+                  whichfit == 2 && getNDaug() == 3 )
+            setProbMax( 700. );
+    } else if ( out_code == 3 ) {
         if ( idVector == EvtPDL::getId( "J/psi" ).getId() && whichfit == 1 &&
              getNDaug() == 4 )
             setProbMax( 42000. );
@@ -139,7 +165,6 @@
         else if ( idVector == EvtPDL::getId( "psi(2S)" ).getId() &&
                   whichfit == 2 && getNDaug() == 4 )
             setProbMax( 2600. );
-
     } else if ( out_code == 5 ) {
         if ( idVector == EvtPDL::getId( "J/psi" ).getId() && whichfit == 1 &&
              getNDaug() == 6 )