Fit to the light curve in X-ray and optical/NIR Vr′Ri′IJHK_S bands including a total of five energy injections. Data before 500 s are excluded from the modelling, as well as the strong X-ray flare. The light curves have been offset by constant factors as given in the figure legend for better legibility. The Bg′ bands are affected by additional Lyman forest absorption and were not included in the fit. For the Bg′ bands, the model predictions have been multiplied by 0.50 and 0.55 respectively to match them with the data. The residuals are clearly improved in comparison to those plotted in Fig. 3, but short-timescale variations like the one at 0.7 days are still not fitted satisfactorily.