No OneTemporary
Actions

Size

151 KB

Subscribers

None

View Options

	diff --git a/Tests/Makefile.am b/Tests/Makefile.am
	--- a/Tests/Makefile.am
	+++ b/Tests/Makefile.am
	@@ -1,387 +1,387 @@
	AM_LDFLAGS += -module -avoid-version -rpath /dummy/path/not/used

	EXTRA_DIST = Inputs python Rivet

	EXTRA_LTLIBRARIES = LeptonTest.la GammaTest.la HadronTest.la DISTest.la

	if WANT_LIBFASTJET
	EXTRA_LTLIBRARIES += HadronJetTest.la LeptonJetTest.la
	HadronJetTest_la_SOURCES = \
	Hadron/VHTest.h Hadron/VHTest.cc\
	Hadron/VTest.h Hadron/VTest.cc\
	Hadron/HTest.h Hadron/HTest.cc
	HadronJetTest_la_CPPFLAGS = $(AM_CPPFLAGS) $(FASTJETINCLUDE) \
	-I$(FASTJETPATH)
	HadronJetTest_la_LIBADD = $(FASTJETLIBS)
	LeptonJetTest_la_SOURCES = \
	Lepton/TopDecay.h Lepton/TopDecay.cc
	LeptonJetTest_la_CPPFLAGS = $(AM_CPPFLAGS) $(FASTJETINCLUDE) \
	-I$(FASTJETPATH)
	LeptonJetTest_la_LIBADD = $(FASTJETLIBS)
	endif

	LeptonTest_la_SOURCES = \
	Lepton/VVTest.h Lepton/VVTest.cc \
	Lepton/VBFTest.h Lepton/VBFTest.cc \
	Lepton/VHTest.h Lepton/VHTest.cc \
	Lepton/FermionTest.h Lepton/FermionTest.cc

	GammaTest_la_SOURCES = \
	Gamma/GammaMETest.h Gamma/GammaMETest.cc \
	Gamma/GammaPMETest.h Gamma/GammaPMETest.cc

	DISTest_la_SOURCES = \
	DIS/DISTest.h DIS/DISTest.cc

	HadronTest_la_SOURCES = \
	Hadron/HadronVVTest.h Hadron/HadronVVTest.cc\
	Hadron/HadronVBFTest.h Hadron/HadronVBFTest.cc\
	Hadron/WHTest.h Hadron/WHTest.cc\
	Hadron/ZHTest.h Hadron/ZHTest.cc\
	Hadron/VGammaTest.h Hadron/VGammaTest.cc\
	Hadron/ZJetTest.h Hadron/ZJetTest.cc\
	Hadron/WJetTest.h Hadron/WJetTest.cc\
	Hadron/QQHTest.h Hadron/QQHTest.cc


	REPO = $(top_builddir)/src/HerwigDefaults.rpo
	HERWIG = $(top_builddir)/src/Herwig
	HWREAD = $(HERWIG) read --repo $(REPO) -L $(builddir)/.libs -i $(top_builddir)/src
	HWBUILD = $(HERWIG) build --repo $(REPO) -L $(builddir)/.libs -i $(top_builddir)/src
	HWINTEGRATE = $(HERWIG) integrate
	HWRUN = $(HERWIG) run -N $${NUMEVENTS:-10000}


	tests : tests-LEP tests-DIS tests-LHC tests-Gamma


	LEPDEPS = \
	test-LEP-VV \
	test-LEP-VH \
	test-LEP-VBF \
	test-LEP-BB \
	test-LEP-Quarks \
	test-LEP-Leptons

	if WANT_LIBFASTJET
	LEPDEPS += test-LEP-TopDecay
	endif

	tests-LEP : $(LEPDEPS)

	tests-DIS : test-DIS-Charged test-DIS-Neutral


	LHCDEPS = \
	test-LHC-WW test-LHC-WZ test-LHC-ZZ \
	test-LHC-ZGamma test-LHC-WGamma \
	test-LHC-ZH test-LHC-WH \
	test-LHC-ZJet test-LHC-WJet \
	test-LHC-Z test-LHC-W \
	test-LHC-ZZVBF test-LHC-VBF \
	test-LHC-WWVBF \
	test-LHC-bbH test-LHC-ttH \
	test-LHC-GammaGamma test-LHC-GammaJet \
	test-LHC-Higgs test-LHC-HiggsJet \
	test-LHC-QCDFast test-LHC-QCD \
	test-LHC-Top


	if WANT_LIBFASTJET
	LHCDEPS += \
	test-LHC-Bottom \
	test-LHC-WHJet test-LHC-ZHJet test-LHC-HJet \
	test-LHC-ZShower test-LHC-WShower \
	test-LHC-WHJet-Powheg test-LHC-ZHJet-Powheg test-LHC-HJet-Powheg \
	test-LHC-ZShower-Powheg test-LHC-WShower-Powheg
	endif

	tests-LHC : $(LHCDEPS)

	tests-Gamma : test-Gamma-FF test-Gamma-WW test-Gamma-P



	LEPLIBS = LeptonTest.la
	HADLIBS = HadronTest.la

	if WANT_LIBFASTJET
	LEPLIBS += LeptonJetTest.la
	HADLIBS += HadronJetTest.la
	endif


	test-LEP-% : Inputs/LEP-%.in $(LEPLIBS)
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<))

	test-Gamma-% : Inputs/Gamma-%.in GammaTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<))

	test-DIS-% : Inputs/DIS-%.in DISTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<))

	test-LHC-% : Inputs/LHC-%.in GammaTest.la $(HADLIBS)
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<))



	tests-Rivet : Rivet-EE Rivet-DIS Rivet-Fixed \
	Rivet-TVT-WZ Rivet-TVT-Photon Rivet-TVT-Jets \
	Rivet-LHC-Jets Rivet-LHC-EW Rivet-LHC-Photon Rivet-LHC-Higgs

	Rivet-%-UE-Cent.yoda : Rivet-%-UE-Cent.run Rivet-%-Cent.run
	rm -rf $(subst .yoda,,$@)
	mkdir $(subst .yoda,,$@)
	mv $(subst .yoda,.run,$@) $(subst UE-,,$(subst .yoda,.run,$@)) $(subst .yoda,,$@)
	cd $(subst .yoda,,$@); ../$(HWRUN) $(subst UE-,,$(subst .yoda,.run,$@)); mv $(subst UE-,,$@) ALICE_2015_PPCentrality.yoda
	cd $(subst .yoda,,$@); export RIVET_ANALYSIS_PATH=`pwd`; ../$(HWRUN) $(subst .yoda,.run,$@)
	mv $(subst .yoda,,$@)/$@ .
	rm -rf $(subst .yoda,,$@)

	Rivet-%.run : Rivet/%.in
	$(HWBUILD) -c .cache/$(subst .run,,$@) $<

	Rivet-Matchbox-%.yoda : Rivet-Matchbox-%.run
	$(HWINTEGRATE) -c .cache/$(subst .run,,$<) $<
	$(HWRUN) -c .cache/$(subst .run,,$<) $<

	Rivet-%.yoda : Rivet-%.run
	$(HWRUN) $<

	Rivet/%.in :
	python/make_input_files.py $(notdir $(subst .in,,$@))

	Rivet-inputfiles: $(shell echo Rivet/EE{,-Powheg,-Matchbox,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox-Powheg,-Merging}-{7.7,9.4,12,13,17,27.6,27.7,29,30.2,30.7,30,31.3,31.6,34,34.8,41,42.1,42.6,43.6,45,50,52,53.3,55,56,57,58,59.5,60.8,60,61.4,66,76,82,85,10,12.8,21.5,22,25,26.8,34.5,35,36.2,44,48.0,91,93.0,130,133,136,161,172,177,183,189,192,196,197,200,202,205,206,207,91-nopi}.in) \
	$(shell echo Rivet/EE{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Powheg,-Matchbox-Powheg}-{14,14.8}.in) \
	$(shell echo Rivet/EE{,-Dipole}-{10.5,11.96,12.8,13.96,16.86,21.84,26.8,28.48,35.44,48.0,97.0}-gg.in) \
	$(shell echo Rivet/EE{,-Powheg,-Matchbox,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox-Powheg}-{2.2,2.6,3.0,3.2,4.17,4.3,4.41,5.0,5.2,4.6,4.8,5.8,6.2,6.6,7.0,7.4,3.63,4.03,4.5,8.8,9.27,9.46,9.51,10.52,10.52-sym,10.54,10.58,10.45,10.47,10.6}.in) \
	$(shell echo Rivet/EE-{Upsilon,Upsilon2,Upsilon3,Upsilon4,Upsilon4-asym,JPsi,Psi2S,Psi2S-All,Psi3770,Tau,Phi,Lambdac,Omega-Meson,Omega-Baryon,Eta,Xi0,Xic0,Xicp,Omegac0,Xim}.in) \
	$(shell echo Rivet/DIS{,-NoME,-Powheg,-Matchbox,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox-Powheg,-Merging}-{e--LowQ2,e+-LowQ2,e+-HighQ2}.in) \
	$(shell echo Rivet/TVT{,-Powheg,-Matchbox,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox-Powheg,-Merging}-{Run-I-Z,Run-I-W,Run-I-WZ,Run-II-Z-e,Run-II-Z-{,LowMass-,HighMass-}mu,Run-II-W}.in) \
	$(shell echo Rivet/TVT{,-Dipole}-Run-II-{DiPhoton-GammaGamma,DiPhoton-GammaJet,PromptPhoton}.in) \
	$(shell echo Rivet/TVT-Powheg-Run-II-{DiPhoton-GammaGamma,DiPhoton-GammaJet}.in) \
	$(shell echo Rivet/TVT{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-{Run-II-Jets-{0..11},Run-I-Jets-{1..8}}.in ) \
	$(shell echo Rivet/TVT{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-{630-Jets-{1..3},300-Jets-1,900-Jets-1}.in ) \
	$(shell echo Rivet/TVT{,-Dipole}-{Run-I,Run-II,300,630,900}-UE.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-7-DiJets-{1..7}-{A,B,C}.in ) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-13-DiJets-{{1..11}-A,{6..11}-B}.in ) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-{7,8,13}-Jets-{0..10}.in ) \
	$(shell echo Rivet/LHC{,-Dipole}-{900,2360,2760,7,8,13}-UE.in ) \
	$(shell echo Rivet/LHC{,-Dipole}-{7,13}{,-UE}-Cent.in ) \
	$(shell echo Rivet/LHC{,-Dipole}-2760-Jets-{1..3}.in ) \
	$(shell echo Rivet/LHC{,-Dipole}-{900,7,13}-UE-Long.in ) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-7-Charm-{1..5}.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-7-Bottom-{0..9}.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-7-Top-{L,SL}.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Matchbox,-Matchbox-Powheg,-Merging}-{8,13}-Top-{All,L,SL}.in) \
	$(shell echo Rivet/Star{,-Dipole}-{UE,Jets-{1..4}}.in ) \
	$(shell echo Rivet/SppS{,-Dipole}-{53,63,200,500,546,900}-UE.in ) \
	- $(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-{W-{e,mu},13-Z-{e,mu},Z-HighMass{1,2}-e,{8,13}-W-mu,{8,13}-Z-Mass{1..4}-{e,mu},Z-{e,mu,mu-SOPHTY},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},{8,13}-WZ,8-ZZ-lv,8-WW-ll,Z-mu-Short}.in) \
	+ $(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-{{,8-}W-{e,mu},13-Z-{e,mu},Z-HighMass{1,2}-e,13-W-mu,{8,13}-Z-Mass{1..4}-{e,mu},Z-{e,mu,mu-SOPHTY},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},{8,13}-WZ,8-ZZ-lv,8-WW-ll,Z-mu-Short}.in) \
	$(shell echo Rivet/LHC{,-Dipole}-7-{W,Z}Gamma-{e,mu}.in) \
	$(shell echo Rivet/LHC{,-Dipole}-8-ZGamma-{e,mu}.in) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-{7-W-Jet-{1..3}-e,7-Z-Jet-{0..3}-e,7-Z-Jet-0-mu}.in) \
	$(shell echo Rivet/LHC{-Matchbox,-Matchbox-Powheg,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-{Z-b,Z-bb,8-Z-b,8-Z-bb,W-b,8-Z-jj}.in) \
	$(shell echo Rivet/LHC{,-Dipole}-{7,8,13}-PromptPhoton-{1..4}.in) Rivet/LHC-GammaGamma-7.in \
	$(shell echo Rivet/LHC{,-Powheg,-Dipole}-{7,8}-{DiPhoton-GammaGamma,DiPhoton-GammaJet}.in) \
	$(shell echo Rivet/LHC{,-Powheg,-Matchbox,-Matchbox-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-{ggH,VBF,WH,ZH}.in) \
	$(shell echo Rivet/LHC{,-Powheg,-Matchbox,-Matchbox-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-8-{{ggH,VBF,WH,ZH}{,-GammaGamma},ggH-WW}.in) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Dipole,-Dipole-MCatNLO,-Dipole-Matchbox-Powheg,-Merging}-ggHJet.in) \
	$(shell echo Rivet/ISR{,-Dipole}-{{30,44,53,62}-UE,{44,62}-Z-mu}.in Rivet/EHS{,-Dipole}-UE.in) \
	$(shell echo Rivet/SPS{,-Dipole}-{17.4-UE,200-Z-mu}.in ) \
	$(shell echo Rivet/Fermilab{,-Dipole}-{27.4,38.8}-Z-mu.in )


	Rivet-GammaGamma: Rivet-GammaGamma/done
	touch $@

	Rivet-GammaGamma/done: $(shell echo Rivet-GammaGamma-mumu-{3.5,4.5,5.5,6.5,7.5,9.0,12.5,17.5,30.0}.yoda )
	rm -rf Rivet-GammaGamma
	python/merge-GammaGamma GammaGamma
	rivet-mkhtml -o Rivet-GammaGamma GammaGamma.yoda:Hw
	touch $@

	Rivet-EE-Gamma: Rivet-EE-Gamma/done
	touch $@

	Rivet-EE-Gamma/done: $(shell echo Rivet-EE-Gamma-Direct-mumu-{161,172,183,189,196,206}.yoda ) \
	$(shell echo Rivet-EE-Gamma-Direct-tautau-{189,196,206}.yoda ) \
	$(shell echo Rivet-EE-Gamma-{Direct,Single-Resolved,Double-Resolved}-Jets-{198,206}.yoda )
	rm -rf Rivet-EE-Gamma
	python/merge-EE-Gamma EE-Gamma
	rivet-mkhtml -o Rivet-EE-Gamma EE-Gamma.yoda:Hw
	touch $@

	Rivet-EE : Rivet-EE/done
	touch $@

	Rivet-EE/done : $(shell echo Rivet{,-Powheg}-EE-{7.7,9.4,12,13,14,14.8,17,27.6,27.7,29,30.2,30.7,30,31.3,31.6,34,34.8,43.6,45,50,52,53.3,55,56,57,58,59.5,60.8,60,61.4,66,76,10,12.8,21.5,22,25,26.8,34.5,35,36.2,41,42.1,42.6,44,48.0,82,85,91,93.0,130,133,136,161,172,177,183,189,192,196,197,200,202,205,206,207,91-nopi}.yoda) \
	$(shell echo Rivet-EE-{10.5,11.96,12.8,13.96,16.86,21.84,26.8,28.48,35.44,48.0,97.0}-gg.yoda) \
	$(shell echo Rivet-EE-{10.52,10.52-sym,10.6,2.2,2.6,3.0,3.2,4.6,4.8,5.8,6.2,6.6,7.0,7.4,3.63,4.03,4.17,4.3,4.41,5.0,5.2,4.5,8.8,9.27,9.46,9.51,10.54,10.58,10.45,10.47,Upsilon,Upsilon2,Upsilon3,Upsilon4,Upsilon4-asym,Tau,Phi,Lambdac,Omega-Meson,Omega-Baryon,Eta,Xi0,Xic0,Xicp,Omegac0,Xim,JPsi,Psi2S,Psi2S-All,Psi3770}.yoda)
	rm -rf Rivet-EE
	python/merge-EE --with-gg --with-decay EE
	python/merge-EE Powheg-EE
	rivet-mkhtml -o Rivet-EE EE.yoda:Hw Powheg-EE.yoda:Hw-Powheg
	python/plot-EE Rivet-EE
	touch $@

	Rivet-LowEnergy-%.yoda:
	$(HWBUILD) -c .cache/$(subst .yoda,,$@) Rivet/$(subst .yoda,.in,$@)
	$(HWRUN) $(subst .yoda,.run,$@)

	Rivet-LowEnergy-EE-%:
	args="--process "$(word 1,$(subst -, ,$(subst Rivet-LowEnergy-EE-,,$@))); if [ -n "$(strip $(word 2,$(subst -, ,$(subst Rivet-LowEnergy-EE-,,$@))))" ]; then args+=" --flavour "$(word 2,$(subst -, ,$(subst Rivet-LowEnergy-EE-,,$@))); fi; OUTPUT=`python/LowEnergy-EE.py $$args --non-perturbative --perturbative`; $(MAKE) $$OUTPUT NUMEVENTS=$${NUMEVENTS:-10000};
	args="--process "$(word 1,$(subst -, ,$(subst Rivet-LowEnergy-EE-,,$@))); plots=`python/LowEnergy-EE.py $$args --plots`; python/mergeLowEnergy.py $(subst Rivet-LowEnergy-,,$@) $$plots; if [ -e LowEnergy-NonPerturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda ] && [ -e LowEnergy-Perturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-EE-$(subst Rivet-LowEnergy-EE-,,$@) LowEnergy-NonPerturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda:"Non-Pert" LowEnergy-Perturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda:"Pert" $$plots; elif [ -e LowEnergy-NonPerturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-EE-$(subst Rivet-LowEnergy-EE-,,$@) LowEnergy-NonPerturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda:"Non-Pert" $$plots; elif [ -e LowEnergy-Perturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-EE-$(subst Rivet-LowEnergy-EE-,,$@) LowEnergy-Perturbative-EE-$(subst Rivet-LowEnergy-EE-,,$@).yoda:"Pert" $$plots; fi

	Rivet-LowEnergy-Photon-%:
	args="--process "$(word 1,$(subst -, ,$(subst Rivet-LowEnergy-Photon-,,$@))); OUTPUT=`python/LowEnergy-Photon.py $$args --non-perturbative`; $(MAKE) $$OUTPUT NUMEVENTS=$${NUMEVENTS:-10000};
	args="--process "$(word 1,$(subst -, ,$(subst Rivet-LowEnergy-Photon-,,$@))); plots=`python/LowEnergy-Photon.py $$args --plots`; python/mergeLowEnergy.py $(subst Rivet-LowEnergy-Photon-,GammaGamma-,$@) $$plots; if [ -e LowEnergy-NonPerturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda ] && [ -e LowEnergy-Perturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-Photon-$(subst Rivet-LowEnergy-Photon-,,$@) LowEnergy-NonPerturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda:"Non-Pert" LowEnergy-Perturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda:"Pert" $$plots; elif [ -e LowEnergy-NonPerturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-Photon-$(subst Rivet-LowEnergy-Photon-,,$@) LowEnergy-NonPerturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda:"Non-Pert" $$plots; elif [ -e LowEnergy-Perturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda ]; then rivet-mkhtml -o Rivet-LowEnergy-Photon-$(subst Rivet-LowEnergy-Photon-,,$@) LowEnergy-Perturbative-GammaGamma-$(subst Rivet-LowEnergy-Photon-,,$@).yoda:"Pert" $$plots; fi

	Rivet-R:
	OUTPUT=`python/R.py --perturbative --non-perturbative`; $(MAKE) $$OUTPUT NUMEVENTS=$${NUMEVENTS:-10000};
	plots=`python/R.py --perturbative --non-perturbative --plots`; python/mergeLowEnergy.py R $$plots; rivet-mkhtml -o Rivet-R LowEnergy-EE-Perturbative-R.yoda:"Pert" LowEnergy-EE-NonPerturbative-R.yoda:"Non-Pert" $$plots

	Rivet-DIS : Rivet-DIS/done
	touch $@

	Rivet-DIS/done: $(shell echo Rivet{-DIS,-DIS-NoME,-Powheg-DIS}-{e--LowQ2,e+-LowQ2,e+-HighQ2}.yoda)
	rm -rf Rivet-DIS
	python/merge-DIS DIS
	python/merge-DIS Powheg-DIS
	python/merge-DIS DIS-NoME
	rivet-mkhtml -o Rivet-DIS DIS.yoda:Hw Powheg-DIS.yoda:Hw-Powheg DIS-NoME.yoda:Hw-NoME
	touch $@



	Rivet-TVT-EW : Rivet-TVT-EW/done
	touch $@

	Rivet-TVT-EW/done: $(shell echo Rivet{,-Powheg}-TVT-{Run-I-Z,Run-I-W,Run-I-WZ,Run-II-Z-{e,{,LowMass-,HighMass-}mu},Run-II-W}.yoda)
	rm -rf Rivet-TVT-EW
	python/merge-TVT-EW TVT
	python/merge-TVT-EW Powheg-TVT
	rivet-mkhtml -o Rivet-TVT-EW TVT-EW.yoda:Hw Powheg-TVT-EW.yoda:Hw-Powheg
	touch $@

	Rivet-TVT-Photon : Rivet-TVT-Photon/done
	touch $@

	Rivet-TVT-Photon/done: $(shell echo Rivet{,-Powheg}-TVT-Run-II-{DiPhoton-GammaGamma,DiPhoton-GammaJet}.yoda Rivet-TVT-Run-II-PromptPhoton.yoda)
	rm -rf Rivet-TVT-Photon
	python/merge-TVT-Photon TVT
	python/merge-TVT-Photon Powheg-TVT
	rivet-mkhtml -o Rivet-TVT-Photon TVT-Photon.yoda:Hw Powheg-TVT-Photon.yoda:Hw-Powheg
	touch $@



	Rivet-TVT-Jets : Rivet-TVT-Jets/done
	touch $@

	Rivet-TVT-Jets/done: $(shell echo Rivet-TVT-{Run-II-Jets-{0..11},Run-I-Jets-{1..8}}.yoda ) \
	$(shell echo Rivet-TVT-{630-Jets-{1..3},300-Jets-1,900-Jets-1}.yoda ) \
	$(shell echo Rivet-TVT-{Run-I,Run-II,300,630,900}-UE.yoda)
	rm -rf Rivet-TVT-Jets
	python/merge-TVT-Jets TVT
	rivet-mkhtml -o Rivet-TVT-Jets TVT-Jets.yoda:Hw
	touch $@

	Rivet-Fixed : Rivet-Fixed/done
	touch $@


	Rivet-Fixed/done : $(shell echo Rivet-SppS-{53,63,200,500,546,900}-UE.yoda ) \
	$(shell echo Rivet-ISR-{{30,44,53,62}-UE,{44,62}-Z-mu}.yoda ) Rivet-EHS-UE.yoda \
	$(shell echo Rivet-Star-{UE,Jets-{1..4}}.yoda ) \
	$(shell echo Rivet-SPS-{17.4-UE,200-Z-mu}.yoda ) \
	$(shell echo Rivet-Fermilab-{27.4,38.8}-Z-mu.yoda )
	rm -rf Rivet-Fixed
	python/merge-Fixed Fixed
	rivet-mkhtml -o Rivet-Fixed Fixed.yoda:Hw
	touch $@

	Rivet-LHC-Jets : Rivet-LHC-Jets/done
	touch $@

	Rivet-LHC-Jets/done : \
	$(shell echo Rivet-LHC-7-DiJets-{1..7}-{A,B,C}.yoda ) \
	$(shell echo Rivet-LHC-13-DiJets-{{1..11}-A,{6..11}-B}.yoda ) \
	$(shell echo Rivet-LHC-{7,8,13}-Jets-{0..10}.yoda ) \
	$(shell echo Rivet-LHC-2760-Jets-{1..3}.yoda ) \
	$(shell echo Rivet-LHC-{900,2360,2760,7,8,13}-UE.yoda ) \
	$(shell echo Rivet-LHC-{7,13}-UE-Cent.yoda ) \
	$(shell echo Rivet-LHC-{900,7,13}-UE-Long.yoda ) \
	$(shell echo Rivet-LHC-7-Charm-{1..5}.yoda ) \
	$(shell echo Rivet-LHC-7-Bottom-{0..9}.yoda ) \
	$(shell echo Rivet-LHC-{7,8,13}-Top-{L,SL}.yoda ) \
	$(shell echo Rivet-LHC-{8,13}-Top-All.yoda )
	rm -rf Rivet-LHC-Jets
	python/merge-LHC-Jets LHC
	rivet-mkhtml -o Rivet-LHC-Jets LHC-Jets.yoda:Hw
	touch $@


	Rivet-LHC-EW : Rivet-LHC-EW/done
	touch $@

	Rivet-LHC-EW/done: \
	- $(shell echo Rivet{,-Powheg}-LHC-{13-Z-{e,mu},{8,13}-W-mu,Z-HighMass{1,2}-e,{8,13}-Z-Mass{1..4}-{e,mu},W-{e,mu},Z-{e,mu,mu-SOPHTY},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},{8,13}-WZ,8-ZZ-lv,8-WW-ll,Z-mu-Short}.yoda) \
	+ $(shell echo Rivet{,-Powheg}-LHC-{13-Z-{e,mu},13-W-mu,Z-HighMass{1,2}-e,{8,13}-Z-Mass{1..4}-{e,mu},{,8-}W-{e,mu},Z-{e,mu,mu-SOPHTY},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},{8,13}-WZ,8-ZZ-lv,8-WW-ll,Z-mu-Short}.yoda) \
	$(shell echo Rivet-LHC-{7-W-Jet-{1..3}-e,7-Z-Jet-{0..3}-e,7-Z-Jet-0-mu}.yoda) \
	$(shell echo Rivet-LHC-7-{W,Z}Gamma-{e,mu}.yoda) \
	$(shell echo Rivet-LHC-8-ZGamma-{e,mu}.yoda)
	rm -rf Rivet-LHC-EW;
	python/merge-LHC-EW LHC
	python/merge-LHC-EW Powheg-LHC
	rivet-mkhtml -o Rivet-LHC-EW LHC-EW.yoda:Hw Powheg-LHC-EW.yoda:Hw-Powheg \
	Rivet-LHC-Z-mu-SOPHTY.yoda:Hw Rivet-Powheg-LHC-Z-mu-SOPHTY.yoda:Hw-Powheg
	touch $@




	Rivet-LHC-Photon : Rivet-LHC-Photon/done
	touch $@

	Rivet-LHC-Photon/done: \
	$(shell echo Rivet-LHC-{7,8,13}-PromptPhoton-{1..4}.yoda) \
	Rivet-LHC-GammaGamma-7.yoda \
	$(shell echo Rivet{,-Powheg}-LHC-{7,8}-{DiPhoton-GammaGamma,DiPhoton-GammaJet}.yoda)
	rm -rf Rivet-LHC-Photon
	python/merge-LHC-Photon LHC
	python/merge-LHC-Photon Powheg-LHC
	rivet-mkhtml -o Rivet-LHC-Photon LHC-Photon.yoda:Hw Powheg-LHC-Photon.yoda:Hw-Powheg
	touch $@




	Rivet-LHC-Higgs : Rivet-LHC-Higgs/done
	touch $@

	Rivet-LHC-Higgs/done: \
	$(shell echo Rivet{,-Powheg}-LHC-{ggH,VBF,WH,ZH}.yoda) \
	$(shell echo Rivet{,-Powheg}-LHC-8-{{ggH,VBF,WH,ZH}{,-GammaGamma},ggH-WW}.yoda) \
	Rivet-LHC-ggHJet.yoda
	yodamerge --add Rivet-Powheg-LHC-8-{ggH{-GammaGamma,-WW,},{VBF,ZH,WH}{,-GammaGamma}}.yoda -o Powheg-LHC-Higgs.yoda
	yodamerge --add Rivet-LHC-8-{ggH{-GammaGamma,-WW,},{VBF,ZH,WH}{,-GammaGamma}}.yoda -o LHC-Higgs.yoda
	rm -rf Rivet-LHC-Higgs
	rivet-mkhtml -o Rivet-LHC-Higgs Powheg-LHC-Higgs.yoda:Hw-Powheg LHC-Higgs.yoda:Hw\
	Rivet-Powheg-LHC-ggH.yoda:gg-Powheg Rivet-LHC-ggH.yoda:gg Rivet-LHC-ggHJet.yoda:HJet \
	Rivet-Powheg-LHC-VBF.yoda:VBF-Powheg Rivet-LHC-VBF.yoda:VBF Rivet-LHC-WH.yoda:WH Rivet-LHC-ZH.yoda:ZH \
	Rivet-Powheg-LHC-WH.yoda:WH-Powheg Rivet-Powheg-LHC-ZH.yoda:ZH-Powheg
	touch $@





	clean-local:
	rm -f .out .log .tex .top .run .dump .mult .Bmult .yoda Rivet/.in anatohepmc.txt hepmctoana.txt
	rm -rf Rivet-*

	distclean-local:
	rm -rf .cache
	diff --git a/Tests/Rivet/LHC/LHC-8-W-mu.in b/Tests/Rivet/LHC/LHC-8-W-e.in
	copy from Tests/Rivet/LHC/LHC-8-W-mu.in
	copy to Tests/Rivet/LHC/LHC-8-W-e.in
	--- a/Tests/Rivet/LHC/LHC-8-W-mu.in
	+++ b/Tests/Rivet/LHC/LHC-8-W-e.in
	@@ -1,6 +1,6 @@
	# -- ThePEG-repository --
	##################################################
	# select the analyses
	##################################################
	# CMS W+jets
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2016_I1491953
	+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2018_I1635273
	diff --git a/Tests/Rivet/LHC/LHC-8-Z-Mass3-e.in b/Tests/Rivet/LHC/LHC-8-Z-Mass3-e.in
	--- a/Tests/Rivet/LHC/LHC-8-Z-Mass3-e.in
	+++ b/Tests/Rivet/LHC/LHC-8-Z-Mass3-e.in
	@@ -1,12 +1,14 @@
	# -- ThePEG-repository --
	##################################################
	# select the analyses
	##################################################
	# ATLAS Z pT and phi*
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2015_I1408516:LMODE=EL
	# ATLAS high mass drell-yan
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2016_I1467454:LMODE=EL
	# ATLAS splittings
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2017_I1589844:LMODE=EL
	# CMS Z+b
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2017_I1499471
	\ No newline at end of file
	+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2017_I1499471
	+# ATLAS Z+jets
	+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2019_I1744201
	diff --git a/Tests/python/make_input_files.py b/Tests/python/make_input_files.py
	--- a/Tests/python/make_input_files.py
	+++ b/Tests/python/make_input_files.py
	@@ -1,2005 +1,2005 @@
	#! /usr/bin/env python
	from __future__ import print_function
	import logging,sys,os
	from string import Template
	import sys
	if sys.version_info[:3] < (2,4,0):
	print ("rivet scripts require Python version >= 2.4.0... exiting")
	sys.exit(1)

	if __name__ == "__main__":
	import logging
	from optparse import OptionParser, OptionGroup
	parser = OptionParser(usage="%prog name [...]")


	simulation=""

	numberOfAddedProcesses=0
	def addProcess(thefactory,theProcess,Oas,Oew,scale,mergedlegs,NLOprocesses):
	global numberOfAddedProcesses
	global simulation
	numberOfAddedProcesses+=1
	res ="set "+thefactory+":OrderInAlphaS "+Oas+"\n"
	res+="set "+thefactory+":OrderInAlphaEW "+Oew+"\n"
	res+="do "+thefactory+":Process "+theProcess+" "
	if ( mergedlegs != 0 ):
	if simulation!="Merging":
	print ("simulation is not Merging, trying to add merged legs.")
	sys.exit(1)
	res+="["
	for j in range(mergedlegs):
	res+=" j "
	res+="]"
	res+="\n"
	if (NLOprocesses!=0):
	if simulation!="Merging":
	print ("simulation is not Merging, trying to add NLOProcesses.")
	sys.exit(1)
	res+="set MergingFactory:NLOProcesses %s \n" % NLOprocesses
	if ( scale != "" ):
	res+="set "+thefactory+":ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/"+scale+"\n"
	return res


	def addLeptonPairCut(minmass,maxmass):
	return "set /Herwig/Cuts/LeptonPairMassCut:MinMass "+minmass+"GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass "+maxmass+"GeV\n"

	didaddfirstjet=False
	def addFirstJet(ptcut):
	global didaddfirstjet
	if(didaddfirstjet):
	logging.error("Can only add jetcut once.")
	sys.exit(1)

	res="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	res+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	res+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	if(ptcut!=""):
	res+="set /Herwig/Cuts/FirstJet:PtMin "+ptcut+"*GeV\n"
	didaddfirstjet=True
	return res



	didaddsecondjet=False
	def addSecondJet(ptcut):
	global didaddsecondjet
	if(didaddsecondjet):
	logging.error("Can only add second jetcut once.")
	sys.exit(1)
	res="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	res+="set /Herwig/Cuts/SecondJet:PtMin "+ptcut+"*GeV\n"
	didaddsecondjet=True
	return res


	didaddjetpair=False
	def addJetPairCut(minmass):
	global didaddjetpair
	if(didaddjetpair):
	logging.error("Can only add second jetcut once.")
	sys.exit(1)
	res="""\
	create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so
	set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet
	set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet
	insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass
	set /Herwig/Cuts/JetPairMass:MassMin {mm}*GeV
	""".format(mm=minmass)
	didaddjetpair=True
	return res

	addedBRReweighter=False
	def addBRReweighter():
	global addedBRReweighter
	if(addedBRReweighter):
	logging.error("Can only add BRReweighter once.")
	sys.exit(1)
	res="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	res+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	addedBRReweighter=True
	return res

	def setHardProcessWidthToZero(list1):
	res=""
	for i in list1:
	res+="set /Herwig/Particles/"+i+":HardProcessWidth 0.\n"
	return res

	selecteddecaymode=False
	def selectDecayMode(particle,decaymodes):
	global selecteddecaymode
	res="do /Herwig/Particles/"+particle+":SelectDecayModes"
	for decay in decaymodes:
	res+=" /Herwig/Particles/"+particle+"/"+decay
	res+="\n"
	selecteddecaymode=True
	return res

	def jet_kt_cut(energy):
	return "set /Herwig/Cuts/JetKtCut:MinKT {E}*GeV\n".format(E=energy)

	def mhatmin_cut(energy):
	return "set /Herwig/Cuts/Cuts:MHatMin {E}*GeV\n".format(E=energy)

	def mhat_minm_maxm(e1,e2,e3):
	return """\
	set /Herwig/Cuts/Cuts:MHatMin {e1}*GeV
	set /Herwig/Cuts/MassCut:MinM {e2}*GeV
	set /Herwig/Cuts/MassCut:MaxM {e3}*GeV
	""".format(**locals())

	def collider_lumi(energy):
	return "set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy {E}*GeV\n".format(E=energy)

	def insert_ME(me,process=None,ifname='Process',subprocess="SubProcess"):
	result = "insert /Herwig/MatrixElements/{subprocess}:MatrixElements 0 /Herwig/MatrixElements/{me}\n".format(**locals())
	if process is not None:
	result += "set /Herwig/MatrixElements/{me}:{ifname} {process}".format(**locals())
	return result

	def particlegroup(factory,name,*particles):
	directory="MatrixElements/Matchbox"
	if(factory!="Factory") : directory="Merging"
	result = ["do /Herwig/{dir}/{fact}:StartParticleGroup {n}".format(n=name,fact=factory,dir=directory)]
	for p in particles:
	result.append(
	"insert /Herwig/{dir}/{fact}:ParticleGroup 0 /Herwig/Particles/{p}".format(p=p,fact=factory,dir=directory)
	)
	result.append("do /Herwig/{dir}/{fact}:EndParticleGroup".format(fact=factory,dir=directory))
	return '\n'.join(result)


	ME_Upsilon = """\
	create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so
	set /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)
	set /Herwig/MatrixElements/MEUpsilon:Coupling 96.72794
	""" + insert_ME("MEUpsilon")


	(opts, args) = parser.parse_args()
	## Check args
	if len(args) != 1:
	logging.error("Must specify at least input file")
	sys.exit(1)

	name = args[0]
	print (name)


	# select the template to load
	# collider
	KNOWN_COLLIDERS = [
	"EE-Gamma",
	"BFactory",
	"EE",
	"DIS",
	"TVT",
	"LHC-GammaGamma",
	"LHC",
	"ISR",
	"Fermilab",
	"SppS",
	"Star",
	"SPS",
	"GammaGamma",
	]
	collider = ""
	for cand_collider in KNOWN_COLLIDERS:
	if cand_collider in name:
	collider = cand_collider
	break
	del cand_collider
	if "EHS" in name : collider="SPS"
	assert collider
	have_hadronic_collider = collider in ["TVT","LHC","ISR","SppS","Star","SPS","Fermilab"]
	print (collider)
	thefactory="Factory"

	parameters = {
	'shower' : '',
	'bscheme' : '',
	}
	# istart determines how many name parts need to be skipped
	istart = 1
	# Dipole shower with Matchbox Powheg
	if "Dipole-Matchbox-Powheg" in name :
	istart = 4
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/Powheg-DipoleShower.in\n"

	# Dipole shower with internal Powheg - Todo: Finish modifying template files.
	'''
	elif "Dipole-Powheg" in name :
	istart = 3
	simulation="Powheg"
	parameters["shower"] = "set /Herwig/EventHandlers/EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler\nread snippets/Dipole_AutoTunes_gss.in\n"
	'''

	# Dipole shower with MCatNLO
	elif "Dipole-MCatNLO" in name :
	istart = 3
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/MCatNLO-DipoleShower.in\n"

	# Dipole shower with Matchbox LO
	elif "Dipole-Matchbox-LO" in name :
	istart = 4
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/LO-DipoleShower.in\n"

	# Dipole shower with internal LO
	elif "Dipole" in name :
	istart = 2
	simulation=""
	parameters["shower"] = "set /Herwig/EventHandlers/EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler\nread snippets/Dipole_AutoTunes_gss.in\n"

	# AO shower with Matchbox Powheg
	elif "Matchbox-Powheg" in name :
	istart = 3
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/Powheg-DefaultShower.in\n"

	# AO shower with MCatNLO
	elif "Matchbox" in name :
	istart = 2
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/MCatNLO-DefaultShower.in\n"

	# AO shower with inernal Powheg
	elif "Powheg" in name :
	istart = 2
	simulation="Powheg"

	# Dipole shower with merging
	elif "Merging" in name :
	istart = 2
	simulation="Merging"
	thefactory="MergingFactory"

	# Flavour settings for Matchbox
	if simulation=="Matchbox" :
	parameters["bscheme"] = "read Matchbox/FiveFlavourScheme.in\n"

	if "Dipole" in parameters["shower"] :
	parameters["bscheme"] += "read Matchbox/FiveFlavourNoBMassScheme.in\n"

	if collider not in ['DIS','EE'] :
	parameters["nlo"] = "read Matchbox/MadGraph-OpenLoops.in\n"

	# Flavour settings for dipole shower with internal ME
	if simulation=="" and "Dipole" in parameters["shower"] :
	parameters["bscheme"] = "read snippets/DipoleShowerFiveFlavours.in"



	# find the template
	if simulation=="" :
	if collider=="LHC-GammaGamma" :
	istart += 1
	templateName="Hadron-Gamma.in"
	elif have_hadronic_collider :
	templateName="Hadron.in"
	elif collider=="EE-Gamma" :
	istart+=1
	if("Direct" in name) :
	templateName="EE-Gamma-Direct.in"
	elif("Single-Resolved" in name) :
	templateName="EE-Gamma-Single-Resolved.in"
	elif("Double-Resolved" in name) :
	templateName="EE-Gamma-Double-Resolved.in"
	else :
	print ("Unknown type of EE-Gamma event ",name)
	quit()
	elif collider=="GammaGamma" :
	templateName="GammaGamma.in"
	elif collider != "BFactory" :
	templateName= "%s.in" % collider
	else :
	templateName= "EE.in"
	else :
	if have_hadronic_collider :
	templateName= "Hadron-%s.in" % simulation
	elif collider != "BFactory" :
	templateName= "%s-%s.in" % (collider,simulation)
	else :
	templateName= "EE-%s.in" % simulation

	# settings for four flavour scheme
	fourFlavour="""
	read Matchbox/FourFlavourScheme.in
	{bjetgroup}
	set /Herwig/Cuts/MatchboxJetMatcher:Group bjet
	""".format(bjetgroup=particlegroup(thefactory,'bjet','b','bbar','c', 'cbar',
	's','sbar','d','dbar','u','ubar','g'))

	# work out the name of the parameter file
	parameterName="-".join(name.split("-")[istart:])
	del istart

	class StringBuilder(object):
	"""
	Avoid expensive string additions until the end
	by building up a list first.

	This helper class avoids rewriting all the += lower down
	to list operations.
	"""
	def __init__(self, init = None):
	self.lines = [] if init is None else [init]

	def __iadd__(self, line):
	self.lines.append(line)
	return self

	def __str__(self):
	return '\n'.join(self.lines)

	# work out the process and parameters
	process=StringBuilder()

	# DIS
	if(collider=="DIS") :
	if(simulation=="") :
	if "NoME" in parameterName :
	process = StringBuilder("set /Herwig/Shower/ShowerHandler:HardEmission None")
	parameterName=parameterName.replace("NoME-","")
	else :
	process = StringBuilder("")
	elif(simulation=="Powheg") :
	process = StringBuilder("")
	elif(simulation=="Matchbox" ) :
	if "e-" in parameterName :
	process = StringBuilder(addProcess(thefactory,"e- p -> e- j","0","2","",0,0))
	else :
	process = StringBuilder(addProcess(thefactory,"e+ p -> e+ j","0","2","",0,0))
	elif(simulation=="Merging" ) :
	if "e-" in parameterName :
	process = StringBuilder(addProcess(thefactory,"e- p -> e- j","0","2","",2,2))
	else :
	process = StringBuilder(addProcess(thefactory,"e+ p -> e+ j","0","2","",2,2))

	# EE
	elif(collider=="EE") :
	if(simulation=="") :
	if "gg" in parameterName :
	process = StringBuilder("create Herwig::MEee2Higgs2SM /Herwig/MatrixElements/MEee2Higgs2SM\n")
	process+=insert_ME("MEee2Higgs2SM","Gluon","Allowed")
	else :
	process = StringBuilder(insert_ME("MEee2gZ2qq"))
	try :
	ecms = float(parameterName)
	if(ecms<=3.75) :
	process+= "set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 3\n"
	elif(ecms<=10.6) :
	process+= "set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 4\n"
	except :
	pass
	elif(simulation=="Powheg") :
	process = StringBuilder()
	try :
	ecms = float(parameterName)
	if(ecms<=3.75) :
	process+= "set /Herwig/MatrixElements/PowhegMEee2gZ2qq:MaximumFlavour 3\n"
	elif(ecms<=10.6) :
	process+= "set /Herwig/MatrixElements/PowhegMEee2gZ2qq:MaximumFlavour 4\n"
	except :
	pass
	elif(simulation=="Matchbox" ) :
	try :
	ecms = float(parameterName)
	if(ecms<=3.75) :
	process = StringBuilder(addProcess(thefactory,"e- e+ -> u ubar","0","2","",0,0))
	process+=addProcess(thefactory,"e- e+ -> d dbar","0","2","",0,0)
	process+=addProcess(thefactory,"e- e+ -> s sbar","0","2","",0,0)
	elif(ecms<=10.6) :
	process = StringBuilder(addProcess(thefactory,"e- e+ -> u ubar","0","2","",0,0))
	process+=addProcess(thefactory,"e- e+ -> d dbar","0","2","",0,0)
	process+=addProcess(thefactory,"e- e+ -> c cbar","0","2","",0,0)
	process+=addProcess(thefactory,"e- e+ -> s sbar","0","2","",0,0)
	else :
	process = StringBuilder(addProcess(thefactory,"e- e+ -> j j","0","2","",0,0))
	except:
	process = StringBuilder(addProcess(thefactory,"e- e+ -> j j","0","2","",0,0))
	elif(simulation=="Merging" ) :
	try :
	ecms = float(parameterName)
	if(ecms<=10.1) :
	process = StringBuilder(addProcess(thefactory,"e- e+ -> j j","0","2","",2,2))
	process+="read Matchbox/FourFlavourScheme.in"
	else :
	process = StringBuilder(addProcess(thefactory,"e- e+ -> j j","0","2","",2,2))
	except:
	process = StringBuilder(addProcess(thefactory,"e- e+ -> j j","0","2","",2,2))
	# EE-Gamma
	elif(collider=="EE-Gamma") :
	if(simulation=="") :
	if("mumu" in parameterName) :
	process = StringBuilder(insert_ME("MEgg2ff","Muon"))
	process +="set /Herwig/Cuts/Cuts:MHatMin 3.\n"
	elif( "tautau" in parameterName) :
	process = StringBuilder(insert_ME("MEgg2ff","Tau"))
	process +="set /Herwig/Cuts/Cuts:MHatMin 3.\n"
	elif( "Jets" in parameterName) :
	if("Direct" in parameterName ) :
	process = StringBuilder(insert_ME("MEgg2ff","Quarks"))
	elif("Single-Resolved" in parameterName ) :
	process = StringBuilder(insert_ME("MEGammaP2Jets",None,"Process","SubProcess"))
	process+= insert_ME("MEGammaP2Jets",None,"Process","SubProcess2")
	else :
	process = StringBuilder(insert_ME("MEQCD2to2"))
	process+="insert /Herwig/Cuts/Cuts:OneCuts[0] /Herwig/Cuts/JetKtCut"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 3."
	else :
	print ("process not supported for Gamma Gamma processes at EE")
	quit()
	else :
	print ("Only internal matrix elements currently supported for Gamma Gamma processes at EE")
	quit()
	elif(collider=="GammaGamma") :
	if(simulation=="") :
	if("mumu" in parameterName) :
	process = StringBuilder(insert_ME("MEgg2ff"))
	process +="set /Herwig/MatrixElements/MEgg2ff:Process Muon\n"
	process +="set /Herwig/Cuts/Cuts:MHatMin 3.\n"
	else :
	print ("process not supported for Gamma Gamma processes at EE")
	quit()
	else :
	print ("Only internal matrix elements currently supported for Gamma Gamma processes at EE")
	quit()
	# TVT
	elif(collider=="TVT") :
	process = StringBuilder("set /Herwig/Generators/EventGenerator:EventHandler:BeamB /Herwig/Particles/pbar-\n")
	if "Run-II" in parameterName : process+=collider_lumi(1960.0)
	elif "Run-I" in parameterName : process+=collider_lumi(1800.0)
	elif "900" in parameterName : process+=collider_lumi(900.0)
	elif "630" in parameterName : process+=collider_lumi(630.0)
	elif "300" in parameterName : process+=collider_lumi(300.0)

	if(simulation=="") :
	if "PromptPhoton" in parameterName :
	process+=insert_ME("MEGammaJet")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 15.\n"
	elif "DiPhoton-GammaGamma" in parameterName :
	process+=insert_ME("MEGammaGamma")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif "DiPhoton-GammaJet" in parameterName :
	process+=insert_ME("MEGammaJet")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	elif "UE" in parameterName :
	if "Dipole" in parameters["shower"]:
	process+="read snippets/MB-DipoleShower.in\n"
	else:
	process+="read snippets/MB.in\n"
	process+="read snippets/Diffraction.in\n"

	process += "set /Herwig/Decays/DecayHandler:LifeTimeOption 0\n"
	process += "set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm\n"
	elif "Jets" in parameterName :
	process+=insert_ME("MEQCD2to2")
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "Run-II-Jets-10" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(500.)
	elif "Run-II-Jets-11" in parameterName: process+=jet_kt_cut( 30.)+mhatmin_cut(900.)
	elif "Run-I-Jets-1" in parameterName : process+=jet_kt_cut( 20.)
	elif "Run-I-Jets-2" in parameterName : process+=jet_kt_cut( 40.)
	elif "Run-I-Jets-3" in parameterName : process+=jet_kt_cut( 65.)
	elif "Run-I-Jets-4" in parameterName : process+=jet_kt_cut( 90.)
	elif "Run-I-Jets-5" in parameterName : process+=jet_kt_cut(160.)
	elif "Run-I-Jets-6" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(100.)
	elif "Run-I-Jets-7" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(400.)
	elif "Run-I-Jets-8" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(700.)
	elif "Run-II-Jets-0" in parameterName : process+=jet_kt_cut( 15.)
	elif "Run-II-Jets-1" in parameterName : process+=jet_kt_cut( 25.)
	elif "Run-II-Jets-2" in parameterName : process+=jet_kt_cut( 40.)
	elif "Run-II-Jets-3" in parameterName : process+=jet_kt_cut( 60.)
	elif "Run-II-Jets-4" in parameterName : process+=jet_kt_cut( 85.)
	elif "Run-II-Jets-5" in parameterName : process+=jet_kt_cut(110.)
	elif "Run-II-Jets-6" in parameterName : process+=jet_kt_cut(160.)
	elif "Run-II-Jets-7" in parameterName : process+=jet_kt_cut(250.)
	elif "Run-II-Jets-8" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(100.)
	elif "Run-II-Jets-9" in parameterName : process+=jet_kt_cut( 30.)+mhatmin_cut(300.)
	elif "900-Jets-1" in parameterName : process+=jet_kt_cut( 10.)
	elif "300-Jets-1" in parameterName : process+=jet_kt_cut( 6.)
	elif "630-Jets-1" in parameterName : process+=jet_kt_cut( 20.)
	elif "630-Jets-2" in parameterName : process+=jet_kt_cut( 40.)
	elif "630-Jets-3" in parameterName : process+=jet_kt_cut( 75.)
	elif "900-Jets-1" in parameterName : process+=jet_kt_cut( 10.)
	elif "Run-I-WZ" in parameterName :
	process+=insert_ME("MEqq2W2ff","Electron")
	process+=insert_ME("MEqq2gZ2ff","Electron")
	elif "Run-II-W" in parameterName or "Run-I-W" in parameterName :
	process+=insert_ME("MEqq2W2ff","Electron")
	elif "Run-II-Z-e" in parameterName or "Run-I-Z" in parameterName :
	process +=insert_ME("MEqq2gZ2ff","Electron")
	elif "Run-II-Z-LowMass-mu" in parameterName :
	process +=insert_ME("MEqq2gZ2ff","Muon")
	process+=addLeptonPairCut("25","70")
	elif "Run-II-Z-HighMass-mu" in parameterName :
	process +=insert_ME("MEqq2gZ2ff","Muon")
	process+=addLeptonPairCut("150","600")
	elif "Run-II-Z-mu" in parameterName :
	process +=insert_ME("MEqq2gZ2ff","Muon")
	elif(simulation=="Powheg") :
	if "Run-I-WZ" in parameterName :
	process+=insert_ME("PowhegMEqq2W2ff","Electron")
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	elif "Run-II-W" in parameterName or "Run-I-W" in parameterName :
	process+=insert_ME("PowhegMEqq2W2ff","Electron")
	elif "Run-II-Z-e" in parameterName or "Run-I-Z" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	elif "Run-II-Z-LowMass-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	process+=addLeptonPairCut("25","70")
	elif "Run-II-Z-HighMass-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	process+=addLeptonPairCut("150","600")
	elif "Run-II-Z-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "DiPhoton-GammaGamma" in parameterName :
	process+=insert_ME("MEGammaGammaPowheg","GammaGamma")
	process+=insert_ME("MEGammaGamma","gg")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+=jet_kt_cut(5.)
	parameterName=parameterName.replace("-GammaGamma","")
	elif "DiPhoton-GammaJet" in parameterName :
	process+=insert_ME("MEGammaGammaPowheg","VJet")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+=jet_kt_cut(5.)
	parameterName=parameterName.replace("-GammaJet","")
	elif(simulation=="Matchbox" or simulation=="Merging" ) :
	if "Jets" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p -> j j","2","0","MaxJetPtScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p -> j j","2","0","MaxJetPtScale",1,0)
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "Run-II-Jets-10" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("500")
	elif "Run-II-Jets-11" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("900")
	elif "Run-II-Jets-12" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("300")
	elif "Run-I-Jets-1" in parameterName : process+=addFirstJet("20")
	elif "Run-I-Jets-2" in parameterName : process+=addFirstJet("40")
	elif "Run-I-Jets-3" in parameterName : process+=addFirstJet("65")
	elif "Run-I-Jets-4" in parameterName : process+=addFirstJet("90")
	elif "Run-I-Jets-5" in parameterName : process+=addFirstJet("160")
	elif "Run-I-Jets-6" in parameterName : process+=addFirstJet("30")+addSecondJet("25")+addJetPairCut("100")
	elif "Run-I-Jets-7" in parameterName : process+=addFirstJet("30")+addSecondJet("25")+addJetPairCut("400")
	elif "Run-I-Jets-8" in parameterName : process+=addFirstJet("30")+addSecondJet("25")+addJetPairCut("700")
	elif "Run-II-Jets-0" in parameterName : process+=addFirstJet("15")
	elif "Run-II-Jets-1" in parameterName : process+=addFirstJet("25")
	elif "Run-II-Jets-2" in parameterName : process+=addFirstJet("40")
	elif "Run-II-Jets-3" in parameterName : process+=addFirstJet("60")
	elif "Run-II-Jets-4" in parameterName : process+=addFirstJet("85")
	elif "Run-II-Jets-5" in parameterName : process+=addFirstJet("110")
	elif "Run-II-Jets-6" in parameterName : process+=addFirstJet("160")
	elif "Run-II-Jets-7" in parameterName : process+=addFirstJet("250")
	elif "Run-II-Jets-8" in parameterName : process+=addFirstJet("30")+addSecondJet("25")+addJetPairCut("100")
	elif "Run-II-Jets-9" in parameterName : process+=addFirstJet("30")+addSecondJet("25")+addJetPairCut("300")
	elif "900-Jets-1" in parameterName : process+=addFirstJet("10")
	elif "300-Jets-1" in parameterName : process+=addFirstJet("6")
	elif "630-Jets-1" in parameterName : process+=addFirstJet("20")
	elif "630-Jets-2" in parameterName : process+=addFirstJet("40")
	elif "630-Jets-3" in parameterName : process+=addFirstJet("75")
	elif "900-Jets-1" in parameterName : process+=addFirstJet("10")
	else :
	logging.error("Exit 00007")
	sys.exit(1)
	elif "Run-I-WZ" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar e+ e-","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p pbar e+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p pbar e- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'epm','e+','e-')
	process+=particlegroup(thefactory,'epmnu','e+','e-','nu_e','nu_ebar')
	process+=addProcess(thefactory,"p pbar epm epmnu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "Run-II-W" in parameterName or "Run-I-W" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar e+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p pbar e- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'epm','e+','e-')
	process+=addProcess(thefactory,"p pbar epm nu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "Run-II-Z-e" in parameterName or "Run-I-Z" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p pbar e+ e-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "Run-II-Z-LowMass-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("25","70")
	elif "Run-II-Z-HighMass-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("150","600")
	elif "Run-II-Z-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p pbar mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	# Star
	elif(collider=="Star" ) :
	process = StringBuilder("set /Herwig/Decays/DecayHandler:LifeTimeOption 0\n")
	process+= "set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm\n"
	process+= "set /Herwig/Generators/EventGenerator:EventHandler:BeamB /Herwig/Particles/p+\n"
	process+= collider_lumi(200.0)
	process+= "set /Herwig/Cuts/Cuts:X2Min 0.01\n"
	if(simulation=="") :
	if "UE" in parameterName :
	if "Dipole" in parameters["shower"]:
	process+="read snippets/MB-DipoleShower.in\n"
	else:
	process+="read snippets/MB.in\n"
	process+="read snippets/Diffraction.in\n"


	else :
	process+=insert_ME("MEQCD2to2")
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "Jets-1" in parameterName : process+=jet_kt_cut(2.)
	elif "Jets-2" in parameterName : process+=jet_kt_cut(5.)
	elif "Jets-3" in parameterName : process+=jet_kt_cut(20.)
	elif "Jets-4" in parameterName : process+=jet_kt_cut(25.)
	else :
	logging.error("Star not supported for %s " % simulation)
	sys.exit(1)
	# ISR and SppS
	elif ( collider=="ISR" or collider =="SppS" or collider == "SPS" or collider == "Fermilab" ) :
	process = StringBuilder("set /Herwig/Decays/DecayHandler:LifeTimeOption 0\n")
	process+="set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm\n"
	if(collider=="SppS") :
	process = StringBuilder("set /Herwig/Generators/EventGenerator:EventHandler:BeamB /Herwig/Particles/pbar-\n")
	if "17.4" in parameterName : process+=collider_lumi( 17.4)
	elif "27.4" in parameterName : process+=collider_lumi( 27.4)
	elif "30" in parameterName : process+=collider_lumi( 30.4)
	elif "38.8" in parameterName : process+=collider_lumi( 38.8)
	elif "44" in parameterName : process+=collider_lumi( 44.4)
	elif "53" in parameterName : process+=collider_lumi( 53.0)
	elif "62" in parameterName : process+=collider_lumi( 62.2)
	elif "63" in parameterName : process+=collider_lumi( 63.0)
	elif "200" in parameterName : process+=collider_lumi(200.0)
	elif "500" in parameterName : process+=collider_lumi(500.0)
	elif "546" in parameterName : process+=collider_lumi(546.0)
	elif "900" in parameterName : process+=collider_lumi(900.0)
	if "UE" in parameterName :
	if(simulation=="") :
	if "Dipole" in parameters["shower"]:
	process+="read snippets/MB-DipoleShower.in\n"
	else:
	process+="read snippets/MB.in\n"
	process+="read snippets/Diffraction.in\n"
	else :
	logging.error(" SppS and ISR not supported for %s " % simulation)
	sys.exit(1)
	elif "Z-mu" in parameterName :
	if simulation == "" :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	process+=mhat_minm_maxm(2,2,20)
	elif simulation == "Powheg" :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	process+=mhat_minm_maxm(2,2,20)
	elif(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	process+=addLeptonPairCut("2","20")
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("2","20")
	else :
	logging.error(" SppS and ISR not supported for %s " % simulation)
	sys.exit(1)
	else :
	logging.error(" Process not supported for SppS and ISR %s " % parameterName )
	sys.exit(1)

	# LHC
	elif(collider=="LHC") :
	if parameterName.startswith("7-") : process = StringBuilder(collider_lumi(7000.0))
	elif parameterName.startswith("8-") : process = StringBuilder(collider_lumi(8000.0))
	elif parameterName.startswith("13-") : process = StringBuilder(collider_lumi(13000.0))
	elif parameterName.startswith("900") : process = StringBuilder(collider_lumi(900.0))
	elif parameterName.startswith("2360") : process = StringBuilder(collider_lumi(2360.0))
	elif parameterName.startswith("2760") : process = StringBuilder(collider_lumi(2760.0))
	else : process = StringBuilder(collider_lumi(7000.0))
	if(simulation=="") :
	if "VBF" in parameterName :
	process+=insert_ME("MEPP2HiggsVBF")
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	addedBRReweighter = True
	process+="set /Herwig/Particles/tau-:Stable Stable\n"

	elif "ggHJet" in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	addedBRReweighter = True
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+=insert_ME("MEHiggsJet")
	process+=jet_kt_cut(20.)
	elif "ggH" in parameterName :
	process+=insert_ME("MEHiggs")
	process+=insert_ME("MEHiggsJet","qqbar")
	process+=jet_kt_cut(0.0)
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	addedBRReweighter = True
	process+="set /Herwig/Particles/tau-:Stable Stable\n"

	elif "PromptPhoton" in parameterName :
	process+=insert_ME("MEGammaJet")
	if "PromptPhoton-1" in parameterName :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	elif "PromptPhoton-2" in parameterName :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 25.\n"
	elif "PromptPhoton-3" in parameterName :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 80.\n"
	elif "PromptPhoton-4" in parameterName :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 150.\n"
	elif "DiPhoton-GammaGamma" in parameterName :
	process+=insert_ME("MEGammaGamma")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif "DiPhoton-GammaJet" in parameterName :
	process+=insert_ME("MEGammaJet")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	elif "8-WH" in parameterName :
	process+=insert_ME("MEPP2WH")
	process+=jet_kt_cut(0.0)
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-ZH" in parameterName :
	process+=insert_ME("MEPP2ZH")
	process+=jet_kt_cut(0.0)
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "WH" in parameterName :
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;"])
	addedBRReweighter = True
	process+=insert_ME("MEPP2WH")
	process+=jet_kt_cut(0.0)
	elif "ZH" in parameterName :
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;"])
	addedBRReweighter = True
	process+=insert_ME("MEPP2ZH")
	process+=jet_kt_cut(0.0)
	- elif "UE" or "Cent" in parameterName :
	+ elif "UE" in parameterName or "Cent" in parameterName :
	if "Dipole" in parameters["shower"]:
	process+="read snippets/MB-DipoleShower.in\n"
	else:
	process+="set /Herwig/Shower/ShowerHandler:IntrinsicPtGaussian 2.2*GeV\n"
	process+="read snippets/MB.in\n"
	process+="read snippets/Diffraction.in\n"
	if "Long" in parameterName :
	process += "set /Herwig/Decays/DecayHandler:MaxLifeTime 100*mm\n"
	elif "8-DiJets" in parameterName or "7-DiJets" in parameterName or "13-DiJets" in parameterName :
	process+=insert_ME("MEQCD2to2")
	process+="set MEQCD2to2:MaximumFlavour 5\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "13-DiJets" not in parameterName :
	if "-A" in parameterName :
	process+=jet_kt_cut(45.)
	process+="set /Herwig/Cuts/JetKtCut:MinEta -3.\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 3.\n"
	elif "-B" in parameterName :
	process+=jet_kt_cut(20.)
	process+="set /Herwig/Cuts/JetKtCut:MinEta -2.7\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 2.7\n"
	elif "-C" in parameterName :
	process+=jet_kt_cut(20.)
	process+="set /Herwig/Cuts/JetKtCut:MinEta -4.8\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 4.8\n"
	else :
	if "-A" in parameterName :
	process+=jet_kt_cut(60.)
	process+="set /Herwig/Cuts/JetKtCut:MinEta -3.\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 3.\n"
	elif "-B" in parameterName :
	process+=jet_kt_cut(180.)
	process+="set /Herwig/Cuts/JetKtCut:MinEta -3.\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 3.\n"

	if "DiJets-1" in parameterName : process+=mhatmin_cut(90.)
	elif "DiJets-2" in parameterName : process+=mhatmin_cut(200.)
	elif "DiJets-3" in parameterName : process+=mhatmin_cut(450.)
	elif "DiJets-4" in parameterName : process+=mhatmin_cut(750.)
	elif "DiJets-5" in parameterName : process+=mhatmin_cut(950.)
	elif "DiJets-6" in parameterName : process+=mhatmin_cut(1550.)
	elif "DiJets-7" in parameterName : process+=mhatmin_cut(2150.)
	elif "DiJets-8" in parameterName : process+=mhatmin_cut(2750.)
	elif "DiJets-9" in parameterName : process+=mhatmin_cut(3750.)
	elif "DiJets-10" in parameterName : process+=mhatmin_cut(4750.)
	elif "DiJets-11" in parameterName : process+=mhatmin_cut(5750.)
	elif( "7-Jets" in parameterName
	or "8-Jets" in parameterName
	or "13-Jets" in parameterName
	or "2760-Jets" in parameterName
	) :
	process+=insert_ME("MEQCD2to2")
	process+="set MEQCD2to2:MaximumFlavour 5\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "Jets-10" in parameterName : process+=jet_kt_cut(1800.)
	elif "Jets-0" in parameterName : process+=jet_kt_cut(5.)
	elif "Jets-1" in parameterName : process+=jet_kt_cut(10.)
	elif "Jets-2" in parameterName : process+=jet_kt_cut(20.)
	elif "Jets-3" in parameterName : process+=jet_kt_cut(40.)
	elif "Jets-4" in parameterName : process+=jet_kt_cut(70.)
	elif "Jets-5" in parameterName : process+=jet_kt_cut(150.)
	elif "Jets-6" in parameterName : process+=jet_kt_cut(200.)
	elif "Jets-7" in parameterName : process+=jet_kt_cut(300.)
	elif "Jets-8" in parameterName : process+=jet_kt_cut(500.)
	elif "Jets-9" in parameterName : process+=jet_kt_cut(800.)
	elif( "7-Charm" in parameterName or "7-Bottom" in parameterName
	or "8-Bottom" in parameterName) :

	if("8-Bottom" in parameterName) :
	addBRReweighter()
	process+=selectDecayMode("Jpsi",["Jpsi->mu-,mu+;"])

	if "Bottom" in parameterName :
	process+="cp MEHeavyQuark MEBottom\n"
	process+="set MEBottom:QuarkType Bottom\n"
	process+=insert_ME("MEBottom")
	else :
	process+="cp MEHeavyQuark MECharm\n"
	process+="set MECharm:QuarkType Charm\n"
	process+=insert_ME("MECharm")
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "-0" in parameterName :
	if "Bottom" in parameterName :
	process+="set MEBottom:Process Pair\n"
	process+=jet_kt_cut(0.)
	elif "-1" in parameterName : process+=jet_kt_cut(5.)
	elif "-2" in parameterName : process+=jet_kt_cut(15.)
	elif "-3" in parameterName : process+=jet_kt_cut(20.)
	elif "-4" in parameterName : process+=jet_kt_cut(50.)
	elif "-5" in parameterName : process+=jet_kt_cut(80.)
	elif "-6" in parameterName : process+=jet_kt_cut(110.)
	elif "-7" in parameterName : process+=jet_kt_cut(30.)+mhatmin_cut(90.)
	elif "-8" in parameterName : process+=jet_kt_cut(30.)+mhatmin_cut(340.)
	elif "-9" in parameterName : process+=jet_kt_cut(30.)+mhatmin_cut(500.)
	elif "Top-L" in parameterName :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+=insert_ME("MEHeavyQuark")
	process+=selectDecayMode("t",["t->nu_e,e+,b;",
	"t->nu_mu,mu+,b;"])
	process+=addBRReweighter()

	elif "Top-SL" in parameterName :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+=insert_ME("MEHeavyQuark")
	process+="set /Herwig/Particles/t:Synchronized Not_synchronized\n"
	process+="set /Herwig/Particles/tbar:Synchronized Not_synchronized\n"
	process+=selectDecayMode("t",["t->nu_e,e+,b;","t->nu_mu,mu+,b;"])
	process+=selectDecayMode("tbar",["tbar->b,bbar,cbar;",
	"tbar->bbar,cbar,d;",
	"tbar->bbar,cbar,s;",
	"tbar->bbar,s,ubar;",
	"tbar->bbar,ubar,d;"])
	process+=addBRReweighter()

	elif "Top-All" in parameterName :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+=insert_ME("MEHeavyQuark")
	elif "WZ" in parameterName :
	process+=insert_ME("MEPP2VV","WZ")
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;"])
	process+=selectDecayMode("W-",["W-->nu_ebar,e-;",
	"W-->nu_mubar,mu-;"])
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;"])
	addedBRReweighter = True

	elif "WW-emu" in parameterName :
	process+=insert_ME("MEPP2VV","WW")
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+=selectDecayMode("W+",["W+->nu_e,e+;"])
	process+=selectDecayMode("W-",["W-->nu_mubar,mu-;"])
	addedBRReweighter = True

	elif "WW-ll" in parameterName :
	process+=insert_ME("MEPP2VV","WW")
	process+=selectDecayMode("W+",["W+->nu_e,e+;","W+->nu_mu,mu+;","W+->nu_tau,tau+;"])
	addedBRReweighter = True

	elif "ZZ-ll" in parameterName :
	process+=insert_ME("MEPP2VV","ZZ")
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;"])
	addedBRReweighter = True

	elif "ZZ-lv" in parameterName :
	process+=insert_ME("MEPP2VV","ZZ")
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;",
	"Z0->nu_e,nu_ebar;",
	"Z0->nu_mu,nu_mubar;",
	"Z0->nu_tau,nu_taubar;"])
	addedBRReweighter = True

	elif "W-Z-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	process+=insert_ME("MEqq2W2ff","Electron")
	elif "W-Z-mu" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	process+=insert_ME("MEqq2W2ff","Muon")
	elif "W-e" in parameterName :
	process+=insert_ME("MEqq2W2ff","Electron")
	elif "W-mu" in parameterName :
	process+=insert_ME("MEqq2W2ff","Muon")
	elif "Z-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	elif "Z-mu" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-LowMass-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	process+=mhat_minm_maxm(20,20,70)
	elif "Z-MedMass-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	process+=mhat_minm_maxm(40,40,130)
	elif "Z-LowMass-mu" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	process+=mhat_minm_maxm(10,10,70)
	elif "Z-Mass1" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(10,10,35)
	else :
	process+=mhat_minm_maxm(10,10,30)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-Mass2" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(25,25,70)
	else :
	process+=mhat_minm_maxm(30,30,65)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-Mass3" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(60,60,120)
	else :
	process+=mhat_minm_maxm(65,65,115)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-Mass4" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(110,110,8000)
	else :
	process+=mhat_minm_maxm(115,115,8000)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-HighMass1" in parameterName :
	process+=mhat_minm_maxm(116,116,400)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "Z-HighMass2" in parameterName :
	process+=mhat_minm_maxm(400,400,7000)
	if "-e" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	elif "W-Jet" in parameterName :
	process+=insert_ME("MEWJet","Electron","WDecay")
	if "W-Jet-1-e" in parameterName :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 100.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-1-e","W-Jet-e")
	elif "W-Jet-2-e" in parameterName :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 190.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-2-e","W-Jet-e")
	elif "W-Jet-3-e" in parameterName :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 270.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-3-e","W-Jet-e")
	elif "Z-Jet" in parameterName :
	if "-e" in parameterName :
	process+=insert_ME("MEZJet","Electron","ZDecay")
	if "Z-Jet-0-e" in parameterName :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 35.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-e","Z-Jet-e")
	elif "Z-Jet-1-e" in parameterName :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 100.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-1-e","Z-Jet-e")
	elif "Z-Jet-2-e" in parameterName :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 190.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-2-e","Z-Jet-e")
	elif "Z-Jet-3-e" in parameterName :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 270.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-3-e","Z-Jet-e")
	else :
	process+=insert_ME("MEZJet","Muon","ZDecay")
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 35.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-mu","Z-Jet-mu")
	elif "WGamma" in parameterName :
	process+=insert_ME("MEPP2VGamma","1")
	process+="set MEPP2VGamma:MassOption 1"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 10.\n"


	if "-e" in parameterName :
	process+=selectDecayMode("W+",["W+->nu_e,e+;"])
	addedBRReweighter=True
	else :
	process+=selectDecayMode("W+",["W+->nu_mu,mu+;"])
	addedBRReweighter=True
	elif "ZGamma" in parameterName :
	process+=insert_ME("MEPP2VGamma","2")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 10.\n"
	if "-e" in parameterName :
	process+=selectDecayMode("Z0",["Z0->e-,e+;"])
	addedBRReweighter=True
	else :
	process+=selectDecayMode("Z0",["Z0->mu-,mu+;"])
	addedBRReweighter=True
	else :
	logging.error(" Process %s not supported for internal matrix elements" % name)
	sys.exit(1)
	elif(simulation=="Powheg") :
	if "VBF" in parameterName :
	process+=insert_ME("PowhegMEPP2HiggsVBF")
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	addedBRReweighter = True
	process+="set /Herwig/Particles/tau-:Stable Stable\n"

	elif "ggHJet" in parameterName :
	logging.error(" Process %s not supported for POWHEG matrix elements" % name)
	sys.exit(1)
	elif "ggH" in parameterName :
	process+=insert_ME("PowhegMEHiggs")
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	addedBRReweighter = True
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	elif "8-WH" in parameterName :
	process+=insert_ME("PowhegMEPP2WH")
	process+=jet_kt_cut(0.0)
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "8-ZH" in parameterName :
	process+=insert_ME("PowhegMEPP2ZH")
	process+=jet_kt_cut(0.0)
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	addedBRReweighter = True
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	addedBRReweighter = True
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	addedBRReweighter = True
	elif "WH" in parameterName :
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;"])
	addedBRReweighter = True
	process+=insert_ME("PowhegMEPP2WH")
	process+=jet_kt_cut(0.0)
	elif "ZH" in parameterName :
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;"])
	addedBRReweighter = True
	process+=insert_ME("PowhegMEPP2ZH")
	process+=jet_kt_cut(0.0)
	elif "UE" in parameterName :
	logging.error(" Process %s not supported for powheg matrix elements" % name)
	sys.exit(1)
	elif "WZ" in parameterName :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\n"
	process+="set /Herwig/Shower/ShowerHandler:SplitHardProcess No\n";
	process+="set /Herwig/Decays/ZDecayer:PhotonGenerator NULL\n";
	process+="set /Herwig/Decays/WDecayer:PhotonGenerator NULL\n";
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 0 /Herwig/Particles/tau-\n"
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 1 /Herwig/Particles/tau+\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PreCascadeHandlers 0 /Herwig/NewPhysics/DecayHandler\n"
	process+=insert_ME("PowhegMEPP2VV","WZ")
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;"])
	process+=selectDecayMode("W-",["W-->nu_ebar,e-;",
	"W-->nu_mubar,mu-;"])
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;"])
	addedBRReweighter = True

	elif "WW-emu" in parameterName :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\n"
	process+="set /Herwig/Shower/ShowerHandler:SplitHardProcess No\n";
	process+="set /Herwig/Decays/ZDecayer:PhotonGenerator NULL\n";
	process+="set /Herwig/Decays/WDecayer:PhotonGenerator NULL\n";
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 0 /Herwig/Particles/tau-\n"
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 1 /Herwig/Particles/tau+\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PreCascadeHandlers 0 /Herwig/NewPhysics/DecayHandler\n"
	process+=insert_ME("PowhegMEPP2VV","WW")
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+=selectDecayMode("W+",["W+->nu_e,e+;"])
	process+=selectDecayMode("W-",["W-->nu_mubar,mu-;"])
	addedBRReweighter = True

	elif "WW-ll" in parameterName :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\n"
	process+="set /Herwig/Shower/ShowerHandler:SplitHardProcess No\n";
	process+="set /Herwig/Decays/ZDecayer:PhotonGenerator NULL\n";
	process+="set /Herwig/Decays/WDecayer:PhotonGenerator NULL\n";
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 0 /Herwig/Particles/tau-\n"
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 1 /Herwig/Particles/tau+\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PreCascadeHandlers 0 /Herwig/NewPhysics/DecayHandler\n"
	process+=insert_ME("PowhegMEPP2VV","WW")
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;",
	"W+->nu_tau,tau+;"])
	addedBRReweighter = True

	elif "ZZ-ll" in parameterName :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\n"
	process+="set /Herwig/Shower/ShowerHandler:SplitHardProcess No\n";
	process+="set /Herwig/Decays/ZDecayer:PhotonGenerator NULL\n";
	process+="set /Herwig/Decays/WDecayer:PhotonGenerator NULL\n";
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 0 /Herwig/Particles/tau-\n"
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 1 /Herwig/Particles/tau+\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PreCascadeHandlers 0 /Herwig/NewPhysics/DecayHandler\n"
	process+=insert_ME("PowhegMEPP2VV","ZZ")
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;"])
	addedBRReweighter = True

	elif "ZZ-lv" in parameterName :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\n"
	process+="set /Herwig/Shower/ShowerHandler:SplitHardProcess No\n";
	process+="set /Herwig/Decays/ZDecayer:PhotonGenerator NULL\n";
	process+="set /Herwig/Decays/WDecayer:PhotonGenerator NULL\n";
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 0 /Herwig/Particles/tau-\n"
	process+="insert /Herwig/NewPhysics/DecayHandler:Excluded 1 /Herwig/Particles/tau+\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PreCascadeHandlers 0 /Herwig/NewPhysics/DecayHandler\n"
	process+=insert_ME("PowhegMEPP2VV","ZZ")
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;",
	"Z0->nu_e,nu_ebar;",
	"Z0->nu_mu,nu_mubar;",
	"Z0->nu_tau,nu_taubar;"])
	addedBRReweighter = True

	elif "W-Z-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	process+=insert_ME("PowhegMEqq2W2ff","Electron")
	elif "W-Z-mu" in parameterName :
	process+=insert_ME("MEqq2gZ2ff","Muon")
	process+=insert_ME("MEqq2W2ff","Muon")
	elif "W-e" in parameterName :
	process+=insert_ME("PowhegMEqq2W2ff","Electron")
	elif "W-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2W2ff","Muon")
	elif "Z-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	elif "Z-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-LowMass-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	process+=mhat_minm_maxm(20,20,70)
	elif "Z-MedMass-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	process+=mhat_minm_maxm(40,40,130)
	elif "Z-LowMass-mu" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	process+=mhat_minm_maxm(10,10,70)
	elif "Z-Mass1" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(10,10,35)
	else :
	process+=mhat_minm_maxm(10,10,30)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-Mass2" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(25,25,70)
	else :
	process+=mhat_minm_maxm(30,30,65)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-Mass3" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(60,60,120)
	else :
	process+=mhat_minm_maxm(65,65,115)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-Mass4" in parameterName :
	if parameterName.startswith("13-") :
	process+=mhat_minm_maxm(110,110,8000)
	else:
	process+=mhat_minm_maxm(115,115,8000)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-HighMass1" in parameterName :
	process+=mhat_minm_maxm(116,116,400)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "Z-HighMass2" in parameterName :
	process+=mhat_minm_maxm(400,400,7000)
	if "-e" in parameterName :
	process+=insert_ME("PowhegMEqq2gZ2ff","Electron")
	else :
	process+=insert_ME("PowhegMEqq2gZ2ff","Muon")
	elif "DiPhoton-GammaGamma" in parameterName :
	process+=insert_ME("MEGammaGammaPowheg","GammaGamma")
	process+=insert_ME("MEGammaGamma","gg")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+=jet_kt_cut(5.)
	parameterName=parameterName.replace("-GammaGamma","")
	elif "DiPhoton-GammaJet" in parameterName :
	process+=insert_ME("MEGammaGammaPowheg","VJet")
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+=jet_kt_cut(5.)
	parameterName=parameterName.replace("-GammaJet","")
	else :
	logging.error(" Process %s not supported for internal POWHEG matrix elements" % name)
	sys.exit(1)

	elif( simulation=="Matchbox" or simulation=="Merging" ) :
	if "VBF" in parameterName :
	parameters["nlo"] = "read Matchbox/VBFNLO.in\n"
	if(simulation=="Merging"):
	process+="cd /Herwig/Merging/\n"
	process+="insert "+thefactory+":DiagramGenerator:RestrictLines 0 /Herwig/Particles/Z0\n"
	process+="insert "+thefactory+":DiagramGenerator:RestrictLines 0 /Herwig/Particles/W+\n"
	process+="insert "+thefactory+":DiagramGenerator:RestrictLines 0 /Herwig/Particles/W-\n"
	process+="insert "+thefactory+":DiagramGenerator:RestrictLines 0 /Herwig/Particles/gamma\n"
	process+="do "+thefactory+":DiagramGenerator:TimeLikeRange 0 0\n"
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p h0 j j","0","3","FixedScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p h0 j j","0","3","FixedScale",1,1)
	process+=setHardProcessWidthToZero(["h0"])
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	process+=addBRReweighter()
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	process+=addBRReweighter()
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	process+=addBRReweighter()
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	process+=addBRReweighter()
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	elif "ggHJet" in parameterName :
	if(simulation=="Merging"):
	logging.warning("ggHJet not explicitly tested for %s " % simulation)
	sys.exit(0)
	parameters["nlo"] = "read Matchbox/MadGraph-GoSam.in\nread Matchbox/HiggsEffective.in\n"
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	process+=addBRReweighter()
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+=setHardProcessWidthToZero(["h0"])
	process+=addProcess(thefactory,"p p h0 j","3","1","FixedScale",0,0)
	process+=addFirstJet("20")
	process+="set "+thefactory+":ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif "ggH" in parameterName :
	parameters["nlo"] = "read Matchbox/MadGraph-GoSam.in\nread Matchbox/HiggsEffective.in\n"
	if(simulation=="Merging"):
	process+= "cd /Herwig/MatrixElements/Matchbox/Amplitudes\nset OpenLoops:HiggsEff Yes\nset MadGraph:Model heft\n"
	process+="cd /Herwig/Merging/\n"
	process+=setHardProcessWidthToZero(["h0"])
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p h0","2","1","FixedScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p h0","2","1","FixedScale",2,2)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	process+=addBRReweighter()
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	process+=addBRReweighter()
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	process+=addBRReweighter()
	elif "8-" not in parameterName :
	process+=selectDecayMode("h0",["h0->tau-,tau+;"])
	process+=addBRReweighter()
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	elif "8-WH" in parameterName :
	if(simulation=="Merging"):
	logging.warning("8-WH not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["h0","W+","W-"])
	process+=addProcess(thefactory,"p p W+ h0","0","2","FixedScale",0,0)
	process+=addProcess(thefactory,"p p W- h0","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	process+=addBRReweighter()
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	process+=addBRReweighter()
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	process+=addBRReweighter()

	elif "8-ZH" in parameterName :
	if(simulation=="Merging"):
	logging.warning("8-ZH not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["h0","Z0"])
	process+=addProcess(thefactory,"p p Z0 h0","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	if "GammaGamma" in parameterName :
	process+=selectDecayMode("h0",["h0->gamma,gamma;"])
	process+=addBRReweighter()
	elif "WW" in parameterName :
	process+=selectDecayMode("h0",["h0->W+,W-;"])
	process+=addBRReweighter()
	elif "ZZ" in parameterName :
	process+=selectDecayMode("h0",["h0->Z0,Z0;"])
	process+=addBRReweighter()

	elif "WH" in parameterName :
	if(simulation=="Merging"):
	logging.warning("WH not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=addBRReweighter()
	process+=setHardProcessWidthToZero(["h0"])
	process+=addProcess(thefactory,"p p e+ nu h0","0","3","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p e- nu h0","0","3","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu+ nu h0","0","3","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu- nu h0","0","3","LeptonPairMassScale",0,0)
	process+=addLeptonPairCut("60","120")
	elif "ZH" in parameterName :
	if(simulation=="Merging"):
	logging.warning("ZH not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=selectDecayMode("h0",["h0->b,bbar;"])
	process+=addBRReweighter()
	process+=setHardProcessWidthToZero(["h0"])
	process+=addProcess(thefactory,"p p e+ e- h0","0","3","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu+ mu- h0","0","3","LeptonPairMassScale",0,0)
	process+=addLeptonPairCut("60","120")
	elif "UE" in parameterName :
	logging.error(" Process %s not supported for Matchbox matrix elements" % name)
	sys.exit(1)
	elif "8-DiJets" in parameterName or "7-DiJets" in parameterName or "13-DiJets" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p j j","2","0","MaxJetPtScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p j j","2","0","MaxJetPtScale",1,1)
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "13-DiJets" not in parameterName :
	if "-A" in parameterName :
	process+=addFirstJet("45")
	process+=addSecondJet("25")
	process+="set /Herwig/Cuts/FirstJet:YRange -3. 3.\n"
	process+="set /Herwig/Cuts/SecondJet:YRange -3. 3.\n"
	elif "-B" in parameterName :
	process+=addFirstJet("20")
	process+=addSecondJet("15")
	process+="set /Herwig/Cuts/FirstJet:YRange -2.7 2.7\n"
	process+="set /Herwig/Cuts/SecondJet:YRange -2.7 2.7\n"
	elif "-C" in parameterName :
	process+=addFirstJet("20")
	process+=addSecondJet("15")
	process+="set /Herwig/Cuts/FirstJet:YRange -4.8 4.8\n"
	process+="set /Herwig/Cuts/SecondJet:YRange -4.8 4.8\n"
	else :
	logging.error("Exit 00001")
	sys.exit(1)
	else :
	if "-A" in parameterName :
	process+= addFirstJet("75.")
	process+=addSecondJet("60.")
	process+="set /Herwig/Cuts/JetKtCut:MinEta -3.\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 3.\n"
	elif "-B" in parameterName :
	process+= addFirstJet("220.")
	process+=addSecondJet("180.")
	process+="set /Herwig/Cuts/JetKtCut:MinEta -3.\n"
	process+="set /Herwig/Cuts/JetKtCut:MaxEta 3.\n"
	else :
	logging.error("Exit 00001")
	sys.exit(1)


	if "DiJets-1" in parameterName : process+=addJetPairCut("90")
	elif "DiJets-2" in parameterName : process+=addJetPairCut("200")
	elif "DiJets-3" in parameterName : process+=addJetPairCut("450")
	elif "DiJets-4" in parameterName : process+=addJetPairCut("750")
	elif "DiJets-5" in parameterName : process+=addJetPairCut("950")
	elif "DiJets-6" in parameterName : process+=addJetPairCut("1550")
	elif "DiJets-7" in parameterName : process+=addJetPairCut("2150")
	elif "DiJets-8" in parameterName : process+=addJetPairCut("2750")
	elif "DiJets-9" in parameterName : process+=mhatmin_cut(3750.)
	elif "DiJets-10" in parameterName : process+=mhatmin_cut(4750.)
	elif "DiJets-11" in parameterName : process+=mhatmin_cut(5750.)
	else :
	logging.error("Exit 00002")
	sys.exit(1)


	elif( "7-Jets" in parameterName
	or "8-Jets" in parameterName
	or "13-Jets" in parameterName
	or "2760-Jets" in parameterName
	) :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p j j","2","0","MaxJetPtScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p j j","2","0","MaxJetPtScale",1,1)
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "Jets-10" in parameterName : process+=addFirstJet("1800")
	elif "Jets-0" in parameterName : process+=addFirstJet("5")
	elif "Jets-1" in parameterName : process+=addFirstJet("10")
	elif "Jets-2" in parameterName : process+=addFirstJet("20")
	elif "Jets-3" in parameterName : process+=addFirstJet("40")
	elif "Jets-4" in parameterName : process+=addFirstJet("70")
	elif "Jets-5" in parameterName : process+=addFirstJet("150")
	elif "Jets-6" in parameterName : process+=addFirstJet("200")
	elif "Jets-7" in parameterName : process+=addFirstJet("300")
	elif "Jets-8" in parameterName : process+=addFirstJet("500")
	elif "Jets-9" in parameterName : process+=addFirstJet("800")
	else :
	logging.error("Exit 00003")
	sys.exit(1)
	elif( "7-Charm" in parameterName or "7-Bottom" in parameterName
	or "8-Bottom" in parameterName) :
	parameters["bscheme"]=fourFlavour
	process+="set /Herwig/Particles/b:HardProcessMass 4.2*GeV\n"
	process+="set /Herwig/Particles/bbar:HardProcessMass 4.2*GeV\n"

	if("8-Bottom" in parameterName) :
	addBRReweighter()
	process+=selectDecayMode("Jpsi",["Jpsi->mu-,mu+;"])

	if "Bottom" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p b bbar","2","0","MaxJetPtScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p b bbar","2","0","MaxJetPtScale",1,0)
	else:
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p c cbar","2","0","MaxJetPtScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p c cbar","2","0","MaxJetPtScale",1,0)

	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if "-0" in parameterName : process+=addFirstJet("0")
	elif "-1" in parameterName : process+=addFirstJet("5")
	elif "-2" in parameterName : process+=addFirstJet("15")
	elif "-3" in parameterName : process+=addFirstJet("20")
	elif "-4" in parameterName : process+=addFirstJet("50")
	elif "-5" in parameterName : process+=addFirstJet("80")
	elif "-6" in parameterName : process+=addFirstJet("110")
	elif "-7" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("90")
	elif "-8" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("340")
	elif "-9" in parameterName :
	process+=addFirstJet("30")
	process+=addSecondJet("25")
	process+=addJetPairCut("500")
	else :
	logging.error("Exit 00004")
	sys.exit(1)

	elif "Top-L" in parameterName :
	process+=setHardProcessWidthToZero(["t","tbar"])
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",2,2)
	process+=selectDecayMode("t",["t->nu_e,e+,b;",
	"t->nu_mu,mu+,b;"])
	process+=addBRReweighter()

	elif "Top-SL" in parameterName :
	process+=setHardProcessWidthToZero(["t","tbar"])
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",2,2)
	process+="set /Herwig/Particles/t:Synchronized Not_synchronized\n"
	process+="set /Herwig/Particles/tbar:Synchronized Not_synchronized\n"
	process+=selectDecayMode("t",["t->nu_e,e+,b;",
	"t->nu_mu,mu+,b;"])
	process+=selectDecayMode("tbar",["tbar->b,bbar,cbar;",
	"tbar->bbar,cbar,d;",
	"tbar->bbar,cbar,s;",
	"tbar->bbar,s,ubar;",
	"tbar->bbar,ubar,d;"])
	process+=addBRReweighter()

	elif "Top-All" in parameterName :
	process+=setHardProcessWidthToZero(["t","tbar"])
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p t tbar","2","0","TopPairMTScale",2,2)
	elif "WZ" in parameterName :
	if(simulation=="Merging"):
	logging.warning("WZ not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["W+","W-","Z0"])
	process+=addProcess(thefactory,"p p W+ Z0","0","2","FixedScale",0,0)
	process+=addProcess(thefactory,"p p W- Z0","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 171.6*GeV\n\n"
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;"])
	process+=selectDecayMode("W-",["W-->nu_ebar,e-;",
	"W-->nu_mubar,mu-;"])
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;"])
	process+=addBRReweighter()
	process+=addLeptonPairCut("60","120")
	elif "WW-emu" in parameterName :
	if(simulation=="Merging"):
	logging.warning("WW-emu not explicitly tested for %s " % simulation)
	sys.exit(0)

	process+=setHardProcessWidthToZero(["W+","W-","Z0"])
	process+=addProcess(thefactory,"p p W+ W-","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 160.8*GeV\n"
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+=selectDecayMode("W+",["W+->nu_e,e+;"])
	process+=selectDecayMode("W-",["W-->nu_mubar,mu-;"])
	process+=addBRReweighter()
	parameters["bscheme"] = "read Matchbox/FourFlavourScheme.in\n"

	process+=addLeptonPairCut("60","120")
	elif "WW-ll" in parameterName :
	if(simulation=="Merging"):
	logging.warning("WW-ll not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["W+","W-","Z0"])
	process+=addProcess(thefactory,"p p W+ W-","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 160.8*GeV\n"
	process+=selectDecayMode("W+",["W+->nu_e,e+;",
	"W+->nu_mu,mu+;",
	"W+->nu_tau,tau+;"])
	process+=addBRReweighter()
	process+=addLeptonPairCut("60","120")
	parameters["bscheme"] = "read Matchbox/FourFlavourScheme.in\n"

	elif "ZZ-ll" in parameterName :
	if(simulation=="Merging"):
	logging.warning("ZZ-ll not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["W+","W-","Z0"])
	process+=addProcess(thefactory,"p p Z0 Z0","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 182.2*GeV\n"
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;"])
	process+=addBRReweighter()
	process+=addLeptonPairCut("60","120")
	elif "ZZ-lv" in parameterName :
	if(simulation=="Merging"):
	logging.warning("ZZ-lv not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=setHardProcessWidthToZero(["W+","W-","Z0"])
	process+=addProcess(thefactory,"p p Z0 Z0","0","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 182.2*GeV\n"
	process+=selectDecayMode("Z0",["Z0->e-,e+;",
	"Z0->mu-,mu+;",
	"Z0->tau-,tau+;",
	"Z0->nu_e,nu_ebar;",
	"Z0->nu_mu,nu_mubar;",
	"Z0->nu_tau,nu_taubar;"])
	process+=addBRReweighter()
	process+=addLeptonPairCut("60","120")
	elif "W-Z-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p e+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p e- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'epm','e+','e-')
	process+=particlegroup(thefactory,'epmnu','e+','e-','nu_e','nu_ebar')
	process+=addProcess(thefactory,"p p epm epmnu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "W-Z-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'mupm','mu+','mu-')
	process+=particlegroup(thefactory,'mupmnu','mu+','mu-','nu_mu','nu_mubar')
	process+=addProcess(thefactory,"p p mupm mupmnu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "W-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p e- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'epm','e+','e-')
	process+=addProcess(thefactory,"p p epm nu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")

	elif "W-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ nu","0","2","LeptonPairMassScale",0,0)
	process+=addProcess(thefactory,"p p mu- nu","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=particlegroup(thefactory,'mupm','mu+','mu-')
	process+=addProcess(thefactory,"p p mupm nu","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")

	elif "Z-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "Z-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("60","120")
	elif "Z-jj" in parameterName :
	if(simulation=="Merging"):
	logging.warning("Z-jj not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=addProcess(thefactory,"p p e+ e- j j","2","2","LeptonPairMassScale",0,0)
	process+=addFirstJet("40")
	process+=addSecondJet("30")
	process+=addLeptonPairCut("60","120")
	elif "Z-LowMass-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("20","70")
	elif "Z-MedMass-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("40","130")
	elif "Z-LowMass-mu" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	process+=addLeptonPairCut("10","70")
	elif "Z-Mass1" in parameterName :
	if parameterName.startswith("13-") :
	process+=addLeptonPairCut("10","35")
	else :
	process+=addLeptonPairCut("10","30")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "Z-Mass2" in parameterName :
	if parameterName.startswith("13-") :
	process+=addLeptonPairCut("25","70")
	else :
	process+=addLeptonPairCut("30","75")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "Z-Mass3" in parameterName :
	if parameterName.startswith("13-") :
	process+=addLeptonPairCut("60","125")
	else :
	process+=addLeptonPairCut("65","120")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "Z-Mass4" in parameterName :
	if parameterName.startswith("13-") :
	process+=addLeptonPairCut("115","8000")
	else :
	process+=addLeptonPairCut("120","8000")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "Z-HighMass1" in parameterName :
	process+=addLeptonPairCut("116","400")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "Z-HighMass2" in parameterName :
	process+=addLeptonPairCut("400","7000")
	if "-e" in parameterName :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p e+ e-","0","2","LeptonPairMassScale",2,2)
	else :
	if(simulation=="Matchbox"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",0,0)
	elif(simulation=="Merging"):
	process+=addProcess(thefactory,"p p mu+ mu-","0","2","LeptonPairMassScale",2,2)
	elif "W-Jet" in parameterName :
	if(simulation=="Merging"):
	logging.warning("W-Jet not explicitly tested for %s " % simulation)
	sys.exit(0)

	process+=addProcess(thefactory,"p p e+ nu j","1","2","HTScale",0,0)
	process+=addProcess(thefactory,"p p e- nu j","1","2","HTScale",0,0)

	process+=addLeptonPairCut("60","120")
	if "W-Jet-1-e" in parameterName :
	process+=addFirstJet("100")
	parameterName=parameterName.replace("W-Jet-1-e","W-Jet-e")
	elif "W-Jet-2-e" in parameterName :
	process+=addFirstJet("190")
	parameterName=parameterName.replace("W-Jet-2-e","W-Jet-e")
	elif "W-Jet-3-e" in parameterName :
	process+=addFirstJet("270")
	parameterName=parameterName.replace("W-Jet-3-e","W-Jet-e")
	else :
	logging.error("Exit 00005")
	sys.exit(1)
	elif "Z-Jet" in parameterName :
	if(simulation=="Merging"):
	logging.warning("Z-Jet not explicitly tested for %s " % simulation)
	sys.exit(0)


	if "-e" in parameterName :
	process+=addProcess(thefactory,"p p e+ e- j","1","2","HTScale",0,0)
	if "Z-Jet-0-e" in parameterName :
	process+=addFirstJet("35")
	parameterName=parameterName.replace("Z-Jet-0-e","Z-Jet-e")
	elif "Z-Jet-1-e" in parameterName :
	process+=addFirstJet("100")
	parameterName=parameterName.replace("Z-Jet-1-e","Z-Jet-e")
	elif "Z-Jet-2-e" in parameterName :
	process+=addFirstJet("190")
	parameterName=parameterName.replace("Z-Jet-2-e","Z-Jet-e")
	elif "Z-Jet-3-e" in parameterName :
	process+=addFirstJet("270")
	parameterName=parameterName.replace("Z-Jet-3-e","Z-Jet-e")
	else :
	logging.error("Exit 00006")
	sys.exit(1)
	else :
	process+=addProcess(thefactory,"p p mu+ mu- j","1","2","HTScale",0,0)
	process+=addFirstJet("35")
	parameterName=parameterName.replace("Z-Jet-0-mu","Z-Jet-mu")
	process+=addLeptonPairCut("60","120")
	elif "Z-bb" in parameterName :
	if(simulation=="Merging"):
	logging.warning("Z-bb not explicitly tested for %s " % simulation)
	sys.exit(0)
	parameters["bscheme"]=fourFlavour
	process+="set /Herwig/Particles/b:HardProcessMass 4.2GeV\nset /Herwig/Particles/bbar:HardProcessMass 4.2GeV\n"
	process+=addProcess(thefactory,"p p e+ e- b bbar","2","2","FixedScale",0,0)
	process+=addLeptonPairCut("66","116")
	process+=addFirstJet("18")
	process+=addSecondJet("15")
	process+=addLeptonPairCut("60","120")
	elif "Z-b" in parameterName :
	if(simulation=="Merging"):
	logging.warning("Z-b not explicitly tested for %s " % simulation)
	sys.exit(0)
	process+=particlegroup(thefactory,'bjet','b','bbar')
	process+=addProcess(thefactory,"p p e+ e- bjet","1","2","FixedScale",0,0)
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 91.2*GeV\n"
	process+=addLeptonPairCut("60","120")
	process+=addFirstJet("15")
	elif "W-b" in parameterName :
	if(simulation=="Merging"):
	logging.warning("W-b not explicitly tested for %s " % simulation)
	sys.exit(0)
	parameters["bscheme"]=fourFlavour
	process += "set /Herwig/Particles/b:HardProcessMass 4.2GeV\nset /Herwig/Particles/bbar:HardProcessMass 4.2GeV\n"
	process+=addProcess(thefactory,"p p e- nu b bbar","2","2","FixedScale",0,0)
	process+=addProcess(thefactory,"p p mu+ nu b bbar","2","2","FixedScale",0,0)
	process += "set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 80.4*GeV\n"
	process+=addFirstJet("30")
	process+=addLeptonPairCut("60","120")
	else :
	logging.error(" Process %s not supported for Matchbox matrix elements" % name)
	sys.exit(1)
	# LHC-GammaGamma
	elif(collider=="LHC-GammaGamma" ) :
	if "-7-" in parameterName : process = StringBuilder(collider_lumi(7000.0))
	elif "-8-" in parameterName : process = StringBuilder(collider_lumi(8000.0))
	else : process = StringBuilder(collider_lumi(7000.0))
	if(simulation=="") :
	if "7" in parameterName : process += insert_ME("MEgg2ff","Muon")
	else :
	logging.error(" Process %s not supported for default matrix elements" % name)
	sys.exit(1)
	else :
	logging.error("LHC-GammaGamma not supported for %s " % simulation)
	sys.exit(1)

	if "EHS" in name :
	parameters['parameterFile'] = os.path.join(collider,"{c}-{pn}.in".format(c="EHS", pn=parameterName))
	else :
	parameters['parameterFile'] = os.path.join(collider,"{c}-{pn}.in".format(c=collider, pn=parameterName))
	parameters['runname'] = 'Rivet-%s' % name
	parameters['process'] = str(process)
	if have_hadronic_collider :
	if "EHS" in name :
	parameters['collider'] = "PPCollider.in\nread snippets/FixedTarget-PP.in"
	else :
	parameters['collider'] = "PPCollider.in"

	#check if selecteddecaymode and addedBRReweighter is consistent

	if selecteddecaymode and not addedBRReweighter:
	logging.error("Decaymode was selected but no BRReweighter was added.")
	sys.exit(1)

	if addedBRReweighter and not selecteddecaymode:
	logging.error("BRReweighter was added but no Decaymode was selected.")
	sys.exit(1)

	# check that we only add one process if in merging mode:

	if numberOfAddedProcesses > 1 and simulation =="Merging":
	logging.error("In Merging only one process is allowed at the moment. See ticket #403.")
	sys.exit(1)

	# Check if a process was added for Merging or Matchbox:

	if numberOfAddedProcesses == 0 and (simulation =="Merging" or simulation =="Matchbox"):
	logging.error("No process was selected.")
	sys.exit(1)

	# get template and write the file

	with open(os.path.join("Rivet/Templates",templateName), 'r') as f:
	templateText = f.read()

	template = Template( templateText )

	with open(os.path.join("Rivet",name+".in"), 'w') as f:
	f.write( template.substitute(parameters) )
	diff --git a/Tests/python/merge-LHC-EW b/Tests/python/merge-LHC-EW
	--- a/Tests/python/merge-LHC-EW
	+++ b/Tests/python/merge-LHC-EW
	@@ -1,469 +1,469 @@
	#! /usr/bin/env python
	from __future__ import print_function
	import logging, sys, os, yoda

	"""%prog

	Script for merging aida files

	"""

	def fillAbove(scale,desthisto, sourcehistosbyptmin) :
	pthigh= 1e100
	ptlow =-1e100
	for pt, h in sorted(sourcehistosbyptmin.items(),reverse=True):
	ptlow=pt
	if(type(desthisto)==yoda.core.Scatter2D) :
	for i in range(0,h.numPoints()) :
	xMin = h.points()[i].x()-h.points()[i].xErrs().minus
	if( xMin*scale >= ptlow and
	xMin*scale < pthigh ) :
	desthisto.addPoint(h.points()[i])
	elif(type(desthisto)==yoda.core.Profile1D) :
	for i in range(0,h.numBins()) :
	if(h.bins()[i].xMin()*scale >= ptlow and
	h.bins()[i].xMin()*scale < pthigh ) :
	desthisto.bins()[i] += h.bins()[i]
	elif(type(desthisto)==yoda.core.Histo1D) :
	for i in range(0,h.numBins()) :
	if(h.bins()[i].xMin()*scale >= ptlow and
	h.bins()[i].xMin()*scale < pthigh ) :
	desthisto.bins()[i] += h.bins()[i]
	else :
	logging.error("Can't merge %s, unknown type" % desthisto.path())
	sys.exit(1)
	pthigh=pt

	def mergeByMass(hpath, sqrts, scale=1.):
	global inhistos_mass
	global outhistos
	try:
	fillAbove(scale,outhistos[hpath], inhistos_mass[hpath][float(sqrts)])
	except:
	pass

	def useOneMass(hpath, sqrts, ptmin):
	global inhistos_mass
	global outhistos
	try:
	## Find best pT_min match
	ptmins = inhistos_mass[hpath][float(sqrts)].keys()
	closest_ptmin = None
	for ptm in ptmins:
	if closest_ptmin is None or \
	abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
	closest_ptmin = ptm
	if closest_ptmin != float(ptmin):
	logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
	"using pTmin=%e instead" % closest_ptmin)
	outhistos[hpath] = inhistos_mass[hpath][float(sqrts)][closest_ptmin]
	except:
	pass

	import sys
	if sys.version_info[:3] < (2,4,0):
	print ("rivet scripts require Python version >= 2.4.0... exiting")
	sys.exit(1)

	if __name__ == "__main__":
	import logging
	from optparse import OptionParser, OptionGroup
	parser = OptionParser(usage="%prog base")
	verbgroup = OptionGroup(parser, "Verbosity control")
	verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
	default=logging.INFO, help="print debug (very verbose) messages")
	verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
	default=logging.INFO, help="be very quiet")
	parser.add_option_group(verbgroup)
	(opts, args) = parser.parse_args()
	logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")

	## Check args
	if len(args) < 1:
	logging.error("Must specify at least the name of the files")
	sys.exit(1)

	yodafiles=["-13-Z-e","-13-Z-mu","-Z-HighMass1-e","-Z-HighMass2-e",
	"-8-Z-Mass1-e" ,"-8-Z-Mass1-mu" ,"-8-Z-Mass2-e" ,"-8-Z-Mass2-mu" ,"-8-Z-Mass3-e" ,"-8-Z-Mass3-mu" ,"-8-Z-Mass4-e" ,"-8-Z-Mass4-mu" ,
	"-13-Z-Mass1-e","-13-Z-Mass1-mu","-13-Z-Mass2-e","-13-Z-Mass2-mu","-13-Z-Mass3-e","-13-Z-Mass3-mu","-13-Z-Mass4-e","-13-Z-Mass4-mu",
	"-W-e","-W-mu","-Z-e","-Z-mu","-Z-mu-Short","-Z-LowMass-e","-Z-LowMass-mu",
	- "-Z-MedMass-e",
	+ "-Z-MedMass-e","8-W-e","8-W-mu",
	"-WW-emu","-WW-ll","-WZ","-ZZ-ll","-ZZ-lv","-8-WZ","-13-WZ","-8-ZZ-lv","-8-WW-ll",
	"-7-W-Jet-1-e","-7-W-Jet-2-e","-7-W-Jet-3-e","-7-Z-Jet-1-e","-7-Z-Jet-2-e","-7-Z-Jet-3-e",
	"-7-WGamma-e","-7-WGamma-mu","-7-ZGamma-e","-7-ZGamma-mu"]
	## Get histos
	outhistos={}

	inhistos_mass = {}
	for f in yodafiles:
	file='Rivet-'+args[0]+f+".yoda"
	if not os.access(file, os.R_OK):
	logging.error("%s can not be read" % file)
	continue
	try:
	aos = yoda.read(file)
	except:
	logging.error("%s can not be parsed as yoda" % file)
	continue
	mass=66
	if(file.find("HighMass1")>=0) :
	mass = 116
	elif(file.find("HighMass2")>=0) :
	mass = 400
	elif(file.find("Mass1")>=0) :
	mass = 12
	elif(file.find("Mass2")>=0) :
	mass = 30
	elif(file.find("Mass3")>=0) :
	mass = 66
	elif(file.find("Mass4")>=0) :
	mass = 116
	## Get histos from this YODA file
	for aopath, ao in aos.items() :
	if("RAW" in aopath) :continue
	if(aopath.find("ATLAS_2010_S8919674")>0) :
	if((aopath.find("d01")>0 or aopath.find("d05")>0 or
	aopath.find("d07")>0) and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif((aopath.find("d02")>0 or aopath.find("d06")>0 or
	aopath.find("d08")>0) and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2011_S9131140")>0) :
	if(aopath.find("d01")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("d02")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2011_I925932")>0) :
	if(aopath.find("d01")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("d02")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2011_I945498")>0) :
	if(aopath.find("y01")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("y02")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("y03")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2013_I1217867")>0) :
	if(aopath.find("y01")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("y02")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2012_I941555")>0) :
	if((aopath.find("y01")>0 or aopath.find("y03")>0 ) and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("y02")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2014_I1300647" )>0) :
	if(aopath.find("y01")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif((not aopath.find("y01")>0) and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2014_I1288706" )>0) :
	if(aopath.find("y02")>0 and file.find("-e")>0) :
	outhistos[aopath] = ao
	elif(aopath.find("y01")>0 and file.find("-mu")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2012_I1204784" )>0) :
	if( file.find("-e")>0 and
	( aopath.find("d03")>0 or
	((aopath.find("d01")>0 or aopath.find("d02")>0) and aopath.find("y01")>0))) :
	outhistos[aopath] = ao
	elif(file.find("-mu")>0 and
	( aopath.find("d04")>0 or
	((aopath.find("d01")>0 or aopath.find("d02")>0) and aopath.find("y02")>0))) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2014_I1312627")>0) :
	if("EL" in aopath) :
	if(file.find("W")>0 and aopath.find("x02")>0) :
	outhistos[aopath] = ao
	elif(file.find("Z")>0 and aopath.find("x03")>0) :
	outhistos[aopath] = ao
	elif ("MU" in aopath) :
	if(file.find("W")>0 and aopath.find("x02")>0) :
	outhistos[aopath] = ao
	elif(file.find("Z")>0 and aopath.find("x03")>0) :
	outhistos[aopath] = ao
	else :
	if(file.find("W")>0 and aopath.find("x02")>0) :
	outhistos[aopath] = ao
	elif(file.find("Z")>0 and aopath.find("x03")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2018_I1711625")>0) :
	if( file.find("-e")>0 and aopath.find("d06-x01-y01") > 0 ):
	if aopath in outhistos :
	outhistos[aopath] += ao
	else :
	outhistos[aopath] = ao
	elif( file.find("-mu")>0 and (aopath.find("d03-x01-y01") > 0 or aopath.find("d05-x01-y01") > 0) ):
	if aopath in outhistos :
	outhistos[aopath] += ao
	else :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2013_I1224539_WJET" )>0) :
	if(file.find("-1-e")>0 and (aopath.find("d52")>0 or aopath.find("d53")>0 or aopath.find("d56")>0 or aopath.find("d57")>0 or aopath.find("d60")>0 or aopath.find("d61")>0 or aopath.find("d64")>0 or aopath.find("d65")>0 or aopath.find("d68")>0 or aopath.find("d69")>0 or aopath.find("d72")>0)) :
	outhistos[aopath] = ao
	elif(file.find("-2-e")>0 and (aopath.find("d54")>0 or aopath.find("d58")>0 or aopath.find("d62")>0 or aopath.find("d66")>0 or aopath.find("d70")>0 or aopath.find("d73")>0)) :
	outhistos[aopath] = ao
	elif(file.find("-3-e")>0 and (aopath.find("d55")>0 or aopath.find("d59")>0 or aopath.find("d63")>0 or aopath.find("d67")>0 or aopath.find("d71")>0 or aopath.find("d74")>0)) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2013_I1224539_ZJET" )>0) :
	if(file.find("-1-e")>0 and (aopath.find("d29")>0 or aopath.find("d30")>0 or aopath.find("d33")>0 or aopath.find("d34")>0 or aopath.find("d37")>0 or aopath.find("d38")>0 or aopath.find("d41")>0 or aopath.find("d42")>0 or aopath.find("d45")>0 or aopath.find("d46")>0 or aopath.find("d49")>0)) :
	outhistos[aopath] = ao
	elif(file.find("-2-e")>0 and (aopath.find("d31")>0 or aopath.find("d35")>0 or aopath.find("d39")>0 or aopath.find("d43")>0 or aopath.find("d47")>0 or aopath.find("d50")>0)) :
	outhistos[aopath] = ao
	elif(file.find("-3-e")>0 and (aopath.find("d32")>0 or aopath.find("d36")>0 or aopath.find("d40")>0 or aopath.find("d44")>0 or aopath.find("d48")>0 or aopath.find("d51")>0)) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2013_I1258128")>0) :
	if(aopath.find("d01")>0 or aopath.find("d02")>0 or
	aopath.find("d03")>0 or aopath.find("d04")>0) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2013_I1209721" )>0 and file.find("-0")>0 ) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2011_I928289")>0) :
	if(file.find("-e")>=0 and (aopath.find("y01")>=0 or aopath.find("y02")>=0)) :
	outhistos[aopath] = ao
	elif(file.find("-mu")>=0 and (aopath.find("y03")>=0 or aopath.find("y04")>=0)) :
	outhistos[aopath] = ao
	elif (aopath.find("CMS_2013_I1122847")>0) :
	if(file.find("-mu")>=0 and aopath.find("d01")>=0 ) :
	outhistos[aopath] = ao
	elif(file.find("-e")>=0 and (aopath.find("d02")>=0 or aopath.find("d03")>=0)) :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2016_I1424838")>0) :
	ihist = int(aopath.split("/")[2].split("-")[0].replace("d",""))
	if(file.find("-mu")>=0 and ihist>30 ) :
	outhistos[aopath] = ao
	elif(file.find("-e")>=0 and ihist<=30) :
	outhistos[aopath] = ao

	elif "ATLAS_2016_I1502620" in aopath :
	newPath = aopath.split(":")[0]+"/"+aopath.split("/")[-1]
	if newPath in outhistos :
	continue
	else :
	outhistos[newPath] = ao
	ao.setPath(newPath)
	elif ("CMS_2015_I1310737" in aopath or "ATLAS_2013_I1219109" in aopath or
	("ATLAS_2014_I1319490" in aopath and "y01" in aopath )) :
	if aopath in outhistos :
	outhistos[aopath] += ao
	else :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2015_I1351916")>=0) :
	if(aopath.find("-y02")>=0) :
	pathBase = "/%s/" % aopath.split("/")[1]
	if("d03" in aopath) :
	hp = aos[pathBase+"_ncf_pos"]
	hm = aos[pathBase+"_ncf_neg"]
	else :
	hp = aos[pathBase+"_ncc_pos"]
	hm = aos[pathBase+"_ncc_neg"]
	ratio = (hp-hm)/(hp+hm)
	title=""
	if hasattr(ao, 'title'):
	title=ao.title()
	hnew = yoda.core.Scatter2D(aopath,title)
	hnew.combineWith(ratio)
	outhistos[aopath] = hnew
	else :
	continue
	elif (aopath.find("ATLAS_2014_I1282447")>=0) :
	if((aopath.find("/ATLAS_2014_I1282447/d02-x01-y01")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d08-x01-y01")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d02-x01-y02")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d02-x01-y01")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d05-x01-y02")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d05-x01-y03")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d06-x01-y01")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d06-x01-y02")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d06-x01-y03")>=0 or
	aopath.find("/ATLAS_2014_I1282447/d06-x01-y04")>=0) and
	not (aopath.find("plus")>=0 or aopath.find("minus")>=0 or
	aopath.find("inc")>=0)) :
	continue
	if aopath in outhistos :
	outhistos[aopath] += ao
	else :
	outhistos[aopath] = ao
	elif (aopath.find("ATLAS_2015_I1408516")>=0) :
	if aopath not in inhistos_mass:
	inhistos_mass[aopath] = {}
	tmpE = inhistos_mass[aopath]
	sqrts=8000
	if sqrts not in tmpE:
	tmpE[sqrts] = {}
	tmpP = tmpE[sqrts]
	if mass not in tmpP:
	tmpP[mass] = ao
	else:
	raise Exception("A set with mass = %s already exists" % ( mass))
	elif (aopath.find("ATLAS_2013_I1234228")>=0) :
	if aopath not in inhistos_mass:
	inhistos_mass[aopath] = {}
	tmpE = inhistos_mass[aopath]
	sqrts=7000
	if sqrts not in tmpE:
	tmpE[sqrts] = {}
	tmpP = tmpE[sqrts]
	if mass not in tmpP:
	tmpP[mass] = ao
	else:
	raise Exception("A set with mass = %s already exists" % ( mass))
	elif (aopath.find("ATLAS_2017_I1514251")>=0) :
	outhistos[aopath] = ao
	if(type(ao)==yoda.core.Scatter2D) :
	ao.scaleY(2.)
	elif(type(ao)==yoda.core.Histo1D) :
	ao.scaleW(2.)
	elif ( "2015_041_" in aopath) :
	aopath=aopath.replace("_EL","").replace("_MU","")
	ao.setPath(aopath)
	outhistos[aopath] = ao
	else :
	outhistos[aopath] = ao
	for hpath,hsets in inhistos_mass.items():
	if(hpath!="/ATLAS_2015_I1408516:LMODE=EL/d41-x01-y01" and
	hpath!="/ATLAS_2015_I1408516:LMODE=MU/d41-x01-y02" and
	hpath!="/ATLAS_2013_I1234228/d01-x01-y02" ) :
	continue
	title=""
	path=""
	histo = list(list(hsets.values())[0].values())[0]
	if hasattr(histo, 'title'):
	title=histo.title()
	if hasattr(histo, 'path'):
	path=histo.path()
	if(type(histo)==yoda.core.Scatter2D) :
	outhistos[hpath] = yoda.core.Scatter2D(path,title)
	elif(type(histo)==yoda.core.Profile1D) :
	outhistos[hpath] = yoda.core.Profile1D(path,title)
	for i in range(0,histo.numBins()) :
	outhistos[hpath].addBin(histo.bins()[i].xMin(),
	histo.bins()[i].xMax())
	elif(type(histo)==yoda.core.Histo1D) :
	outhistos[hpath] = yoda.core.Histo1D(path,title)
	for i in range(0,histo.numBins()) :
	outhistos[hpath].addBin(histo.bins()[i].xMin(),
	histo.bins()[i].xMax())
	else :
	logging.error("Histogram %s is of unknown type" % hpath)
	sys.exit(1)
	# sort out mass bins for ATLAS Z-> e,mu at 8 TeV
	for ltype in ["EL","MU"] :
	if(ltype=="EL") :
	y = "y01"
	mergeByMass("/ATLAS_2015_I1408516:LMODE=EL/d41-x01-y01", "8000")
	else :
	y = "y04"
	mergeByMass("/ATLAS_2015_I1408516:LMODE=MU/d41-x01-y02", "8000")
	for d in [2,3,4,14,26,38]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "30" )
	for d in [5,6,7,8,9,10,15,17,18,19,20,21,22,27,29,30,31,32,33,34,39]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "66" )
	for d in [11,12,13,16,28,40]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "116" )
	for d in [23,35]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "12" )
	for d in [24,36]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "12" )
	for d in [25,37]:
	useOneMass("/ATLAS_2015_I1408516:LMODE=%s/d%02d-x01-%s" % (ltype,d,y), "8000", "30" )
	# sort out ratios for ATLAS W+c
	if("/ATLAS_2014_I1282447/d02-x01-y01_plus" in outhistos and
	"/ATLAS_2014_I1282447/d02-x01-y01_minus" in outhistos) :
	d02y01_plus = outhistos["/ATLAS_2014_I1282447/d02-x01-y01_plus"]
	d02y01_minus = outhistos["/ATLAS_2014_I1282447/d02-x01-y01_minus"]
	ratio_d02y01 = d02y01_plus/d02y01_minus
	ratio_d02y01.setPath("/ATLAS_2014_I1282447/d02-x01-y01")
	del outhistos["/ATLAS_2014_I1282447/d02-x01-y01_plus"]
	del outhistos["/ATLAS_2014_I1282447/d02-x01-y01_minus"]
	outhistos["/ATLAS_2014_I1282447/d02-x01-y01"] = ratio_d02y01
	if("/ATLAS_2014_I1282447/d02-x01-y02_plus" in outhistos and
	"/ATLAS_2014_I1282447/d02-x01-y02_minus" in outhistos) :
	d02y02_plus = outhistos["/ATLAS_2014_I1282447/d02-x01-y02_plus"]
	d02y02_minus = outhistos["/ATLAS_2014_I1282447/d02-x01-y02_minus"]
	ratio_d02y02 = d02y02_plus/d02y02_minus
	ratio_d02y02.setPath("/ATLAS_2014_I1282447/d02-x01-y02")
	del outhistos["/ATLAS_2014_I1282447/d02-x01-y02_plus"]
	del outhistos["/ATLAS_2014_I1282447/d02-x01-y02_minus"]
	outhistos["/ATLAS_2014_I1282447/d02-x01-y02"] = ratio_d02y02
	if("/ATLAS_2014_I1282447/d08-x01-y01_plus" in outhistos and
	"/ATLAS_2014_I1282447/d08-x01-y01_minus" in outhistos) :
	d08y01_plus = outhistos["/ATLAS_2014_I1282447/d08-x01-y01_plus"]
	d08y01_minus = outhistos["/ATLAS_2014_I1282447/d08-x01-y01_minus"]
	ratio_d08y01 = d08y01_plus/d08y01_minus
	ratio_d08y01.setPath("/ATLAS_2014_I1282447/d08-x01-y01")
	del outhistos["/ATLAS_2014_I1282447/d08-x01-y01_plus"]
	del outhistos["/ATLAS_2014_I1282447/d08-x01-y01_minus"]
	outhistos["/ATLAS_2014_I1282447/d08-x01-y01"] = ratio_d08y01
	if ("/ATLAS_2014_I1282447/d05-x01-y01" in outhistos and
	"/ATLAS_2014_I1282447/d01-x01-y02" in outhistos) :
	h_winc = outhistos["/ATLAS_2014_I1282447/d05-x01-y01"]
	h_d = outhistos["/ATLAS_2014_I1282447/d01-x01-y02"]
	ratio_wd = h_d/h_winc
	ratio_wd.setPath("/ATLAS_2014_I1282447/d05-x01-y02")
	outhistos["/ATLAS_2014_I1282447/d05-x01-y02"] = ratio_wd
	if ("/ATLAS_2014_I1282447/d05-x01-y01" in outhistos and
	"/ATLAS_2014_I1282447/d01-x01-y03" in outhistos) :
	h_winc = outhistos["/ATLAS_2014_I1282447/d05-x01-y01"]
	h_dstar= outhistos["/ATLAS_2014_I1282447/d01-x01-y03"]
	ratio_wdstar = h_dstar/h_winc
	ratio_wdstar.setPath("/ATLAS_2014_I1282447/d05-x01-y03")
	outhistos["/ATLAS_2014_I1282447/d05-x01-y03"] = ratio_wdstar
	if("/ATLAS_2014_I1282447/d06-x01-y01_winc" in outhistos and
	"/ATLAS_2014_I1282447/d06-x01-y02_winc" in outhistos) :
	h_winc_plus = outhistos["/ATLAS_2014_I1282447/d06-x01-y01_winc"]
	h_winc_minus = outhistos["/ATLAS_2014_I1282447/d06-x01-y02_winc"]
	if( "/ATLAS_2014_I1282447/d06-x01-y01_wplus" in outhistos ) :
	h_wd_plus = outhistos["/ATLAS_2014_I1282447/d06-x01-y01_wplus"]
	ratio_wd_plus = h_wd_plus/h_winc_plus
	ratio_wd_plus.setPath("/ATLAS_2014_I1282447/d06-x01-y01")
	outhistos["/ATLAS_2014_I1282447/d06-x01-y01"] = ratio_wd_plus
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y01_wplus"]
	if( "/ATLAS_2014_I1282447/d06-x01-y02_wminus" in outhistos ) :
	h_wd_minus = outhistos["/ATLAS_2014_I1282447/d06-x01-y02_wminus"]
	ratio_wd_minus = h_wd_minus/h_winc_minus
	ratio_wd_minus.setPath("/ATLAS_2014_I1282447/d06-x01-y02")
	outhistos["/ATLAS_2014_I1282447/d06-x01-y02"] = ratio_wd_minus
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y02_wminus"]
	if ( "/ATLAS_2014_I1282447/d06-x01-y03_wplus" in outhistos) :
	h_wdstar_plus = outhistos["/ATLAS_2014_I1282447/d06-x01-y03_wplus"]
	ratio_wdstar_plus = h_wdstar_plus/h_winc_plus
	ratio_wdstar_plus.setPath("/ATLAS_2014_I1282447/d06-x01-y03")
	outhistos["/ATLAS_2014_I1282447/d06-x01-y03"] = ratio_wdstar_plus
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y03_wplus"]
	if ( "/ATLAS_2014_I1282447/d06-x01-y04_wminus" in outhistos) :
	h_wdstar_minus = outhistos["/ATLAS_2014_I1282447/d06-x01-y04_wminus"]
	ratio_wdstar_minus = h_wdstar_minus/h_winc_minus
	ratio_wdstar_minus.setPath("/ATLAS_2014_I1282447/d06-x01-y04")
	outhistos["/ATLAS_2014_I1282447/d06-x01-y04"] = ratio_wdstar_minus
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y04_wminus"]
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y01_winc"]
	del outhistos["/ATLAS_2014_I1282447/d06-x01-y02_winc"]


	mergeByMass("/ATLAS_2013_I1234228/d01-x01-y02", "7000")
	# ATLAS W/Z ratios
	iy=0
	for suffix in ["",":LMODE=EL",":LMODE=MU"] :
	iy+=1
	basePath = "/ATLAS_2014_I1312627%s/" % suffix
	for ihist in range(1,21) :
	num = "%sd%02d-x02-y%02d" % (basePath,ihist,iy)
	den = "%sd%02d-x03-y%02d" % (basePath,ihist,iy)
	if(den in outhistos and num in outhistos) :
	ratio = outhistos[num]/outhistos[den]
	ratio.setPath("%sd%02d-x01-y%02d" % (basePath,ihist,iy))
	del outhistos[num]
	del outhistos[den]
	outhistos[ "%sd%02d-x01-y%02d" % (basePath,ihist,iy)]=ratio
	elif den in outhistos :
	del outhistos[den]
	elif num in outhistos :
	del outhistos[num]
	# Choose output file
	name = args[0]+"-EW.yoda"
	yoda.writeYODA(outhistos,name)
	sys.exit(0)
	diff --git a/src/DM.in b/src/DM.in
	--- a/src/DM.in
	+++ b/src/DM.in
	@@ -1,89 +1,96 @@
	# -- ThePEG-repository --


	##################################################
	# Read the model parameters
	##################################################
	read DM.model

	##################################################
	# Set the beams
	##################################################

	read snippets/EECollider.in
	cd /Herwig/EventHandlers
	set EventHandler:BeamA /Herwig/Particles/chi
	set EventHandler:BeamB /Herwig/Particles/chi
	set /Herwig/Particles/chi:PDF /Herwig/Partons/NoPDF

	##################################################
	# Set long-lived hadrons/muon to be unstable
	##################################################
	set /Herwig/Particles/pi0:Stable Unstable
	set /Herwig/Particles/pi+:Stable Unstable
	set /Herwig/Particles/mu-:Stable Unstable
	set /Herwig/Particles/n0:Stable Unstable
	set /Herwig/Particles/K_L0:Stable Unstable
	set /Herwig/Particles/K+:Stable Unstable

	##################################################
	# Selected the hard processes
	##################################################
	cd /Herwig/MatrixElements

	+# using currents
	insert SubProcess:MatrixElements 0 MEDM2Kaons
	insert SubProcess:MatrixElements 0 MEDM2Pions
	insert SubProcess:MatrixElements 0 MEDM2Kaons
	insert SubProcess:MatrixElements 0 MEDM3Pions
	insert SubProcess:MatrixElements 0 MEDM4Pions
	insert SubProcess:MatrixElements 0 MEDM2EtaPiPi
	insert SubProcess:MatrixElements 0 MEDM2EtaPrimePiPi
	insert SubProcess:MatrixElements 0 MEDM2OmegaPiPi
	insert SubProcess:MatrixElements 0 MEDM2OmegaPi
	insert SubProcess:MatrixElements 0 MEDM2PiGamma
	insert SubProcess:MatrixElements 0 MEDM2EtaPhoton
	insert SubProcess:MatrixElements 0 MEDM2EtaPhi
	insert SubProcess:MatrixElements 0 MEDM2EtaOmega
	insert SubProcess:MatrixElements 0 MEDM2ppbar
	insert SubProcess:MatrixElements 0 MEDM2KKPi
	insert SubProcess:MatrixElements 0 MEDM2PhiPi
	-
	+# using standard matrix elements and hadronization
	+# cd /Herwig/NewPhysics
	+# insert HPConstructor:Incoming 0 /Herwig/Particles/chi
	+# insert HPConstructor:Incoming 1 /Herwig/Particles/chi
	+# insert HPConstructor:Outgoing 0 /Herwig/Particles/u
	+# insert HPConstructor:Outgoing 1 /Herwig/Particles/d
	+# insert HPConstructor:Outgoing 2 /Herwig/Particles/s

	##################################################
	# LEP physics parameters (override defaults)
	##################################################
	cd /Herwig/Generators
	# must be slightly greater than 2 x DM mass
	set EventGenerator:EventHandler:LuminosityFunction:Energy 2.0001
	# mass of the DM particle
	set /Herwig/Particles/chi:NominalMass 1.0
	# mass of the mediator
	set /Herwig/Particles/Zp:NominalMass 2.0
	# mediator couplings
	# DM
	set /Herwig/NewPhysics/DM/Model:cDMmed 1.
	# quark-anti-quark
	set /Herwig/NewPhysics/DM/Model:cSMmed 0 1.0
	set /Herwig/NewPhysics/DM/Model:cSMmed 1 1.0
	set /Herwig/NewPhysics/DM/Model:cSMmed 2 1.0
	# cut-off for partonic mediator decays
	set /Herwig/NewPhysics/VectorDecayConstructor:MassCut 2.
	# limit for isr
	set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 0.2
	cd /Herwig/Generators
	##################################################
	## prepare for Rivet analysis or HepMC output
	## when running with parton shower
	##################################################
	# create ThePEG::RivetAnalysis /Herwig/Analysis/Rivet RivetAnalysis.so
	# insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/Rivet
	# insert /Herwig/Analysis/Rivet:Analyses 0 ANALYSIS


	###################################################
	# Save run for later usage with 'Herwig run'
	##################################################
	set EventGenerator:MaxErrors 10000
	set EventGenerator:EventHandler:StatLevel Full
	set EventGenerator:EventHandler:CascadeHandler NULL
	saverun DM EventGenerator

File Metadata

Mime Type: text/x-diff
Expires: Sat, Dec 21, 12:30 PM (1 d, 20 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 4022780
Default Alt Text: (151 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions