No OneTemporary
Actions

Size

113 KB

Subscribers

None

View Options

	diff --git a/AnalysisTools/GridSetup/BatchSubmitHad.py b/AnalysisTools/GridSetup/BatchSubmitHad.py
	--- a/AnalysisTools/GridSetup/BatchSubmitHad.py
	+++ b/AnalysisTools/GridSetup/BatchSubmitHad.py
	@@ -1,78 +1,80 @@
	#!/usr/bin/python

	import sys
	import os
	import string
	import tarfile
	import random
	import errno
	import time
	import subprocess

	def mkdir_p(path):
	try:
	os.makedirs(path)
	except OSError as exc: # Python >2.5
	if exc.errno == errno.EEXIST and os.path.isdir(path):
	pass
	else:
	raise

	import shutil

	def copytree(src, dst, symlinks=False, ignore=None):
	if not os.path.exists(dst):
	os.makedirs(dst)
	for item in os.listdir(src):
	s = os.path.join(src, item)
	d = os.path.join(dst, item)
	if os.path.isdir(s):
	copytree(s, d, symlinks, ignore)
	else:
	if not os.path.exists(d) or os.stat(s).st_mtime - os.stat(d).st_mtime > 1:
	shutil.copy2(s, d)

	pwd = os.getcwd()
	i=0
	k=0

	+# You should check these are executing the correct scripts. Especially the setupContur script
	+# should probably be one in your own area.
	HerwigSetup="source /unix/cedar/software/sl6/Herwig-7.1.0/setupEnv.sh"
	-ConturSetup="source /home/jmb/svn-managed/contur/setupContur.sh"
	+ConturSetup="source $HOME/contur/setupContur.sh"

	for i in range(100,3100,100):
	for j in range(100,2100,100):
	modelpath = 'mY_'+ str(i) + '_mX_' + str(j)
	copytree(pwd + '/GridPack',modelpath)
	#mkdir_p(str(modelpath))
	HerwigString = ''
	HC=open('HerwigCommandHad', 'r')

	HerwigString += 'read FRModel.model \n'
	HerwigString += 'set /Herwig/FRModel/Particles/Y1:NominalMass ' + str(i) + '.*GeV \n'
	HerwigString += 'set /Herwig/FRModel/Particles/Xm:NominalMass ' + str(j) + '.*GeV \n'
	HerwigString += str(HC.read())
	HC.close()
	RunCard = open(str(modelpath + '/LHC.in'), 'w')
	RunCard.write(str(HerwigString))
	RunCard.close()

	subprocess.call([HerwigSetup], shell=True)
	subprocess.call([ConturSetup], shell=True)
	os.chdir(modelpath)
	subprocess.call(['Herwig read LHC.in'], shell=True)
	batch_command = ''
	batch_command += HerwigSetup + '; '
	batch_command += 'cd ' + pwd + '/' + modelpath +'; '
	batch_command += ConturSetup + "; "
	if i < 600:
	numEv=30000
	else:
	numEv=15000
	batch_command += 'Herwig run --seed='+str(i)+str(j)+' --tag='+str(modelpath)+' --jobs=2 --numevents='+ str(numEv) +' LHC.run;'
	batch_filename = str(modelpath)+'.sh'
	batch_submit = open(batch_filename, 'w')
	batch_submit.write(batch_command)
	batch_submit.close()
	subprocess.call([ "qsub -q medium " + batch_filename],shell=True )
	os.chdir(pwd)

	diff --git a/AnalysisTools/GridSetup/BatchSubmitWeak.py b/AnalysisTools/GridSetup/BatchSubmitWeak.py
	--- a/AnalysisTools/GridSetup/BatchSubmitWeak.py
	+++ b/AnalysisTools/GridSetup/BatchSubmitWeak.py
	@@ -1,81 +1,83 @@
	#!/usr/bin/python

	import sys
	import os
	import string
	import tarfile
	import random
	import errno
	import time
	import subprocess
	import numpy as np

	def mkdir_p(path):
	try:
	os.makedirs(path)
	except OSError as exc: # Python >2.5
	if exc.errno == errno.EEXIST and os.path.isdir(path):
	pass
	else:
	raise

	import shutil

	def copytree(src, dst, symlinks=False, ignore=None):
	if not os.path.exists(dst):
	os.makedirs(dst)
	for item in os.listdir(src):
	s = os.path.join(src, item)
	d = os.path.join(dst, item)
	if os.path.isdir(s):
	copytree(s, d, symlinks, ignore)
	else:
	if not os.path.exists(d) or os.stat(s).st_mtime - os.stat(d).st_mtime > 1:
	shutil.copy2(s, d)

	pwd = os.getcwd()
	i=0
	k=0

	+# You should check these are executing the correct scripts. Especially the setupContur script
	+# should probably be one in your own area.
	HerwigSetup="source /unix/cedar/software/sl6/Herwig-7.1.0/setupEnv.sh"
	-ConturSetup="source /home/jmb/svn-managed/contur/setupContur.sh"
	+ConturSetup="source $HOME/contur/setupContur.sh"

	for i in range(100,3100,100):
	for j in range(100,2100,100):
	i = int(i)
	j = int(j)
	modelpath = 'mY_'+ str(i) + '_mX_' + str(j)
	copytree(pwd + '/GridPack',modelpath)
	#mkdir_p(str(modelpath))
	HerwigString = ''
	HC=open('HerwigCommandWeak', 'r')

	HerwigString += 'read FRModel.model \n'
	HerwigString += 'set /Herwig/FRModel/Particles/Y1:NominalMass ' + str(i) + '.*GeV \n'
	HerwigString += 'set /Herwig/FRModel/Particles/Xm:NominalMass ' + str(j) + '.*GeV \n'
	HerwigString += str(HC.read())
	HC.close()
	RunCard = open(str(modelpath + '/LHC.in'), 'w')
	RunCard.write(str(HerwigString))
	RunCard.close()

	subprocess.call([HerwigSetup], shell=True)
	subprocess.call([ConturSetup], shell=True)
	os.chdir(modelpath)
	subprocess.call(['Herwig read LHC.in'], shell=True)
	batch_command = ''
	batch_command += HerwigSetup + '; '
	batch_command += 'cd ' + pwd + '/' + modelpath +'; '
	batch_command += ConturSetup + "; "
	if i < 1500:
	numEv=30000
	else:
	numEv=5000
	batch_command += 'Herwig run --seed='+str(i)+str(j)+' --tag='+str(modelpath)+' --jobs=5 --numevents='+ str(numEv) +' LHC.run;'
	batch_filename = str(modelpath)+'.sh'
	batch_submit = open(batch_filename, 'w')
	batch_submit.write(batch_command)
	batch_submit.close()

	subprocess.call([ "qsub -q medium " + batch_filename],shell=True )
	os.chdir(pwd)
	diff --git a/AnalysisTools/GridSetup/GridPack/13_WEAK.ana b/AnalysisTools/GridSetup/GridPack/13_WEAK.ana
	--- a/AnalysisTools/GridSetup/GridPack/13_WEAK.ana
	+++ b/AnalysisTools/GridSetup/GridPack/13_WEAK.ana
	@@ -1,3 +1,5 @@
	- insert Rivet:Analyses 0 ATLAS_2017_I1514251_MU # Z+jet
	- insert Rivet:Analyses 0 ATLAS_2017_I1514251_EL # Z+jet
	+insert Rivet:Analyses 0 ATLAS_2017_I1514251_MU # Z+jet
	+insert Rivet:Analyses 0 ATLAS_2017_I1514251_EL # Z+jet
	+insert Rivet:Analyses 0 CMS_2017_I1610623 # W+jet
	+
	# insert Rivet:Analyses 0 ATLAS_2017_I1609448 # MET-JET analysis
	diff --git a/AnalysisTools/GridSetup/GridPack/8_HAD.ana b/AnalysisTools/GridSetup/GridPack/8_HAD.ana
	--- a/AnalysisTools/GridSetup/GridPack/8_HAD.ana
	+++ b/AnalysisTools/GridSetup/GridPack/8_HAD.ana
	@@ -1,9 +1,13 @@
	insert Rivet:Analyses 0 ATLAS_2015_I1394679 # four jet
	insert Rivet:Analyses 0 ATLAS_2017_I1598613 # b hadron (muon final state)
	+insert Rivet:Analyses 0 CMS_2017_I1499471 # b jets
	+
	+
	+
	# insert Rivet:Analyses 0 ATLAS_2017_I1598613_BB # b hadron (hadron final state. can't use at present, see analysis db.



	#insert Rivet:Analyses 0 CMS_2016_I1421646 # Dijet decorrelations. All normalised and the total is not reported, so can't be used until we have the SM cross section prediction too.


	diff --git a/AnalysisTools/contur/contur/TestingFunctions/analyses.sql b/AnalysisTools/contur/contur/TestingFunctions/analyses.sql
	--- a/AnalysisTools/contur/contur/TestingFunctions/analyses.sql
	+++ b/AnalysisTools/contur/contur/TestingFunctions/analyses.sql
	@@ -1,467 +1,474 @@
	PRAGMA foreign_keys=ON;
	BEGIN TRANSACTION;
	CREATE TABLE analysis_pool (
	pool TEXT NOT NULL UNIQUE,
	PRIMARY KEY(pool)
	);
	INSERT INTO analysis_pool VALUES('ATLAS_7_JETS');
	INSERT INTO analysis_pool VALUES('CMS_7_JETS');

	INSERT INTO analysis_pool VALUES('ATLAS_7_GAMMA');
	INSERT INTO analysis_pool VALUES('CMS_7_GAMMA');

	INSERT INTO analysis_pool VALUES('ATLAS_7_LMDY');

	INSERT INTO analysis_pool VALUES('ATLAS_7_Zjj_EL');
	INSERT INTO analysis_pool VALUES('ATLAS_7_Zjj_MU');
	INSERT INTO analysis_pool VALUES('ATLAS_7_Wjj_EL');
	INSERT INTO analysis_pool VALUES('ATLAS_7_Wjj_MU');
	INSERT INTO analysis_pool VALUES('ATLAS_7_Wjj');
	INSERT INTO analysis_pool VALUES('CMS_7_Zjj');
	INSERT INTO analysis_pool VALUES('CMS_7_Wjj');
	INSERT INTO analysis_pool VALUES('ATLAS_7_Z_GAMMA');
	INSERT INTO analysis_pool VALUES('ATLAS_7_W_GAMMA_MU');
	INSERT INTO analysis_pool VALUES('ATLAS_7_W_GAMMA_EL');
	INSERT INTO analysis_pool VALUES('ATLAS_7_ZZ');
	INSERT INTO analysis_pool VALUES('ATLAS_7_WW');

	INSERT INTO analysis_pool VALUES('ATLAS_8_JETS');
	INSERT INTO analysis_pool VALUES('CMS_8_JETS');

	INSERT INTO analysis_pool VALUES('ATLAS_8_GAMMA');
	INSERT INTO analysis_pool VALUES('CMS_8_GAMMA');

	INSERT INTO analysis_pool VALUES('ATLAS_8_HMDY');

	INSERT INTO analysis_pool VALUES('ATLAS_8_Zjj');
	INSERT INTO analysis_pool VALUES('ATLAS_8_Wjj');
	INSERT INTO analysis_pool VALUES('CMS_8_Wjj');

	INSERT INTO analysis_pool VALUES('ATLAS_8_ZZ');
	INSERT INTO analysis_pool VALUES('ATLAS_8_WW');
	INSERT INTO analysis_pool VALUES('CMS_8_WW');
	INSERT INTO analysis_pool VALUES('ATLAS_8_Z_GAMMA');
	INSERT INTO analysis_pool VALUES('ATLAS_8_GAMMA_MET');

	INSERT INTO analysis_pool VALUES('ATLAS_13_JETS');
	INSERT INTO analysis_pool VALUES('CMS_13_JETS');
	INSERT INTO analysis_pool VALUES('ATLAS_13_Zjj');
	+INSERT INTO analysis_pool VALUES('CMS_13_Wjj');
	INSERT INTO analysis_pool VALUES('ATLAS_13_METJET');

	CREATE TABLE analysis (
	id TEXT NOT NULL UNIQUE,
	lumi REAL NOT NULL,
	pool TEXT,
	PRIMARY KEY(id),
	FOREIGN KEY(pool) REFERENCES analysis_pool(pool)
	);
	-- Note, the LUMI value given here should be in units which match the normalisation
	-- of the cross section plots. Most in pb, some in fb.

	-- Superseded/deprecated analyses
	INSERT INTO analysis VALUES('ATLAS_2012_I1083318',36.0,NULL);
	INSERT INTO analysis VALUES('ATLAS_2011_I945498',36.0,NULL);
	INSERT INTO analysis VALUES('ATLAS_2011_I921594',35.0,NULL);
	INSERT INTO analysis VALUES('ATLAS_2011_S9128077',2.4,NULL);

	-- 7 TeV fully hadronic
	INSERT INTO analysis VALUES('ATLAS_2014_I1325553',4500.0,'ATLAS_7_JETS');
	INSERT INTO analysis VALUES('ATLAS_2014_I1268975',4500.0,'ATLAS_7_JETS');
	INSERT INTO analysis VALUES('ATLAS_2014_I1326641',4510.0,'ATLAS_7_JETS');
	INSERT INTO analysis VALUES('ATLAS_2014_I1307243',4500.0,'ATLAS_7_JETS');
	INSERT INTO analysis VALUES('CMS_2014_I1298810',5000.0,'CMS_7_JETS');
	INSERT INTO analysis VALUES('CMS_2014_I1273574',5000.0,'CMS_7_JETS');
	INSERT INTO analysis VALUES('CMS_2013_I1224539_DIJET',5000.0,'CMS_7_JETS');

	-- 7 TeV photons
	INSERT INTO analysis VALUES('ATLAS_2012_I1093738',37.0,'ATLAS_7_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2013_I1244522',37.0,'ATLAS_7_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2013_I1263495',4600.0,'ATLAS_7_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2012_I1199269',4900.0,'ATLAS_7_GAMMA');
	INSERT INTO analysis VALUES('CMS_2014_I1266056',4500.0,'CMS_7_GAMMA');

	-- 7 TeV Z+jets
	INSERT INTO analysis VALUES('ATLAS_2013_I1230812_EL',4600.0,'ATLAS_7_Zjj_EL');
	INSERT INTO analysis VALUES('ATLAS_2013_I1230812_MU',4600.0,'ATLAS_7_Zjj_MU');
	INSERT INTO analysis VALUES('ATLAS_2014_I1306294_EL',4600.0,'ATLAS_7_Zjj_EL');
	INSERT INTO analysis VALUES('ATLAS_2014_I1306294_MU',4600.0,'ATLAS_7_Zjj_MU');
	INSERT INTO analysis VALUES('CMS_2015_I1310737',4900.0,'CMS_7_Zjj');

	-- 7 TeV ttbar
	INSERT INTO analysis VALUES('ATLAS_2015_I1345452',4600.0,'ATLAS_7_Wjj');

	-- 7 TeV Low mass DY
	INSERT INTO analysis VALUES('ATLAS_2014_I1288706',1600.0,'ATLAS_7_LMDY');

	-- 7 TeV single jet masses
	INSERT INTO analysis VALUES('CMS_2013_I1224539_WJET',5.0,'CMS_7_Wjj');
	INSERT INTO analysis VALUES('CMS_2013_I1224539_ZJET',5.0,'CMS_7_Zjj');

	-- 7 TeV W+jets.
	INSERT INTO analysis VALUES('ATLAS_2014_I1319490_MU',4600.0,'ATLAS_7_Wjj_MU');
	INSERT INTO analysis VALUES('ATLAS_2014_I1319490_EL',4600.0,'ATLAS_7_Wjj_EL');
	INSERT INTO analysis VALUES('CMS_2014_I1303894',5000.0,'CMS_7_Wjj');

	-- 7 TeV W+charm
	INSERT INTO analysis VALUES('ATLAS_2014_I1282447',4600.0,'ATLAS_7_Wjj');

	-- 7 TeV W+b
	INSERT INTO analysis VALUES('ATLAS_2013_I1219109_EL',4600.0,'ATLAS_7_Wjj');
	INSERT INTO analysis VALUES('ATLAS_2013_I1219109_MU',4600.0,'ATLAS_7_Wjj');

	-- 7 TeV Z+bb
	INSERT INTO analysis VALUES('CMS_2013_I1256943',5000.0,'CMS_7_Zjj');

	-- 7 TeV WW. Plots in fb.
	INSERT INTO analysis VALUES('ATLAS_2013_I1190187',4.6,'ATLAS_7_WW');

	-- 7 TeV dibosons, in fb-1
	INSERT INTO analysis VALUES('ATLAS_2013_I1217863_Z',4.6,'ATLAS_7_Z_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2013_I1217863_W_MU',4.6,'ATLAS_7_W_GAMMA_MU');
	INSERT INTO analysis VALUES('ATLAS_2013_I1217863_W_EL',4.6,'ATLAS_7_W_GAMMA_EL');
	INSERT INTO analysis VALUES('ATLAS_2012_I1203852',4.6,'ATLAS_7_ZZ');

	-- 8 TeV fully hadronic
	-- plots in fb
	INSERT INTO analysis VALUES('ATLAS_2015_I1394679',20.3,'ATLAS_8_JETS');
	INSERT INTO analysis VALUES('ATLAS_2017_I1598613',11.4,'ATLAS_8_JETS');
	INSERT INTO analysis VALUES('ATLAS_2017_I1598613_BB',11.4,'ATLAS_8_JETS');
	+INSERT INTO analysis VALUES('CMS_2017_I1499471',11.4,'CMS_8_JETS');
	+
	+
	+
	-- normalised, no total xsec yet INSERT INTO analysis VALUES('CMS_2016_I1421646',19700,'CMS_8_JETS');


	-- 8 TeV photons
	-- in pb
	INSERT INTO analysis VALUES('ATLAS_2016_I1457605',20200.0,'ATLAS_8_GAMMA');
	-- in fb
	INSERT INTO analysis VALUES('ATLAS_2017_I1591327',20.3,'ATLAS_8_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2014_I1306615',20.3,'ATLAS_8_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2017_I1632756',20.2,'ATLAS_8_GAMMA');

	-- 8 TeV High mass DY
	-- in pb
	INSERT INTO analysis VALUES('ATLAS_2016_I1467454_EL',20300.0,'ATLAS_8_HMDY');
	INSERT INTO analysis VALUES('ATLAS_2016_I1467454_MU',20300.0,'ATLAS_8_HMDY');

	-- 8 TeV leptons+MET
	-- normalised
	INSERT INTO analysis VALUES('ATLAS_2014_I1279489',20300.0,'ATLAS_8_Zjj');
	INSERT INTO analysis VALUES('ATLAS_2015_I1408516_EL',20300.0,'ATLAS_8_Zjj');
	INSERT INTO analysis VALUES('ATLAS_2015_I1408516_MU',20300.0,'ATLAS_8_Zjj');
	-- plots in pb
	-- INSERT INTO analysis VALUES('ATLAS_2017_I1589844',20200.0,'ATLAS_8_Zjj');

	-- plots in fb
	INSERT INTO analysis VALUES('CMS_2016_I1454211',19.7,'CMS_8_Wjj');
	INSERT INTO analysis VALUES('ATLAS_2015_I1397637',20.3,'ATLAS_8_Wjj');
	INSERT INTO analysis VALUES('CMS_2016_I1491953',19.6,'CMS_8_Wjj');


	-- plots in fb
	INSERT INTO analysis VALUES('CMS_2017_I1518399',19.7,'CMS_8_Wjj');
	INSERT INTO analysis VALUES('ATLAS_2016_I1426515',20.3,'ATLAS_8_WW');
	INSERT INTO analysis VALUES('CMS_2017_I1467451',19.4,'CMS_8_WW');

	-- fb
	INSERT INTO analysis VALUES('ATLAS_2015_I1394865',20.3,'ATLAS_8_ZZ');
	INSERT INTO analysis VALUES('ATLAS_2016_I1448301_EL',20.3,'ATLAS_8_Z_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2016_I1448301_MU',20.3,'ATLAS_8_Z_GAMMA');
	INSERT INTO analysis VALUES('ATLAS_2016_I1448301_NU',20.3,'ATLAS_8_GAMMA_MET');

	-- 13 TeV fully hadronic
	-- plots in pb
	INSERT INTO analysis VALUES('CMS_2016_I1459051',71.0,'CMS_13_JETS');

	-- TBD normalised cross sections INSERT INTO analysis VALUES('CMS_2017_I1519995',2600,'CMS_13_JETS');

	-- 13 TeV photons

	-- 13 TeV leptons+MET
	-- plots in pb
	INSERT INTO analysis VALUES('ATLAS_2017_I1514251_EL',3160.0,'ATLAS_13_Zjj');
	INSERT INTO analysis VALUES('ATLAS_2017_I1514251_MU',3160.0,'ATLAS_13_Zjj');
	INSERT INTO analysis VALUES('ATLAS_2017_I1609448_Zee',3200.0,'ATLAS_13_Zjj');
	INSERT INTO analysis VALUES('ATLAS_2017_I1609448_Zmumu',3200.0,'ATLAS_13_Zjj');
	+INSERT INTO analysis VALUES('CMS_2017_I1610623',2200.0,'CMS_13_Wjj');


	-- 13 TeV MET+JET
	INSERT INTO analysis VALUES('ATLAS_2017_I1609448',3200.0,'ATLAS_13_METJET');
	INSERT INTO analysis VALUES('ATLAS_2017_I1609448_Znunu',3200.0,'ATLAS_13_METJET');



	-- CMS_2016_I1421646 – Dijet azimuthal decorrelations in $pp$ collisions at $\sqrt{s} = 8$ TeV
	-- commented out in python 'CMS_2013_I1256943',5200.0 CMS Z+b



	CREATE TABLE blacklist (
	id TEXT NOT NULL,
	pattern TEXT NOT NULL,
	UNIQUE(id,pattern),
	FOREIGN KEY(id) REFERENCES analysis(id)
	);
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_EL','d02');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_EL','d04');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_EL','d06');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_EL','d08');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_MU','d02');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_MU','d04');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_MU','d06');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1230812_MU','d08');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d13');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d14');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d15');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d16');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d17');
	INSERT INTO blacklist VALUES('CMS_2014_I1298810','d18');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_Z','d17');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_Z','d18');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_Z','d20');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_MU','d15');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_MU','d16');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_MU','d19');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_EL','d15');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_EL','d16');
	INSERT INTO blacklist VALUES('ATLAS_2013_I1217863_W_EL','d19');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d03');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d04');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d05');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d06');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d07');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1203852','d08');
	INSERT INTO blacklist VALUES('ATLAS_2012_I1199269','d04');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d10-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d11-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d12-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d13-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d14-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2016_I1426515','d15-x01-y02');
	INSERT INTO blacklist VALUES('ATLAS_2014_I1306615','d29');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_EL','d01-x02-y01');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_MU','d01-x02-y01');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_EL','d01-x02-y02');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_MU','d01-x02-y02');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_EL','d01-x02-y03');
	INSERT INTO blacklist VALUES('ATLAS_2017_I1514251_MU','d01-x02-y03');
	INSERT INTO blacklist VALUES('CMS_2017_I1518399','d02');
	INSERT INTO blacklist VALUES('CMS_2016_I1491953','d36');
	INSERT INTO blacklist VALUES('CMS_2016_I1491953','d37');
	INSERT INTO blacklist VALUES('CMS_2016_I1491953','d38');
	INSERT INTO blacklist VALUES('CMS_2016_I1491953','d39');




	-- have to veto all these at the moment because contur doesn't know
	--- how to handle weight differential xsecs presented as a 2D scatter.
	+-- how to handle weighted differential xsecs presented as a 2D scatter.
	INSERT INTO blacklist VALUES('ATLAS_2017_I1598613_BB','d01');

	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d01-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d02-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d03-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d04-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d01-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d02-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d03-x01-y01');
	--- INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d04-x01-y01');
	+
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d01-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d02-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d03-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448','d04-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d01-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d02-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d03-x01-y01');
	+INSERT INTO blacklist VALUES('ATLAS_2017_I1609448_Znunu','d04-x01-y01');



	CREATE TABLE whitelist (
	id TEXT NOT NULL,
	pattern TEXT NOT NULL,
	UNIQUE(id,pattern),
	FOREIGN KEY(id) REFERENCES analysis(id)
	);
	-- W charm
	INSERT INTO whitelist VALUES('ATLAS_2014_I1282447','d03-x01-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1282447','d04-x01-y03');
	--
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d01-x01-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d01-x02-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d02-x01-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d03-x01-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d03-x02-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d04-x01-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d04-x02-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d05-x03-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d05-x04-y01');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1279489','d05-x05-y01');
	--
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d13');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d14');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d15');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d16');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d17');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d18');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d19');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d20');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d21');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d22');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d23');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d24');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d25');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d26');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d27');
	INSERT INTO whitelist VALUES('ATLAS_2014_I1307243','d28');
	--
	INSERT INTO whitelist VALUES('CMS_2016_I1454211','d10');
	INSERT INTO whitelist VALUES('CMS_2016_I1454211','d12');

	CREATE TABLE subpool (
	id TEXT NOT NULL,
	pattern TEXT NOT NULL,
	subid INTEGER NOT NULL,
	UNIQUE(id,pattern),
	FOREIGN KEY(id) REFERENCES analysis(id)
	);
	INSERT INTO subpool VALUES('ATLAS_2014_I1325553','d01-x01-y0[1-6]',0);
	INSERT INTO subpool VALUES('ATLAS_2014_I1325553','d02-x01-y0[1-6]',1);
	INSERT INTO subpool VALUES('ATLAS_2014_I1268975','d01-x01-y0[1-6]',0);
	INSERT INTO subpool VALUES('ATLAS_2014_I1268975','d02-x01-y0[1-6]',1);
	INSERT INTO subpool VALUES('ATLAS_2014_I1326641','d0[1-5]-x01-y01',0);
	INSERT INTO subpool VALUES('ATLAS_2014_I1326641','(d10\|d0[6-9])-x01-y01',1);
	INSERT INTO subpool VALUES('CMS_2014_I1298810','d0[1-6]-x01-y01',0);
	-- check this next one, one of the .py entries had y02 listed
	INSERT INTO subpool VALUES('CMS_2014_I1298810','(d1[0-2]\|d0[7-9])-x01-y01',1);
	INSERT INTO subpool VALUES('ATLAS_2014_I1307243','(d20\|d1[3-9])-x01-y01',0);
	INSERT INTO subpool VALUES('ATLAS_2014_I1307243','d2[1-8]-x01-y01',1);
	INSERT INTO subpool VALUES('ATLAS_2013_I1263495','d01-x01-y0[13]',0);


	INSERT INTO subpool VALUES('ATLAS_2016_I1457605','.',0);

	CREATE TABLE normalization (
	-- analysis id, plot pattern, norm factor for ref data, was mc plot area normalised by rivet? (0=no, 1=yes)
	id TEXT NOT NULL,
	pattern TEXT NOT NULL,
	norm REAL NOT NULL,
	scalemc INT NOT NULL,
	UNIQUE(id,pattern),
	FOREIGN KEY(id) REFERENCES analysis(id)
	);
	-- this is BR to a single charged lepton, needed when the xsec is report as a W
	-- and the generator reports the final state.
	-- INSERT INTO normalization VALUES('ATLAS_2014_I1319490_MU','d',0.108059,0);
	-- INSERT INTO normalization VALUES('ATLAS_2014_I1319490_EL','d',0.108059,0);
	-- INSERT INTO normalization VALUES('CMS_2014_I1303894','d',0.108059,0);

	-- these are the integrated cross section of the plot, in pb
	INSERT INTO normalization VALUES('ATLAS_2014_I1279489','d01',5.88,1);
	INSERT INTO normalization VALUES('ATLAS_2014_I1279489','d02',1.82,1);
	INSERT INTO normalization VALUES('ATLAS_2014_I1279489','d05',0.066,1);
	INSERT INTO normalization VALUES('ATLAS_2014_I1279489','d03',1.10,1);
	INSERT INTO normalization VALUES('ATLAS_2014_I1279489','d04',0.447,1);
	-- these are the integrated cross section of the plot, in pb
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d23',1.45,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d24',1.03,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d25',0.97,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d02',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d03',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d04',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d14',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d26',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d05',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d06',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d07',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d08',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d09',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d10',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d15',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d17',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d18',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d19',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d20',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d21',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d22',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d27',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d11',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d12',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d13',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d16',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_EL','d28',5.59,1);
	-- these are the integrated cross section of the plot, in pb
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d23',1.45,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d24',1.03,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d25',0.97,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d02',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d03',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d04',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d14',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d26',14.96,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d05',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d06',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d07',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d08',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d09',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d10',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d15',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d17',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d18',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d19',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d20',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d21',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d22',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d27',537.10,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d11',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d12',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d13',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d16',5.59,1);
	INSERT INTO normalization VALUES('ATLAS_2015_I1408516_MU','d28',5.59,1);

	-- CMS single jet mass stuff here
	-- WJETS. Cross section in pb
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d52',1.06,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d56',1.06,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d60',1.06,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d64',1.06,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d68',1.06,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d53',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d57',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d61',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d65',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d69',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d72',2.3,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d54',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d58',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d62',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d66',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d70',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d73',0.962,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d55',0.43,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d59',0.43,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d63',0.43,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d67',0.43,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d71',0.43,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_WJET','d74',0.43,1);
	-- ZJETS. Cross section in pb
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d29',0.852,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d33',0.852,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d37',0.852,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d41',0.852,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d45',0.852,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d30',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d34',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d38',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d42',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d46',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d49',1.22,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d31',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d35',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d39',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d43',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d47',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d50',0.377,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d32',0.141,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d36',0.141,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d40',0.141,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d44',0.141,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d48',0.141,1);
	INSERT INTO normalization VALUES('CMS_2013_I1224539_ZJET','d51',0.141,1);
	-- 7 TeV WW cross section in fb
	INSERT INTO normalization VALUES('ATLAS_2013_I1190187','d04-x01-y01',392.6,1);
	-- 8 TeV b hadron cross section in fb
	INSERT INTO normalization VALUES('ATLAS_2017_I1598613','d',17700000,1);
	INSERT INTO normalization VALUES('ATLAS_2017_I1598613_BB','d',17700000,1);
	--
	COMMIT;
	diff --git a/modified_analyses/Analyses/CMS_2017_I1499471.cc b/modified_analyses/Analyses/CMS_2017_I1499471.cc
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/Analyses/CMS_2017_I1499471.cc
	@@ -0,0 +1,295 @@
	+#include "Rivet/Analysis.hh"
	+#include "Rivet/Tools/Cuts.hh"
	+#include "Rivet/Projections/FinalState.hh"
	+#include "Rivet/Projections/VetoedFinalState.hh"
	+#include "Rivet/Projections/DressedLeptons.hh"
	+#include "Rivet/Projections/FastJets.hh"
	+#include "Rivet/Projections/ZFinder.hh"
	+
	+//#define DebugLog
	+
	+namespace Rivet {
	+
	+ class CMS_2017_I1499471 : public Analysis {
	+ public:
	+
	+ /// Constructor
	+ DEFAULT_RIVET_ANALYSIS_CTOR(CMS_2017_I1499471);
	+
	+ /// Book histograms and initialise projections before the run
	+ void init() {
	+
	+#ifdef DebugLog
	+ // set optionally the verbosity for the internal Rivet message system
	+ getLog().setLevel(0);
	+#endif
	+
	+ FinalState fs; ///< @todo No cuts?
	+ VisibleFinalState visfs(fs);
	+
	+ ZFinder zeeFinder(fs, Cuts::abseta < 2.4 && Cuts::pT > 20GeV, PID::ELECTRON, 71.0GeV, 111.0*GeV, 0.1 );
	+ addProjection(zeeFinder, "ZeeFinder");
	+
	+ ZFinder zmumuFinder(fs, Cuts::abseta < 2.4 && Cuts::pT > 20GeV, PID::MUON, 71.0GeV, 111.0*GeV, 0.1 );
	+ addProjection(zmumuFinder, "ZmumuFinder");
	+
	+ VetoedFinalState jetConstits(visfs);
	+ jetConstits.addVetoOnThisFinalState(zeeFinder);
	+ jetConstits.addVetoOnThisFinalState(zmumuFinder);
	+
	+ FastJets akt05Jets(jetConstits, FastJets::ANTIKT, 0.5);
	+ addProjection(akt05Jets, "AntiKt05Jets");
	+
	+ //Histograms booking
	+
	+ _h_first_bjet_pt_b = bookHisto1D(1,1,1);
	+ _h_first_bjet_abseta_b = bookHisto1D(3,1,1);
	+ _h_Z_pt_b = bookHisto1D(5,1,1);
	+ _h_HT_b = bookHisto1D(7,1,1);
	+ _h_Dphi_Zb_b = bookHisto1D(9,1,1);
	+
	+ _h_first_jet_pt_ratio = bookScatter2D(2,1,1);
	+ _h_first_jet_abseta_ratio = bookScatter2D(4,1,1);
	+ _h_Z_pt_ratio = bookScatter2D(6,1,1);
	+ _h_HT_ratio = bookScatter2D(8,1,1);
	+ _h_Dphi_Zj_ratio = bookScatter2D(10,1,1);
	+
	+ _h_first_jet_pt = std::make_shared<YODA::Histo1D>(*_h_first_bjet_pt_b);
	+ _h_first_jet_abseta = std::make_shared<YODA::Histo1D>(*_h_first_bjet_abseta_b);
	+ _h_Z_pt = std::make_shared<YODA::Histo1D>(*_h_Z_pt_b);
	+ _h_HT = std::make_shared<YODA::Histo1D>(*_h_HT_b);
	+ _h_Dphi_Zj = std::make_shared<YODA::Histo1D>(*_h_Dphi_Zb_b);
	+
	+ _h_first_bjet_pt_bb = bookHisto1D(11,1,1);
	+ _h_second_bjet_pt_bb = bookHisto1D(12,1,1);
	+ _h_Z_pt_bb = bookHisto1D(13,1,1);
	+ _h_bb_mass_bb = bookHisto1D(14,1,1);
	+ _h_Zbb_mass_bb = bookHisto1D(15,1,1);
	+ _h_Dphi_bb = bookHisto1D(16,1,1);
	+ _h_DR_bb = bookHisto1D(17,1,1);
	+ _h_DR_Zbmin_bb = bookHisto1D(18,1,1);
	+ _h_A_DR_Zb_bb = bookHisto1D(19,1,1);
	+
	+ _h_bjet_multiplicity = bookHisto1D(20,1,1);
	+
	+ }
	+
	+
	+ /// Perform the per-event analysis
	+ void analyze(const Event& event) {
	+
	+ const ZFinder& zeeFS = applyProjection<ZFinder>(event, "ZeeFinder");
	+ const ZFinder& zmumuFS = applyProjection<ZFinder>(event, "ZmumuFinder");
	+
	+ const Particles& zees = zeeFS.bosons();
	+ const Particles& zmumus = zmumuFS.bosons();
	+
	+ // We did not find exactly one Z. No good.
	+ if (zees.size() + zmumus.size() != 1) {
	+ MSG_DEBUG("Did not find exactly one good Z candidate");
	+ vetoEvent;
	+ }
	+
	+ //event identification depending on mass window
	+ bool ee_event=false;
	+ bool mm_event=false;
	+
	+ if (zees.size() == 1) { ee_event = true; }
	+ if (zmumus.size() == 1) { mm_event = true; }
	+ const Particles& theLeptons = zees.size() ? zeeFS.constituents() : zmumuFS.constituents();
	+
	+ // Cluster jets
	+ // NB. Veto has already been applied on leptons and photons used for dressing
	+ const FastJets& fj = applyProjection<FastJets>(event, "AntiKt05Jets");
	+ const Jets& jets = fj.jetsByPt(Cuts::abseta < 2.4 && Cuts::pT > 30*GeV);
	+
	+ // Perform lepton-jet overlap and HT calculation
	+ double Ht = 0;
	+ Jets goodjets;
	+ foreach (const Jet& j, jets) {
	+ // Decide if this jet is "good", i.e. isolated from the leptons
	+ /// @todo Nice use-case for any() and a C++11 lambda
	+ bool overlap = false;
	+ foreach (const Particle& l, theLeptons) {
	+ if (Rivet::deltaR(j, l) < 0.5) {
	+ overlap = true;
	+ break;
	+ }
	+ }
	+
	+ // Fill HT and good-jets collection
	+ if (overlap) continue;
	+ goodjets.push_back(j);
	+ Ht += j.pT();
	+ }
	+
	+ // We don't care about events with no isolated jets
	+ if (goodjets.empty()) {
	+ MSG_DEBUG("No jets in event");
	+ vetoEvent;
	+ }
	+
	+ Jets jb_final;
	+
	+ //identification of bjets
	+
	+ foreach (const Jet& j, goodjets) {
	+ if ( j.bTagged() ) { jb_final.push_back(j); }
	+ }
	+
	+ //Event weight
	+ const double w = 0.5*event.weight();
	+
	+ //histogram filling
	+
	+ if ((ee_event \|\| mm_event) && goodjets.size() > 0) {
	+
	+ FourMomentum j1(goodjets[0].momentum());
	+
	+ _h_first_jet_pt->fill(j1.pt(),w);
	+ _h_first_jet_abseta->fill(fabs(j1.eta()),w);
	+ if ( ee_event ) _h_Z_pt->fill(zees[0].pt(),w);
	+ if ( mm_event ) _h_Z_pt->fill(zmumus[0].pt(),w);
	+ _h_HT->fill(Ht,w);
	+ if ( ee_event ) _h_Dphi_Zj->fill(deltaPhi(zees[0], j1),w);
	+ if ( mm_event ) _h_Dphi_Zj->fill(deltaPhi(zmumus[0], j1),w);
	+
	+ if ( jb_final.size() > 0 ) {
	+
	+ FourMomentum b1(jb_final[0].momentum());
	+
	+ _h_bjet_multiplicity->fill(1.,w);
	+
	+ _h_first_bjet_pt_b->fill(b1.pt(),w);
	+ _h_first_bjet_abseta_b->fill(fabs(b1.eta()),w);
	+ if ( ee_event ) _h_Z_pt_b->fill(zees[0].pt(),w);
	+ if ( mm_event ) _h_Z_pt_b->fill(zmumus[0].pt(),w);
	+ _h_HT_b->fill(Ht,w);
	+ if ( ee_event ) _h_Dphi_Zb_b->fill(deltaPhi(zees[0], b1.phi()),w);
	+ if ( mm_event ) _h_Dphi_Zb_b->fill(deltaPhi(zmumus[0], b1.phi()),w);
	+
	+ if ( jb_final.size() > 1 ) {
	+
	+ FourMomentum b2(jb_final[1].momentum());
	+
	+ _h_bjet_multiplicity->fill(2.,w);
	+
	+ _h_first_bjet_pt_bb->fill(b1.pt(),w);
	+ _h_second_bjet_pt_bb->fill(b2.pt(),w);
	+ if ( ee_event ) _h_Z_pt_bb->fill(zees[0].pt(),w);
	+ if ( mm_event ) _h_Z_pt_bb->fill(zmumus[0].pt(),w);
	+
	+ FourMomentum bb = add(b1,b2);
	+ FourMomentum Zbb;
	+ if (ee_event) Zbb = add(zees[0],bb);
	+ if (mm_event) Zbb = add(zmumus[0],bb);
	+
	+ _h_bb_mass_bb->fill(bb.mass(),w);
	+ _h_Zbb_mass_bb->fill(Zbb.mass(),w);
	+
	+ _h_Dphi_bb->fill(deltaPhi(b1,b2),w);
	+ _h_DR_bb->fill(deltaR(b1,b2),w);
	+
	+ double DR_Z_b1(0.), DR_Z_b2(0.);
	+ if ( ee_event ) {
	+ DR_Z_b1 = deltaR(zees[0],b1);
	+ DR_Z_b2 = deltaR(zees[0],b2);
	+ }
	+ if ( mm_event ) {
	+ DR_Z_b1 = deltaR(zmumus[0],b1);
	+ DR_Z_b2 = deltaR(zmumus[0],b2);
	+ }
	+
	+ double DR_Zb_min = DR_Z_b1;
	+ double DR_Zb_max = DR_Z_b2;
	+ if ( DR_Zb_min > DR_Zb_max ) {
	+ DR_Zb_min = DR_Z_b2;
	+ DR_Zb_max = DR_Z_b1;
	+ }
	+ double A_Zbb = (DR_Zb_max - DR_Zb_min)/(DR_Zb_max + DR_Zb_min);
	+
	+ _h_DR_Zbmin_bb->fill(DR_Zb_min,w);
	+ _h_A_DR_Zb_bb->fill(A_Zbb,w);
	+
	+ }
	+
	+ }
	+
	+ }
	+
	+ }
	+
	+
	+ /// Normalise histograms etc., after the run
	+ void finalize() {
	+
	+ const double norm = (sumOfWeights() != 0) ? crossSection()/picobarn/sumOfWeights() : 1.0;
	+
	+ MSG_INFO("Cross section = " << std::setfill(' ') << std::setw(14) << std::fixed << std::setprecision(3) << crossSection() << " pb");
	+ MSG_INFO("# Events = " << std::setfill(' ') << std::setw(14) << std::fixed << std::setprecision(3) << numEvents() );
	+ MSG_INFO("SumW = " << std::setfill(' ') << std::setw(14) << std::fixed << std::setprecision(3) << sumOfWeights());
	+ MSG_INFO("Norm factor = " << std::setfill(' ') << std::setw(14) << std::fixed << std::setprecision(6) << norm);
	+
	+ scale( _h_first_bjet_pt_b, 100. );
	+ scale( _h_first_bjet_abseta_b, 100. );
	+ scale( _h_Z_pt_b, 100. );
	+ scale( _h_HT_b, 100. );
	+ scale( _h_Dphi_Zb_b, 100. );
	+
	+ divide( _h_first_bjet_pt_b , _h_first_jet_pt , _h_first_jet_pt_ratio );
	+ divide( _h_first_bjet_abseta_b , _h_first_jet_abseta , _h_first_jet_abseta_ratio );
	+ divide( _h_Z_pt_b , _h_Z_pt , _h_Z_pt_ratio );
	+ divide( _h_HT_b , _h_HT , _h_HT_ratio );
	+ divide( _h_Dphi_Zb_b , _h_Dphi_Zj , _h_Dphi_Zj_ratio );
	+
	+ scale( _h_first_bjet_pt_b, norm/100. );
	+ scale( _h_first_bjet_abseta_b, norm/100. );
	+ scale( _h_Z_pt_b, norm/100. );
	+ scale( _h_HT_b, norm/100. );
	+ scale( _h_Dphi_Zb_b, norm/100. );
	+
	+ scale( _h_first_bjet_pt_bb, norm);
	+ scale( _h_second_bjet_pt_bb, norm);
	+ scale( _h_Z_pt_bb, norm);
	+ scale( _h_bb_mass_bb, norm);
	+ scale( _h_Zbb_mass_bb, norm);
	+ scale( _h_Dphi_bb, norm);
	+ scale( _h_DR_bb, norm);
	+ scale( _h_DR_Zbmin_bb, norm);
	+ scale( _h_A_DR_Zb_bb, norm);
	+
	+ scale( _h_bjet_multiplicity, norm );
	+
	+ }
	+
	+
	+ private:
	+
	+ /// @name Histograms
	+
	+ Histo1DPtr _h_first_jet_pt, _h_first_bjet_pt_b;
	+ Histo1DPtr _h_first_jet_abseta, _h_first_bjet_abseta_b;
	+ Histo1DPtr _h_Z_pt, _h_Z_pt_b;
	+ Histo1DPtr _h_HT, _h_HT_b;
	+ Histo1DPtr _h_Dphi_Zj, _h_Dphi_Zb_b;
	+
	+ Scatter2DPtr _h_first_jet_pt_ratio;
	+ Scatter2DPtr _h_first_jet_abseta_ratio;
	+ Scatter2DPtr _h_Z_pt_ratio;
	+ Scatter2DPtr _h_HT_ratio;
	+ Scatter2DPtr _h_Dphi_Zj_ratio;
	+
	+ Histo1DPtr _h_first_bjet_pt_bb, _h_second_bjet_pt_bb;
	+ Histo1DPtr _h_Z_pt_bb;
	+ Histo1DPtr _h_bb_mass_bb, _h_Zbb_mass_bb;
	+ Histo1DPtr _h_Dphi_bb, _h_DR_bb, _h_DR_Zbmin_bb, _h_A_DR_Zb_bb;
	+
	+ Histo1DPtr _h_bjet_multiplicity;
	+
	+ };
	+
	+
	+ // The hook for the plugin system
	+ DECLARE_RIVET_PLUGIN(CMS_2017_I1499471);
	+
	+}
	diff --git a/modified_analyses/Analyses/CMS_2017_I1610623.cc b/modified_analyses/Analyses/CMS_2017_I1610623.cc
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/Analyses/CMS_2017_I1610623.cc
	@@ -0,0 +1,262 @@
	+#include "Rivet/Analysis.hh"
	+#include "Rivet/Projections/FinalState.hh"
	+
	+#include "Rivet/Tools/Logging.hh"
	+#include "Rivet/Projections/FastJets.hh"
	+#include "Rivet/Projections/VetoedFinalState.hh"
	+#include "Rivet/Projections/WFinder.hh"
	+
	+#include "Rivet/AnalysisLoader.hh"
	+#include "Rivet/AnalysisInfo.hh"
	+#include "Rivet/Tools/RivetYODA.hh"
	+
	+#include <iostream>
	+
	+namespace Rivet {
	+
	+
	+ /// @brief Add a short analysis description here
	+ class CMS_2017_I1610623 : public Analysis {
	+ public:
	+
	+ /// Constructor
	+ DEFAULT_RIVET_ANALYSIS_CTOR(CMS_2017_I1610623);
	+
	+
	+ /// @name Analysis methods
	+ //@{
	+
	+ /// Book histograms and initialise projections before the run
	+ void init() {
	+
	+ // Initialise and register projections
	+ FinalState fs;
	+ WFinder wfinder_mu(fs, Cuts::abseta < 2.4 && Cuts::pT > 0GeV, PID::MUON, 0GeV, 1000000GeV, 0GeV, 0.1, WFinder::CLUSTERNODECAY, WFinder::TRACK, WFinder::TRANSMASS);
	+ //WFinder wfinder_mu(fs, Cuts::abseta < 2.4 && Cuts::pT > 0GeV, PID::MUON, 0GeV, 1000000GeV, 0GeV, 0.1, WFinder::CLUSTERNODECAY, WFinder::NOTRACK, WFinder::TRANSMASS);
	+ addProjection(wfinder_mu, "WFinder_mu");
	+
	+ // Define veto FS
	+ VetoedFinalState vfs;
	+ vfs.addVetoOnThisFinalState(wfinder_mu);
	+ vfs.addVetoPairId(PID::MUON);
	+ vfs.vetoNeutrinos();
	+
	+ FastJets fastjets(vfs, FastJets::ANTIKT, 0.4);
	+ addProjection(fastjets, "Jets");
	+
	+ //-------------
	+ _hist_Mult_exc = bookHisto1D("d01-x01-y01");
	+ _hist_inc_WJetMult = bookHisto1D("d02-x01-y01");
	+
	+ //-------------
	+ _hist_JetPt1j = bookHisto1D("d03-x01-y01");
	+ _hist_JetPt2j = bookHisto1D("d04-x01-y01");
	+ _hist_JetPt3j = bookHisto1D("d05-x01-y01");
	+ _hist_JetPt4j = bookHisto1D("d06-x01-y01");
	+
	+ //-------------
	+ _hist_JetRap1j = bookHisto1D("d07-x01-y01");
	+ _hist_JetRap2j = bookHisto1D("d08-x01-y01");
	+ _hist_JetRap3j = bookHisto1D("d09-x01-y01");
	+ _hist_JetRap4j = bookHisto1D("d10-x01-y01");
	+
	+ //-------------
	+ _hist_Ht_1j = bookHisto1D("d11-x01-y01");
	+ _hist_Ht_2j = bookHisto1D("d12-x01-y01");
	+ _hist_Ht_3j = bookHisto1D("d13-x01-y01");
	+ _hist_Ht_4j = bookHisto1D("d14-x01-y01");
	+
	+ //-------------
	+ _hist_dphij1mu_1j =bookHisto1D("d15-x01-y01");
	+ _hist_dphij2mu_2j =bookHisto1D("d16-x01-y01");
	+ _hist_dphij3mu_3j =bookHisto1D("d17-x01-y01");
	+ _hist_dphij4mu_4j =bookHisto1D("d18-x01-y01");
	+
	+ //-------------
	+ _hist_dRmuj_1j =bookHisto1D("d19-x01-y01");
	+
	+ }
	+
	+ // define function used for filiing inc Njets histo
	+ void Fill(Histo1DPtr& _histJetMult, const double& weight, std::vector<FourMomentum>& finaljet_list){
	+ _histJetMult->fill(0, weight);
	+ for (size_t i=0 ; i<finaljet_list.size() ; ++i) {
	+ if (i==6) break;
	+ _histJetMult->fill(i+1, weight); // inclusive multiplicity
	+ }
	+ }
	+
	+ /// Perform the per-event analysis
	+ void analyze(const Event& event) {
	+
	+ /// @todo Do the event by event analysis here
	+ const double weight = event.weight();
	+ const WFinder& wfinder_mu = applyProjection<WFinder>(event, "WFinder_mu");
	+
	+ if (wfinder_mu.bosons().size() != 1) {
	+ vetoEvent;
	+ }
	+
	+ if (wfinder_mu.bosons().size() == 1) {
	+
	+ const FourMomentum lepton0 = wfinder_mu.constituentLepton().momentum();
	+ const FourMomentum neutrino = wfinder_mu.constituentNeutrino().momentum();
	+ double WmT = wfinder_mu.mT();
	+
	+ if (WmT < 50.0*GeV) vetoEvent;
	+
	+ double pt0 = lepton0.pT();
	+ double eta0 = lepton0.eta();
	+
	+ if ( (fabs(eta0) > 2.4) \|\| (pt0 < 25.0*GeV) ) vetoEvent;
	+
	+ // Obtain the jets.
	+ vector<FourMomentum> finaljet_list;
	+ vector<FourMomentum> jet100_list;
	+ double HT = 0.0;
	+
	+ // loop over jets in an event, pushback in finaljet_list collection
	+ foreach (const Jet& j, applyProjection<FastJets>(event, "Jets").jetsByPt(30.0*GeV)) {
	+ const double jrap = j.momentum().rap();
	+ const double jpt = j.momentum().pT();
	+ if ( (fabs(jrap) < 2.4) && (deltaR(lepton0, j.momentum()) > 0.4) ) {
	+ if(jpt > 30.0*GeV) {
	+ finaljet_list.push_back(j.momentum());
	+ HT += j.momentum().pT();
	+ }
	+ if(jpt > 100.0*GeV) {
	+ jet100_list.push_back(j.momentum());
	+ }
	+ }
	+ } // end looping over jets
	+
	+ //---------------------- FILL HISTOGRAMS ------------------
	+
	+ // Multiplicity exc plot.
	+ _hist_Mult_exc->fill(finaljet_list.size(), weight);
	+
	+ // Multiplicity inc plot.
	+ Fill(_hist_inc_WJetMult, weight, finaljet_list);
	+
	+ // dRmuj plot.
	+ double mindR(99999);
	+ if(jet100_list.size()>=1) {
	+ for (unsigned ji = 0; ji < jet100_list.size(); ji++){
	+ double dr_(9999);
	+ dr_ = fabs(deltaR(lepton0, jet100_list[ji]));
	+ if (dr_ < mindR){
	+ mindR = dr_;
	+ }
	+ }
	+ if(jet100_list[0].pT() > 300.0*GeV){
	+ _hist_dRmuj_1j->fill(mindR, weight);
	+ }
	+ }
	+
	+ if(finaljet_list.size()>=1) {
	+ _hist_JetPt1j->fill(finaljet_list[0].pT(), weight);
	+ _hist_JetRap1j->fill(fabs(finaljet_list[0].rap()), weight);
	+ _hist_Ht_1j->fill(HT, weight);
	+ _hist_dphij1mu_1j->fill(deltaPhi(finaljet_list[0].phi(), lepton0.phi()), weight);
	+ }
	+
	+ if(finaljet_list.size()>=2) {
	+ _hist_JetPt2j->fill(finaljet_list[1].pT(), weight);
	+ _hist_JetRap2j->fill(fabs(finaljet_list[1].rap()), weight);
	+ _hist_Ht_2j->fill(HT, weight);
	+ _hist_dphij2mu_2j->fill(deltaPhi(finaljet_list[1].phi(), lepton0.phi()), weight);
	+ }
	+
	+ if(finaljet_list.size()>=3) {
	+ _hist_JetPt3j->fill(finaljet_list[2].pT(), weight);
	+ _hist_JetRap3j->fill(fabs(finaljet_list[2].rap()), weight);
	+ _hist_Ht_3j->fill(HT, weight);
	+ _hist_dphij3mu_3j->fill(deltaPhi(finaljet_list[2].phi(), lepton0.phi()), weight);
	+ }
	+
	+ if(finaljet_list.size()>=4) {
	+ _hist_JetPt4j->fill(finaljet_list[3].pT(), weight);
	+ _hist_JetRap4j->fill(fabs(finaljet_list[3].rap()), weight);
	+ _hist_Ht_4j->fill(HT, weight);
	+ _hist_dphij4mu_4j->fill(deltaPhi(finaljet_list[3].phi(), lepton0.phi()), weight);
	+ }
	+ } // close the Wboson loop
	+
	+ } //void loop
	+
	+
	+ /// Normalise histograms etc., after the run
	+ void finalize() {
	+ const double crossec = !std::isnan(crossSectionPerEvent()) ? crossSection() : 61526.7*picobarn;
	+ if (std::isnan(crossSectionPerEvent())){
	+ MSG_INFO("No valid cross-section given, using NNLO xsec calculated by FEWZ " << crossec/picobarn << " pb");
	+ }
	+
	+ scale(_hist_Mult_exc, crossec/picobarn/sumOfWeights());
	+ scale(_hist_inc_WJetMult, crossec/picobarn/sumOfWeights());
	+
	+ scale(_hist_JetPt1j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetPt2j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetPt3j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetPt4j, crossec/picobarn/sumOfWeights());
	+
	+ scale(_hist_JetRap1j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetRap2j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetRap3j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_JetRap4j, crossec/picobarn/sumOfWeights());
	+
	+ scale(_hist_Ht_1j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_Ht_2j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_Ht_3j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_Ht_4j, crossec/picobarn/sumOfWeights());
	+
	+ scale(_hist_dphij1mu_1j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_dphij2mu_2j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_dphij3mu_3j, crossec/picobarn/sumOfWeights());
	+ scale(_hist_dphij4mu_4j, crossec/picobarn/sumOfWeights());
	+
	+ scale(_hist_dRmuj_1j, crossec/picobarn/sumOfWeights());
	+
	+ }
	+
	+ //@}
	+
	+
	+ private:
	+
	+ /// @name Histograms
	+ //@{
	+ Histo1DPtr _hist_Mult_exc;
	+ Histo1DPtr _hist_inc_WJetMult;
	+
	+ Histo1DPtr _hist_JetPt1j;
	+ Histo1DPtr _hist_JetPt2j;
	+ Histo1DPtr _hist_JetPt3j;
	+ Histo1DPtr _hist_JetPt4j;
	+
	+ Histo1DPtr _hist_JetRap1j;
	+ Histo1DPtr _hist_JetRap2j;
	+ Histo1DPtr _hist_JetRap3j;
	+ Histo1DPtr _hist_JetRap4j;
	+
	+ Histo1DPtr _hist_Ht_1j;
	+ Histo1DPtr _hist_Ht_2j;
	+ Histo1DPtr _hist_Ht_3j;
	+ Histo1DPtr _hist_Ht_4j;
	+
	+ Histo1DPtr _hist_dphij1mu_1j;
	+ Histo1DPtr _hist_dphij2mu_2j;
	+ Histo1DPtr _hist_dphij3mu_3j;
	+ Histo1DPtr _hist_dphij4mu_4j;
	+
	+ Histo1DPtr _hist_dRmuj_1j;
	+ //@}
	+
	+ };
	+
	+
	+ // The hook for the plugin system
	+ DECLARE_RIVET_PLUGIN(CMS_2017_I1610623);
	+
	+
	+}
	diff --git a/modified_analyses/Analyses/buildrivet.sh b/modified_analyses/Analyses/buildrivet.sh
	old mode 100644
	new mode 100755
	diff --git a/modified_analyses/anainfo/CMS_2017_I1499471.info b/modified_analyses/anainfo/CMS_2017_I1499471.info
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/anainfo/CMS_2017_I1499471.info
	@@ -0,0 +1,56 @@
	+Name: CMS_2017_I1499471
	+Year: 2017
	+Summary: Measurements of the associated production of a Z boson and b jets in pp collisions at $\sqrt{s} = 8$ TeV
	+Experiment: CMS
	+Collider: LHC
	+SpiresID: I1499471
	+Status: VALIDATED
	+Authors:
	+ - Fabio Cossutti (fabio.cossutti@ts.infn.it)
	+References:
	+ - DOI:10.1140/epjc/s10052-017-5140-y
	+ - arXiv:1611.06507
	+ - CMS-SMP-14-010
	+
	+RunInfo: 'Run MC generators with Z decaying leptonically into both electrons and muons at 8 TeV CoM energy.
	+For jet-flavour inclusive samples, order of 20 million unweighted events can give a reasonable global comparison,
	+but precision in the high jet multiplicity region/high jet pt may require substantially larger samples or statistical
	+enhancement of high jet multiplicities.'
	+NumEvents: 20000000
	+NeedCrossSection: yes
	+Beams: [p+, p+]
	+Energies: [8000]
	+Description:
	+'Measurements of the associated production of a Z boson with at least
	+one jet originating from a b quark in proton-proton collisions
	+at $\sqrt{s} = 8$ TeV are presented. Differential cross sections are
	+measured with data collected by the CMS experiment corresponding to an
	+integrated luminosity of 19.8 $\text{fb}^{-1}$. Z bosons are reconstructed
	+through their decays to electrons and muons. Cross sections are
	+measured as a function of observables characterizing the kinematics of
	+the b jet and the Z boson. Ratios of differential cross
	+sections for the associated production with at least one b jet
	+to the associated production with any jet are also presented. The
	+production of a Z boson with at least two b jets is investigated,
	+and differential cross sections are measured for the dijet system.
	+Results are compared to theoretical predictions, testing two different
	+flavour schemes for the choice of initial-state partons. Note of the authors:
	+the Z boson selection and b-tagging algorithms used in the code are based on
	+standard Rivet code and do not exactly coincide with those used in the paper.
	+They anyway reproduce the paper results at the percent level or better for most
	+of the bins.'
	+BibKey: Khachatryan:2016iob
	+BibTeX: '@article{Khachatryan:2016iob,
	+ author = "Khachatryan, Vardan and others",
	+ title = "{Measurements of the associated production of a Z boson
	+ and b jets in pp collisions at $\sqrt{s}$ = 8 TeV}",
	+ collaboration = "CMS",
	+ journal = "Submitted to: Eur. Phys. J. C",
	+ year = "2016",
	+ eprint = "1611.06507",
	+ archivePrefix = "arXiv",
	+ primaryClass = "hep-ex",
	+ reportNumber = "CMS-SMP-14-010, CERN-EP-2016-254",
	+ SLACcitation = "%%CITATION = ARXIV:1611.06507;%%"
	+}'
	+
	diff --git a/modified_analyses/anainfo/CMS_2017_I1610623.info b/modified_analyses/anainfo/CMS_2017_I1610623.info
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/anainfo/CMS_2017_I1610623.info
	@@ -0,0 +1,42 @@
	+Name: CMS_2017_I1610623
	+Year: 2017
	+Summary: Measurements of differential cross sections for the associated production of a W boson and jets in proton-proton collisions at sqrt(s) = 13 TeV
	+Experiment: CMS
	+Collider: LHC
	+SpiresID: 1610623
	+Status: VALIDATED
	+Authors:
	+ - Emanuela.Barberis@cern.ch
	+ - apichart.hortiangtham@cern.ch
	+ - Kadir.Ocalan@cern.ch
	+References:
	+ - DOI:10.1103/PhysRevD.96.072005
	+ - arXiv:1707.05979
	+ - CMS-SMP-16-005
	+RunInfo: 13 TeV $pp \to W+jets$.
	+NumEvents: 1000000
	+NeedCrossSection: True
	+Beams: [p+, p+]
	+Energies: [13000]
	+Description:
	+ 'Differential cross sections for a W boson produced in association with jets are measured in a data sample of proton-proton collisions at a center-of-mass energy of 13 TeV recorded with the CMS detector and corresponding to an integrated luminosity of 2.2 inverse femtobarns. The W bosons are identified through their decay into a muon and a neutrino. The cross sections are reported as functions of jet multiplicity, jet transverse momenta, jet rapidity $\|y\|$, the scalar sum of jet transverse momenta ($H_T$), and angular correlations between the muon and each jet for different jet multiplicities. The cross sections are measured in the fiducial region defined by a muon with $p_T > 25$ GeV and pseudorapidity $\|\eta\| < 2.4$, and by a transverse mass between the muon and the missing transverse energy $M_T > 50$ GeV. Jets are reconstructed using the anti-kT algorithm with a distance parameter R = 0.4, and only jets with $p_T > 30$ GeV, $\|y\| < 2.4$, and a separation of $\Delta R > 0.4$ from the muon are considered. In addition, the differential cross section is measured as a function of the angular distance (DeltaR) between the muon and the closest jet for events with one or more jets, requiring jets with $p_T > 100$ GeV and the leading jet with $p_T > 300$ GeV.'
	+BibKey: Sirunyan:2017wgx
	+BibTeX: '@article{Sirunyan:2017wgx,
	+ author = "Sirunyan, Albert M and others",
	+ title = "{Measurement of the differential cross sections for the
	+ associated production of a $W$ boson and jets in
	+ proton-proton collisions at $\sqrt{s}=13$ TeV}",
	+ collaboration = "CMS",
	+ journal = "Phys. Rev.",
	+ volume = "D96",
	+ year = "2017",
	+ number = "7",
	+ pages = "072005",
	+ doi = "10.1103/PhysRevD.96.072005",
	+ eprint = "1707.05979",
	+ archivePrefix = "arXiv",
	+ primaryClass = "hep-ex",
	+ reportNumber = "CMS-SMP-16-005, CERN-EP-2017-142",
	+ SLACcitation = "%%CITATION = ARXIV:1707.05979;%%"
	+}'
	+
	diff --git a/modified_analyses/anainfo/analysis-info-files-go-here b/modified_analyses/anainfo/analysis-info-files-go-here
	new file mode 100644
	diff --git a/modified_analyses/plotinfo/CMS_2017_I1499471.plot b/modified_analyses/plotinfo/CMS_2017_I1499471.plot
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/plotinfo/CMS_2017_I1499471.plot
	@@ -0,0 +1,177 @@
	+# BEGIN PLOT /CMS_2017_I1499471/d01-x01-y01
	+Title=$\text{CMS, 8 TeV, Leading b jet transverse momentum, at least one b jet}$
	+XLabel=$p_{t}(b)$ GeV
	+YLabel=$d\sigma_{fiducial}/dp_{t}(b)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d02-x01-y01
	+Title=$\text{CMS, 8 TeV, Leading b\/inclusive jet transverse momentum ratio, at least one b jet over inclusive}$
	+XLabel=$p_{t}(b/j)$ GeV
	+YLabel=b jet/inclusive jets $d\sigma_{fiducial}/dp_{t}(b/j)$ ratio \%
	+LogY = 0
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d03-x01-y01
	+Title=$\text{CMS, 8 TeV, Leading b jet absolute pseudorapidity, at least one b jet}$
	+XLabel=$\eta(b)$
	+YLabel=$d\sigma_{fiducial}/d\eta(b)$ pb
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d04-x01-y01
	+Title=$\text{CMS, 8 TeV, Leading b/inclusive jet absolute pseudorapidity ratio, at least one b jet over inclusive}$
	+XLabel=$\eta(b/j)$
	+YLabel=b jet/inclusive jets $d\sigma_{fiducial}/d\eta(b/j)$ ratio \%
	+LogY = 0
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d05-x01-y01
	+Title=$\text{CMS, 8 TeV, Z boson pt, at least one b jet}$
	+XLabel=$p_{t}(Z)$ GeV
	+YLabel=$d\sigma_{fiducial}/dp_{t}(Z)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d06-x01-y01
	+Title=$\text{CMS, 8 TeV, Z boson pt ratio, at least one b jet over inclusive}$
	+XLabel=$p_{t}(Z)$ GeV
	+YLabel=b jet/inclusive jets $d\sigma_{fiducial}/dp_{t}(Z)$ ratio \%
	+LogY = 0
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d07-x01-y01
	+Title=$\text{CMS, 8 TeV, HT, at least one b jet}$
	+XLabel=$H_{t}$ GeV
	+YLabel=$d\sigma_{fiducial}/dH_{t}$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d08-x01-y01
	+Title=$\text{CMS, 8 TeV, HT ratio, at least one b jet over inclusive}$
	+XLabel=$H_{t}$ GeV
	+YLabel=b jet/inclusive jets $d\sigma_{fiducial}/dH_{t}$ ratio \%
	+LogY = 0
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d09-x01-y01
	+Title=$\text{CMS, 8 TeV, DeltaPhi_Zb, at least one b jet}$
	+XLabel=$d\phi(Zb)$ rad
	+YLabel=$d\sigma_{fiducial}/d\phi(Zb)$ pb/rad
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d10-x01-y01
	+Title=$\text{CMS, 8 TeV, DeltaPhi_Z(b/j) ratio, at least one b jet over inclusive}$
	+XLabel=$d\phi(Zb/j)$
	+YLabel=b jet/inclusive jets $d\sigma_{fiducial}/d\phi(Zb/j)$ ratio \%
	+LogY = 0
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d11-x01-y01
	+Title=$\text{CMS, 8 TeV, Leading b jet transverse momentum, at least two b jets}$
	+XLabel=$p_{t}(b)$ GeV
	+YLabel=$d\sigma_{fiducial}/dp_{t}(b)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d12-x01-y01
	+Title=$\text{CMS, 8 TeV, Subleading b jet transverse momentum, at least two b jets}$
	+XLabel=$p_{t}(2nd b)$ GeV
	+YLabel=$d\sigma_{fiducial}/dp_{t}(2nd b)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d13-x01-y01
	+Title=$\text{CMS, 8 TeV, Z boson pt, at least two b jets}$
	+XLabel=$p_{t}(Z)$ GeV
	+YLabel=$d\sigma_{fiducial}/dp_{t}(Z)$ pb/GeV
	+LogY = 1
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d14-x01-y01
	+Title=$\text{CMS, 8 TeV, bb system mass, at least two b jets}$
	+XLabel=$M(bb)$ GeV
	+YLabel=$d\sigma_{fiducial}/dM(bb)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d15-x01-y01
	+Title=$\text{CMS, 8 TeV, Zbb system mass, at least two b jets}$
	+XLabel=$M(Zbb)$ GeV
	+YLabel=$d\sigma_{fiducial}/dM(Zbb)$ pb/GeV
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d16-x01-y01
	+Title=$\text{CMS, 8 TeV, DeltaPhi bb, at least two b jets}$
	+XLabel=$d\phi(bb)$ rad
	+YLabel=$d\sigma_{fiducial}/d\phi(bb)$ pb/rad
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d17-x01-y01
	+Title=$\text{CMS, 8 TeV, DeltaR bb, at least two b jets}$
	+XLabel=$dR(bb)$
	+YLabel=$d\sigma_{fiducial}/dR(bb)$ pb
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d18-x01-y01
	+Title=$\text{CMS, 8 TeV, DeltaR Zb minimum, at least two b jets}$
	+XLabel=$min dR(Zb)$
	+YLabel=$d\sigma_{fiducial}/d(min\, dR(Zb))$ pb
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d19-x01-y01
	+Title=$\text{CMS, 8 TeV, Asymmetry dR Zb max dR Zb min , at least two b jets}$
	+XLabel=$A_{Zb}(max\, dR(Zb)\, min\, dR(Zb))$
	+YLabel=$d\sigma_{fiducial}/dA_{Zb}$ pb
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1499471/d20-x01-y01
	+Title=$\text{CMS, 8 TeV, Inclusive b jet multiplicity}$
	+XLabel=$n_{jet}$
	+YLabel=$d\sigma_{fiducial}/dn_{jet}$ pb
	+LogY = 1
	+LegendXPos=0.7
	+LegendYPos=0.9
	+# END PLOT
	diff --git a/modified_analyses/plotinfo/CMS_2017_I1610623.plot b/modified_analyses/plotinfo/CMS_2017_I1610623.plot
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/plotinfo/CMS_2017_I1610623.plot
	@@ -0,0 +1,113 @@
	+# BEGIN PLOT /CMS_2017_I1610623/d01-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=EXCLUSIVE JET MULTIPLICITY
	+YLabel=$\sigma(W\rightarrow \mu \nu + N \text{ jets})$ [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d02-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=INCLUSIVE JET MULTIPLICITY
	+YLabel=$\sigma(W\rightarrow \mu \nu + \ge N \text{ jets})$ [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d03-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=FIRST LEADING JET PT [GEV]
	+YLabel=D(SIG)/D(PT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d04-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=SECOND LEADING JET PT [GEV]
	+YLabel=D(SIG)/D(PT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d05-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=THIRD LEADING JET PT [GEV]
	+YLabel=D(SIG)/D(PT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d06-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=FOURTH LEADING JET PT [GEV]
	+YLabel=D(SIG)/D(PT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d07-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=FIRST LEADING JET ABS(YRAP)
	+YLabel=D(SIG)/D(ABS(YRAP)) [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d08-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=SECOND LEADING JET ABS(YRAP)
	+YLabel=D(SIG)/D(ABS(YRAP)) [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d09-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=THIRD LEADING JET ABS(YRAP)
	+YLabel=D(SIG)/D(ABS(YRAP)) [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d10-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=FOURTH LEADING JET ABS(YRAP)
	+YLabel=D(SIG)/D(ABS(YRAP)) [PB]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d11-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=HT ($\ge 1$ JET) [GEV]
	+YLabel=D(SIG)/D(HT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d12-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=HT ($\ge 2$ JET) [GEV]
	+YLabel=D(SIG)/D(HT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d13-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=HT ($\ge 3$ JET) [GEV]
	+YLabel=D(SIG)/D(HT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d14-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=HT ($\ge 4$ JET) [GEV]
	+YLabel=D(SIG)/D(HT) [PB/GEV]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d15-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=DELTA PHI (FIRST LEADING JET, MUON) [RAD]
	+YLabel=D(SIG)/D(DELTAPHI) [PB/RAD]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d16-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=DELTA PHI (SECOND LEADING JET, MUON) [RAD]
	+YLabel=D(SIG)/D(DELTAPHI) [PB/RAD]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d17-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=DELTA PHI (THIRD LEADING JET, MUON) [RAD]
	+YLabel=D(SIG)/D(DELTAPHI) [PB/RAD]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d18-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=DELTA PHI (FOURTH LEADING JET, MUON) [RAD]
	+YLabel=D(SIG)/D(DELTAPHI) [PB/RAD]
	+# END PLOT
	+
	+# BEGIN PLOT /CMS_2017_I1610623/d19-x01-y01
	+Title=CMS, sqrt(s) = 13 TeV, W+jets
	+XLabel=DELTA R (CLOSEST JET, MUON) ($\ge 1$ JET)
	+YLabel=D(SIG)/D(DELTAR) [PB]
	+# END PLOT
	diff --git a/modified_analyses/refdata/CMS_2017_I1499471.yoda b/modified_analyses/refdata/CMS_2017_I1499471.yoda
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/refdata/CMS_2017_I1499471.yoda
	@@ -0,0 +1,401 @@
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d01-x01-y01
	+Path=/REF/CMS_2017_I1499471/d01-x01-y01
	+Title=CMS, 8 TeV, Leading b jet transverse momentum, at least one b jet
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 36.700000 6.700000 6.700000 0.1075 0.0085 0.0085
	+ 50.100000 6.700000 6.700000 0.0596 0.0043 0.0043
	+ 63.499999 6.700000 6.700000 0.0347 0.0021 0.0021
	+ 76.899999 6.700000 6.700000 0.0210 0.0022 0.0022
	+ 90.299999 6.700000 6.700000 0.0125 0.0012 0.0012
	+ 103.699998 6.700000 6.700000 0.00805 0.00073 0.00073
	+ 117.099998 6.700000 6.700000 0.00516 0.00056 0.00056
	+ 130.499998 6.700000 6.700000 0.00347 0.00035 0.00035
	+ 143.899998 6.700000 6.700000 0.00256 0.00023 0.00023
	+ 157.299997 6.700000 6.700000 0.00162 0.00016 0.00016
	+ 170.699997 6.700000 6.700000 0.00125 0.00027 0.00027
	+ 190.799997 13.400000 13.400000 0.000806 0.000089 0.000089
	+ 224.299996 20.100000 20.100000 0.000357 0.000074 0.000074
	+ 277.899995 33.499999 33.499999 0.000131 0.000031 0.000031
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d02-x01-y01
	+Path=/REF/CMS_2017_I1499471/d02-x01-y01
	+Title=CMS, 8 TeV, Leading b/inclusive jet transverse momentum ratio, at least one b jet over inclusive
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 36.700000 6.700000 6.700000 4.72 0.46 0.46
	+ 50.100000 6.700000 6.700000 5.00 0.36 0.36
	+ 63.499999 6.700000 6.700000 5.05 0.26 0.26
	+ 76.899999 6.700000 6.700000 5.11 0.43 0.43
	+ 90.299999 6.700000 6.700000 4.80 0.48 0.48
	+ 103.699998 6.700000 6.700000 4.70 0.33 0.33
	+ 117.099998 6.700000 6.700000 4.52 0.48 0.48
	+ 130.499998 6.700000 6.700000 4.41 0.40 0.40
	+ 143.899998 6.700000 6.700000 4.54 0.36 0.36
	+ 157.299997 6.700000 6.700000 4.05 0.39 0.39
	+ 170.699997 6.700000 6.700000 4.23 0.89 0.89
	+ 190.799997 13.400000 13.400000 4.14 0.41 0.41
	+ 224.299996 20.100000 20.100000 3.48 0.68 0.68
	+ 277.899995 33.499999 33.499999 3.09 0.80 0.80
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d03-x01-y01
	+Path=/REF/CMS_2017_I1499471/d03-x01-y01
	+Title=CMS, 8 TeV, Leading b jet absolute pseudorapidity, at least one b jet
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.156250 0.156250 0.156250 1.96 0.18 0.18
	+ 0.468750 0.156250 0.156250 1.89 0.21 0.21
	+ 0.781250 0.156250 0.156250 1.85 0.14 0.14
	+ 1.093750 0.156250 0.156250 1.66 0.10 0.10
	+ 1.406250 0.156250 0.156250 1.428 0.097 0.097
	+ 1.718750 0.156250 0.156250 1.197 0.096 0.096
	+ 2.031250 0.156250 0.156250 0.948 0.097 0.097
	+ 2.343750 0.156250 0.156250 0.488 0.058 0.058
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d04-x01-y01
	+Path=/REF/CMS_2017_I1499471/d04-x01-y01
	+Title=CMS, 8 TeV, Leading b/inclusive jet absolute pseudorapidity ratio, at least one b jet over inclusive
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.156250 0.156250 0.156250 4.91 0.44 0.44
	+ 0.468750 0.156250 0.156250 4.86 0.50 0.50
	+ 0.781250 0.156250 0.156250 5.04 0.36 0.36
	+ 1.093750 0.156250 0.156250 4.95 0.30 0.30
	+ 1.406250 0.156250 0.156250 4.95 0.37 0.37
	+ 1.718750 0.156250 0.156250 4.93 0.40 0.40
	+ 2.031250 0.156250 0.156250 4.82 0.50 0.50
	+ 2.343750 0.156250 0.156250 4.55 0.54 0.54
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d05-x01-y01
	+Path=/REF/CMS_2017_I1499471/d05-x01-y01
	+Title=CMS, 8 TeV, Z boson pt, at least one b jet
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 5.000000 5.000000 5.000000 0.0095 0.0023 0.0023
	+ 15.000000 5.000000 5.000000 0.0226 0.0037 0.0037
	+ 25.000000 5.000000 5.000000 0.0365 0.0047 0.0047
	+ 35.000000 5.000000 5.000000 0.0508 0.0046 0.0046
	+ 45.000000 5.000000 5.000000 0.0530 0.0042 0.0042
	+ 55.000000 5.000000 5.000000 0.0448 0.0049 0.0049
	+ 65.000000 5.000000 5.000000 0.0337 0.0024 0.0024
	+ 75.000000 5.000000 5.000000 0.0253 0.0020 0.0020
	+ 85.000000 5.000000 5.000000 0.0180 0.0012 0.0012
	+ 95.000000 5.000000 5.000000 0.01330 0.00086 0.00086
	+ 105.000000 5.000000 5.000000 0.00985 0.00062 0.00062
	+ 115.000000 5.000000 5.000000 0.0073 0.0010 0.0010
	+ 125.000000 5.000000 5.000000 0.00559 0.00057 0.00057
	+ 135.000000 5.000000 5.000000 0.00395 0.00030 0.00030
	+ 150.000000 10.000000 10.000000 0.00270 0.00027 0.00027
	+ 170.000000 10.000000 10.000000 0.00154 0.00026 0.00026
	+ 195.000000 15.000000 15.000000 0.00091 0.00017 0.00017
	+ 230.000000 20.000000 20.000000 0.000416 0.000088 0.000088
	+ 280.000000 30.000000 30.000000 0.000179 0.000033 0.000033
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d06-x01-y01
	+Path=/REF/CMS_2017_I1499471/d06-x01-y01
	+Title=CMS, 8 TeV, Z boson pt ratio, at least one b jet over inclusive
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 5.000000 5.000000 5.000000 4.9 1.5 1.5
	+ 15.000000 5.000000 5.000000 4.96 0.90 0.90
	+ 25.000000 5.000000 5.000000 4.29 0.48 0.48
	+ 35.000000 5.000000 5.000000 3.95 0.35 0.35
	+ 45.000000 5.000000 5.000000 4.40 0.34 0.34
	+ 55.000000 5.000000 5.000000 4.84 0.50 0.50
	+ 65.000000 5.000000 5.000000 5.11 0.34 0.34
	+ 75.000000 5.000000 5.000000 5.49 0.40 0.40
	+ 85.000000 5.000000 5.000000 5.66 0.34 0.34
	+ 95.000000 5.000000 5.000000 5.98 0.34 0.34
	+ 105.000000 5.000000 5.000000 6.16 0.35 0.35
	+ 115.000000 5.000000 5.000000 6.11 0.63 0.63
	+ 125.000000 5.000000 5.000000 6.41 0.62 0.62
	+ 135.000000 5.000000 5.000000 5.95 0.43 0.43
	+ 150.000000 10.000000 10.000000 6.00 0.58 0.58
	+ 170.000000 10.000000 10.000000 5.6 1.3 1.3
	+ 195.000000 15.000000 15.000000 5.7 1.1 1.1
	+ 230.000000 20.000000 20.000000 5.6 1.2 1.2
	+ 280.000000 30.000000 30.000000 5.8 1.3 1.3
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d07-x01-y01
	+Path=/REF/CMS_2017_I1499471/d07-x01-y01
	+Title=CMS, 8 TeV, HT, at least one b jet
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 39.700000 9.700000 9.700000 0.0680 0.0074 0.0074
	+ 59.100000 9.700000 9.700000 0.0281 0.0031 0.0031
	+ 78.500000 9.700000 9.700000 0.0222 0.0020 0.0020
	+ 97.900000 9.700000 9.700000 0.0152 0.0014 0.0014
	+ 117.300000 9.700000 9.700000 0.0106 0.0010 0.0010
	+ 136.700000 9.700000 9.700000 0.00786 0.00069 0.00069
	+ 156.100000 9.700000 9.700000 0.00599 0.00085 0.00085
	+ 175.500000 9.700000 9.700000 0.00447 0.00045 0.00045
	+ 194.900000 9.700000 9.700000 0.00345 0.00041 0.00041
	+ 214.300000 9.700000 9.700000 0.00256 0.00024 0.00024
	+ 233.700000 9.700000 9.700000 0.00190 0.00018 0.00018
	+ 262.800000 19.400000 19.400000 0.00138 0.00019 0.00019
	+ 311.300000 29.100000 29.100000 0.000788 0.000083 0.000083
	+ 379.200000 38.800000 38.800000 0.000396 0.000052 0.000052
	+ 466.500000 48.500000 48.500000 0.000174 0.000022 0.000022
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d08-x01-y01
	+Path=/REF/CMS_2017_I1499471/d08-x01-y01
	+Title=CMS, 8 TeV, HT ratio, at least one b jet over inclusive
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 39.700000 9.700000 9.700000 3.76 0.49 0.49
	+ 59.100000 9.700000 9.700000 4.24 0.46 0.46
	+ 78.500000 9.700000 9.700000 5.38 0.44 0.44
	+ 97.900000 9.700000 9.700000 5.97 0.50 0.50
	+ 117.300000 9.700000 9.700000 6.36 0.53 0.53
	+ 136.700000 9.700000 9.700000 6.98 0.55 0.55
	+ 156.100000 9.700000 9.700000 7.66 0.94 0.94
	+ 175.500000 9.700000 9.700000 7.90 0.57 0.57
	+ 194.900000 9.700000 9.700000 8.37 0.85 0.85
	+ 214.300000 9.700000 9.700000 8.23 0.62 0.62
	+ 233.700000 9.700000 9.700000 8.13 0.62 0.62
	+ 262.800000 19.400000 19.400000 8.7 1.0 1.0
	+ 311.300000 29.100000 29.100000 8.77 0.78 0.78
	+ 379.200000 38.800000 38.800000 9.0 1.1 1.1
	+ 466.500000 48.500000 48.500000 8.9 1.1 1.1
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d09-x01-y01
	+Path=/REF/CMS_2017_I1499471/d09-x01-y01
	+Title=CMS, 8 TeV, DeltaPhi_Zb, at least one b jet
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.130900 0.130900 0.130900 0.172 0.022 0.022
	+ 0.392699 0.130900 0.130900 0.205 0.029 0.029
	+ 0.654498 0.130900 0.130900 0.228 0.025 0.025
	+ 0.916298 0.130900 0.130900 0.223 0.029 0.029
	+ 1.178097 0.130900 0.130900 0.266 0.038 0.038
	+ 1.439897 0.130900 0.130900 0.380 0.041 0.041
	+ 1.701696 0.130900 0.130900 0.478 0.045 0.045
	+ 1.963495 0.130900 0.130900 0.645 0.077 0.077
	+ 2.225295 0.130900 0.130900 0.98 0.10 0.10
	+ 2.487094 0.130900 0.130900 1.67 0.16 0.16
	+ 2.748894 0.130900 0.130900 2.96 0.20 0.20
	+ 3.010693 0.130900 0.130900 5.33 0.66 0.66
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d10-x01-y01
	+Path=/REF/CMS_2017_I1499471/d10-x01-y01
	+Title=CMS, 8 TeV, DeltaPhi_Z(b/j) ratio, at least one b jet over inclusive
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.130900 0.130900 0.130900 9.2 1.5 1.5
	+ 0.392699 0.130900 0.130900 10.2 2.1 2.1
	+ 0.654498 0.130900 0.130900 10.3 1.6 1.6
	+ 0.916298 0.130900 0.130900 8.8 1.7 1.7
	+ 1.178097 0.130900 0.130900 8.4 1.5 1.5
	+ 1.439897 0.130900 0.130900 9.1 1.2 1.2
	+ 1.701696 0.130900 0.130900 7.94 0.87 0.87
	+ 1.963495 0.130900 0.130900 6.61 0.82 0.82
	+ 2.225295 0.130900 0.130900 6.02 0.55 0.55
	+ 2.487094 0.130900 0.130900 5.81 0.71 0.71
	+ 2.748894 0.130900 0.130900 5.17 0.34 0.34
	+ 3.010693 0.130900 0.130900 3.71 0.56 0.56
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d11-x01-y01
	+Path=/REF/CMS_2017_I1499471/d11-x01-y01
	+Title=CMS, 8 TeV, Leading b jet transverse momentum, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 36.700000 6.700000 6.700000 0.0043 0.0020 0.0020
	+ 50.100000 6.700000 6.700000 0.0072 0.0016 0.0016
	+ 63.499999 6.700000 6.700000 0.0043 0.0012 0.0012
	+ 76.899999 6.700000 6.700000 0.00253 0.00067 0.00067
	+ 90.299999 6.700000 6.700000 0.0017 0.0010 0.0010
	+ 103.699999 6.700000 6.700000 0.00145 0.00031 0.00031
	+ 117.099998 6.700000 6.700000 0.00083 0.00029 0.00029
	+ 130.499998 6.700000 6.700000 0.00056 0.00025 0.00025
	+ 143.899998 6.700000 6.700000 0.00046 0.00020 0.00020
	+ 157.299998 6.700000 6.700000 0.00027 0.00013 0.00013
	+ 170.699998 6.700000 6.700000 0.00024 0.00018 0.00018
	+ 190.799997 13.400000 13.400000 0.00016 0.00015 0.00015
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d12-x01-y01
	+Path=/REF/CMS_2017_I1499471/d12-x01-y01
	+Title=CMS, 8 TeV, Subleading b jet transverse momentum, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 34.700000 4.700000 4.700000 0.0219 0.0046 0.0046
	+ 44.100000 4.700000 4.700000 0.0070 0.0016 0.0016
	+ 53.500000 4.700000 4.700000 0.00329 0.00077 0.00077
	+ 62.900000 4.700000 4.700000 0.00199 0.00078 0.00078
	+ 72.300000 4.700000 4.700000 0.00101 0.00055 0.00055
	+ 81.700000 4.700000 4.700000 0.00068 0.00027 0.00027
	+ 91.100000 4.700000 4.700000 0.00041 0.00021 0.00021
	+ 100.500000 4.700000 4.700000 0.00031 0.00021 0.00021
	+ 109.900000 4.700000 4.700000 0.00026 0.00012 0.00012
	+ 119.300000 4.700000 4.700000 0.00021 0.00015 0.00015
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d13-x01-y01
	+Path=/REF/CMS_2017_I1499471/d13-x01-y01
	+Title=CMS, 8 TeV, Z boson pt, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 5.000000 5.000000 5.000000 0.00083 0.00037 0.00037
	+ 15.000000 5.000000 5.000000 0.00172 0.00053 0.00053
	+ 25.000000 5.000000 5.000000 0.00251 0.00074 0.00074
	+ 35.000000 5.000000 5.000000 0.00318 0.00067 0.00067
	+ 45.000000 5.000000 5.000000 0.00319 0.00091 0.00091
	+ 55.000000 5.000000 5.000000 0.0033 0.0011 0.0011
	+ 65.000000 5.000000 5.000000 0.00307 0.00093 0.00093
	+ 75.000000 5.000000 5.000000 0.00259 0.00064 0.00064
	+ 85.000000 5.000000 5.000000 0.00224 0.00054 0.00054
	+ 95.000000 5.000000 5.000000 0.00188 0.00071 0.00071
	+ 105.000000 5.000000 5.000000 0.00143 0.00042 0.00042
	+ 115.000000 5.000000 5.000000 0.00115 0.00066 0.00066
	+ 125.000000 5.000000 5.000000 0.00112 0.00033 0.00033
	+ 135.000000 5.000000 5.000000 0.00057 0.00014 0.00014
	+ 145.000000 5.000000 5.000000 0.00052 0.00015 0.00015
	+ 155.000000 5.000000 5.000000 0.00038 0.00014 0.00014
	+ 175.000000 15.000000 15.000000 0.000265 0.000061 0.000061
	+ 215.000000 25.000000 25.000000 0.000126 0.000042 0.000042
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d14-x01-y01
	+Path=/REF/CMS_2017_I1499471/d14-x01-y01
	+Title=CMS, 8 TeV, bb system mass, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 44.230769 14.230769 14.230769 0.00106 0.00032 0.00032
	+ 72.692308 14.230769 14.230769 0.00233 0.00062 0.00062
	+ 101.153846 14.230769 14.230769 0.00269 0.00044 0.00044
	+ 129.615385 14.230769 14.230769 0.00153 0.00032 0.00032
	+ 158.076923 14.230769 14.230769 0.00091 0.00037 0.00037
	+ 186.538462 14.230769 14.230769 0.00081 0.00021 0.00021
	+ 215.000000 14.230769 14.230769 0.00058 0.00012 0.00012
	+ 243.461538 14.230769 14.230769 0.00042 0.00010 0.00010
	+ 271.923077 14.230769 14.230769 0.00024 0.00010 0.00010
	+ 300.384615 14.230769 14.230769 0.000098 0.000056 0.000056
	+ 357.307692 42.692308 42.692308 0.000105 0.000052 0.000052
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d15-x01-y01
	+Path=/REF/CMS_2017_I1499471/d15-x01-y01
	+Title=CMS, 8 TeV, Zbb system mass, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 212.500000 12.500000 12.500000 0.00144 0.00064 0.00064
	+ 237.500000 12.500000 12.500000 0.00165 0.00050 0.00050
	+ 262.500000 12.500000 12.500000 0.00180 0.00051 0.00051
	+ 287.500000 12.500000 12.500000 0.00145 0.00037 0.00037
	+ 312.500000 12.500000 12.500000 0.00112 0.00024 0.00024
	+ 337.500000 12.500000 12.500000 0.00080 0.00019 0.00019
	+ 362.500000 12.500000 12.500000 0.00058 0.00039 0.00039
	+ 400.000000 25.000000 25.000000 0.00039 0.00025 0.00025
	+ 462.500000 37.500000 37.500000 0.000276 0.000083 0.000083
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d16-x01-y01
	+Path=/REF/CMS_2017_I1499471/d16-x01-y01
	+Title=CMS, 8 TeV, DeltaPhi bb, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.130900 0.130900 0.130900 0.072 0.037 0.037
	+ 0.392699 0.130900 0.130900 0.085 0.022 0.022
	+ 0.654498 0.130900 0.130900 0.113 0.036 0.036
	+ 0.916298 0.130900 0.130900 0.085 0.015 0.015
	+ 1.178097 0.130900 0.130900 0.083 0.015 0.015
	+ 1.439897 0.130900 0.130900 0.075 0.019 0.019
	+ 1.701696 0.130900 0.130900 0.088 0.013 0.013
	+ 1.963495 0.130900 0.130900 0.108 0.027 0.027
	+ 2.225295 0.130900 0.130900 0.112 0.018 0.018
	+ 2.487094 0.130900 0.130900 0.136 0.020 0.020
	+ 2.748894 0.130900 0.130900 0.133 0.030 0.030
	+ 3.010693 0.130900 0.130900 0.143 0.044 0.044
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d17-x01-y01
	+Path=/REF/CMS_2017_I1499471/d17-x01-y01
	+Title=CMS, 8 TeV, DeltaR bb, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.250000 0.250000 0.250000 0 0 0
	+ 0.750000 0.250000 0.250000 0.058 0.015 0.015
	+ 1.250000 0.250000 0.250000 0.067 0.011 0.011
	+ 1.750000 0.250000 0.250000 0.092 0.018 0.018
	+ 2.250000 0.250000 0.250000 0.108 0.019 0.019
	+ 2.750000 0.250000 0.250000 0.150 0.024 0.024
	+ 3.250000 0.250000 0.250000 0.119 0.017 0.017
	+ 3.750000 0.250000 0.250000 0.0399 0.0083 0.0083
	+ 4.250000 0.250000 0.250000 0.017 0.012 0.012
	+ 4.750000 0.250000 0.250000 0.0032 0.0067 0.0067
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d18-x01-y01
	+Path=/REF/CMS_2017_I1499471/d18-x01-y01
	+Title=CMS, 8 TeV, DeltaR Zb minimum, at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.133333 0.133333 0.133333 0.0078 0.0053 0.0053
	+ 0.400000 0.133333 0.133333 0.0212 0.0088 0.0088
	+ 0.666667 0.133333 0.133333 0.039 0.012 0.012
	+ 0.933333 0.133333 0.133333 0.047 0.012 0.012
	+ 1.200000 0.133333 0.133333 0.064 0.016 0.016
	+ 1.466667 0.133333 0.133333 0.089 0.019 0.019
	+ 1.733333 0.133333 0.133333 0.094 0.019 0.019
	+ 2.000000 0.133333 0.133333 0.147 0.031 0.031
	+ 2.266667 0.133333 0.133333 0.169 0.028 0.028
	+ 2.533333 0.133333 0.133333 0.175 0.029 0.029
	+ 2.800000 0.133333 0.133333 0.153 0.036 0.036
	+ 3.066667 0.133333 0.133333 0.097 0.025 0.025
	+ 3.333333 0.133333 0.133333 0.047 0.012 0.012
	+ 3.600000 0.133333 0.133333 0.0254 0.0072 0.0072
	+ 3.866667 0.133333 0.133333 0.0173 0.0067 0.0067
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d19-x01-y01
	+Path=/REF/CMS_2017_I1499471/d19-x01-y01
	+Title=CMS, 8 TeV, Asymmetry DR Zb max DR Zb min , at least two b jets
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 0.050000 0.050000 0.050000 1.31 0.20 0.20
	+ 0.150000 0.050000 0.050000 0.84 0.14 0.14
	+ 0.250000 0.050000 0.050000 0.45 0.10 0.10
	+ 0.350000 0.050000 0.050000 0.244 0.069 0.069
	+ 0.450000 0.050000 0.050000 0.150 0.058 0.058
	+ 0.550000 0.050000 0.050000 0.120 0.041 0.041
	+ 0.650000 0.050000 0.050000 0.082 0.030 0.030
	+ 0.850000 0.150000 0.150000 0.032 0.016 0.016
	+# END YODA_SCATTER2D
	+
	+
	+# BEGIN YODA_SCATTER2D /REF/CMS_2017_I1499471/d20-x01-y01
	+Path=/REF/CMS_2017_I1499471/d20-x01-y01
	+Title=CMS, 8 TeV, inclusive b jet multiplicity
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+ 1.000000 0.500000 0.500000 3.55 0.24 0.24
	+ 2.000000 0.500000 0.500000 0.331 0.037 0.037
	+# END YODA_SCATTER2D
	diff --git a/modified_analyses/refdata/CMS_2017_I1610623.yoda b/modified_analyses/refdata/CMS_2017_I1610623.yoda
	new file mode 100644
	--- /dev/null
	+++ b/modified_analyses/refdata/CMS_2017_I1610623.yoda
	@@ -0,0 +1,375 @@
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d01-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d01-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t1
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.000000e+00 5.000000e-01 5.000000e-01 8.326000e+02 2.412000e+01 2.412000e+01
	+2.000000e+00 5.000000e-01 5.000000e-01 1.894000e+02 1.608000e+01 1.608000e+01
	+3.000000e+00 5.000000e-01 5.000000e-01 3.907000e+01 4.482000e+00 4.482000e+00
	+4.000000e+00 5.000000e-01 5.000000e-01 8.587000e+00 1.376000e+00 1.376000e+00
	+5.000000e+00 5.000000e-01 5.000000e-01 2.210000e+00 4.704000e-01 4.704000e-01
	+6.000000e+00 5.000000e-01 5.000000e-01 4.923000e-01 1.516000e-01 1.516000e-01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d02-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d02-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t2
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.000000e+00 5.000000e-01 5.000000e-01 1.072000e+03 2.937000e+01 2.937000e+01
	+2.000000e+00 5.000000e-01 5.000000e-01 2.399000e+02 1.676000e+01 1.676000e+01
	+3.000000e+00 5.000000e-01 5.000000e-01 5.046000e+01 4.715000e+00 4.715000e+00
	+4.000000e+00 5.000000e-01 5.000000e-01 1.139000e+01 1.463000e+00 1.463000e+00
	+5.000000e+00 5.000000e-01 5.000000e-01 2.801000e+00 4.967000e-01 4.967000e-01
	+6.000000e+00 5.000000e-01 5.000000e-01 5.903000e-01 1.596000e-01 1.596000e-01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d03-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d03-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t3
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+3.450000e+01 4.500000e+00 4.500000e+00 3.947000e+01 2.404000e+00 2.404000e+00
	+4.400000e+01 5.000000e+00 5.000000e+00 2.319000e+01 1.560000e+00 1.560000e+00
	+5.550000e+01 6.500000e+00 6.500000e+00 1.359000e+01 7.748000e-01 7.748000e-01
	+7.050000e+01 8.500000e+00 8.500000e+00 7.383000e+00 3.438000e-01 3.438000e-01
	+9.200000e+01 1.300000e+01 1.300000e+01 3.531000e+00 1.635000e-01 1.635000e-01
	+1.215000e+02 1.650000e+01 1.650000e+01 1.518000e+00 7.479000e-02 7.479000e-02
	+1.595000e+02 2.150000e+01 2.150000e+01 5.849000e-01 2.982000e-02 2.982000e-02
	+2.060000e+02 2.500000e+01 2.500000e+01 2.292000e-01 1.535000e-02 1.535000e-02
	+2.625000e+02 3.150000e+01 3.150000e+01 8.556000e-02 5.253000e-03 5.253000e-03
	+3.345000e+02 4.050000e+01 4.050000e+01 2.928000e-02 1.921000e-03 1.921000e-03
	+4.625000e+02 8.750000e+01 8.750000e+01 6.415000e-03 6.080000e-04 6.080000e-04
	+7.250000e+02 1.750000e+02 1.750000e+02 7.925000e-04 1.164000e-04 1.164000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d04-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d04-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t4
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+3.450000e+01 4.500000e+00 4.500000e+00 1.039000e+01 7.847000e-01 7.847000e-01
	+4.400000e+01 5.000000e+00 5.000000e+00 5.272000e+00 3.950000e-01 3.950000e-01
	+5.550000e+01 6.500000e+00 6.500000e+00 2.722000e+00 1.751000e-01 1.751000e-01
	+7.000000e+01 8.000000e+00 8.000000e+00 1.372000e+00 7.473000e-02 7.473000e-02
	+9.150000e+01 1.350000e+01 1.350000e+01 6.209000e-01 3.604000e-02 3.604000e-02
	+1.235000e+02 1.850000e+01 1.850000e+01 2.552000e-01 1.401000e-02 1.401000e-02
	+1.635000e+02 2.150000e+01 2.150000e+01 9.426000e-02 6.058000e-03 6.058000e-03
	+2.100000e+02 2.500000e+01 2.500000e+01 4.108000e-02 3.239000e-03 3.239000e-03
	+2.675000e+02 3.250000e+01 3.250000e+01 1.459000e-02 1.501000e-03 1.501000e-03
	+3.400000e+02 4.000000e+01 4.000000e+01 4.704000e-03 6.895000e-04 6.895000e-04
	+4.400000e+02 6.000000e+01 6.000000e+01 1.619000e-03 2.709000e-04 2.709000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d05-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d05-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t5
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+3.550000e+01 5.500000e+00 5.500000e+00 2.407000e+00 2.271000e-01 2.271000e-01
	+5.000000e+01 9.000000e+00 9.000000e+00 7.765000e-01 6.669000e-02 6.669000e-02
	+7.000000e+01 1.100000e+01 1.100000e+01 2.319000e-01 2.398000e-02 2.398000e-02
	+9.550000e+01 1.450000e+01 1.450000e+01 7.904000e-02 6.557000e-03 6.557000e-03
	+1.310000e+02 2.100000e+01 2.100000e+01 1.898000e-02 2.288000e-03 2.288000e-03
	+1.760000e+02 2.400000e+01 2.400000e+01 6.092000e-03 9.320000e-04 9.320000e-04
	+2.500000e+02 5.000000e+01 5.000000e+01 9.463000e-04 2.958000e-04 2.958000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d06-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d06-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t6
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+3.450000e+01 4.500000e+00 4.500000e+00 6.514000e-01 7.782000e-02 7.782000e-02
	+4.400000e+01 5.000000e+00 5.000000e+00 2.412000e-01 3.533000e-02 3.533000e-02
	+5.550000e+01 6.500000e+00 6.500000e+00 9.452000e-02 1.427000e-02 1.427000e-02
	+7.000000e+01 8.000000e+00 8.000000e+00 4.072000e-02 6.103000e-03 6.103000e-03
	+8.700000e+01 9.000000e+00 9.000000e+00 1.680000e-02 2.994000e-03 2.994000e-03
	+1.380000e+02 4.200000e+01 4.200000e+01 1.989000e-03 4.843000e-04 4.843000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d07-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d07-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t7
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.000000e-01 1.000000e-01 1.000000e-01 5.772000e+02 3.437000e+01 3.437000e+01
	+3.000000e-01 1.000000e-01 1.000000e-01 5.710000e+02 3.295000e+01 3.295000e+01
	+5.000000e-01 1.000000e-01 1.000000e-01 5.571000e+02 3.245000e+01 3.245000e+01
	+7.000000e-01 1.000000e-01 1.000000e-01 5.328000e+02 3.246000e+01 3.246000e+01
	+9.000000e-01 1.000000e-01 1.000000e-01 5.238000e+02 3.073000e+01 3.073000e+01
	+1.100000e+00 1.000000e-01 1.000000e-01 4.911000e+02 3.121000e+01 3.121000e+01
	+1.300000e+00 1.000000e-01 1.000000e-01 4.406000e+02 2.763000e+01 2.763000e+01
	+1.500000e+00 1.000000e-01 1.000000e-01 4.151000e+02 2.699000e+01 2.699000e+01
	+1.700000e+00 1.000000e-01 1.000000e-01 3.921000e+02 2.457000e+01 2.457000e+01
	+1.900000e+00 1.000000e-01 1.000000e-01 3.548000e+02 2.509000e+01 2.509000e+01
	+2.100000e+00 1.000000e-01 1.000000e-01 3.153000e+02 2.318000e+01 2.318000e+01
	+2.300000e+00 1.000000e-01 1.000000e-01 2.690000e+02 2.081000e+01 2.081000e+01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d08-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d08-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t8
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.000000e-01 1.000000e-01 1.000000e-01 1.218000e+02 9.328000e+00 9.328000e+00
	+3.000000e-01 1.000000e-01 1.000000e-01 1.232000e+02 9.264000e+00 9.264000e+00
	+5.000000e-01 1.000000e-01 1.000000e-01 1.196000e+02 8.902000e+00 8.902000e+00
	+7.000000e-01 1.000000e-01 1.000000e-01 1.199000e+02 9.404000e+00 9.404000e+00
	+9.000000e-01 1.000000e-01 1.000000e-01 1.148000e+02 9.040000e+00 9.040000e+00
	+1.100000e+00 1.000000e-01 1.000000e-01 1.057000e+02 8.762000e+00 8.762000e+00
	+1.300000e+00 1.000000e-01 1.000000e-01 9.849000e+01 8.095000e+00 8.095000e+00
	+1.500000e+00 1.000000e-01 1.000000e-01 9.482000e+01 7.842000e+00 7.842000e+00
	+1.700000e+00 1.000000e-01 1.000000e-01 8.956000e+01 7.422000e+00 7.422000e+00
	+1.900000e+00 1.000000e-01 1.000000e-01 8.233000e+01 7.769000e+00 7.769000e+00
	+2.100000e+00 1.000000e-01 1.000000e-01 7.262000e+01 6.951000e+00 6.951000e+00
	+2.300000e+00 1.000000e-01 1.000000e-01 6.426000e+01 6.547000e+00 6.547000e+00
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d09-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d09-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t9
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.500000e-01 1.500000e-01 1.500000e-01 2.564000e+01 2.880000e+00 2.880000e+00
	+4.500000e-01 1.500000e-01 1.500000e-01 2.600000e+01 2.654000e+00 2.654000e+00
	+7.500000e-01 1.500000e-01 1.500000e-01 2.471000e+01 2.579000e+00 2.579000e+00
	+1.050000e+00 1.500000e-01 1.500000e-01 2.289000e+01 2.649000e+00 2.649000e+00
	+1.350000e+00 1.500000e-01 1.500000e-01 2.085000e+01 2.264000e+00 2.264000e+00
	+1.650000e+00 1.500000e-01 1.500000e-01 1.903000e+01 2.111000e+00 2.111000e+00
	+1.950000e+00 1.500000e-01 1.500000e-01 1.647000e+01 2.045000e+00 2.045000e+00
	+2.250000e+00 1.500000e-01 1.500000e-01 1.321000e+01 1.786000e+00 1.786000e+00
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d10-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d10-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t10
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.500000e-01 1.500000e-01 1.500000e-01 6.185000e+00 9.756000e-01 9.756000e-01
	+4.500000e-01 1.500000e-01 1.500000e-01 5.440000e+00 8.302000e-01 8.302000e-01
	+7.500000e-01 1.500000e-01 1.500000e-01 5.504000e+00 9.149000e-01 9.149000e-01
	+1.050000e+00 1.500000e-01 1.500000e-01 5.176000e+00 8.631000e-01 8.631000e-01
	+1.350000e+00 1.500000e-01 1.500000e-01 4.347000e+00 7.842000e-01 7.842000e-01
	+1.650000e+00 1.500000e-01 1.500000e-01 3.630000e+00 6.907000e-01 6.907000e-01
	+1.950000e+00 1.500000e-01 1.500000e-01 3.790000e+00 7.183000e-01 7.183000e-01
	+2.250000e+00 1.500000e-01 1.500000e-01 2.419000e+00 5.672000e-01 5.672000e-01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d11-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d11-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t11
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+3.450000e+01 4.500000e+00 4.500000e+00 3.768000e+01 2.766000e+00 2.766000e+00
	+4.400000e+01 5.000000e+00 5.000000e+00 1.976000e+01 1.219000e+00 1.219000e+00
	+5.550000e+01 6.500000e+00 6.500000e+00 1.041000e+01 5.030000e-01 5.030000e-01
	+7.000000e+01 8.000000e+00 8.000000e+00 7.057000e+00 4.032000e-01 4.032000e-01
	+8.700000e+01 9.000000e+00 9.000000e+00 4.827000e+00 3.086000e-01 3.086000e-01
	+1.070000e+02 1.100000e+01 1.100000e+01 3.023000e+00 1.777000e-01 1.777000e-01
	+1.340000e+02 1.600000e+01 1.600000e+01 1.754000e+00 1.098000e-01 1.098000e-01
	+1.700000e+02 2.000000e+01 2.000000e+01 9.311000e-01 6.115000e-02 6.115000e-02
	+2.150000e+02 2.500000e+01 2.500000e+01 4.759000e-01 3.183000e-02 3.183000e-02
	+2.700000e+02 3.000000e+01 3.000000e+01 2.323000e-01 1.750000e-02 1.750000e-02
	+3.350000e+02 3.500000e+01 3.500000e+01 1.151000e-01 8.458000e-03 8.458000e-03
	+4.100000e+02 4.000000e+01 4.000000e+01 5.771000e-02 4.692000e-03 4.692000e-03
	+4.950000e+02 4.500000e+01 4.500000e+01 3.066000e-02 2.875000e-03 2.875000e-03
	+5.950000e+02 5.500000e+01 5.500000e+01 1.471000e-02 1.529000e-03 1.529000e-03
	+7.350000e+02 8.500000e+01 8.500000e+01 5.222000e-03 6.205000e-04 6.205000e-04
	+9.600000e+02 1.400000e+02 1.400000e+02 1.860000e-03 2.428000e-04 2.428000e-04
	+1.300000e+03 2.000000e+02 2.000000e+02 3.299000e-04 7.843000e-05 7.843000e-05
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d12-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d12-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t12
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+6.900000e+01 9.000000e+00 9.000000e+00 1.924000e+00 3.157000e-01 3.157000e-01
	+8.700000e+01 9.000000e+00 9.000000e+00 2.386000e+00 2.160000e-01 2.160000e-01
	+1.070000e+02 1.100000e+01 1.100000e+01 1.869000e+00 1.271000e-01 1.271000e-01
	+1.340000e+02 1.600000e+01 1.600000e+01 1.278000e+00 8.799000e-02 8.799000e-02
	+1.700000e+02 2.000000e+01 2.000000e+01 7.628000e-01 5.288000e-02 5.288000e-02
	+2.150000e+02 2.500000e+01 2.500000e+01 4.250000e-01 2.848000e-02 2.848000e-02
	+2.700000e+02 3.000000e+01 3.000000e+01 2.179000e-01 1.625000e-02 1.625000e-02
	+3.350000e+02 3.500000e+01 3.500000e+01 1.107000e-01 8.015000e-03 8.015000e-03
	+4.100000e+02 4.000000e+01 4.000000e+01 5.667000e-02 4.496000e-03 4.496000e-03
	+4.950000e+02 4.500000e+01 4.500000e+01 3.033000e-02 2.787000e-03 2.787000e-03
	+5.950000e+02 5.500000e+01 5.500000e+01 1.450000e-02 1.483000e-03 1.483000e-03
	+7.250000e+02 7.500000e+01 7.500000e+01 5.335000e-03 6.661000e-04 6.661000e-04
	+1.000000e+03 2.000000e+02 2.000000e+02 1.656000e-03 1.810000e-04 1.810000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d13-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d13-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t13
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.430000e+02 2.500000e+01 2.500000e+01 2.747000e-01 3.618000e-02 3.618000e-02
	+1.940000e+02 2.600000e+01 2.600000e+01 2.210000e-01 2.156000e-02 2.156000e-02
	+2.600000e+02 4.000000e+01 4.000000e+01 1.278000e-01 1.314000e-02 1.314000e-02
	+3.500000e+02 5.000000e+01 5.000000e+01 6.032000e-02 5.831000e-03 5.831000e-03
	+4.750000e+02 7.500000e+01 7.500000e+01 2.577000e-02 2.475000e-03 2.475000e-03
	+6.650000e+02 1.150000e+02 1.150000e+02 6.927000e-03 6.607000e-04 6.607000e-04
	+9.400000e+02 1.600000e+02 1.600000e+02 1.781000e-03 2.152000e-04 2.152000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d14-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d14-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t14
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+2.280000e+02 2.500000e+01 2.500000e+01 4.037000e-02 7.844000e-03 7.844000e-03
	+2.865000e+02 3.350000e+01 3.350000e+01 3.020000e-02 5.600000e-03 5.600000e-03
	+3.650000e+02 4.500000e+01 4.500000e+01 2.203000e-02 3.116000e-03 3.116000e-03
	+4.700000e+02 6.000000e+01 6.000000e+01 1.219000e-02 1.721000e-03 1.721000e-03
	+6.100000e+02 8.000000e+01 8.000000e+01 5.927000e-03 8.501000e-04 8.501000e-04
	+8.000000e+02 1.100000e+02 1.100000e+02 1.860000e-03 3.304000e-04 3.304000e-04
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d15-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d15-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t15
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+7.854000e-02 7.854000e-02 7.854000e-02 8.503000e+01 8.089000e+00 8.089000e+00
	+2.356195e-01 7.853950e-02 7.853950e-02 9.005000e+01 8.426000e+00 8.426000e+00
	+3.926990e-01 7.854000e-02 7.854000e-02 1.069000e+02 9.810000e+00 9.810000e+00
	+5.497790e-01 7.854000e-02 7.854000e-02 1.220000e+02 1.102000e+01 1.102000e+01
	+7.068585e-01 7.853950e-02 7.853950e-02 1.383000e+02 1.200000e+01 1.200000e+01
	+8.639380e-01 7.854000e-02 7.854000e-02 1.516000e+02 1.235000e+01 1.235000e+01
	+1.021019e+00 7.854100e-02 7.854100e-02 1.724000e+02 1.330000e+01 1.330000e+01
	+1.178100e+00 7.854000e-02 7.854000e-02 1.950000e+02 1.408000e+01 1.408000e+01
	+1.335180e+00 7.854000e-02 7.854000e-02 2.274000e+02 1.577000e+01 1.577000e+01
	+1.492260e+00 7.854000e-02 7.854000e-02 2.647000e+02 1.693000e+01 1.693000e+01
	+1.649340e+00 7.854000e-02 7.854000e-02 3.039000e+02 1.826000e+01 1.826000e+01
	+1.806420e+00 7.854000e-02 7.854000e-02 3.534000e+02 2.067000e+01 2.067000e+01
	+1.963500e+00 7.854000e-02 7.854000e-02 4.073000e+02 2.283000e+01 2.283000e+01
	+2.120575e+00 7.853500e-02 7.853500e-02 4.557000e+02 2.546000e+01 2.546000e+01
	+2.277650e+00 7.854000e-02 7.854000e-02 5.261000e+02 2.872000e+01 2.872000e+01
	+2.434730e+00 7.854000e-02 7.854000e-02 5.726000e+02 3.160000e+01 3.160000e+01
	+2.591810e+00 7.854000e-02 7.854000e-02 6.345000e+02 3.454000e+01 3.454000e+01
	+2.748890e+00 7.854000e-02 7.854000e-02 6.899000e+02 3.909000e+01 3.909000e+01
	+2.905970e+00 7.854000e-02 7.854000e-02 7.137000e+02 4.463000e+01 4.463000e+01
	+3.063050e+00 7.854000e-02 7.854000e-02 7.157000e+02 5.144000e+01 5.144000e+01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d16-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d16-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t16
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+7.854000e-02 7.854000e-02 7.854000e-02 5.548000e+01 5.180000e+00 5.180000e+00
	+2.356195e-01 7.853950e-02 7.853950e-02 5.533000e+01 4.919000e+00 4.919000e+00
	+3.926990e-01 7.854000e-02 7.854000e-02 6.054000e+01 5.651000e+00 5.651000e+00
	+5.497790e-01 7.854000e-02 7.854000e-02 6.954000e+01 6.201000e+00 6.201000e+00
	+7.068585e-01 7.853950e-02 7.853950e-02 6.970000e+01 5.906000e+00 5.906000e+00
	+8.639380e-01 7.854000e-02 7.854000e-02 7.120000e+01 6.073000e+00 6.073000e+00
	+1.021019e+00 7.854100e-02 7.854100e-02 7.427000e+01 6.004000e+00 6.004000e+00
	+1.178100e+00 7.854000e-02 7.854000e-02 7.555000e+01 6.419000e+00 6.419000e+00
	+1.335180e+00 7.854000e-02 7.854000e-02 7.968000e+01 6.686000e+00 6.686000e+00
	+1.492260e+00 7.854000e-02 7.854000e-02 7.800000e+01 6.233000e+00 6.233000e+00
	+1.649340e+00 7.854000e-02 7.854000e-02 8.189000e+01 6.953000e+00 6.953000e+00
	+1.806420e+00 7.854000e-02 7.854000e-02 8.573000e+01 7.020000e+00 7.020000e+00
	+1.963500e+00 7.854000e-02 7.854000e-02 8.243000e+01 6.771000e+00 6.771000e+00
	+2.120575e+00 7.853500e-02 7.853500e-02 8.384000e+01 6.886000e+00 6.886000e+00
	+2.277650e+00 7.854000e-02 7.854000e-02 8.480000e+01 6.903000e+00 6.903000e+00
	+2.434730e+00 7.854000e-02 7.854000e-02 8.608000e+01 6.549000e+00 6.549000e+00
	+2.591810e+00 7.854000e-02 7.854000e-02 8.391000e+01 6.885000e+00 6.885000e+00
	+2.748890e+00 7.854000e-02 7.854000e-02 8.638000e+01 7.251000e+00 7.251000e+00
	+2.905970e+00 7.854000e-02 7.854000e-02 8.718000e+01 7.370000e+00 7.370000e+00
	+3.063050e+00 7.854000e-02 7.854000e-02 8.497000e+01 6.848000e+00 6.848000e+00
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d17-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d17-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t17
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+9.817500e-02 9.817500e-02 9.817500e-02 1.152000e+01 1.617000e+00 1.617000e+00
	+2.945245e-01 9.817450e-02 9.817450e-02 1.175000e+01 1.807000e+00 1.807000e+00
	+4.908740e-01 9.817500e-02 9.817500e-02 1.418000e+01 1.979000e+00 1.979000e+00
	+6.872235e-01 9.817450e-02 9.817450e-02 1.346000e+01 1.873000e+00 1.873000e+00
	+8.835730e-01 9.817500e-02 9.817500e-02 1.544000e+01 2.025000e+00 2.025000e+00
	+1.079924e+00 9.817600e-02 9.817600e-02 1.433000e+01 1.903000e+00 1.903000e+00
	+1.276275e+00 9.817500e-02 9.817500e-02 1.434000e+01 1.540000e+00 1.540000e+00
	+1.472625e+00 9.817500e-02 9.817500e-02 1.598000e+01 2.049000e+00 2.049000e+00
	+1.668975e+00 9.817500e-02 9.817500e-02 1.622000e+01 1.919000e+00 1.919000e+00
	+1.865325e+00 9.817500e-02 9.817500e-02 1.671000e+01 2.037000e+00 2.037000e+00
	+2.061670e+00 9.817000e-02 9.817000e-02 1.679000e+01 1.991000e+00 1.991000e+00
	+2.258015e+00 9.817500e-02 9.817500e-02 1.870000e+01 2.204000e+00 2.204000e+00
	+2.454365e+00 9.817500e-02 9.817500e-02 1.886000e+01 2.002000e+00 2.002000e+00
	+2.650715e+00 9.817500e-02 9.817500e-02 1.975000e+01 2.209000e+00 2.209000e+00
	+2.847065e+00 9.817500e-02 9.817500e-02 2.104000e+01 2.281000e+00 2.281000e+00
	+3.043415e+00 9.817500e-02 9.817500e-02 1.909000e+01 2.411000e+00 2.411000e+00
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d18-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d18-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t18
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+9.817500e-02 9.817500e-02 9.817500e-02 2.350000e+00 5.350000e-01 5.350000e-01
	+2.945245e-01 9.817450e-02 9.817450e-02 2.678000e+00 5.408000e-01 5.408000e-01
	+4.908740e-01 9.817500e-02 9.817500e-02 3.320000e+00 6.792000e-01 6.792000e-01
	+6.872235e-01 9.817450e-02 9.817450e-02 3.155000e+00 7.588000e-01 7.588000e-01
	+8.835730e-01 9.817500e-02 9.817500e-02 4.090000e+00 8.843000e-01 8.843000e-01
	+1.079924e+00 9.817600e-02 9.817600e-02 3.108000e+00 7.162000e-01 7.162000e-01
	+1.276275e+00 9.817500e-02 9.817500e-02 3.053000e+00 6.824000e-01 6.824000e-01
	+1.472625e+00 9.817500e-02 9.817500e-02 2.517000e+00 6.476000e-01 6.476000e-01
	+1.668975e+00 9.817500e-02 9.817500e-02 4.130000e+00 8.105000e-01 8.105000e-01
	+1.865325e+00 9.817500e-02 9.817500e-02 4.226000e+00 9.618000e-01 9.618000e-01
	+2.061670e+00 9.817000e-02 9.817000e-02 3.786000e+00 7.497000e-01 7.497000e-01
	+2.258015e+00 9.817500e-02 9.817500e-02 4.610000e+00 7.465000e-01 7.465000e-01
	+2.454365e+00 9.817500e-02 9.817500e-02 3.755000e+00 7.251000e-01 7.251000e-01
	+2.650715e+00 9.817500e-02 9.817500e-02 3.920000e+00 7.816000e-01 7.816000e-01
	+2.847065e+00 9.817500e-02 9.817500e-02 3.190000e+00 8.885000e-01 8.885000e-01
	+3.043415e+00 9.817500e-02 9.817500e-02 4.092000e+00 6.768000e-01 6.768000e-01
	+END YODA_SCATTER2D
	+
	+BEGIN YODA_SCATTER2D /REF/CMS_2017_I1610623/d19-x01-y01
	+IsRef=1
	+Path=/REF/CMS_2017_I1610623/d19-x01-y01
	+Title=doi:10.17182/hepdata.79859.v1/t19
	+Type=Scatter2D
	+# xval xerr- xerr+ yval yerr- yerr+
	+1.000000e-01 1.000000e-01 1.000000e-01 0.000000e+00 0.000000e+00 0.000000e+00
	+3.000000e-01 1.000000e-01 1.000000e-01 0.000000e+00 0.000000e+00 0.000000e+00
	+5.000000e-01 1.000000e-01 1.000000e-01 7.429000e-01 2.184000e-01 2.184000e-01
	+7.000000e-01 1.000000e-01 1.000000e-01 1.271000e+00 2.197000e-01 2.197000e-01
	+9.000000e-01 1.000000e-01 1.000000e-01 1.512000e+00 2.464000e-01 2.464000e-01
	+1.100000e+00 1.000000e-01 1.000000e-01 1.045000e+00 1.897000e-01 1.897000e-01
	+1.300000e+00 1.000000e-01 1.000000e-01 1.323000e+00 2.192000e-01 2.192000e-01
	+1.500000e+00 1.000000e-01 1.000000e-01 1.299000e+00 1.944000e-01 1.944000e-01
	+1.700000e+00 1.000000e-01 1.000000e-01 1.107000e+00 1.673000e-01 1.673000e-01
	+1.900000e+00 1.000000e-01 1.000000e-01 8.555000e-01 1.408000e-01 1.408000e-01
	+2.100000e+00 1.000000e-01 1.000000e-01 6.746000e-01 1.522000e-01 1.522000e-01
	+2.300000e+00 1.000000e-01 1.000000e-01 7.273000e-01 1.405000e-01 1.405000e-01
	+2.500000e+00 1.000000e-01 1.000000e-01 6.110000e-01 1.550000e-01 1.550000e-01
	+2.700000e+00 1.000000e-01 1.000000e-01 1.093000e+00 1.798000e-01 1.798000e-01
	+2.900000e+00 1.000000e-01 1.000000e-01 2.040000e+00 2.556000e-01 2.556000e-01
	+3.100000e+00 1.000000e-01 1.000000e-01 2.642000e+00 3.167000e-01 3.167000e-01
	+3.300000e+00 1.000000e-01 1.000000e-01 8.321000e-01 1.641000e-01 1.641000e-01
	+3.500000e+00 1.000000e-01 1.000000e-01 3.947000e-01 1.502000e-01 1.502000e-01
	+3.700000e+00 1.000000e-01 1.000000e-01 5.237000e-01 2.717000e-01 2.717000e-01
	+3.900000e+00 1.000000e-01 1.000000e-01 1.550000e-01 7.565000e-02 7.565000e-02
	+END YODA_SCATTER2D
	+

File Metadata

Mime Type: text/x-diff
Expires: Sat, Dec 21, 1:52 PM (15 h, 32 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4023031
Default Alt Text: (113 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions