diff --git a/Changes-API.md b/Changes-API.md
index 2c8667e..a6fae4f 100644
--- a/Changes-API.md
+++ b/Changes-API.md
@@ -1,116 +1,125 @@
 # Changelog for HEJ API
 
 This log lists only changes on the HEJ API. These are primarily code changes
 relevant for calling HEJ as an API. This file should only be read as an addition
 to [`Changes.md`](Changes.md), where the main features are documented.
 
+## Version 2.2
+
+### 2.2.0
+
+#### New and updated functions
+* Updated function `Event::EventData.reconstruct_intermediate()`
+  - Now requires an `EWConstants` argument
+  - Added support for WpWp events.
+  - In the case of WpWp same-flavour the reconstruction will minimise the total
+    off-shell momentum.
+
+
 ## Version 2.1
 
 ### 2.1.0
 
 #### Changes to Event class
 * Restructured `Event` class
   - `Event` can now only be build from a (new) `Event::EventData` class
   - Removed default constructor for `Event`
   - `Event::EventData` replaces the old `UnclusteredEvent` struct.
   - `UnclusteredEvent` is now **deprecated**, and will be removed in HEJ Version
     2.3.0
   - Removed `Event.unclustered()` function
   - Added new member function `Events.parameters()`, to directly access
     (underlying) `Parameters<EventParameters>`
   - New member functions `begin_partons`, `end_partons` (`rbegin_partons`,
     `rend_partons`) with aliases `cbegin_partons`, `cend_partons`
     (`crbegin_partons`, `crend_partons`) for constant (reversed) iterators over
     outgoing partons.
 * New function `Event::EventData.reconstruct_intermediate()` to reconstruct
   bosons from decays, e.g. `positron + nu_e => Wp`
-  - Only supports 1-boson or `WpWp` events.
-  - In the case of WpWp same-flavour the reconstruction will minimise the total
-    off-shell momentum.
 * Added optional Colour charges to particles (`Particle.colour`)
   - Colour connection in the HEJ limit can be generated via
     `Event.generate_colours` (automatically done in the resummation)
   - Colour configuration of input events can be checked with
     `Event.is_leading_colour`
 * Added function `Event.valid_hej_state` to check if event _could have_ been
   produced by `HEJ` according to the `soft pt regulator` cut on the jets
 
 #### New and updated functions
 * Renamed `EventType::nonHEJ` to `EventType::non_resummable` and `is_HEJ()`
   to `is_resummable()` such that run card is consistent with internal workings
 * Made `MatrixElement.tree_kin(...)` and `MatrixElement.tree_param(...)` public
 * New `EventReweighter` member function `treatment` to query the
   treatment with respect to resummation for the various event types.
 * Added auxiliary functions to obtain a `std::string` from `EventDescription`
   (`to_string` for human readable, and `to_simple_string` for easy parsable
   string)
 * New `get_analyses` function to read in multiple `HEJ::Analysis` at once,
   similar to `get_analysis`
 * New `get_ew_parameters` function to extract electroweak parameters from
   YAML configuration.
 
 #### New classes
 * New class `Unweighter` to do unweighting
 * New class `CrossSectionAccumulator` to keep track of Cross Section of the
   different subprocess
 * New template struct `Parameters` similar to old `Weights`
   - `Weights` are now an alias for `Parameters<double>`. Calling `Weights` did
     not change
   - `Weights.hh` was replaced by `Parameters.hh`. The old `Weights.hh` header
     will be removed in HEJ Version 2.2.0
 * Function to multiplication and division of `EventParameters.weight` by double
   - This can be combined with `Parameters`, e.g.
     `Parameters<EventParameters>*Weights`, see also `Events.parameters()`
   - Moved `EventParameters` to `Parameters.hh` header
 * new class `EWConstants` replaces previously hard coded `vev`
   - `EWConstants` have to be set in the general `Config` and the
     `MatrixElementConfig`
 
 #### Input/Output
 * New abstract `EventReader` class, as base for reading events from files
   - Moved LHE file reader to `HEJ::LesHouchesReader`
 * New (optional) function `finish()` in abstract class `EventWriter`. `finish()`
   is called _after_ all events are written.
 * Support reading (`HDF5Reader`) and writing (`HDF5Writer`) `hdf5` files
 * New `BufferedEventReader` class that allows to put events back into
   the reader.
 * New `SherpaLHEReader` to read Sherpa LHE files with correct weights
 * `get_analysis` now requires `YAML::Node` and `LHEF::HEPRUP` as arguments
 * Replaced `HepMCInterface` and `HepMCWriter` by `HepMCInterfaceX` and
   `HepMCWriterX` respectively, with `X` being the major version of HepMC (2 or
   3)
   - Renamed `HepMCInterfaceX::init_kinematics` to `HepMCInterfaceX::init_event`
     and protected it, use `HepMCInterfaceX::operator()` instead
   - Removed redundant `HepMCInterfaceX::add_variation` function
 
 #### Linking
 * Export cmake target `HEJ::HEJ` to link directly against `libHEJ`
 * Preprocessor flags (`HEJ_BUILD_WITH_XYZ`) for enabled/disabled dependencies
   are now written to `ConfigFlags.hh`
 * Provide links to version specific object files, e.g. `libHEJ.so ->
   libHEJ.so.2.1 (soname) -> libHEJ.so.2.1.0`
 * Removed `LHAPDF` from public interface
 
 #### Miscellaneous
 * Capitalisation of `Config.hh` now matches class `Config` (was `config.hh`)
 * Renamed `Config::max_ext_soft_pt_fraction` to `Config::soft_pt_regulator`.
   The new `Config::soft_pt_regulator` parameter is optional and the old
   `Config::max_ext_soft_pt_fraction` is **deprecated**.
 * Replaced redundant member `EventReweighterConfig::jet_param` with getter
   function `EventReweighter.jet_param()` (`JetParameters` are already in
   `PhaseSpacePointConfig`)
 * Avoid storing reference to the Random Number Generator inside classes
   - Constructors of `EventReweighter` now expect `std::shared_ptr<HEJ::RNG>`
     (was reference)
   - Moved reference to `HEJ::RNG` from constructor of `JetSplitter` to
     `JetSplitter.split`
 
 ## Version 2.0
 
 ### 2.0.4
 
 * Fixed wrong path of `HEJ_INCLUDE_DIR` in `hej-config.cmake`
 
 ### 2.0.0
 
 * First release
diff --git a/Changes.md b/Changes.md
index 5d9987b..073079a 100644
--- a/Changes.md
+++ b/Changes.md
@@ -1,147 +1,160 @@
 # Changelog
 
 This is the log for changes to the HEJ program. Further changes to the HEJ API
 are documented in [`Changes-API.md`](Changes-API.md). If you are using HEJ as a
 library, please also read the changes there.
 
+## Version 2.2
+
+This release adds support for a new process:
+
+* Wp+Wp with jets.
+
+This release includes minor changes which affect users of HEJ as a library (see
+[`Changes-API.md`](Changes-API.md)).
+
+### 2.2.0
+
+#### New Processes
+
+* Resummation for WpWp with jets, including interference between configurations.
+
 ## Version 2.1
 
-This release adds support for three new processes:
+This release adds support for two new processes:
 
 * W boson with jets.
-* Wp+Wp with jets.
 * Jets with a charged lepton-antilepton pair via a virtual Z boson or photon.
 
 In addition, the complete set of first subleading processes (unordered
 gluon, central and extremal quark-antiquark pair) is implemented for
 pure jets and W + jets, see
 [arXiv:TODO](https://arxiv.org/abs/TODO). Unordered gluon emission is
 also supported for Higgs boson + jets and Z boson/photon + jets.
 
 This release include many changes to the code, which affect users of
 HEJ as a library (see [`Changes-API.md`](Changes-API.md)).
 
 ### 2.1.0
 
 #### New Processes
 * Resummation for W bosons with jets
   - New subleading processes `extremal qqx` & `central qqx` for a quark and
     anti-quark in the final state, e.g. `g g => u d_bar Wm g` (the other
     subleading processes also work with W's)
   - `HEJFOG` can generate multiple jets together with a (off-shell) W bosons
     decaying into lepton & neutrino
-* Resummation for WpWp with jets, including interference between configurations.
 * Resummation for jets with a charged lepton-antilepton pair via a
   virtual Z boson or photon. Includes the `unordered` subleading process.
 * Resummation can now be performed on all subleading processes within pure
   jets also. This includes `unordered`, `extremal qqbar` and `central
   qqbar` processes.
 
 #### More Physics implementation
 * Partons now have a Colour charge
   - Colours are read from and written to LHE files
   - For reweighted events the colours are created according to leading colour in
     the FKL limit
 * Use relative fraction for soft transverse momentum in tagging jets (`soft pt
   regulator`) as new (optional) parameter.
   - This supersedes `min extparton pt`, which is marked **deprecated** and will
     be removed in future versions.
   - This is a direct replacement for the old `max ext soft pt fraction`, which
     is also **deprecated**.
 * Changed the redistribution of the momenta of soft emissions. Now also bosons
   can take part of the recoil (previously only jets).
 
 #### Updates to Runcard
 * Allow multiplication and division of multiple scale functions e.g.
   `H_T/2*m_j1j2`
 * Grouped `event treatment` for subleading channels together in runcard
   - Rename `non-HEJ` processes to `non-resummable`
 * Read electro-weak constants from input
   - new mandatory setting `vev` to change vacuum expectation value
   - new mandatory settings `particle properties` to specify mass & width of
     bosons
   - `HEJFOG`: decays are now specified in `decays` setting (previously under
     `particle properties`)
 * Allow loading multiple analyses with `analyses`. The old `analysis` (with "i")
   is marked **deprecated**.
 * Optional setting to specify maximal number of Fixed Order events (`max
   events`, default is all)
 * Allow changing the regulator lambda in input (`regulator parameter`, only for
   advanced users)
 
 #### Changes to Input/Output
 * Added support to read & write `hdf5` event files suggested in
   [arXiv:1905.05120](https://arxiv.org/abs/1905.05120) (needs
   [HighFive](https://github.com/BlueBrain/HighFive))
 * Support input with average weight equal to the cross section (`IDWTUP=1 or 4`)
 * Support unmodified Les Houches Event Files written by Sherpa with
   `cross section = sum(weights)/sum(trials)`
 * Analyses now get general run information (`LHEF::HEPRUP`) in the
   constructor. **This might break previously written, external analyses!**
   - external analyses should now be created with
     `make_analysis(YAML::Node const & config, LHEF::HEPRUP const & heprup)`
 * Support `rivet` version 3 with both `HepMC` version 2 and 3
   - Multiple weights with `rivet 3` will only create one `.yoda` file (instead
     of one per weight/scale)
 * Added option to unweight only resummation events
   (`unweight: {type: resummation}`)
 * Added option for partially unweighting resummation events, similar to
   the fixed-order generator.
 * Improved unweighting algorithm.
 * Follow HepMC convention for particle Status codes: incoming = 11,
   decaying = 2, outgoing = 1 (unchanged)
 
 #### Miscellaneous
 * Print cross sections at end of run
 * Added example analysis & scale to `examples/`. Everything in `examples/` will
   be build when the flag `-DBUILD_EXAMPLES=TRUE` is set in `cmake`.
 * Dropped support for HepMC 3.0.0, either HepMC version 2 or >3.1 is required
   - It is now possible to write out both HepMC 2 and HepMC 3 events at the same
     time
 * Require LHADPF version 6. Dropped support for all other versions.
 * Use `git-lfs` for raw data in test (`make test` now requires `git-lfs`)
 * Currents are now generated with [`FORM`](https://github.com/vermaseren/form)
   - `FORM` is included as a `git submodule`, use `git submodule update --init`
     to download `FORM`
 * Create [Sphinx](http://sphinx-doc.org/) and [Doxygen](http://doxygen.org/)
   documentation by `make sphinx` or `make doxygen` in your `build/` folder
 
 ## Version 2.0
 
 First release of HEJ 2. Complete code rewrite compared to HEJ 1. Improved
 matching to Fixed Order ([arXiv:1805.04446](https://arxiv.org/abs/1805.04446)).
 Implemented processes: Higgs boson with jets (FKL and unordered gluon emission,
 with finite quark mass loop,
 [arXiv:1812.08072](https://arxiv.org/abs/1812.08072)), and pure jets (only FKL).
 See [arXiv:1902.08430](https://arxiv.org/abs/1902.08430)
 
 
 ### 2.0.6
 
 * Fixed compiling rivet when YODA headers are _outside_ of rivet directory
 
 ### 2.0.5
 
 * Fixed event classification for input not ordered in rapidity
 
 ### 2.0.4
 
 * Fixed wrong path of `HEJ_INCLUDE_DIR` in `hej-config.cmake`
 
 ### 2.0.3
 
 * Fixed parsing of (numerical factor) * (base scale) in configuration
 * Don't change scale names, but sanitise Rivet output file names instead
 
 ### 2.0.2
 
 * Changed scale names to `"_over_"` and `"_times_"` for proper file names (was
   `"/"` and `"*"` before)
 
 ### 2.0.1
 
 * Fixed name of fixed-order generator in error message.
 
 ### 2.0.0
 
 * First release