97{
98 if(pReaction->GetReactionType() != fRectionType)
99 {
100 std::vector<MolType> molVector;
101 molVector.push_back(pReaction->GetReactant1());
102 molVector.push_back(pReaction->GetReactant2());
103 const G4int nbProducts = pReaction->GetNbProducts();
104 if (nbProducts) {
105 for (
G4int j = 0; j < nbProducts; ++j) {
106 auto product = pReaction->GetProduct(j);
107 molVector.push_back(product);
108 }
109 }
110 for(const auto& it : molVector)
111 {
112 if(it == fReactant1 || it == fReactant2 )
113 {
114 fAddEquilibrium = true;
115 fEquilibriumTime = fGlobalTime;
116 if(fVerbose >1) {
117 G4cout <<
"Reaction type : " << pReaction->GetReactionType() <<
" : "
118 << pReaction->GetReactant1()->GetName() << " + "
119 << pReaction->GetReactant2()->GetName() <<
G4endl;
120 G4cout <<
"SetEquilibrium : on " << fRectionType <<
" fEquilibriumTime : "
122 }
123 break;
124 }
125 }
126 }
127}