diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.cc b/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.cc
@@ -1,128 +1,128 @@
 // -*- C++ -*-
 #include "Rivet/Analysis.hh"
 #include "Rivet/Projections/ChargedLeptons.hh"
 #include "Rivet/Projections/VetoedFinalState.hh"
 #include "Rivet/Projections/MissingMomentum.hh"
 #include "Rivet/Projections/FastJets.hh"
 #include "Rivet/Projections/DressedLeptons.hh"
 
 namespace Rivet {
 
   /// ATLAS same-sign WW production at 8 TeV (PRL)
   class ATLAS_2014_I1298023 : public Analysis {
   public:
 
     /// @name Constructors etc.
     //@{
     ATLAS_2014_I1298023()
       : Analysis("ATLAS_2014_I1298023")
     {
       setNeedsCrossSection(true);
     }
     //@}
 
 
   public:
 
     /// @name Analysis methods
     //@{
 
     void init() {
 
       const FinalState fs;
 
       // bare leptons
       ChargedLeptons bare_leptons(fs);
 
       // dressed leptons
       Cut cuts = (Cuts::abseta < 2.5) & (Cuts::pT > 25*GeV);
       DressedLeptons leptons(fs, bare_leptons, 0.1, cuts);
       declare(leptons, "leptons");
 
       // MET
       declare(MissingMomentum(fs), "MissingET");
 
       // jets
       VetoedFinalState vfs(fs);
       vfs.addVetoPairId(PID::MUON);
       vfs.vetoNeutrinos();
       declare(FastJets(vfs, FastJets::ANTIKT, 0.4), "jets");
 
       // book histogram
       _hist = bookHisto1D(1, 1, 1);
     }
 
 
     /// Perform the per-event analysis
     void analyze(const Event& event) {
       const double weight = event.weight();
 
       const vector<DressedLepton>& leptons = apply<DressedLeptons>(event, "leptons").dressedLeptons();
       if ( leptons.size() < 2 )  vetoEvent;
 
       double minDR_ll = MAXDOUBLE, mll = -1.0;
       for (unsigned int i = 0; i < leptons.size(); ++i) {
         DressedLepton lep1 = leptons[i];
         for (unsigned int j = i + 1; j < leptons.size(); ++j) {
           DressedLepton lep2 = leptons[j];
           double dr = deltaR(lep1, lep2);
           if ( dr < minDR_ll )  minDR_ll = dr;
           double m = (lep1.momentum() + lep2.momentum()).mass();
           if ( mll < 0. || m < mll )  mll = m;
         }
       }
       if ( minDR_ll <= 0.3 || mll <= 20*GeV )  vetoEvent;
 
       if ( leptons[0].charge() * leptons[1].charge() < 0.0 )  vetoEvent;
 
       const MissingMomentum& met = apply<MissingMomentum>(event, "MissingET");
       if ( met.vectorEt().mod() <= 40*GeV )  vetoEvent;
 
       const Jets& all_jets = apply<FastJets>(event, "jets").jetsByPt( (Cuts::abseta < 4.5) && (Cuts::pT > 30*GeV) );
       Jets jets;
       double minDR_overall = MAXDOUBLE;
       foreach (const Jet& jet, all_jets) {
         double minDR_jet = MAXDOUBLE, minDR_electrons = MAXDOUBLE;
         foreach( DressedLepton lep, leptons ) {
           double dr = deltaR(jet, lep);
           if ( dr < minDR_jet )  minDR_jet = dr;
           if ( lep.abspid() == 11 && dr < minDR_electrons )  minDR_electrons = dr;
         }
         if ( minDR_electrons < 0.05 )  continue; // veto jet if it overlaps with electron
         if ( minDR_jet < minDR_overall )  minDR_overall = minDR_jet;
         jets += jet;
       }
       if ( jets.size() < 2 || minDR_overall <= 0.3 )  vetoEvent;
 
       FourMomentum dijet = jets[0].momentum() + jets[1].momentum();
       if ( dijet.mass() <= 500*GeV )  vetoEvent;
 
       // inclusive region
-      _hist->fill(0.5, weight);
+      _hist->fill(1.0, weight);
 
       // VBS region
-      if ( deltaRap(jets[0], jets[1]) > 2.4 )  _hist->fill(1.5, weight);
+      if ( deltaRap(jets[0], jets[1]) > 2.4 )  _hist->fill(2.0, weight);
     }
 
     /// Normalise histograms etc., after the run
     void finalize() {
 
       const double scalefactor( crossSection() / sumOfWeights() );
       scale(_hist, scalefactor);
 
     }
     //@}
 
   private:
 
     /// @name Histograms
     //@{
     Histo1DPtr _hist;
     //@}
 
   };
 
   // The hook for the plugin system
   DECLARE_RIVET_PLUGIN(ATLAS_2014_I1298023);
 
 }
diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.plot b/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.plot
@@ -1,7 +1,7 @@
 # BEGIN PLOT /ATLAS_2014_I1298023/d.*
-XCustomMajorTicks=0.5	inclusive	1.5	VBS
+XCustomMajorTicks=1.0	inclusive	2.0	VBS
 XLabel=phase space
 YLabel=$\sigma^\text{fid}$ [pb]
 LogY=0
 # END PLOT
 
diff --git a/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda b/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
--- a/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
+++ b/analyses/pluginATLAS/ATLAS_2014_I1298023.yoda
@@ -1,10 +1,10 @@
 BEGIN YODA_SCATTER2D_V2 /REF/ATLAS_2014_I1298023/d01-x01-y01
 IsRef: 1
 Path: /REF/ATLAS_2014_I1298023/d01-x01-y01
-Title: ~
+Title: doi:10.17182/hepdata.64182.v1/t1
 Type: Scatter2D
 ---
 # xval	 xerr-	 xerr+	 yval	 yerr-	 yerr+
-5.000000e-01	5.000000e-01	5.000000e-01	2.100000e+00	5.830952e-01	5.830952e-01
-1.500000e+00	5.000000e-01	5.000000e-01	1.300000e+00	4.472136e-01	4.472136e-01
+1.000000e+00	5.000000e-01	5.000000e-01	2.100000e+00	5.830952e-01	5.830952e-01
+2.000000e+00	5.000000e-01	5.000000e-01	1.300000e+00	4.472136e-01	4.472136e-01
 END YODA_SCATTER2D_V2