No OneTemporary
Actions

Size

129 KB

Subscribers

None

View Options

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

	EXTRA_DIST = Inputs python Rivet

	dist-hook:
	rm -rf $(distdir)/Inputs/.svn
	rm -rf $(distdir)/python/.svn
	rm -rf $(distdir)/Rivet/.svn

	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
	HWRUN = $(HERWIG) run

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

	if WANT_LIBFASTJET
	tests-LEP : test-LEP-VV test-LEP-VH test-LEP-VBF test-LEP-BB test-LEP-Quarks test-LEP-Leptons \
	test-LEP-default test-LEP-Powheg test-LEP-TopDecay
	else
	tests-LEP : test-LEP-VV test-LEP-VH test-LEP-VBF test-LEP-BB test-LEP-Quarks test-LEP-Leptons
	endif

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

	if WANT_LIBFASTJET
	tests-LHC : 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 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
	else
	tests-LHC : 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
	endif

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

	if WANT_LIBFASTJET
	test-LEP-% : Inputs/LEP-%.in LeptonTest.la LeptonJetTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
	else
	test-LEP-% : Inputs/LEP-%.in LeptonTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
	endif

	Rivet-LHC-Matchbox-% : Rivet/LHC-Matchbox-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-TVT-Matchbox-% : Rivet/TVT-Matchbox-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-TVT-Dipole-% : Rivet/TVT-Dipole-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-LHC-Dipole-% : Rivet/LHC-Dipole-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

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

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

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

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

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

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

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

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

	Rivet-LEP-Matchbox-% : Rivet/LEP-Matchbox-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-LEP-Dipole-% : Rivet/LEP-Dipole-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-BFactory-Matchbox-% : Rivet/BFactory-Matchbox-%.in
	if [ ! -d Rivet-$(notdir $(subst .in,,$<)) ]; then mkdir Rivet-$(notdir $(subst .in,,$<)); fi;
	cd Rivet-$(notdir $(subst .in,,$<)); echo `pwd`; \
	../$(HERWIG) read --repo ../$(REPO) -L ../$(top_builddir)/lib -i ../$(top_builddir)/src ../$<; \
	../$(HERWIG) run $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}; \
	mv $(notdir $(subst .in,.yoda,$<)) ..; \
	cd ..

	Rivet-LEP-% : Rivet/LEP-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-BFactory-% : Rivet/BFactory-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-TVT-% : Rivet/TVT-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-DIS-% : Rivet/DIS-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-LHC-% : Rivet/LHC-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-Star-% : Rivet/Star-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-SppS-% : Rivet/SppS-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-ISR-% : Rivet/ISR-%.in
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}

	Rivet-inputfiles: $(shell echo Rivet/LEP{,-Powheg,-Matchbox,-Dipole,-Matchbox-Powheg}-{10,14,22,35,44,91,130,133,136,161,172,177,183,189,192,196,197,200,202,206,91-nopi}.in) \
	- $(shell echo Rivet/BFactory{,-Powheg,-Matchbox,-Dipole,-Matchbox-Powheg}-{10.52,10.52-sym,10.54,10.45,10.58}.in) \
	- $(shell echo Rivet/BFactory-{Upsilon,Upsilon2,Upsilon4,Tau}.in) \
	+ $(shell echo Rivet/BFactory{,-Powheg,-Matchbox,-Dipole,-Matchbox-Powheg}-{10.52,10.52-sym,10.54,10.45,10.58-non}.in) \
	+ $(shell echo Rivet/BFactory-{Upsilon,Upsilon2,Upsilon4,Tau,10.58-res}.in) \
	$(shell echo Rivet/DIS{,-NoME,-Powheg,-Matchbox,-Dipole,-Matchbox-Powheg}-{e--LowQ2,e+-LowQ2,e+-HighQ2}.in) \
	$(shell echo Rivet/TVT{,-Powheg,-Matchbox,-Dipole,-Matchbox-Powheg}-{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-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,-Matchbox,-Matchbox-Powheg}-{Run-II-Jets-{0..11},Run-I-Jets-{1..8}}.in ) \
	$(shell echo Rivet/TVT{,-Dipole,-Matchbox,-Matchbox-Powheg}-{630-Jets-{1..3},300-Jets-1,900-Jets-1}.in ) \
	$(shell echo Rivet/TVT-{Run-I,Run-II,300,630,900}-UE.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Matchbox,-Matchbox-Powheg}-7-Jets-{0..15}.in ) \
	$(shell echo Rivet/LHC-{900,2360,2760,7,8,13}-UE.in ) \
	$(shell echo Rivet/LHC-{900,7}-UE-Long.in ) \
	$(shell echo Rivet/LHC{,-Dipole,-Matchbox,-Matchbox-Powheg}-7-Charm-{1..5}.in) \
	$(shell echo Rivet/LHC{,-Dipole,-Matchbox,-Matchbox-Powheg}-7-Bottom-{0..8}.in) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg}-7-Top-{L,SL,All}.in) \
	$(shell echo Rivet/Star-{UE,Jets-{1..4}}.in ) \
	$(shell echo Rivet/ISR-{30,44,53,62}-UE.in ) \
	$(shell echo Rivet/SppS-{53,63,200,500,900,546}-UE.in ) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Powheg,-Dipole}-{W-{e,mu},13-Z-{e,mu},Z-{e,mu},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},W-Z-{e,mu},8-Z-mu}.in) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Dipole}-{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}-{Z-b,Z-bb,W-b,8-Z-jj}.in) \
	$(shell echo Rivet/LHC-7-PromptPhoton-{1..4}.in) Rivet/LHC-GammaGamma-7.in \
	$(shell echo Rivet/LHC{,-Powheg}-7-{DiPhoton-GammaGamma,DiPhoton-GammaJet}.in) \
	$(shell echo Rivet/LHC{,-Powheg,-Matchbox,-Matchbox-Powheg,-Dipole}-{ggH,VBF,WH,ZH}.in) \
	$(shell echo Rivet/LHC{,-Powheg,-Matchbox,-Matchbox-Powheg,-Dipole}-8-{ggH,VBF,WH,ZH}{,-GammaGamma}.in) \
	$(shell echo Rivet/LHC{,-Matchbox,-Matchbox-Powheg,-Dipole}-ggHJet.in)

	Rivet-LEP: $(shell echo Rivet-LEP{,-Powheg,-Matchbox,-Dipole}-{10,14,22,35,44,91,130,133,136,161,172,177,183,189,192,196,197,200,202,206,91-nopi})
	rm -rf Rivet-LEP
	python/merge-LEP LEP
	python/merge-LEP LEP-Powheg
	python/merge-LEP LEP-Matchbox
	python/merge-LEP LEP-Dipole
	rivet-mkhtml -o Rivet-LEP LEP.yoda:Hw++ LEP-Powheg.yoda:Hw++-Powheg LEP-Matchbox.yoda:Hw++-Matchbox LEP-Dipole.yoda:Hw++-Dipole

	-Rivet-BFactory: $(shell echo Rivet-BFactory{,-Powheg,-Matchbox,-Dipole}-{10.52,10.52-sym,10.54,10.45,10.58}) \
	- $(shell echo Rivet-BFactory-{Upsilon,Upsilon2,Upsilon4,Tau})
	+Rivet-BFactory: $(shell echo Rivet-BFactory{,-Powheg,-Matchbox,-Dipole}-{10.52,10.52-sym,10.54,10.45}) \
	+ $(shell echo Rivet-BFactory-{Upsilon,Upsilon2,Upsilon4,Tau,10.58-res,10.58})
	rm -rf Rivet-BFactory
	python/merge-BFactory BFactory
	python/merge-BFactory BFactory-Powheg
	python/merge-BFactory BFactory-Matchbox
	python/merge-BFactory BFactory-Dipole
	rivet-mkhtml -o Rivet-BFactory BFactory.yoda:Hw++ BFactory-Powheg.yoda:Hw++-Powheg BFactory-Matchbox.yoda:Hw++-Matchbox BFactory-Dipole.yoda:Hw++-Dipole

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

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

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

	Rivet-TVT-Jets: $(shell echo Rivet-TVT-{Run-II-Jets-{0..11},Run-I-Jets-{1..8}} ) \
	$(shell echo Rivet-TVT-{630-Jets-{1..3},300-Jets-1,900-Jets-1} ) \
	$(shell echo Rivet-TVT-{Run-I,Run-II,300,630,900}-UE)
	python/merge-TVT-Energy TVT
	rivet-merge-CDF_2012_NOTE10874 TVT-300-Energy.yoda TVT-900-Energy.yoda TVT-1960-Energy.yoda
	flat2yoda RatioPlots.dat -o TVT-RatioPlots.yoda
	rm -rf Rivet-TVT-Jets
	python/merge-TVT-Jets TVT
	rivet-mkhtml -o Rivet-TVT-Jets TVT-Jets.yoda:Hw++

	Rivet-LHC-Jets: $(shell echo Rivet-LHC-7-Jets-{0..15} ) \
	$(shell echo Rivet-LHC-{900,2360,2760,7,8,13}-UE ) \
	$(shell echo Rivet-LHC-{900,7}-UE-Long ) \
	$(shell echo Rivet-LHC-7-Charm-{1..5}) \
	$(shell echo Rivet-LHC-7-Bottom-{0..8}) \
	$(shell echo Rivet-LHC-7-Top-{L,SL,All})
	rm -rf Rivet-LHC-Jets
	python/merge-LHC-Jets
	rivet-mkhtml -o Rivet-LHC-Jets LHC-Jets.yoda:Hw++

	Rivet-Star: $(shell echo Rivet-Star-{UE,Jets-{1..4}} )
	rm -rf Rivet-Star
	python/merge-Star Star
	rivet-mkhtml -o Rivet-Star Star.yoda

	Rivet-SppS: $(shell echo Rivet-ISR-{30,44,53,62}-UE ) \
	$(shell echo Rivet-SppS-{53,63,200,500,900,546}-UE )
	rm -rf Rivet-SppS
	python/merge-SppS SppS
	rivet-mkhtml -o Rivet-SppS SppS.yoda

	Rivet-LHC-EW: $(shell echo Rivet-LHC{,-Matchbox,-Powheg,-Dipole}-{13-Z-{e,mu},W-{e,mu},Z-{e,mu},Z-LowMass-{e,mu},Z-MedMass-e,WZ,WW-{emu,ll},ZZ-{ll,lv},W-Z-{e,mu},8-Z-mu}) \
	$(shell echo Rivet-LHC{,-Matchbox,-Dipole}-{7-W-Jet-{1..3}-e,7-Z-Jet-{0..3}-e,7-Z-Jet-0-mu}) \
	$(shell echo Rivet-LHC{-Matchbox,-Dipole}-{Z-b,Z-bb,W-b,8-Z-jj})
	rm -rf Rivet-LHC-EW;
	python/merge-LHC-EW LHC-{13-Z-{e,mu},W-{e,mu},Z-e,Z-mu,Z-LowMass-{e,mu},Z-MedMass-e,W-Z-{e,mu},WW-{emu,ll},WZ,ZZ-{ll,lv}}.yoda LHC-7-{W,Z}-Jet-{1,2,3}-e.yoda -o LHC-EW.yoda;
	python/merge-LHC-EW LHC-Matchbox-{13-Z-{e,mu},W-{e,mu},Z-{e,mu},Z-LowMass-{e,mu},Z-MedMass-e,W-Z-{e,mu},WW-{emu,ll},WZ,ZZ-{ll,lv}}.yoda LHC-Matchbox-7-{W,Z}-Jet-{1,2,3}-e.yoda -o LHC-Matchbox-EW.yoda;
	python/merge-LHC-EW LHC-Dipole-{13-Z-{e,mu},W-{e,mu},Z-{e,mu},Z-LowMass-{e,mu},Z-MedMass-e,W-Z-{e,mu},WW-{emu,ll},WZ,ZZ-{ll,lv}}.yoda LHC-Dipole-7-{W,Z}-Jet-{1,2,3}-e.yoda -o LHC-Dipole-EW.yoda;
	python/merge-LHC-EW LHC-Powheg-{W-{e,mu},Z-{e,mu},Z-LowMass-{e,mu},Z-MedMass-e,W-Z-{e,mu},WW-{emu,ll},WZ,ZZ-{ll,lv}}.yoda -o LHC-Powheg-EW.yoda;
	rivet-mkhtml -o Rivet-LHC-EW LHC-EW.yoda:Hw++ LHC-Powheg-EW.yoda:Hw++-Powheg LHC-Matchbox-EW.yoda:Hw++-Matchbox LHC-Matchbox-Z-b.yoda:Hw++-Matchbox-Zb \
	LHC-Matchbox-Z-bb.yoda:Hw++-Matchbox-Zbb LHC-Matchbox-W-b.yoda:Hw++-Matchbox-W-bb LHC-Dipole-EW.yoda:Hw++-Dipole \
	LHC-Dipole-Z-b.yoda:Hw++-Dipole-Zb LHC-Dipole-Z-bb.yoda:Hw++-Dipole-Zbb LHC-Dipole-W-b.yoda:Hw++-Dipole-W-bb;

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

	Rivet-LHC-Higgs: $(shell echo Rivet-LHC{,-Powheg}-{ggH,VBF,WH,ZH})\
	$(shell echo Rivet-LHC{,-Powheg}-8-{ggH,VBF,WH,ZH}{,-GammaGamma}) Rivet-LHC-ggHJet
	rm -rf Rivet-LHC-Higgs
	rivet-mkhtml -o Rivet-LHC-Higgs LHC-Powheg-ggH.yoda:gg-Powheg LHC-ggH.yoda:gg LHC-ggHJet.yoda:HJet \
	LHC-Powheg-VBF.yoda:VBF-Powheg LHC-VBF.yoda:VBF LHC-WH.yoda:WH LHC-ZH.yoda:ZH \
	LHC-Powheg-WH.yoda:WH-Powheg LHC-Powheg-ZH.yoda:ZH-Powheg

	tests-Rivet : Rivet-LEP Rivet-BFactory Rivet-DIS Rivet-TVT-WZ Rivet-TVT-Photon Rivet-TVT-Jets Rivet-LHC-Jets Rivet-Star Rivet-SppS Rivet-LHC-EW Rivet-LHC-Photon


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

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

	if WANT_LIBFASTJET
	test-LHC-% : Inputs/LHC-%.in HadronTest.la GammaTest.la HadronJetTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
	else
	test-LHC-% : Inputs/LHC-%.in HadronTest.la GammaTest.la
	$(HWREAD) $<
	$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
	endif

	clean-local:
	rm -f .out .log .tex .top .run .dump .mult .Bmult *.yoda
	diff --git a/Tests/Rivet/BFactory/BFactory-10.58.in b/Tests/Rivet/BFactory/BFactory-10.58-res.in
	copy from Tests/Rivet/BFactory/BFactory-10.58.in
	copy to Tests/Rivet/BFactory/BFactory-10.58-res.in
	--- a/Tests/Rivet/BFactory/BFactory-10.58.in
	+++ b/Tests/Rivet/BFactory/BFactory-10.58-res.in
	@@ -1,12 +1,11 @@
	create ThePEG::LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity
	set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxA 3.5*GeV
	set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxB 8.*GeV
	set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 9.99
	set /Herwig/Particles/e-:PDF /Herwig/Partons/NoPDF
	set /Herwig/Particles/e+:PDF /Herwig/Partons/NoPDF
	set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity
	# BELLE charm hadron production
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2001_S4598261
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2005_S6181155
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2007_S6895344
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2003_I593379
	diff --git a/Tests/Rivet/BFactory/BFactory-10.58.in b/Tests/Rivet/BFactory/BFactory-10.58.in
	--- a/Tests/Rivet/BFactory/BFactory-10.58.in
	+++ b/Tests/Rivet/BFactory/BFactory-10.58.in
	@@ -1,12 +1,9 @@
	create ThePEG::LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity
	set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxA 3.5*GeV
	set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxB 8.*GeV
	set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 9.99
	set /Herwig/Particles/e-:PDF /Herwig/Partons/NoPDF
	set /Herwig/Particles/e+:PDF /Herwig/Partons/NoPDF
	set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity
	# BELLE charm hadron production
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2001_S4598261
	insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2005_S6181155
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2007_S6895344
	-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2003_I593379
	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,1356 +1,1356 @@
	#! /usr/bin/env python
	import logging,sys,os
	from string import strip, 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 [...]")

	(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]

	collider=""
	# select the template to load
	# collider
	parameters = {}
	if(name.find("BFactory")==0) :
	collider="BFactory"
	elif(name.find("LEP")==0) :
	collider="LEP"
	elif(name.find("DIS")==0) :
	collider="DIS"
	elif(name.find("TVT")==0) :
	collider="TVT"
	elif(name.find("LHC-GammaGamma")==0) :
	collider="LHC-GammaGamma"
	elif(name.find("LHC")==0) :
	collider="LHC"
	elif(name.find("ISR")==0) :
	collider="ISR"
	elif(name.find("SppS")==0) :
	collider="SppS"
	elif(name.find("Star")==0) :
	collider="Star"
	simulation=""
	istart = 1
	print name
	if(name.find("Matchbox-Powheg")>0) :
	istart = 3
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/Powheg-DefaultShower.in\n"
	elif(name.find("Matchbox")>0) :
	istart = 2
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/MCatNLO-DefaultShower.in\n"
	elif(name.find("Dipole")>0) :
	istart = 2
	simulation="Matchbox"
	parameters["shower"] = "read Matchbox/MCatNLO-DipoleShower.in\n"
	elif(name.find("Powheg")>0) :
	istart = 2
	simulation="Powheg"

	if(simulation=="Matchbox") :
	parameters["bscheme"] = "read Matchbox/FiveFlavourScheme.in\n"
	if(parameters["shower"].find("Dipole")>=0) :
	parameters["bscheme"] += "read Matchbox/FiveFlavourNoBMassScheme.in\n"
	if(collider.find("DIS")<0) :
	parameters["nlo"] = "read Matchbox/MadGraph-OpenLoops.in\n"

	if(collider=="") :
	logging.error("Can\'t find collider")
	sys.exit(1)

	# find the template
	if(simulation=="") :
	if(collider.find("LHC-GammaGamma") >=0) :
	istart += 1
	templateName="Hadron-Gamma.in"
	elif(collider.find("TVT")>=0 or collider.find("LHC") >=0 or
	collider.find("ISR")>=0 or collider.find("SppS")>=0 or
	collider.find("Star")>=0) :
	templateName="Hadron.in"
	elif(collider.find("BFactory")<0) :
	templateName= "%s.in" % (collider)
	else :
	templateName= "LEP.in"
	else :
	if(collider.find("TVT")>=0 or collider.find("LHC") >=0 or
	collider.find("ISR")>=0 or collider.find("SppS")>=0 or
	collider.find("Star")>=0) :
	templateName= "Hadron-%s.in" % (simulation)
	elif(collider.find("BFactory")<0) :
	templateName= "%s-%s.in" % (collider,simulation)
	else :
	templateName= "LEP-%s.in" % (simulation)
	with open(os.path.join("Rivet/Templates",templateName), 'r') as f:
	templateText = f.read()
	template = Template( templateText )
	# work out the name of the parameter file
	nameSplit=name.split("-")
	parameterName=nameSplit[istart]
	for i in range(istart+1,len(nameSplit)) :
	parameterName += "-%s" % nameSplit[i]

	# work out the process and parameters
	process=""
	# Bfactory
	if(collider=="BFactory") :
	if(simulation=="") :
	- process = "set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 4"
	- if(parameterName=="10.58") :
	- process += "\ncreate Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so\nset /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)\nset /Herwig/MatrixElements/MEUpsilon:Coupling 0.0004151809\ninsert /Herwig/MatrixElements/SimpleEE:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon"
	+ if(parameterName=="10.58-res") :
	+ process += "\ncreate Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so\nset /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)\nset /Herwig/MatrixElements/MEUpsilon:Coupling 0.0004151809\nset /Herwig/MatrixElements/SimpleEE:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon"
	+ elif(parameterName=="10.58") :
	+ process += "\ncreate Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so\nset /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)\nset /Herwig/MatrixElements/MEUpsilon:Coupling 0.0004151809\ninsert /Herwig/MatrixElements/SimpleEE:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon\n"
	+ process += "set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 4\n"
	+ else :
	+ process += "set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 4\n"
	elif(simulation=="Powheg") :
	- process = "set /Herwig/MatrixElements/PowhegMEee2gZ2qq:MaximumFlavour 4"
	- if(parameterName=="10.58") :
	- process += "\ncreate Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so\nset /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)\nset /Herwig/MatrixElements/MEUpsilon:Coupling 0.0004151809\ninsert /Herwig/MatrixElements/SimpleEE:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon"
	+ process = "set /Herwig/MatrixElements/PowhegMEee2gZ2qq:MaximumFlavour 4\n"
	elif(simulation=="Matchbox" ) :
	- process = "do Factory:Process e- e+ -> u ubar\ndo Factory:Process e- e+ -> d dbar\ndo Factory:Process e- e+ -> c cbar\ndo Factory:Process e- e+ -> s sbar"
	- if(parameterName=="10.58") :
	- process += "\ninsert /Herwig/Generators/EventGenerator:EventHandler:SubProcessHandlers 0 /Herwig/MatrixElements/SimpleEE\ncreate Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so\nset /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S)\nset /Herwig/MatrixElements/MEUpsilon:Coupling 0.0004151809\ninsert /Herwig/MatrixElements/SimpleEE:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon\n"
	+ process = "do Factory:Process e- e+ -> u ubar\ndo Factory:Process e- e+ -> d dbar\ndo Factory:Process e- e+ -> c cbar\ndo Factory:Process e- e+ -> s sbar\n"
	# DIS
	elif(collider=="DIS") :
	if(simulation=="") :
	if(parameterName.find("NoME")>=0) :
	process = "set /Herwig/Shower/Evolver:MECorrMode 0"
	parameterName=parameterName.replace("NoME-","")
	else :
	process = ""
	elif(simulation=="Powheg") :
	process = ""
	elif(simulation=="Matchbox" ) :
	if(parameterName.find("e-")>=0) :
	process="do Factory:Process e- p -> e- j"
	else :
	process="do Factory:Process e+ p -> e+ j"
	# LEP
	elif(collider=="LEP") :
	if(simulation=="") :
	process=""
	if(parameterName=="10") :
	process="set /Herwig/MatrixElements/MEee2gZ2qq:MaximumFlavour 4"
	elif(simulation=="Powheg") :
	process=""
	if(parameterName=="10") :
	process="set /Herwig/MatrixElements/PowhegMEee2gZ2qq:MaximumFlavour 4"
	elif(simulation=="Matchbox" ) :
	if(parameterName=="10") :
	process="do Factory:Process e- e+ -> u ubar\ndo Factory:Process e- e+ -> d dbar\ndo Factory:Process e- e+ -> c cbar\ndo Factory:Process e- e+ -> s sbar"
	else :
	process="do Factory:Process e- e+ -> j j"
	# TVT
	elif(collider=="TVT") :
	process="set /Herwig/Generators/EventGenerator:EventHandler:BeamB /Herwig/Particles/pbar-\n"
	if(parameterName.find("Run-II")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 1960.0\n"
	elif(parameterName.find("Run-I")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 1800.0\n"
	elif(parameterName.find("900")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 900.0\n"
	elif(parameterName.find("630")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 630.0\n"
	elif(parameterName.find("300")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 300.0\n"

	if(simulation=="") :
	if(parameterName.find("PromptPhoton")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaJet\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 15.\n"
	elif(parameterName.find("DiPhoton-GammaGamma")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGamma\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif(parameterName.find("DiPhoton-GammaJet")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaJet\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	elif(parameterName.find("UE")>=0) :
	process += "insert SimpleQCD:MatrixElements[0] MEMinBias\n"
	process += "set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	process += "set /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts\n"
	process += "create Herwig::MPIXSecReweighter /Herwig/Generators/MPIXSecReweighter\n"
	process += "insert /Herwig/Generators/EventGenerator:EventHandler:PostSubProcessHandlers 0 /Herwig/Generators/MPIXSecReweighter\n"
	process += "set /Herwig/Decays/DecayHandler:LifeTimeOption 0\n"
	process += "set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm\n"
	elif(parameterName.find("Jets")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEQCD2to2\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("Run-II-Jets-10")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 500.*GeV\n"
	elif(parameterName.find("Run-II-Jets-11")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 900.*GeV\n"
	elif(parameterName.find("Run-I-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	elif(parameterName.find("Run-I-Jets-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 40.\n"
	elif(parameterName.find("Run-I-Jets-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 65.\n"
	elif(parameterName.find("Run-I-Jets-4")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 90.\n"
	elif(parameterName.find("Run-I-Jets-5")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 160.\n"
	elif(parameterName.find("Run-I-Jets-6")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 100.*GeV\n"
	elif(parameterName.find("Run-I-Jets-7")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 400.*GeV\n"
	elif(parameterName.find("Run-I-Jets-8")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 700.*GeV\n"
	elif(parameterName.find("Run-II-Jets-0")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 15.\n"
	elif(parameterName.find("Run-II-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 25.\n"
	elif(parameterName.find("Run-II-Jets-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 40.\n"
	elif(parameterName.find("Run-II-Jets-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 60.\n"
	elif(parameterName.find("Run-II-Jets-4")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 85.\n"
	elif(parameterName.find("Run-II-Jets-5")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 110.\n"
	elif(parameterName.find("Run-II-Jets-6")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 160.\n"
	elif(parameterName.find("Run-II-Jets-7")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 250.\n"
	elif(parameterName.find("Run-II-Jets-8")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 100.*GeV\n"
	elif(parameterName.find("Run-II-Jets-9")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 300.*GeV\n"
	elif(parameterName.find("900-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 10.\n"
	elif(parameterName.find("300-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 6.\n"
	elif(parameterName.find("630-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	elif(parameterName.find("630-Jets-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 40.\n"
	elif(parameterName.find("630-Jets-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 75.\n"
	elif(parameterName.find("900-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 10.\n"
	elif(parameterName.find("Run-I-WZ")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Electron\ninsert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Run-I-W")>=0 or parameterName.find("Run-II-W")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Electron\n"
	elif(parameterName.find("Run-I-Z")>=0 or parameterName.find("Run-II-Z-e")>=0) :
	process +="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Run-II-Z-LowMass-mu")>=0) :
	process +="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 25GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 70GeV\n"
	elif(parameterName.find("Run-II-Z-HighMass-mu")>=0) :
	process +="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 150GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 600GeV\n"
	elif(parameterName.find("Run-II-Z-mu")>=0) :
	process +="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	elif(simulation=="Powheg") :
	if(parameterName.find("Run-I-WZ")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff\nset PowhegMEqq2W2ff:Process Electron\ninsert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Run-I-W")>=0 or parameterName.find("Run-II-W")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff\nset PowhegMEqq2W2ff:Process Electron\n"
	elif(parameterName.find("Run-I-Z")>=0 or parameterName.find("Run-II-Z-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Run-II-Z-LowMass-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 25GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 70GeV\n"
	elif(parameterName.find("Run-II-Z-HighMass-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 150GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 600GeV\n"
	elif(parameterName.find("Run-II-Z-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Muon\n"
	elif(parameterName.find("DiPhoton-GammaGamma")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGammaPowheg\n"
	process+="set MEGammaGammaPowheg:Process GammaGamma\n"
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGamma\n"
	process+="set MEGammaGamma:Process gg\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif(parameterName.find("DiPhoton-GammaJet")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGammaPowheg\n"
	process+="set MEGammaGammaPowheg:Process VJet\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	elif(simulation=="Matchbox" ) :
	if(parameterName.find("Jets")>=0) :
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	process+="do Factory:Process p p j j\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/MaxJetPtScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("Run-II-Jets-10")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 500.*GeV\n"
	elif(parameterName.find("Run-II-Jets-11")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 900.*GeV\n"
	elif(parameterName.find("Run-II-Jets-12")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 300.*GeV\n"
	elif(parameterName.find("Run-I-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.\n"
	elif(parameterName.find("Run-I-Jets-2")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 40.\n"
	elif(parameterName.find("Run-I-Jets-3")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 65.\n"
	elif(parameterName.find("Run-I-Jets-4")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 90.\n"
	elif(parameterName.find("Run-I-Jets-5")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 160.\n"
	elif(parameterName.find("Run-I-Jets-6")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 100.*GeV\n"
	elif(parameterName.find("Run-I-Jets-7")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 400.*GeV\n"
	elif(parameterName.find("Run-I-Jets-8")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 700.*GeV\n"
	elif(parameterName.find("Run-II-Jets-0")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 15.\n"
	elif(parameterName.find("Run-II-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 25.\n"
	elif(parameterName.find("Run-II-Jets-2")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 40.\n"
	elif(parameterName.find("Run-II-Jets-3")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 60.\n"
	elif(parameterName.find("Run-II-Jets-4")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 85.\n"
	elif(parameterName.find("Run-II-Jets-5")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 110.\n"
	elif(parameterName.find("Run-II-Jets-6")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 160.\n"
	elif(parameterName.find("Run-II-Jets-7")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 250.\n"
	elif(parameterName.find("Run-II-Jets-8")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 100.*GeV\n"
	elif(parameterName.find("Run-II-Jets-9")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 300.*GeV\n"
	elif(parameterName.find("900-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 10.\n"
	elif(parameterName.find("300-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 6.\n"
	elif(parameterName.find("630-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.\n"
	elif(parameterName.find("630-Jets-2")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 40.\n"
	elif(parameterName.find("630-Jets-3")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 75.\n"
	elif(parameterName.find("900-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 10.\n"
	elif(parameterName.find("Run-I-WZ")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar e+ e-\ndo Factory:Process p pbar e+ nu\ndo Factory:Process p pbar e- nu\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Run-I-W")>=0 or parameterName.find("Run-II-W")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar e+ nu\ndo Factory:Process p pbar e- nu\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Run-I-Z")>=0 or parameterName.find("Run-II-Z-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar e+ e-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Run-II-Z-LowMass-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar mu+ mu-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 25GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 70GeV\n"
	elif(parameterName.find("Run-II-Z-HighMass-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar mu+ mu-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 150.GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 600GeV\n"
	elif(parameterName.find("Run-II-Z-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p pbar mu+ mu-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	# Star
	elif(collider=="Star" ) :
	process = "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+= "set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 200.0\n"
	process+= "set /Herwig/Cuts/QCDCuts:X2Min 0.01\n"
	if(simulation=="") :
	if(parameterName.find("UE")>=0) :
	process += "insert SimpleQCD:MatrixElements[0] MEMinBias\n"
	process += "set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	process += "set /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts\n"
	process += "create Herwig::MPIXSecReweighter /Herwig/Generators/MPIXSecReweighter\n"
	process += "insert /Herwig/Generators/EventGenerator:EventHandler:PostSubProcessHandlers 0 /Herwig/Generators/MPIXSecReweighter\n"
	else :
	process+="insert SimpleQCD:MatrixElements[0] MEQCD2to2\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 2.\n"
	elif(parameterName.find("Jets-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	elif(parameterName.find("Jets-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	elif(parameterName.find("Jets-4")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 25.\n"
	else :
	logging.error("Star not supported for %s " % simulation)
	sys.exit(1)
	# ISR and SppS
	elif(collider=="ISR" or collider =="SppS" ) :
	process="set /Herwig/Decays/DecayHandler:LifeTimeOption 0\n"
	process+="set /Herwig/Decays/DecayHandler:MaxLifeTime 10*mm\n"
	if(collider=="SppS") :
	process ="set /Herwig/Generators/EventGenerator:EventHandler:BeamB /Herwig/Particles/pbar-\n"
	if(parameterName.find("30")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 30.4\n"
	elif(parameterName.find("44")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 44.4\n"
	elif(parameterName.find("53")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 53.0\n"
	elif(parameterName.find("62")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 62.2\n"
	elif(parameterName.find("63")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 63.0\n"
	elif(parameterName.find("200")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 200.0\n"
	elif(parameterName.find("500")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 500.0\n"
	elif(parameterName.find("546")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 546.0\n"
	elif(parameterName.find("900")>=0) :
	process+="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 900.0\n"
	if(simulation=="") :
	process += "insert SimpleQCD:MatrixElements[0] MEMinBias\n"
	process += "set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	process += "set /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts\n"
	process += "create Herwig::MPIXSecReweighter /Herwig/Generators/MPIXSecReweighter\n"
	process += "insert /Herwig/Generators/EventGenerator:EventHandler:PostSubProcessHandlers 0 /Herwig/Generators/MPIXSecReweighter\n"
	else :
	logging.error(" SppS and ISR not supported for %s " % simulation)
	sys.exit(1)
	# LHC
	elif(collider=="LHC") :
	if(parameterName.find("7-")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0\n"
	elif(parameterName.find("8-")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 8000.0\n"
	elif(parameterName.find("13-")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 13000.0\n"
	elif(parameterName.find("900")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 900.0\n"
	elif(parameterName.find("2360")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 2360.0\n"
	elif(parameterName.find("2760")==0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 2760.0\n"
	else :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0\n"
	if(simulation=="") :
	if(parameterName.find("8-VBF")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2HiggsVBF\n"
	elif(parameterName.find("VBF")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="insert SimpleQCD:MatrixElements[0] MEPP2HiggsVBF\n"
	elif(parameterName.find("ggHJet")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHiggsJet\n"
	elif(parameterName.find("8-ggH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEHiggs\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHiggsJet\n"
	process+="set MEHiggsJet:Process qqbar\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("ggH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHiggs\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHiggsJet\n"
	process+="set MEHiggsJet:Process qqbar\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("PromptPhoton")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaJet\n"
	if(parameterName.find("PromptPhoton-1")>=0) :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	elif(parameterName.find("PromptPhoton-2")>=0) :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 25.\n"
	elif(parameterName.find("PromptPhoton-3")>=0) :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 80.\n"
	elif(parameterName.find("PromptPhoton-4")>=0) :
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 150.\n"
	elif(parameterName.find("DiPhoton-GammaGamma")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGamma\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif(parameterName.find("DiPhoton-GammaJet")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaJet\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	elif(parameterName.find("8-WH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2WH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("8-ZH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2ZH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("WH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+;\n"
	process+="insert SimpleQCD:MatrixElements[0] MEPP2WH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("ZH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes Z0->e-,e+; Z0->mu-,mu+;\n"
	process+="insert SimpleQCD:MatrixElements[0] MEPP2ZH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("UE")>=0) :
	process += "insert SimpleQCD:MatrixElements[0] MEMinBias\n"
	process += "set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	process += "set /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts\n"
	process += "create Herwig::MPIXSecReweighter /Herwig/Generators/MPIXSecReweighter\n"
	process += "insert /Herwig/Generators/EventGenerator:EventHandler:PostSubProcessHandlers 0 /Herwig/Generators/MPIXSecReweighter\n"
	if(parameterName.find("Long")>=0) :
	process += "set /Herwig/Decays/DecayHandler:MaxLifeTime 100*mm\n"
	elif(parameterName.find("7-Jets")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEQCD2to2\n"
	process+="set MEQCD2to2:MaximumFlavour 5\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("7-Jets-0")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	elif(parameterName.find("7-Jets-10")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 200.*GeV\n"
	elif(parameterName.find("7-Jets-11")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 600.*GeV\n"
	elif(parameterName.find("7-Jets-12")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 1000.*GeV\n"
	elif(parameterName.find("7-Jets-13")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 1600.*GeV\n"
	elif(parameterName.find("7-Jets-14")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 2200.*GeV\n"
	elif(parameterName.find("7-Jets-15")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 2800.*GeV\n"
	elif(parameterName.find("7-Jets-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 10.\n"
	elif(parameterName.find("7-Jets-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	elif(parameterName.find("7-Jets-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 40.\n"
	elif(parameterName.find("7-Jets-4")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 70.\n"
	elif(parameterName.find("7-Jets-5")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 150.\n"
	elif(parameterName.find("7-Jets-6")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 200.\n"
	elif(parameterName.find("7-Jets-7")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 300.\n"
	elif(parameterName.find("7-Jets-8")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 500.\n"
	elif(parameterName.find("7-Jets-9")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 90.*GeV\n"
	elif(parameterName.find("7-Charm")>=0 or \
	parameterName.find("7-Bottom")>=0) :
	if(parameterName.find("7-Bottom")>=0) :
	process+="cp MEHeavyQuark MEBottom\n"
	process+="set MEBottom:QuarkType Bottom\n"
	process+="insert SimpleQCD:MatrixElements[0] MEBottom\n"
	else :
	process+="cp MEHeavyQuark MECharm\n"
	process+="set MECharm:QuarkType Charm\n"
	process+="insert SimpleQCD:MatrixElements[0] MECharm\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("7-Heavy-0")>=0) :
	if(parameterName.find("7-Bottom")>=0) :
	process+="set MEBottom:Process Pair\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.\n"
	elif(parameterName.find("-1")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	elif(parameterName.find("-2")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 20.\n"
	elif(parameterName.find("-3")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 50.\n"
	elif(parameterName.find("-4")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 80.\n"
	elif(parameterName.find("-5")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 110.\n"
	elif(parameterName.find("-6")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 90.*GeV\n"
	elif(parameterName.find("-7")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 340.*GeV\n"
	elif(parameterName.find("-8")>=0) :
	process+="set /Herwig/Cuts/JetKtCut:MinKT 30.\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 500.*GeV\n"
	elif(parameterName.find("7-Top-L")>=0) :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHeavyQuark\n"
	process+="do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("7-Top-SL")>=0) :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHeavyQuark\n"
	process+="set /Herwig/Particles/t:Synchronized Not_synchronized\n"
	process+="set /Herwig/Particles/tbar:Synchronized Not_synchronized\n"
	process+="do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;\n"
	process+="do /Herwig/Particles/tbar:SelectDecayModes tbar->b,bbar,cbar; tbar->bbar,cbar,d; tbar->bbar,cbar,s; tbar->bbar,s,ubar; tbar->bbar,ubar,d;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("7-Top-All")>=0) :
	process+="set MEHeavyQuark:QuarkType Top\n"
	process+="insert SimpleQCD:MatrixElements[0] MEHeavyQuark\n"
	elif(parameterName.find("WZ")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2VV\nset MEPP2VV:Process WZ\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_ebar,e-; /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("WW-emu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2VV\nset MEPP2VV:Process WW\n"
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("WW-ll")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2VV\nset MEPP2VV:Process WW\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+; /Herwig/Particles/W+/W+->nu_tau,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("ZZ-ll")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2VV\nset MEPP2VV:Process ZZ\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("ZZ-lv")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEPP2VV\nset MEPP2VV:Process ZZ\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+; /Herwig/Particles/Z0/Z0->nu_e,nu_ebar; /Herwig/Particles/Z0/Z0->nu_mu,nu_mubar; /Herwig/Particles/Z0/Z0->nu_tau,nu_taubar;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("W-Z-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Electron\n"
	elif(parameterName.find("W-Z-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Muon\n"
	elif(parameterName.find("W-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Electron\n"
	elif(parameterName.find("W-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Muon\n"
	elif(parameterName.find("Z-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Z-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	elif(parameterName.find("Z-LowMass-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 20.GeV\nset /Herwig/Cuts/MassCut:MinM 20.GeV\nset /Herwig/Cuts/MassCut:MaxM 70.*GeV\n"
	elif(parameterName.find("Z-MedMass-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Electron\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 40.GeV\nset /Herwig/Cuts/MassCut:MinM 40.GeV\nset /Herwig/Cuts/MassCut:MaxM 130.*GeV\n"
	elif(parameterName.find("Z-LowMass-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 10.GeV\nset /Herwig/Cuts/MassCut:MinM 10.GeV\nset /Herwig/Cuts/MassCut:MaxM 70.*GeV\n"
	elif(parameterName.find("W-Jet")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEWJet\nset MEWJet:WDecay Electron\n"
	if(parameterName.find("W-Jet-1-e")>=0) :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 100.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-1-e","W-Jet-e")
	elif(parameterName.find("W-Jet-2-e")>=0) :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 190.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-2-e","W-Jet-e")
	elif(parameterName.find("W-Jet-3-e")>=0) :
	process+="set /Herwig/Cuts/WBosonKtCut:MinKT 270.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-3-e","W-Jet-e")
	elif(parameterName.find("Z-Jet")>=0) :
	if(parameterName.find("-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEZJet\nset MEZJet:ZDecay Electron\n"
	if(parameterName.find("Z-Jet-0-e")>=0) :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 35.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-1-e")>=0) :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 100.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-1-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-2-e")>=0) :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 190.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-2-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-3-e")>=0) :
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 270.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-3-e","Z-Jet-e")
	else :
	process+="insert SimpleQCD:MatrixElements[0] MEZJet\nset MEZJet:ZDecay Muon\n"
	process+="set /Herwig/Cuts/ZBosonKtCut:MinKT 35.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-mu","Z-Jet-mu")
	else :
	logging.error(" Process %s not supported for internal matrix elements" % name)
	sys.exit(1)
	elif(simulation=="Powheg") :
	if(parameterName.find("8-VBF")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2HiggsVBF\n"
	elif(parameterName.find("VBF")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2HiggsVBF\n"
	elif(parameterName.find("ggHJet")>=0) :
	logging.error(" Process %s not supported for POWHEG matrix elements" % name)
	sys.exit(1)
	elif(parameterName.find("8-ggH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEHiggs\n"
	elif(parameterName.find("ggH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEHiggs\n"
	elif(parameterName.find("8-WH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2WH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("8-ZH")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2ZH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("WH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+;\n"
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2WH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("ZH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes Z0->e-,e+; Z0->mu-,mu+;\n"
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEPP2ZH\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV\n"
	elif(parameterName.find("UE")>=0) :
	logging.error(" Process %s not supported for powheg matrix elements" % name)
	sys.exit(1)
	elif(parameterName.find("WZ")>=0) :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\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 SimpleQCD:MatrixElements[0] PowhegMEPP2VV\nset PowhegMEPP2VV:Process WZ\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_ebar,e-; /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("WW-emu")>=0) :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\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 SimpleQCD:MatrixElements[0] PowhegMEPP2VV\nset PowhegMEPP2VV:Process WW\n"
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("WW-ll")>=0) :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\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 SimpleQCD:MatrixElements[0] PowhegMEPP2VV\nset PowhegMEPP2VV:Process WW\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+; /Herwig/Particles/W+/W+->nu_tau,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("ZZ-ll")>=0) :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\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 SimpleQCD:MatrixElements[0] PowhegMEPP2VV\nset PowhegMEPP2VV:Process ZZ\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("ZZ-lv")>=0) :
	process+="create Herwig::HwDecayHandler /Herwig/NewPhysics/DecayHandler\n"
	process+="set /Herwig/NewPhysics/DecayHandler:NewStep No\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 SimpleQCD:MatrixElements[0] PowhegMEPP2VV\nset PowhegMEPP2VV:Process ZZ\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+; /Herwig/Particles/Z0/Z0->nu_e,nu_ebar; /Herwig/Particles/Z0/Z0->nu_mu,nu_mubar; /Herwig/Particles/Z0/Z0->nu_tau,nu_taubar;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("W-Z-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff\nset PowhegMEqq2W2ff:Process Electron\n"
	elif(parameterName.find("W-Z-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff\nset MEqq2gZ2ff:Process Muon\n"
	process+="insert SimpleQCD:MatrixElements[0] MEqq2W2ff\nset MEqq2W2ff:Process Muon\n"
	elif(parameterName.find("W-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff\nset PowhegMEqq2W2ff:Process Electron\n"
	elif(parameterName.find("W-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2W2ff\nset PowhegMEqq2W2ff:Process Muon\n"
	elif(parameterName.find("Z-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	elif(parameterName.find("Z-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Muon\n"
	elif(parameterName.find("Z-LowMass-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 20.GeV\nset /Herwig/Cuts/MassCut:MinM 20.GeV\nset /Herwig/Cuts/MassCut:MaxM 70.*GeV\n"
	elif(parameterName.find("Z-MedMass-e")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Electron\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 40.GeV\nset /Herwig/Cuts/MassCut:MinM 40.GeV\nset /Herwig/Cuts/MassCut:MaxM 130.*GeV\n"
	elif(parameterName.find("Z-LowMass-mu")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff\nset PowhegMEqq2gZ2ff:Process Muon\n"
	process+="set /Herwig/Cuts/QCDCuts:MHatMin 10.GeV\nset /Herwig/Cuts/MassCut:MinM 10.GeV\nset /Herwig/Cuts/MassCut:MaxM 70.*GeV\n"
	elif(parameterName.find("DiPhoton-GammaGamma")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGammaPowheg\n"
	process+="set MEGammaGammaPowheg:Process GammaGamma\n"
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGamma\n"
	process+="set MEGammaGamma:Process gg\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaGamma","")
	elif(parameterName.find("DiPhoton-GammaJet")>=0) :
	process+="insert SimpleQCD:MatrixElements[0] MEGammaGammaPowheg\n"
	process+="set MEGammaGammaPowheg:Process VJet\n"
	process+="set /Herwig/Cuts/PhotonKtCut:MinKT 5.\n"
	process+="set /Herwig/Cuts/JetKtCut:MinKT 5.\n"
	parameterName=parameterName.replace("-GammaJet","")
	else :
	logging.error(" Process %s not supported for internal POWHEG matrix elements" % name)
	sys.exit(1)

	elif(simulation=="Matchbox" ) :
	if(parameterName.find("8-VBF")>=0) :
	parameters["nlo"] = "read Matchbox/VBFNLO.in\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 3\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/Z0\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/W+\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/W-\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/gamma\n"
	process+="do Factory:DiagramGenerator:TimeLikeRange 0 0\n"
	process+="do Factory:Process p p h0 j j\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("VBF")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	parameters["nlo"] = "read Matchbox/VBFNLO.in\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 3\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/Z0\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/W+\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/W-\n"
	process+="insert Factory:DiagramGenerator:RestrictLines 0 /Herwig/Particles/gamma\n"
	process+="do Factory:DiagramGenerator:TimeLikeRange 0 0\n"
	process+="do Factory:Process p p h0 j j\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("ggHJet")>=0) :
	parameters["nlo"] = "read Matchbox/MadGraph-GoSam.in\nread Matchbox/HiggsEffective.in\n"
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="set Factory:OrderInAlphaS 3\nset Factory:OrderInAlphaEW 1\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p h0 j\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("8-ggH")>=0) :
	parameters["nlo"] = "read Matchbox/MadGraph-GoSam.in\nread Matchbox/HiggsEffective.in\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 1\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p h0\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("ggH")>=0) :
	parameters["nlo"] = "read Matchbox/MadGraph-GoSam.in\nread Matchbox/HiggsEffective.in\n"
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;\n"
	process+="set /Herwig/Particles/tau-:Stable Stable\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 1\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p h0\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("8-WH")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p W+ h0\n"
	process+="do Factory:Process p p W- h0\n"
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("8-ZH")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p Z0 h0\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 125.7\n"
	elif(parameterName.find("WH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 3\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p e+ nu h0\n"
	process+="do Factory:Process p p e- nu h0\n"
	process+="do Factory:Process p p mu+ nu h0\n"
	process+="do Factory:Process p p mu- nu h0\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("ZH")>=0) :
	process+="do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 3\n"
	process+="set /Herwig/Particles/h0:HardProcessWidth 0.\n"
	process+="do Factory:Process p p e+ e- h0\n"
	process+="do Factory:Process p p mu+ mu- h0\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("UE")>=0) :
	logging.error(" Process %s not supported for Matchbox matrix elements" % name)
	sys.exit(1)
	elif(parameterName.find("7-Jets")>=0) :
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	process+="do Factory:Process p p j j\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/MaxJetPtScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("7-Jets-0")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 5.\n"
	elif(parameterName.find("7-Jets-10")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 200.*GeV\n"
	elif(parameterName.find("7-Jets-11")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 600.*GeV\n"
	elif(parameterName.find("7-Jets-12")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 1000.*GeV\n"
	elif(parameterName.find("7-Jets-13")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 1600.*GeV\n"
	elif(parameterName.find("7-Jets-14")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 2200.*GeV\n"
	elif(parameterName.find("7-Jets-15")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 2800.*GeV\n"
	elif(parameterName.find("7-Jets-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 10.\n"
	elif(parameterName.find("7-Jets-2")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.\n"
	elif(parameterName.find("7-Jets-3")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 40.\n"
	elif(parameterName.find("7-Jets-4")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 70.\n"
	elif(parameterName.find("7-Jets-5")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 150.\n"
	elif(parameterName.find("7-Jets-6")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 200.\n"
	elif(parameterName.find("7-Jets-7")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 300.\n"
	elif(parameterName.find("7-Jets-8")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 500.\n"
	elif(parameterName.find("7-Jets-9")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 90.*GeV\n"
	elif(parameterName.find("7-Charm")>=0 or \
	parameterName.find("7-Bottom")>=0) :
	parameters["bscheme"]="read Matchbox/FourFlavourScheme.in"
	process+="set /Herwig/Particles/b:HardProcessMass 4.2*GeV\n"
	process+="set /Herwig/Particles/bbar:HardProcessMass 4.2*GeV\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	if(parameterName.find("7-Bottom")>=0) :
	process+="do Factory:Process p p b bbar\n"
	else:
	process+="do Factory:Process p p c cbar\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/MaxJetPtScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0\n"
	if(parameterName.find("-0")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 0.\n"
	elif(parameterName.find("-1")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 5.\n"
	elif(parameterName.find("-2")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 20.\n"
	elif(parameterName.find("-3")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 50.\n"
	elif(parameterName.find("-4")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 80.\n"
	elif(parameterName.find("-5")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 110.\n"
	elif(parameterName.find("-6")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 90.*GeV\n"
	elif(parameterName.find("-7")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 340.*GeV\n"
	elif(parameterName.find("-8")>=0) :
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 25.\n"
	process+="create ThePEG::JetPairRegion /Herwig/Cuts/JetPairMass JetCuts.so\n"
	process+="set /Herwig/Cuts/JetPairMass:FirstRegion /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/JetPairMass:SecondRegion /Herwig/Cuts/SecondJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetPairRegions 0 /Herwig/Cuts/JetPairMass\n"
	process+="set /Herwig/Cuts/JetPairMass:MassMin 500.*GeV\n"
	elif(parameterName.find("7-Top-L")>=0) :
	process+="set /Herwig/Particles/t:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/tbar:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	process+="do Factory:Process p p t tbar\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/TopPairMTScale\n"
	process+="do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("7-Top-SL")>=0) :
	process+="set /Herwig/Particles/t:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/tbar:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	process+="do Factory:Process p p t tbar\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/TopPairMTScale\n"
	process+="set /Herwig/Particles/t:Synchronized Not_synchronized\n"
	process+="set /Herwig/Particles/tbar:Synchronized Not_synchronized\n"
	process+="do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;\n"
	process+="do /Herwig/Particles/tbar:SelectDecayModes tbar->b,bbar,cbar; tbar->bbar,cbar,d; tbar->bbar,cbar,s; tbar->bbar,s,ubar; tbar->bbar,ubar,d;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	elif(parameterName.find("7-Top-All")>=0) :
	process+="set /Herwig/Particles/t:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/tbar:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 0\n"
	process+="do Factory:Process p p t tbar\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/TopPairMTScale\n"
	elif(parameterName.find("WZ")>=0) :
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/W-:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p W+ Z0\ndo Factory:Process p p W- Z0\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 171.6*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_ebar,e-; /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("WW-emu")>=0) :
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/W-:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p W+ W-\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 160.8*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/Particles/W+:Synchronized 0\n"
	process+="set /Herwig/Particles/W-:Synchronized 0\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+;\n"
	process+="do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_mubar,mu-;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("WW-ll")>=0) :
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/W-:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p W+ W-\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 160.8*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+; /Herwig/Particles/W+/W+->nu_tau,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("ZZ-ll")>=0) :
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/W-:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p Z0 Z0\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 182.2*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("ZZ-lv")>=0) :
	process+="set /Herwig/Particles/W+:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/W-:HardProcessWidth 0.*GeV\n"
	process+="set /Herwig/Particles/Z0:HardProcessWidth 0.*GeV\n"
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p Z0 Z0\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 182.2*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+; /Herwig/Particles/Z0/Z0->tau-,tau+; /Herwig/Particles/Z0/Z0->nu_e,nu_ebar; /Herwig/Particles/Z0/Z0->nu_mu,nu_mubar; /Herwig/Particles/Z0/Z0->nu_tau,nu_taubar;\n"
	process+="create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter\n"
	process+="insert /Herwig/Generators/EventGenerator:EventHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("W-Z-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\n"
	process+="do Factory:Process p p e+ e-\ndo Factory:Process p p e+ nu\ndo Factory:Process p p e- nu\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("W-Z-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\n"
	process+="do Factory:Process p p mu+ mu-\ndo Factory:Process p p mu+ nu\ndo Factory:Process p p mu- nu\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("W-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ nu\ndo Factory:Process p p e- nu\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("W-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p mu+ nu\ndo Factory:Process p p mu- nu\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Z-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ e-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Z-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p mu+ mu-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Z-jj")>=0) :
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 2\n"
	process+="do Factory:Process p p e+ e- j j\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 40.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Z-LowMass-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ e-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 20GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 70GeV\n"
	elif(parameterName.find("Z-MedMass-e")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ e-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 40GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 130GeV\n"
	elif(parameterName.find("Z-LowMass-mu")>=0) :
	process+="set Factory:OrderInAlphaS 0\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p mu+ mu-\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/LeptonPairMassScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 10GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 70GeV\n"
	elif(parameterName.find("W-Jet")>=0) :
	process+="set Factory:OrderInAlphaS 1\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ nu j\ndo Factory:Process p p e- nu j\n\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/HTScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	if(parameterName.find("W-Jet-1-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 100.*GeV\n"
	parameterName=parameterName.replace("W-Jet-1-e","W-Jet-e")
	elif(parameterName.find("W-Jet-2-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 190.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-2-e","W-Jet-e")
	elif(parameterName.find("W-Jet-3-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 270.0*GeV\n"
	parameterName=parameterName.replace("W-Jet-3-e","W-Jet-e")
	elif(parameterName.find("Z-Jet")>=0) :
	process+="set Factory:OrderInAlphaS 1\nset Factory:OrderInAlphaEW 2\n"
	if(parameterName.find("-e")>=0) :
	process+="do Factory:Process p p e+ e- j\n"
	if(parameterName.find("Z-Jet-0-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 35.*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-1-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 100.*GeV\n"
	parameterName=parameterName.replace("Z-Jet-1-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-2-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 190.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-2-e","Z-Jet-e")
	elif(parameterName.find("Z-Jet-3-e")>=0) :
	process+="set /Herwig/Cuts/FirstJet:PtMin 270.0*GeV\n"
	parameterName=parameterName.replace("Z-Jet-3-e","Z-Jet-e")
	else :
	process+="do Factory:Process p p mu+ mu- j\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 35.*GeV\n"
	parameterName=parameterName.replace("Z-Jet-0-mu","Z-Jet-mu")
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	process+="set Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/HTScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	elif(parameterName.find("Z-bb")>=0) :
	parameters["bscheme"]="read Matchbox/FourFlavourScheme.in"
	process+="set /Herwig/Particles/b:HardProcessMass 4.2GeV\nset /Herwig/Particles/bbar:HardProcessMass 4.2GeV\n"
	process+="set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ e- b bbar\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 91.2*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 66GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 116GeV\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/SecondJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 18.*GeV\n"
	process+="set /Herwig/Cuts/SecondJet:PtMin 15.*GeV\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	elif(parameterName.find("Z-b")>=0) :
	process+="do Factory:StartParticleGroup bjet\n"
	process+="insert Factory:ParticleGroup 0 /Herwig/Particles/b\n"
	process+="insert Factory:ParticleGroup 0 /Herwig/Particles/bbar\n"
	process+="do Factory:EndParticleGroup\n"
	process+="set Factory:OrderInAlphaS 1\nset Factory:OrderInAlphaEW 2\ndo Factory:Process p p e+ e- bjet\n"
	process+="set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 91.2*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 15.*GeV\n"
	elif(parameterName.find("W-b")>=0) :
	parameters["bscheme"]="read Matchbox/FourFlavourScheme.in"
	process += "set /Herwig/Particles/b:HardProcessMass 4.2GeV\nset /Herwig/Particles/bbar:HardProcessMass 4.2GeV\n"
	process += "set Factory:OrderInAlphaS 2\nset Factory:OrderInAlphaEW 2\n"
	process += "do Factory:Process p p e+ nu b bbar\ndo Factory:Process p p e- nu b bbar\n"
	process += "do Factory:Process p p mu+ nu b bbar\ndo Factory:Process p p mu- nu b bbar\n"
	process += "set /Herwig/MatrixElements/Matchbox/Scales/FixedScale:FixedScale 80.4*GeV\nset Factory:ScaleChoice /Herwig/MatrixElements/Matchbox/Scales/FixedScale\n"
	process+="set /Herwig/Cuts/Cuts:JetFinder /Herwig/Cuts/JetFinder\n"
	process+="insert /Herwig/Cuts/Cuts:MultiCuts 0 /Herwig/Cuts/JetCuts\n"
	process+="insert /Herwig/Cuts/JetCuts:JetRegions 0 /Herwig/Cuts/FirstJet\n"
	process+="set /Herwig/Cuts/FirstJet:PtMin 30.*GeV\n"
	process+="set /Herwig/Cuts/LeptonPairMassCut:MinMass 60GeV\nset /Herwig/Cuts/LeptonPairMassCut:MaxMass 120GeV\n"
	else :
	logging.error(" Process %s not supported for Matchbox matrix elements" % name)
	sys.exit(1)
	# Star
	elif(collider=="LHC-GammaGamma" ) :
	if(parameterName.find("-7-")>=0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0\n"
	elif(parameterName.find("-8-")>=0) :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 8000.0\n"
	else :
	process="set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0\n"
	if(simulation=="") :
	if(parameterName.find("7")>=0) :
	process += "insert SimpleQCD:MatrixElements 0 /Herwig/MatrixElements/MEgg2ff\n"
	process += "set /Herwig/MatrixElements/MEgg2ff:Process Muon\n"
	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)


	parameters['parameterFile'] = os.path.join(collider,collider+"-"+parameterName+".in")
	parameters['runname'] = name
	parameters['process'] = process

	# write the file
	if(simulation=="Matchbox" ) :
	with open(os.path.join("Rivet",name+".in") ,'w') as f:
	f.write( template.substitute(parameters))
	else :
	with open(os.path.join("Rivet",name+".in") ,'w') as f:
	f.write( template.substitute(parameters))






	diff --git a/Tests/python/merge-BFactory b/Tests/python/merge-BFactory
	--- a/Tests/python/merge-BFactory
	+++ b/Tests/python/merge-BFactory
	@@ -1,211 +1,212 @@
	#! /usr/bin/env python
	import logging
	import sys

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

	import os, yoda

	# #############################################

	def fillAbove(desthisto, sourcehistosbysqrts):
	for sqrts, h in sorted(sourcehistosbysqrts.iteritems()) :
	for i in range(0,h.numPoints) :
	if(h.points[i].x-0.01<sqrts and
	h.points[i].x+0.01>sqrts) :
	desthisto.addPoint(h.points[i])

	def merge(hpath):
	global inhistos
	global outhistos
	try:
	fillAbove(outhistos[hpath], inhistos[hpath])
	except:
	logging.error("Unknown analysis object" + hpath)
	sys.exit(1)

	def useOne(hpath, sqrts):
	global inhistos
	global outhistos
	try:
	outhistos[hpath] = inhistos[hpath][float(sqrts)]
	except:
	pass

	if __name__ == "__main__":
	import logging
	from optparse import OptionParser, OptionGroup
	parser = OptionParser(usage="%prog name")
	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=["-10.52.yoda","-10.52-sym.yoda","-10.54.yoda","-10.58.yoda","-Upsilon.yoda",
	+yodafiles=["-10.52.yoda","-10.52-sym.yoda","-10.54.yoda",
	+ "-10.58-res.yoda","-10.58.yoda","-Upsilon.yoda",
	"-Upsilon2.yoda","-Upsilon4.yoda","-10.45.yoda","-Tau.yoda"]

	## Get histos
	inhistos = {}
	outhistos={}
	weights = {}
	for f in yodafiles:
	file = args[0]+f
	if(file.find("10.45")>0) :
	sqrts=10.45
	elif(file.find("10.52")>0) :
	sqrts=10.52
	elif(file.find("10.54")>0) :
	sqrts=10.54
	elif(file.find("10.58")>0) :
	sqrts=10.58
	elif(file.find("Tau")>0) :
	sqrts=10.58
	elif(file.find("Upsilon4")>0) :
	sqrts=10.58
	elif(file.find("Upsilon2")>0) :
	sqrts=10.02
	elif(file.find("Upsilon")>0) :
	sqrts=9.46
	if not os.access(file, os.R_OK):
	if(file.find("Upsilon")>=0 or file.find("Tau")>=0) : continue
	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
	## Get histos from this YODA file
	for aopath, ao in aos.iteritems() :
	# jet shapes
	if(aopath.find("6265367")>0 or aopath.find("6895344")>0 or
	aopath.find("6181155")>0 or aopath.find("2789213")>0 or
	aopath.find("2669951")>0) :
	if not inhistos.has_key(aopath):
	inhistos[aopath] = {}
	tmpE = inhistos[aopath]
	sqrttemp=sqrts
	if(aopath.find("2669951")>0 and aopath.find("d01")>0 and
	sqrts==10.45) :
	sqrts=9.9

	if not tmpE.has_key(sqrts):
	tmpE[sqrts] = ao
	else:
	- raise Exception("A set with sqrts = %s already exists" % ( sqrts))
	+ raise Exception("A set with sqrts = %s already exists" % ( sqrts))
	sqrts=sqrttemp
	else :
	outhistos[aopath] = ao

	## Make empty output histos if needed
	for hpath,hsets in inhistos.iteritems():
	if( hpath.find("2669951")>0) :
	if(type(hsets.values()[0])==yoda.core.Scatter2D) :
	outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].path,hsets.values()[0].title)
	else :
	outhistos[hpath] = hsets.values()[0]

	# BELLE
	useOne("/BELLE_2006_S6265367/d01-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y03","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y04","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y05","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y06","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y07","10.52")
	useOne("/BELLE_2006_S6265367/d01-x01-y08","10.52")
	useOne("/BELLE_2006_S6265367/d02-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d02-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d03-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d03-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d04-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d04-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d05-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d05-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d06-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d06-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d07-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d07-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d08-x01-y01","10.52")
	useOne("/BELLE_2006_S6265367/d08-x01-y02","10.52")
	useOne("/BELLE_2006_S6265367/d09-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d09-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d10-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d10-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d11-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d11-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d12-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d12-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d13-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d13-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d14-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d14-x01-y02","10.58")
	useOne("/BELLE_2006_S6265367/d15-x01-y01","10.58")
	useOne("/BELLE_2006_S6265367/d15-x01-y02","10.58")
	# BABAR
	useOne("/BABAR_2007_S6895344/d01-x01-y01","10.54")
	useOne("/BABAR_2007_S6895344/d02-x01-y01","10.54")
	useOne("/BABAR_2007_S6895344/d03-x01-y01","10.58")
	useOne("/BABAR_2007_S6895344/d04-x01-y01","10.58")
	# BABAR
	useOne("/BABAR_2005_S6181155/d01-x01-y01","10.58")
	useOne("/BABAR_2005_S6181155/d02-x01-y01","10.58")
	useOne("/BABAR_2005_S6181155/d02-x01-y02","10.54")
	useOne("/BABAR_2005_S6181155/d03-x01-y01","10.54")
	useOne("/BABAR_2005_S6181155/d04-x01-y01","10.58")
	useOne("/BABAR_2005_S6181155/d05-x01-y01","10.58")
	useOne("/BABAR_2005_S6181155/d05-x01-y02","10.54")
	# ARGUS
	useOne("/ARGUS_1993_S2789213/d01-x01-y01","10.45")
	useOne("/ARGUS_1993_S2789213/d01-x01-y02","10.45")
	useOne("/ARGUS_1993_S2789213/d01-x01-y03","10.45")
	useOne("/ARGUS_1993_S2789213/d01-x01-y04","10.45")
	useOne("/ARGUS_1993_S2789213/d01-x01-y05","10.45")
	useOne("/ARGUS_1993_S2789213/d02-x01-y01", "9.46")
	useOne("/ARGUS_1993_S2789213/d02-x01-y02", "9.46")
	useOne("/ARGUS_1993_S2789213/d02-x01-y03", "9.46")
	useOne("/ARGUS_1993_S2789213/d02-x01-y04", "9.46")
	useOne("/ARGUS_1993_S2789213/d02-x01-y05", "9.46")
	useOne("/ARGUS_1993_S2789213/d03-x01-y01","10.58")
	useOne("/ARGUS_1993_S2789213/d03-x01-y02","10.58")
	useOne("/ARGUS_1993_S2789213/d03-x01-y03","10.58")
	useOne("/ARGUS_1993_S2789213/d03-x01-y04","10.58")
	useOne("/ARGUS_1993_S2789213/d03-x01-y05","10.58")
	useOne("/ARGUS_1993_S2789213/d04-x01-y01","10.45")
	useOne("/ARGUS_1993_S2789213/d05-x01-y01", "9.46")
	useOne("/ARGUS_1993_S2789213/d06-x01-y01","10.58")
	useOne("/ARGUS_1993_S2789213/d07-x01-y01","10.45")
	useOne("/ARGUS_1993_S2789213/d08-x01-y01", "9.46")
	useOne("/ARGUS_1993_S2789213/d09-x01-y01","10.58")
	useOne("/ARGUS_1993_S2789213/d10-x01-y01","10.45")
	useOne("/ARGUS_1993_S2789213/d11-x01-y01", "9.46")
	useOne("/ARGUS_1993_S2789213/d12-x01-y01","10.58")
	useOne("/ARGUS_1993_S2789213/d13-x01-y01","10.45")
	useOne("/ARGUS_1993_S2789213/d14-x01-y01", "9.46")
	useOne("/ARGUS_1993_S2789213/d15-x01-y01","10.58")

	useOne("/ARGUS_1993_S2669951/d02-x01-y01","10.45")
	useOne("/ARGUS_1993_S2669951/d03-x01-y01","9.46")
	useOne("/ARGUS_1993_S2669951/d04-x01-y01","10.02")
	merge("/ARGUS_1993_S2669951/d01-x01-y01")
	merge("/ARGUS_1993_S2669951/d01-x01-y02")
	merge("/ARGUS_1993_S2669951/d05-x01-y01")

	# Choose output file
	name = args[0]+".yoda"
	# output the yoda file
	yoda.writeYODA(outhistos,name)
	sys.exit(0)

File Metadata

Mime Type: text/x-diff
Expires: Tue, Nov 19, 4:35 PM (1 d, 10 h)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 3797442
Default Alt Text: (129 KB)

No OneTemporaryActions

View Options

File Metadata

Event Timeline

No OneTemporary
Actions