E = \sum_{p=1}^{N_{\rm prod}} \left[ R_O \cdot q^{\uparrow}_{O, p} + R_G \cdot  q^{\uparrow}_{G, p} \right] 
- \sum_{p=1}^{N_{\rm prod}} C_L  \cdot q^{\uparrow}_{L, p}
- \sum_{p=1}^{N_{\rm prod}} C_O  \cdot q^{\uparrow}_{O, p}
- \sum_{p=1}^{N_{\rm prod}} C_G  \cdot q^{\uparrow}_{G, p}
- \sum_{p=1}^{N_{\rm prod}} C_W  \cdot q^{\uparrow}_{W, p}
- \sum_{i=1}^{N_{W, \rm inj}} C_{W, \rm inj}  \cdot q^{\downarrow}_{W, i}
- \sum_{j=1}^{N_{G, \rm inj}} C_{G, \rm inj}  \cdot q^{\downarrow}_{G, i}