We present direct measurements of cubic bias parameters of dark matter halos from the halo-matter-matter-matter trispectrum. We measure this statistic efficiently by cross-correlating the halo field measured in N-body simulations with specific third-order nonlocal transformations of the initial density field in the same simulation. Together with the recent Ref.[1], these are the first measurements of halo bias using the four-point function that have been reported to date. We also obtain constraints on the quadratic bias parameters. For all individual cubic parameters involving the tidal field K ij , we find broad consistency with the prediction of the Lagrangian local-in-matter-density ansatz, with some indications of a positive Lagrangian coefficient b L td multiplying the time derivative of K ij . For the quadratic tidal bias (b K 2 ), we obtain a significant detection of a negative Lagrangian tidal bias.as well, as we will see). The parameter b K 2 has also been measured from the tree-level bispectrum in [10,16,27,28], the Lagrangian bispectrum [29], the 3-point function [30], and from Lagrangian moments-based measurements [23,31]. Some disagreement has been found between [16] and [31] in the results for b K 2 . Finally, the parameter b 3nl , a certain combination of quadratic and cubic tidal biases, has also been constrained from the 1-loop halo-matter power spectrum in [16], but, following our discussion, it is degenerate with the higher-derivative bias, which was set to zero in that reference.The goal of the present paper is to measure all the cubic bias terms, which are: b 3 = 6b δ 3 , b td , b K 3 , and b δK 2 . The leading statistic to which these contribute is the four-point function (trispectrum). In order to measure them, we generalize a technique proposed by Ref.[26] to measure the relevant trispectrum contributions efficiently. This technique allows us to measure all the cubic bias parameters at once. Together with the recent Ref.[1], these are the first measurements of halo bias using the four-point function that have been reported to date. We further use the analogous technique for the bispectrum to measure b 1 , b 2 and b K 2 , to cross-check our results obtained from the trispectrum, and to compare with previous measurements.The very similar study of [1] came out shortly after this paper appeared on the arXiv preprint server. While they use the same technique to obtain cubic order bias parameters from the trispectrum, some details such as higher-order corrections are treated differently. Overall, our results are in good agreement with theirs.This paper is organised as follows: in section 2 we present our estimator for the trispectrum and show how to obtain the bias parameters from it. Section 3 describes our set of simulations, shortly explains the halo finding procedure, and gives details on the actual procedure to measure the bias parameters. Section 4 reviews previous measurements and theoretical predictions for the parameters. Finally, section 5 presents and discusses our results. We conclu...