diff --git a/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C b/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
index 11b6bae0938b390856da1ec7ca6dfbbe1061f0ce..3de79763ccb07254a74da9dd63f3e879bee1949f 100644
--- a/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
+++ b/src/thermophysicalModels/reactionThermo/psiuReactionThermo/heheuPsiThermo.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2017 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 2017-2019 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
                             | Copyright (C) 2011-2016 OpenFOAM Foundation
@@ -611,7 +611,7 @@ Foam::heheuPsiThermo<BasicPsiThermo, MixtureType>::psib() const
     forAll(psibCells, celli)
     {
         psibCells[celli] =
-            this->cellReactants(celli).psi(pCells[celli], TbCells[celli]);
+            this->cellProducts(celli).psi(pCells[celli], TbCells[celli]);
     }
 
     volScalarField::Boundary& psibBf = psib.boundaryFieldRef();