The data are tab separated and can be used without any specific script. Explenations of the content are as follow:

Fig3:
Field-induced magnetic Bragg intensity at Q = (111), measured at T = 50 mK as a function of the magnetic field applied along the [112] direction.

Fig4
(a) and (b):
 A selection of unprocessed INS data taken at T = 0.05 K at Q = (111) and (110), respectively, in magnetic fields of 7, 9.5, 12, and 13.5 T applied along [112], presented as offset plot. Files (a_1) and (b_1) are the fits of the respective measured curves
(c):
Data collected at the highet magnetic field of 14.5 T at different Q vectors. File (c_1) contains the fits of the respective curves.
(d) and (e):
Magnetic field dependence of the inelastic peak positions extracted from the fits from the curves measured at Q = (111) and Q = (-110).

Fig5:
Dispersion of field-induced magnetic excitations in Ce3Pd20Si6 measured in a magnetic field of 14.5 T applied along the [112] axis. The color map is composed of energy scans such as those in Fig. 4(c), measured along high-symmetry directions following the polygonal path (0.4 0.4 0.4)–(111)–(110)–(-0.5 0.5 0). The fitted peak position are shown as data points Fig5_1; the lines are an empirical fit to the data (Fig5_1).

Important Note on how to make Fig.5: 
Fig5 is a 15x50 matrix with the color coded intensities. (NaN = Not a Number means the point was not measured to save time). Here the data presented can be reproduced as follows:
 
Y-spacing is 0.05meV (range is from 0 to 2.5meV), X-spacing is (2*pi/a)*sqrt(h^2 + k^2 + l^2) and range is (0.4,0.4,0.4)-(1,1,1)-(-1,1,0)-(-0.5,0.5,0), where h,k,l are the Miller indices. 

Fig5_1 contains the fitted peak positions and dispersion as well as the X and Y-spacings (r.l.u. and energy respectively).