This work is supported by NSFC No.11773028, 11633001, 11653002, 11603020, 11903030, 12003029, 11903033, the Fundamental Research Funds for the Central Universities under Grant No.WK2030000036, WK3440000004 and WK2030000044, the Strategic Priority Research Program of the Chinese Academy of Sciences Grant No. XDB23010200, Key Research Program of the Chinese Academy of Sciences, Grant No. XDPB15, and the China Manned Space Program through its Space Application System, and the China Postdoctoral Science Foundation grant No.2019M662168. This research has made use of data, software and/or web tools obtained from the Gravitational Wave Open Science Center (https://www.gw-openscience.org/ ), a service of LIGO Laboratory, the LIGO Scientific Collaboration and the Virgo Collaboration. LIGO Laboratory and Advanced LIGO are funded by the United States National Science Foundation (NSF) as well as the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. Virgo is funded, through the European Gravitational Observatory (EGO), by the French Centre National de Recherche Scientifique (CNRS), the Italian Istituto Nazionale di Fisica Nucleare (INFN) and the Dutch Nikhef, with contributions by institutions from Belgium, Germany, Greece, Hungary, Ireland, Japan, Monaco, Poland, Portugal, Spain.

This data release corresponds to the paper "Constraining scalar-tensor theories by neutron star-balck hole gravitational wave events" (arXiv:2105.13644). In this paper, we consider three specific models of scalar-tensor theories, including the Brans-Dicke theory (BD), the theory with scalarization phenomena proposed by Damour and Esposito-Far\`{e}se (DEF), and Screened Modified Gravity (SMG). From all 4 possible NSBH events so far, we use two of them to place the constraints. The other two are excluded in this work due to the possible unphysical deviations. Four equations of state (EoSs), sly, alf2, H4 and mpa1, are used to derive the scalar charges of neutron stars for BD and DEF. The constraints are obtained by performing the full Bayesian inference with the help of the open source software Bilby. This dataset contains all posterior samples of the runs discussed in the paper. The models and EoSs can be read form the filenames for the runs of BD and DEF. The files of "*_half_dipole.json" correspond to the runs for constraining the dipole radiation without considering specific model parameters. All files are JSON format which is the default output format of Bilby. They are human readable and also can be processed or visualized by Bilby or PESummary conveniently.