Two datasets of quark and gluon jets generated with Herwig 7.1.4, one with all kinematically realizable quark jets and one that excludes charm and bottom quark jets (at the level of the hard process), analogous to this dataset of Pythia jets. Note that the two datasets in this record should not be combined. Generation parameters are listed below: Herwig 7.1.4, \(\sqrt{s}=14\,\text{TeV} \) Quarks from \(gq\to Z(\to\nu\bar\nu)q\), gluons from \(q\bar q\to Z(\to\nu\bar\nu)g\) FastJet 3.3.0, anti-kT jets with R=0.4 \(p_T^\text{jet}\in[500,550]\,\text{GeV},\,|y^\text{jet} |<1.7\) There are 20 files in each dataset, each in compressed NumPy format. Files including charm and bottom jets have 'withbc' in their filename. There are two arrays in each file X: (100000,M,4), exactly 50k quark and 50k gluon jets, randomly sorted, where M is the max multiplicity of the jets in that file (other jets have been padded with zero-particles), and the features of each particle are its pt, rapidity, azimuthal angle, and pdgid. y: (100000,), an array of labels for the jets where gluon is 0 and quark is 1. If you use this dataset, please cite this Zenodo record and, optionally, the Pythia dataset which inspired it. The datasets can be downloaded and read into python automatically using the EnergyFlow Python package. Changes: v1 - Renamed files from v0 to include 'withbc' (events were also shuffled around), added files without b and c quark jets.