diff --git a/Tests/python/LowEnergy.py b/Tests/python/LowEnergy.py --- a/Tests/python/LowEnergy.py +++ b/Tests/python/LowEnergy.py @@ -1,395 +1,395 @@ #! /usr/bin/env python import yoda,os,math,subprocess,optparse from string import Template # get the path for the rivet data p = subprocess.Popen(["rivet-config", "--datadir"],stdout=subprocess.PIPE) path=p.communicate()[0].strip() #Define the arguments op = optparse.OptionParser(usage=__doc__) op.add_option("--process" , dest="processes" , default=[], action="append") op.add_option("--path" , dest="path" , default=path) op.add_option("--non-perturbative", dest="nonPerturbative" , default=False, action="store_true") op.add_option("--perturbative" , dest="perturbative" , default=False, action="store_true") op.add_option("--dest" , dest="dest" , default="Rivet") op.add_option("--list" , dest="list" , default=False, action="store_true") op.add_option("--flavour" , dest="flavour" , default="All" ) op.add_option("--plots" , dest="plot" , default=False, action="store_true") opts, args = op.parse_args() path=opts.path thresholds = [0.7,2.*.5,2.*1.87,2.*5.28] # the list of analyses and processes -analyses={ 'KK' : {} , 'pipi' : {}, 'ppbar' : {}, "3pi" : {}, "etapipi" : {}, - "etaprimepipi" : {} , "4pi" : {}, "etaPhi" : {}, "etaOmega" : {}, - "2K1pi" : {} ,"2K2pi" : {} , "4K" : {}, "6m" : {}, "omegapi" : {} , "pigamma" : {}, "etagamma" : {}, - "phipi" : {}, "omegapipi" : {} , "DD" : {} , "BB" : {}, "5pi" :{} , +analyses={ 'KK' : {} , 'PiPi' : {}, 'PPbar' : {}, "3Pi" : {}, "EtaPiPi" : {}, + "EtaprimePiPi" : {} , "4Pi" : {}, "EtaPhi" : {}, "EtaOmega" : {}, + "2K1Pi" : {} ,"2K2Pi" : {} , "4K" : {}, "6m" : {}, "OmegaPi" : {} , "PiGamma" : {}, "EtaGamma" : {}, + "PhiPi" : {}, "OmegaPiPi" : {} , "DD" : {} , "BB" : {}, "5Pi" :{} , "LL" : {} } # pi+pi- -analyses["pipi"]["KLOE_2009_I797438"] = ["d02-x01-y01"] -analyses["pipi"]["KLOE_2005_I655225"] = ["d02-x01-y01"] -analyses["pipi"]["KLOE2_2017_I1634981"] = ["d01-x01-y01"] -analyses["pipi"]["BABAR_2009_I829441"] = ["d01-x01-y01"] -analyses["pipi"]["DM1_1978_I134061"] = ["d01-x01-y01"] -analyses["pipi"]["DM2_1989_I267118"] = ["d01-x01-y01"] -analyses["pipi"]["CMD2_2007_I728302"] = ["d02-x01-y01"] -analyses["pipi"]["CMD2_2006_I728191"] = ["d03-x01-y01"] -analyses["pipi"]["BESIII_2016_I1385603"] = ["d01-x01-y01"] -analyses["pipi"]["SND_2005_I686349"] = ["d01-x01-y01"] -analyses["pipi"]["CMD_1985_I221309"] = ["d01-x01-y01","d02-x01-y01"] -analyses["pipi"]["CMD2_2002_I568807"] = ["d01-x01-y02"] -analyses["pipi"]["CMD2_1999_I498859"] = ["d01-x01-y01"] -analyses['pipi']["CLEOC_2005_I693873"] = ["d01-x01-y01"] -analyses['pipi']["ND_1991_I321108"] = ["d11-x01-y01"] -analyses['pipi']["OLYA_1984_I208231"] = ["d01-x01-y01"] +analyses["PiPi"]["KLOE_2009_I797438"] = ["d02-x01-y01"] +analyses["PiPi"]["KLOE_2005_I655225"] = ["d02-x01-y01"] +analyses["PiPi"]["KLOE2_2017_I1634981"] = ["d01-x01-y01"] +analyses["PiPi"]["BABAR_2009_I829441"] = ["d01-x01-y01"] +analyses["PiPi"]["DM1_1978_I134061"] = ["d01-x01-y01"] +analyses["PiPi"]["DM2_1989_I267118"] = ["d01-x01-y01"] +analyses["PiPi"]["CMD2_2007_I728302"] = ["d02-x01-y01"] +analyses["PiPi"]["CMD2_2006_I728191"] = ["d03-x01-y01"] +analyses["PiPi"]["BESIII_2016_I1385603"] = ["d01-x01-y01"] +analyses["PiPi"]["SND_2005_I686349"] = ["d01-x01-y01"] +analyses["PiPi"]["CMD_1985_I221309"] = ["d01-x01-y01","d02-x01-y01"] +analyses["PiPi"]["CMD2_2002_I568807"] = ["d01-x01-y02"] +analyses["PiPi"]["CMD2_1999_I498859"] = ["d01-x01-y01"] +analyses['PiPi']["CLEOC_2005_I693873"] = ["d01-x01-y01"] +analyses['PiPi']["ND_1991_I321108"] = ["d11-x01-y01"] +analyses['PiPi']["OLYA_1984_I208231"] = ["d01-x01-y01"] # K+K- and K_S^0 K_L^0 analyses['KK']["BESIII_2018_I1704558"] = ["d01-x01-y01"] analyses['KK']["BABAR_2013_I1238807"] = ["d01-x01-y01"] analyses['KK']["DM1_1981_I156053"] = ["d01-x01-y01"] analyses['KK']["DM1_1981_I156054"] = ["d01-x01-y01"] analyses['KK']["CLEOC_2005_I693873"] = ["d01-x01-y02"] analyses['KK']["BABAR_2015_I1383130"] = ["d01-x01-y04"] analyses['KK']["DM2_1988_I262690" ] = ["d01-x01-y01"] analyses['KK']["SND_2007_I755881"] = ["d01-x01-y01"] analyses['KK']["SND_2001_I533574"] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03", "d02-x01-y01","d02-x01-y02","d02-x01-y03"] analyses['KK']["SND_2006_I720035"] = ["d01-x01-y01"] analyses['KK']["BABAR_2014_I1287920"] = ["d09-x01-y01"] analyses['KK']["CMD2_2003_I601222"] = ["d01-x01-y01"] analyses['KK']["CMD3_2016_I1444990"] = ["d01-x01-y06"] analyses['KK']["CMD2_1995_I406880"] = ["d01-x01-y01","d01-x01-y02"] analyses['KK']["CMD2_1999_I502164"] = ["d01-x01-y01","d02-x01-y01", "d03-x01-y01","d04-x01-y01"] analyses['KK']["CMD2_2008_I782516"] = ["d01-x01-y01","d02-x01-y01"] analyses['KK']["ND_1991_I321108"] = ["d12-x01-y01","d13-x01-y01"] analyses['KK']["OLYA_1981_I173076"] = ["d01-x01-y01"] analyses['KK']["SND_2016_I1484677"] = ["d01-x01-y01","d02-x01-y01"] # proton-antiproton -analyses['ppbar']["BESIII_2019_I1736235"] = ["d01-x01-y01"] -analyses['ppbar']["BESIII_2019_I1718337"] = ["d01-x01-y01"] -analyses['ppbar']["BESIII_2015_I1358937"] = ["d01-x01-y05"] -analyses['ppbar']["BABAR_2013_I1217421"] = ["d01-x01-y01"] -analyses['ppbar']["SND_2014_I1321689"] = ["d01-x01-y01","d02-x01-y01"] -analyses['ppbar']["CMD3_2016_I1385598"] = ["d01-x01-y06"] -analyses['ppbar']["CLEOC_2005_I693873"] = ["d01-x01-y03"] -analyses['ppbar']["BABAR_2006_I700020"] = ["d01-x01-y01","d02-x01-y01"] -analyses['ppbar']["DM2_1983_I190558"] = ["d01-x01-y01"] -analyses["ppbar"]["DM2_1990_I297706"] = ["d01-x01-y01"] -analyses["ppbar"]["DM1_1979_I141565"] = ["d01-x01-y01"] -analyses["ppbar"]["FENICE_1998_I471263"] = ["d01-x01-y01"] -analyses["ppbar"]["FENICE_1994_I377833"] = ["d01-x01-y01"] -analyses['ppbar']["BESII_2005_I685906"] = ["d01-x01-y01"] -analyses['ppbar']["BESIII_2014_I1286898"] = ["d01-x01-y06"] +analyses['PPbar']["BESIII_2019_I1736235"] = ["d01-x01-y01"] +analyses['PPbar']["BESIII_2019_I1718337"] = ["d01-x01-y01"] +analyses['PPbar']["BESIII_2015_I1358937"] = ["d01-x01-y05"] +analyses['PPbar']["BABAR_2013_I1217421"] = ["d01-x01-y01"] +analyses['PPbar']["SND_2014_I1321689"] = ["d01-x01-y01","d02-x01-y01"] +analyses['PPbar']["CMD3_2016_I1385598"] = ["d01-x01-y06"] +analyses['PPbar']["CLEOC_2005_I693873"] = ["d01-x01-y03"] +analyses['PPbar']["BABAR_2006_I700020"] = ["d01-x01-y01","d02-x01-y01"] +analyses['PPbar']["DM2_1983_I190558"] = ["d01-x01-y01"] +analyses["PPbar"]["DM2_1990_I297706"] = ["d01-x01-y01"] +analyses["PPbar"]["DM1_1979_I141565"] = ["d01-x01-y01"] +analyses["PPbar"]["FENICE_1998_I471263"] = ["d01-x01-y01"] +analyses["PPbar"]["FENICE_1994_I377833"] = ["d01-x01-y01"] +analyses['PPbar']["BESII_2005_I685906"] = ["d01-x01-y01"] +analyses['PPbar']["BESIII_2014_I1286898"] = ["d01-x01-y06"] # pi0 gamma -analyses["pigamma"]["SND_2018_I1694988"] = ["d01-x01-y01"] -analyses["pigamma"]["SND_2016_I1418483"] = ["d01-x01-y05"] -analyses["pigamma"]["SND_2003_I612867"] = ["d01-x01-y01"] -analyses["pigamma"]["CMD2_2005_I658856"] = ["d02-x01-y01"] -analyses["pigamma"]["SND_2000_I524221"] = ["d01-x01-y02"] -# eta gamma -analyses["etagamma"]["CMD2_2005_I658856"] = ["d01-x01-y01"] -analyses["etagamma"]["SND_2006_I717778" ] = ["d01-x01-y01","d02-x01-y01"] -analyses["etagamma"]["SND_2014_I1275333"] = ["d01-x01-y01"] -analyses["etagamma"]["SND_2000_I524221"] = ["d01-x01-y01"] -analyses["etagamma"]["CMD2_1999_I503154"] = ["d01-x01-y01"] -analyses["etagamma"]["CMD2_2001_I554522"] = ["d01-x01-y01"] -analyses['etagamma']["CMD2_1995_I406880"] = ["d01-x01-y04"] -analyses['etagamma']["BABAR_2006_I716277"] = ["d01-x01-y01"] +analyses["PiGamma"]["SND_2018_I1694988"] = ["d01-x01-y01"] +analyses["PiGamma"]["SND_2016_I1418483"] = ["d01-x01-y05"] +analyses["PiGamma"]["SND_2003_I612867"] = ["d01-x01-y01"] +analyses["PiGamma"]["CMD2_2005_I658856"] = ["d02-x01-y01"] +analyses["PiGamma"]["SND_2000_I524221"] = ["d01-x01-y02"] +# Eta gamma +analyses["EtaGamma"]["CMD2_2005_I658856"] = ["d01-x01-y01"] +analyses["EtaGamma"]["SND_2006_I717778" ] = ["d01-x01-y01","d02-x01-y01"] +analyses["EtaGamma"]["SND_2014_I1275333"] = ["d01-x01-y01"] +analyses["EtaGamma"]["SND_2000_I524221"] = ["d01-x01-y01"] +analyses["EtaGamma"]["CMD2_1999_I503154"] = ["d01-x01-y01"] +analyses["EtaGamma"]["CMD2_2001_I554522"] = ["d01-x01-y01"] +analyses['EtaGamma']["CMD2_1995_I406880"] = ["d01-x01-y04"] +analyses['EtaGamma']["BABAR_2006_I716277"] = ["d01-x01-y01"] # 3 meson -analyses["3pi"]["BABAR_2004_I656680"] = ["d01-x01-y01"] -analyses["3pi"]["SND_2002_I582183"] = ["d01-x01-y01"] -analyses["3pi"]["SND_2003_I619011"] = ["d01-x01-y01"] -analyses["3pi"]["SND_1999_I508003"] = ["d01-x01-y01"] -analyses["3pi"]["SND_2001_I533574"] = ["d01-x01-y04","d02-x01-y04"] -analyses["3pi"]["CMD2_2000_I523691"] = ["d01-x01-y01"] -analyses["3pi"]["CMD2_1998_I480170"] = ["d01-x01-y01"] -analyses['3pi']["CMD2_1995_I406880"] = ["d01-x01-y03"] -analyses['3pi']["DM2_1992_I339265" ] = ["d01-x01-y01"] -analyses['3pi']["DM1_1980_I140174" ] = ["d01-x01-y01"] -analyses['3pi']["ND_1991_I321108"] = ["d05-x01-y01","d10-x01-y04"] -analyses['3pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y01"] -analyses["3pi"]["CLEO_2006_I691720"] = ["d01-x01-y01"] -analyses["3pi"]["SND_2015_I1389908"] = ["d01-x01-y01"] -# etapipi -analyses["etapipi"]["BABAR_2018_I1700745"] = ["d02-x01-y01"] -analyses["etapipi"]["BABAR_2018_I1647139"] = ["d01-x01-y01"] -analyses["etapipi"]["SND_2015_I1332929"] = ["d01-x01-y01"] -analyses["etapipi"]["SND_2018_I1638368"] = ["d01-x01-y01"] -analyses["etapipi"]["BABAR_2007_I758568"] = ["d01-x01-y01","d02-x01-y01"] -analyses["etapipi"]["CMD2_2000_I532970"] = ["d02-x01-y01"] -analyses["etapipi"]["DM2_1988_I264144"] = ["d01-x01-y01"] -analyses['etapipi']["ND_1991_I321108"] = ["d06-x01-y01","d14-x01-y01"] -analyses['etapipi']["CMD3_2019_I1744510"] = ["d02-x01-y01"] -# eta Phi -analyses["etaPhi"]["BABAR_2008_I765258"] = ["d04-x01-y01","d05-x01-y01"] -analyses["etaPhi"]["SND_2018_I1693737"] = ["d01-x01-y01"] -analyses["etaPhi"]["BABAR_2017_I1511276"] = ["d03-x01-y01"] -analyses["etaPhi"]["SND_2019_I1726419"] = ["d01-x01-y03"] -analyses["etaPhi"]["CMD3_2019_I1740541"] = ["d01-x01-y06","d02-x01-y06","d03-x01-y06"] -# eta Omega -analyses["etaOmega"]["SND_2016_I1473343"] = ["d01-x01-y01"] -analyses["etaOmega"]["BABAR_2006_I709730"] = ["d02-x01-y01"] -analyses["etaOmega"]["SND_2019_I1726419"] = ["d01-x01-y01","d01-x01-y02"] -analyses["etaOmega"]["CMD3_2017_I1606078"] = ["d01-x01-y01","d01-x01-y02"] +analyses["3Pi"]["BABAR_2004_I656680"] = ["d01-x01-y01"] +analyses["3Pi"]["SND_2002_I582183"] = ["d01-x01-y01"] +analyses["3Pi"]["SND_2003_I619011"] = ["d01-x01-y01"] +analyses["3Pi"]["SND_1999_I508003"] = ["d01-x01-y01"] +analyses["3Pi"]["SND_2001_I533574"] = ["d01-x01-y04","d02-x01-y04"] +analyses["3Pi"]["CMD2_2000_I523691"] = ["d01-x01-y01"] +analyses["3Pi"]["CMD2_1998_I480170"] = ["d01-x01-y01"] +analyses['3Pi']["CMD2_1995_I406880"] = ["d01-x01-y03"] +analyses['3Pi']["DM2_1992_I339265" ] = ["d01-x01-y01"] +analyses['3Pi']["DM1_1980_I140174" ] = ["d01-x01-y01"] +analyses['3Pi']["ND_1991_I321108"] = ["d05-x01-y01","d10-x01-y04"] +analyses['3Pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y01"] +analyses["3Pi"]["CLEO_2006_I691720"] = ["d01-x01-y01"] +analyses["3Pi"]["SND_2015_I1389908"] = ["d01-x01-y01"] +# EtaPiPi +analyses["EtaPiPi"]["BABAR_2018_I1700745"] = ["d02-x01-y01"] +analyses["EtaPiPi"]["BABAR_2018_I1647139"] = ["d01-x01-y01"] +analyses["EtaPiPi"]["SND_2015_I1332929"] = ["d01-x01-y01"] +analyses["EtaPiPi"]["SND_2018_I1638368"] = ["d01-x01-y01"] +analyses["EtaPiPi"]["BABAR_2007_I758568"] = ["d01-x01-y01","d02-x01-y01"] +analyses["EtaPiPi"]["CMD2_2000_I532970"] = ["d02-x01-y01"] +analyses["EtaPiPi"]["DM2_1988_I264144"] = ["d01-x01-y01"] +analyses['EtaPiPi']["ND_1991_I321108"] = ["d06-x01-y01","d14-x01-y01"] +analyses['EtaPiPi']["CMD3_2019_I1744510"] = ["d02-x01-y01"] +# Eta Phi +analyses["EtaPhi"]["BABAR_2008_I765258"] = ["d04-x01-y01","d05-x01-y01"] +analyses["EtaPhi"]["SND_2018_I1693737"] = ["d01-x01-y01"] +analyses["EtaPhi"]["BABAR_2017_I1511276"] = ["d03-x01-y01"] +analyses["EtaPhi"]["SND_2019_I1726419"] = ["d01-x01-y03"] +analyses["EtaPhi"]["CMD3_2019_I1740541"] = ["d01-x01-y06","d02-x01-y06","d03-x01-y06"] +# Eta Omega +analyses["EtaOmega"]["SND_2016_I1473343"] = ["d01-x01-y01"] +analyses["EtaOmega"]["BABAR_2006_I709730"] = ["d02-x01-y01"] +analyses["EtaOmega"]["SND_2019_I1726419"] = ["d01-x01-y01","d01-x01-y02"] +analyses["EtaOmega"]["CMD3_2017_I1606078"] = ["d01-x01-y01","d01-x01-y02"] # 4 pions -analyses["4pi"]["BABAR_2017_I1621593"] = ["d01-x01-y01","d02-x01-y01"] -analyses["4pi"]["BABAR_2012_I1086164"] = ["d01-x01-y01"] -analyses["4pi"]["CMD2_2000_I531667"] = ["d01-x01-y01"] -analyses["4pi"]["CMD2_2004_I648023"] = ["d01-x01-y01"] -analyses["4pi"]["BABAR_2005_I676691"] = ["d01-x01-y01"] -analyses["4pi"]["CMD2_2000_I511375"] = ["d01-x01-y01"] -analyses["4pi"]["CMD2_1999_I483994"] = ["d01-x01-y01","d02-x01-y01"] -analyses["4pi"]["BESII_2008_I801210"] = ["d01-x01-y01"] -analyses["4pi"]["KLOE_2008_I791841"] = ["d01-x01-y01"] -analyses['4pi']["ND_1991_I321108"] = ["d07-x01-y01","d08-x01-y01","d10-x01-y01","d10-x01-y02"] -analyses['4pi']["BESII_2007_I750713"] = ["d01-x01-y03"] -analyses['4pi']["SND_2001_I579319"] = ["d01-x01-y01","d02-x01-y01"] -analyses['4pi']["DM1_1982_I168552"] = ["d01-x01-y01"] -analyses['4pi']["DM1_1979_I132828"] = ["d01-x01-y01"] -analyses['4pi']["GAMMAGAMMA_1980_I153382"] = ["d01-x01-y01"] -analyses['4pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y02"] -# (these are omega(-> pi0 gamma) pi0) -analyses["omegapi"]["SND_2016_I1489182" ] = ["d01-x01-y01"] -analyses["omegapi"]["SND_2000_I527752" ] = ["d01-x01-y01"] -analyses["omegapi"]["SND_2000_I503946" ] = ["d01-x01-y01"] -analyses["omegapi"]["CMD2_2003_I616446" ] = ["d01-x01-y01"] -# non omega -analyses["omegapi"]["SND_2002_I587084" ] = ["d01-x01-y01"] -analyses["omegapi"]["CMD2_2004_I630009" ] = ["d01-x01-y01"] -analyses["omegapi"]["KLOE_2008_I791841"] = ["d02-x01-y01"] -# from 4pion -analyses["omegapi"]["CMD2_1999_I483994" ] = ["d03-x01-y01"] -analyses['omegapi']["ND_1991_I321108"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", +analyses["4Pi"]["BABAR_2017_I1621593"] = ["d01-x01-y01","d02-x01-y01"] +analyses["4Pi"]["BABAR_2012_I1086164"] = ["d01-x01-y01"] +analyses["4Pi"]["CMD2_2000_I531667"] = ["d01-x01-y01"] +analyses["4Pi"]["CMD2_2004_I648023"] = ["d01-x01-y01"] +analyses["4Pi"]["BABAR_2005_I676691"] = ["d01-x01-y01"] +analyses["4Pi"]["CMD2_2000_I511375"] = ["d01-x01-y01"] +analyses["4Pi"]["CMD2_1999_I483994"] = ["d01-x01-y01","d02-x01-y01"] +analyses["4Pi"]["BESII_2008_I801210"] = ["d01-x01-y01"] +analyses["4Pi"]["KLOE_2008_I791841"] = ["d01-x01-y01"] +analyses['4Pi']["ND_1991_I321108"] = ["d07-x01-y01","d08-x01-y01","d10-x01-y01","d10-x01-y02"] +analyses['4Pi']["BESII_2007_I750713"] = ["d01-x01-y03"] +analyses['4Pi']["SND_2001_I579319"] = ["d01-x01-y01","d02-x01-y01"] +analyses['4Pi']["DM1_1982_I168552"] = ["d01-x01-y01"] +analyses['4Pi']["DM1_1979_I132828"] = ["d01-x01-y01"] +analyses['4Pi']["GAMMAGAMMA_1980_I153382"] = ["d01-x01-y01"] +analyses['4Pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y02"] +# (these are Omega(-> pi0 gamma) pi0) +analyses["OmegaPi"]["SND_2016_I1489182" ] = ["d01-x01-y01"] +analyses["OmegaPi"]["SND_2000_I527752" ] = ["d01-x01-y01"] +analyses["OmegaPi"]["SND_2000_I503946" ] = ["d01-x01-y01"] +analyses["OmegaPi"]["CMD2_2003_I616446" ] = ["d01-x01-y01"] +# non Omega +analyses["OmegaPi"]["SND_2002_I587084" ] = ["d01-x01-y01"] +analyses["OmegaPi"]["CMD2_2004_I630009" ] = ["d01-x01-y01"] +analyses["OmegaPi"]["KLOE_2008_I791841"] = ["d02-x01-y01"] +# from 4 Pion +analyses["OmegaPi"]["CMD2_1999_I483994" ] = ["d03-x01-y01"] +analyses['OmegaPi']["ND_1991_I321108"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", "d04-x01-y01","d10-x01-y03"] # 5 pion and related -analyses["omegapipi"]["DM1_1981_I166964" ] = ["d01-x01-y01"] -analyses["omegapipi"]["DM2_1992_I339265" ] = ["d02-x01-y01"] -analyses["omegapipi"]["CMD2_2000_I532970" ] = ["d01-x01-y01"] -analyses["omegapipi"]["BABAR_2018_I1700745"] = ["d03-x01-y01"] -analyses["omegapipi"]["BABAR_2007_I758568" ] = ["d03-x01-y01","d04-x01-y01"] -analyses["5pi"]["CMD2_2000_I532970"] = ["d03-x01-y01"] -analyses["5pi"]["BABAR_2007_I758568"] = ["d01-x01-y01"] -analyses['5pi']["ND_1991_I321108"] = ["d14-x01-y01"] -analyses['5pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y03"] -analyses["5pi"]["BABAR_2018_I1700745"] = ["d01-x01-y01"] +analyses["OmegaPiPi"]["DM1_1981_I166964" ] = ["d01-x01-y01"] +analyses["OmegaPiPi"]["DM2_1992_I339265" ] = ["d02-x01-y01"] +analyses["OmegaPiPi"]["CMD2_2000_I532970" ] = ["d01-x01-y01"] +analyses["OmegaPiPi"]["BABAR_2018_I1700745"] = ["d03-x01-y01"] +analyses["OmegaPiPi"]["BABAR_2007_I758568" ] = ["d03-x01-y01","d04-x01-y01"] +analyses["5Pi"]["CMD2_2000_I532970"] = ["d03-x01-y01"] +analyses["5Pi"]["BABAR_2007_I758568"] = ["d01-x01-y01"] +analyses['5Pi']["ND_1991_I321108"] = ["d14-x01-y01"] +analyses['5Pi']["GAMMAGAMMA_1981_I158474"] = ["d01-x01-y03"] +analyses["5Pi"]["BABAR_2018_I1700745"] = ["d01-x01-y01"] # 2K 1 pi -analyses["2K1pi"]["BABAR_2008_I765258"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] -analyses["2K1pi"]["DM1_1982_I176801"] = ["d01-x01-y01"] -analyses["2K1pi"]["DM2_1991_I318558"] = ["d01-x01-y01","d02-x01-y01"] -analyses["2K1pi"]["BESII_2008_I801208"] = ["d01-x01-y01"] -analyses["2K1pi"]["SND_2018_I1637194"] = ["d01-x01-y01"] -analyses["2K1pi"]["BESIII_2018_I1691798"] = ["d01-x01-y01"] -analyses["2K1pi"]["BABAR_2017_I1511276"] = ["d01-x01-y01"] -analyses["phipi"]["BABAR_2017_I1511276"] = ["d01-x01-y01","d02-x01-y01"] -analyses["phipi"]["BABAR_2008_I765258"] = ["d02-x01-y01","d03-x01-y01"] +analyses["2K1Pi"]["BABAR_2008_I765258"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] +analyses["2K1Pi"]["DM1_1982_I176801"] = ["d01-x01-y01"] +analyses["2K1Pi"]["DM2_1991_I318558"] = ["d01-x01-y01","d02-x01-y01"] +analyses["2K1Pi"]["BESII_2008_I801208"] = ["d01-x01-y01"] +analyses["2K1Pi"]["SND_2018_I1637194"] = ["d01-x01-y01"] +analyses["2K1Pi"]["BESIII_2018_I1691798"] = ["d01-x01-y01"] +analyses["2K1Pi"]["BABAR_2017_I1511276"] = ["d01-x01-y01"] +analyses["PhiPi"]["BABAR_2017_I1511276"] = ["d01-x01-y01","d02-x01-y01"] +analyses["PhiPi"]["BABAR_2008_I765258"] = ["d02-x01-y01","d03-x01-y01"] # 2K 2 pi -analyses["2K2pi"]["DM1_1982_I169382"] = ["d01-x01-y01"] -analyses["2K2pi"]["BABAR_2005_I676691"] = ["d02-x01-y01"] -analyses["2K2pi"]["BABAR_2014_I1287920"] = ["d09-x01-y01","d10-x01-y01","d11-x01-y01"] -analyses["2K2pi"]["BABAR_2012_I892684"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", +analyses["2K2Pi"]["DM1_1982_I169382"] = ["d01-x01-y01"] +analyses["2K2Pi"]["BABAR_2005_I676691"] = ["d02-x01-y01"] +analyses["2K2Pi"]["BABAR_2014_I1287920"] = ["d09-x01-y01","d10-x01-y01","d11-x01-y01"] +analyses["2K2Pi"]["BABAR_2012_I892684"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", "d04-x01-y01","d05-x01-y01", "d06-x01-y01","d07-x01-y01"] -analyses["2K2pi"]["BABAR_2007_I747875"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", +analyses["2K2Pi"]["BABAR_2007_I747875"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", "d04-x01-y01","d05-x01-y01","d07-x01-y01"] -analyses["2K2pi"]["BESII_2008_I801210"] = ["d01-x01-y02"] -analyses["2K2pi"]["BESII_2008_I801208"] = ["d01-x01-y02"] -analyses["2K2pi"]["BELLE_2009_I809630"] = ["d01-x01-y01"] -analyses["2K2pi"]["CMD3_2016_I1395968"] = ["d01-x01-y01"] -analyses['2K2pi']["BESII_2007_I750713"] = ["d01-x01-y04"] -analyses["2K2pi"]["BABAR_2017_I1511276"] = ["d03-x01-y01","d04-x01-y01"] -analyses["2K2pi"]["BABAR_2017_I1591716"] = ["d01-x01-y01","d02-x01-y01"] -analyses['2K2pi']["BESIII_2018_I1699641"] = ["d01-x01-y01","d02-x01-y01"] +analyses["2K2Pi"]["BESII_2008_I801210"] = ["d01-x01-y02"] +analyses["2K2Pi"]["BESII_2008_I801208"] = ["d01-x01-y02"] +analyses["2K2Pi"]["BELLE_2009_I809630"] = ["d01-x01-y01"] +analyses["2K2Pi"]["CMD3_2016_I1395968"] = ["d01-x01-y01"] +analyses['2K2Pi']["BESII_2007_I750713"] = ["d01-x01-y04"] +analyses["2K2Pi"]["BABAR_2017_I1511276"] = ["d03-x01-y01","d04-x01-y01"] +analyses["2K2Pi"]["BABAR_2017_I1591716"] = ["d01-x01-y01","d02-x01-y01"] +analyses['2K2Pi']["BESIII_2018_I1699641"] = ["d01-x01-y01","d02-x01-y01"] # 4K analyses["4K"]["BESIII_2019_I1743841"] = ["d01-x01-y01","d02-x01-y01"] analyses["4K"]["BABAR_2005_I676691"] = ["d03-x01-y01"] analyses["4K"]["BABAR_2014_I1287920"] = ["d12-x01-y01"] analyses["4K"]["BABAR_2012_I892684"] = ["d08-x01-y01"] analyses["4K"]["BABAR_2007_I747875"] = ["d07-x01-y01"] analyses['4K']["BESII_2007_I750713"] = ["d01-x01-y06","d01-x01-y07"] # 6 mesons analyses["6m"]["CMD3_2013_I1217420"] = ["d01-x01-y01"] analyses["6m"]["SND_2019_I1726419"] = ["d01-x01-y01","d01-x01-y04"] analyses["6m"]["CMD3_2017_I1606078"] = ["d01-x01-y03","d01-x01-y04"] analyses["6m"]["CMD3_2019_I1720610"] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03"] analyses["6m"]["BABAR_2018_I1700745"] = ["d04-x01-y01","d05-x01-y01"] analyses["6m"]["SND_2016_I1471515"] = ["d01-x01-y06"] analyses["6m"]["DM1_1981_I166353"] = ["d01-x01-y01"] analyses["6m"]["BABAR_2006_I709730"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] analyses["6m"]["BABAR_2007_I758568"] = ["d05-x01-y01","d07-x01-y01", "d08-x01-y01","d09-x01-y01","d10-x01-y01","d11-x01-y01"] -analyses["etaprimepipi"]["BABAR_2007_I758568"] = ["d05-x01-y01","d06-x01-y01"] +analyses["EtaprimePiPi"]["BABAR_2007_I758568"] = ["d05-x01-y01","d06-x01-y01"] analyses["6m"]["BESII_2007_I763880"] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03","d01-x01-y04", "d01-x01-y05","d01-x01-y06","d01-x01-y07"] analyses["6m"]["BESII_2007_I762901"] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03","d01-x01-y04", "d01-x01-y06","d01-x01-y07","d01-x01-y08","d01-x01-y09","d01-x01-y10"] analyses["6m"]["CLEO_2006_I691720"] = ["d01-x01-y02","d01-x01-y03","d01-x01-y04","d01-x01-y05", "d01-x01-y07","d01-x01-y08","d01-x01-y09","d01-x01-y10","d01-x01-y11", "d01-x01-y12","d01-x01-y13","d01-x01-y14","d01-x01-y15","d01-x01-y17"] analyses["6m"]["BESII_2008_I801210"] = ["d01-x01-y03","d01-x01-y04","d01-x01-y05"] analyses["6m"]["BESII_2008_I801208"] = ["d01-x01-y03","d01-x01-y04","d01-x01-y05","d01-x01-y06"] analyses["6m"]["JADE_1979_I142874"] = ["d02-x01-y01"] analyses["6m"]["MARKI_1982_I169326"] = ["d06-x01-y01"] analyses["6m"]["MARKI_1975_I100592"] = ["d01-x01-y01","d02-x01-y01"] analyses['6m']["BESII_2007_I750713"] = ["d01-x01-y08","d01-x01-y09","d01-x01-y11", "d01-x01-y12","d01-x01-y13","d01-x01-y14", "d01-x01-y15","d01-x01-y16","d01-x01-y17","d01-x01-y18"] analyses['6m']["SND_2016_I1473343"] = ["d01-x01-y01"] # DD analyses["DD"]["BELLE_2007_I723333" ] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["BELLE_2007_I756012" ] = ["d01-x01-y01"] analyses["DD"]["BELLE_2007_I756643" ] = ["d01-x01-y01"] analyses["DD"]["BELLE_2008_I757220" ] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["BELLE_2008_I759073" ] = ["d01-x01-y01"] analyses["DD"]["BABAR_2008_I776519" ] = ["d01-x01-y01","d01-x01-y02"] analyses["DD"]["BELLE_2008_I791660" ] = ["d01-x01-y01"] analyses["DD"]["BELLE_2013_I1225975"] = ["d01-x01-y01"] analyses["DD"]["BELLE_2014_I1282602"] = ["d01-x01-y01"] analyses["DD"]["BELLE_2015_I1324785"] = ["d01-x01-y01"] analyses["DD"]["BESIII_2016_I1457597"] = ["d01-x01-y07"] analyses["DD"]["BESIII_2015_I1355215"] = ["d01-x01-y10"] analyses["DD"]["BESIII_2015_I1377204"] = ["d01-x01-y10"] analyses["DD"]["BESIII_2016_I1495838"] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["CRYSTAL_BALL_1986_I238081"] = ["d02-x01-y01"] analyses["DD"]["CLEOC_2008_I777917"] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03", "d02-x01-y01","d02-x01-y02","d02-x01-y03", "d03-x01-y01","d03-x01-y02","d03-x01-y03", "d04-x01-y01","d04-x01-y02", "d05-x01-y01","d05-x01-y02"] analyses["DD"]["BELLE_2017_I1613517"] = ["d01-x01-y01","d01-x01-y02"] analyses["DD"]["BESIII_2014_I1323621"] = ["d01-x01-y01"] analyses["DD"]["BESIII_2015_I1406939"] = ["d02-x01-y06","d03-x01-y06"] analyses["DD"]["BESIII_2017_I1628093"] = ["d01-x01-y01"] analyses["DD"]["BESIII_2019_I1723934"] = ["d01-x01-y01"] analyses["DD"]["BESIII_2019_I1756876"] = ["d01-x01-y09","d01-x01-y10"] # BB analyses["BB"]["BELLE_2016_I1389855"] = ["d01-x01-y02","d01-x01-y03"] analyses["BB"]["BELLE_2008_I764099"] = ["d01-x01-y01","d02-x01-y01", "d03-x01-y01","d04-x01-y01"] analyses["BB"]["CLEO_1999_I474676"] = ["d01-x01-y01","d01-x01-y02"] analyses["BB"]["CUSB_1991_I325661"] = ["d01-x01-y01"] analyses["BB"]["BABAR_2001_I558091"] = ["d01-x01-y01"] analyses["BB"]["CLEO_1991_I29927"] = ["d01-x01-y01"] analyses["LL"]["BESIII_2018_I1627871"] = ["d01-x01-y01"] analyses["LL"]["DM2_1990_I297706"] = ["d02-x01-y01"] analyses["LL"]["BESIII_2019_I1758883"] = ["d01-x01-y05"] analyses["LL"]["BABAR_2007_I760730"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] # list the analysis if required and quit() if "all" in opts.processes : processes = sorted(list(analyses.keys())) else : processes = sorted(list(set(opts.processes))) if(opts.list) : for process in processes : print " ".join(analyses[process]) quit() if(opts.plot) : output="" for process in processes: for analysis in analyses[process] : for plot in analyses[process][analysis]: output+= " -m/%s/%s" % (analysis,plot) print output quit() # mapping of process to me to use -me = { "pipi" : "MEee2Pions", +me = { "PiPi" : "MEee2Pions", "KK" : "MEee2Kaons", - "3pi" : "MEee3Pions", - "4pi" : "MEee4Pions", - "etapipi" : "MEee2EtaPiPi", - "etaprimepipi" : "MEee2EtaPrimePiPi", - "etaPhi" : "MEee2EtaPhi", - "etaOmega" : "MEee2EtaOmega", - "omegapi" : "MEee2OmegaPi", - "omegapipi" : "MEee2OmegaPiPi", - "phipi" : "MEee2PhiPi", - "pigamma" : "MEee2PiGamma", - "etagamma" : "MEee2EtaGamma", - "ppbar" : "MEee2ppbar", + "3Pi" : "MEee3Pions", + "4Pi" : "MEee4Pions", + "EtaPiPi" : "MEee2EtaPiPi", + "EtaprimePiPi" : "MEee2EtaPrimePiPi", + "EtaPhi" : "MEee2EtaPhi", + "EtaOmega" : "MEee2EtaOmega", + "OmegaPi" : "MEee2OmegaPi", + "OmegaPiPi" : "MEee2OmegaPiPi", + "PhiPi" : "MEee2PhiPi", + "PiGamma" : "MEee2PiGamma", + "EtaGamma" : "MEee2EtaGamma", + "PPbar" : "MEee2PPbar", "LL" : "MEee2LL" , - "2K1pi" : "MEee2KKPi" } + "2K1Pi" : "MEee2KKPi" } # energies we need energies={} def nearestEnergy(en) : Emin=0 delta=1e30 anals=[] for val in energies : if(abs(val-en)200) : energy *= 0.001 emin,delta,anals = nearestEnergy(energy) if(energy in energies) : if(analysis not in energies[energy][1]) : energies[energy][1].append(analysis) if(matrix!="" and matrix not in energies[energy][0]) : energies[energy][0].append(matrix) elif(delta<1e-7) : if(analysis not in anals[1]) : anals[1].append(analysis) if(matrix!="" and matrix not in anals[0]) : anals[0].append(matrix) else : if(matrix=="") : energies[energy]=[[],[analysis]] else : energies[energy]=[[matrix],[analysis]] with open("python/LowEnergy-EE-Perturbative.in", 'r') as f: templateText = f.read() perturbative=Template( templateText ) with open("python/LowEnergy-EE-NonPerturbative.in", 'r') as f: templateText = f.read() nonPerturbative=Template( templateText ) targets="" for energy in sorted(energies) : anal="" for analysis in energies[energy][1]: anal+= "insert /Herwig/Analysis/Rivet:Analyses 0 %s\n" %analysis proc="" for matrix in energies[energy][0] : proc+="insert SubProcess:MatrixElements 0 %s\n" % matrix proc+="set %s:Flavour %s\n" % (matrix,opts.flavour) maxflavour =5 if(energy thresholds[0]) : inputPerturbative = perturbative.substitute({"ECMS" : "%8.6f" % energy, "ANALYSES" : anal, "lepton" : "", "maxflavour" : maxflavour}) with open(opts.dest+"/Rivet-LowEnergy-EE-Perturbative-%8.6f.in" % energy ,'w') as f: f.write(inputPerturbative) targets += "Rivet-LowEnergy-EE-Perturbative-%8.6f.yoda " % energy # input file for currents if(opts.nonPerturbative and proc!="") : inputNonPerturbative = nonPerturbative.substitute({"ECMS" : "%8.6f" % energy, "ANALYSES" : anal, "processes" : proc}) with open(opts.dest+"/Rivet-LowEnergy-EE-NonPerturbative-%8.6f.in" % energy ,'w') as f: f.write(inputNonPerturbative) targets += "Rivet-LowEnergy-EE-NonPerturbative-%8.6f.yoda " % energy print targets diff --git a/src/defaults/MatrixElements.in b/src/defaults/MatrixElements.in --- a/src/defaults/MatrixElements.in +++ b/src/defaults/MatrixElements.in @@ -1,323 +1,323 @@ # -*- ThePEG-repository -*- ############################################################################## # Setup of default matrix elements. # # Only one ME is activated by default, but this file lists # some alternatives. All available MEs can be found in the # 'include/Herwig/MatrixElements' subdirectory of your Herwig # installation. # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # Instead of editing this file directly, you should reset # the matrix elements in your own input files: # # - create your custom SubProcessHandler # - insert the MEs you need # - set your SubProcessHandler instead of the default (see HerwigDefaults.in) ############################################################################## mkdir /Herwig/MatrixElements cd /Herwig/MatrixElements library HwMELepton.so library HwMEHadron.so library HwMEDIS.so ############################################################ # e+e- matrix elements ############################################################ # e+e- > q qbar create Herwig::MEee2gZ2qq MEee2gZ2qq newdef MEee2gZ2qq:MinimumFlavour 1 newdef MEee2gZ2qq:MaximumFlavour 5 newdef MEee2gZ2qq:AlphaQCD /Herwig/Shower/AlphaQCD newdef MEee2gZ2qq:AlphaQED /Herwig/Shower/AlphaQED # e+e- -> l+l- create Herwig::MEee2gZ2ll MEee2gZ2ll newdef MEee2gZ2ll:Allowed Charged set MEee2gZ2ll:AlphaQED /Herwig/Shower/AlphaQED # e+e- -> W+W- ZZ create Herwig::MEee2VV MEee2VV # e+e- -> ZH create Herwig::MEee2ZH MEee2ZH newdef MEee2ZH:Coupling /Herwig/Shower/AlphaQCD # e+e- -> e+e-H/nu_enu_ebarH create Herwig::MEee2HiggsVBF MEee2HiggsVBF ############################################################ # Low energy matrix elements ############################################################ # e+ e- -> pi+pi- create Herwig::MEee2Mesons MEee2Pions HwMELeptonLowEnergy.so create Herwig::TwoPionCzyzCurrent /Herwig/Decays/TwoPionCzyzCurrent HwWeakCurrents.so set MEee2Pions:WeakCurrent /Herwig/Decays/TwoPionCzyzCurrent # e+ e- -> K+K-/ K0K0 create Herwig::MEee2Mesons MEee2Kaons HwMELeptonLowEnergy.so create Herwig::TwoKaonCzyzCurrent /Herwig/Decays/TwoKaonCzyzCurrent set MEee2Kaons:WeakCurrent /Herwig/Decays/TwoKaonCzyzCurrent # e+ e- -> pi+ pi- pi0 create Herwig::MEee2Mesons MEee3Pions HwMELeptonLowEnergy.so create Herwig::ThreePionCzyzCurrent /Herwig/Decays/ThreePionCzyzCurrent set MEee3Pions:WeakCurrent /Herwig/Decays/ThreePionCzyzCurrent # e+ e- -> 2pi+ 2pi-, 2pi0, pi+ pi- create Herwig::MEee2Mesons MEee4Pions HwMELeptonLowEnergy.so create Herwig::FourPionCzyzCurrent /Herwig/Decays/FourPionCzyzCurrent set MEee4Pions:WeakCurrent /Herwig/Decays/FourPionCzyzCurrent # e+ e- -> eta pi+ pi- create Herwig::MEee2Mesons MEee2EtaPiPi HwMELeptonLowEnergy.so create Herwig::EtaPiPiCurrent /Herwig/Decays/EtaPiPiCurrent set MEee2EtaPiPi:WeakCurrent /Herwig/Decays/EtaPiPiCurrent # e+ e- -> eta' pi+ pi- create Herwig::MEee2Mesons MEee2EtaPrimePiPi HwMELeptonLowEnergy.so create Herwig::EtaPrimePiPiCurrent /Herwig/Decays/EtaPrimePiPiCurrent set MEee2EtaPrimePiPi:WeakCurrent /Herwig/Decays/EtaPrimePiPiCurrent # e+ e- -> omega pi (omega -> pi0 gamma) create Herwig::MEee2Mesons MEee2OmegaPi HwMELeptonLowEnergy.so create Herwig::TwoPionPhotonSNDCurrent /Herwig/Decays/OmegaPiCurrent set MEee2OmegaPi:WeakCurrent /Herwig/Decays/OmegaPiCurrent # e+ e- > pi0 gamma create Herwig::MEee2Mesons MEee2PiGamma HwMELeptonLowEnergy.so create Herwig::PionPhotonCurrent /Herwig/Decays/PiGammaCurrent set MEee2PiGamma:WeakCurrent /Herwig/Decays/PiGammaCurrent # e+e- -> eta gamma create Herwig::MEee2Mesons MEee2EtaGamma HwMELeptonLowEnergy.so create Herwig::EtaPhotonCurrent /Herwig/Decays/EtaGammaCurrent set MEee2EtaGamma:WeakCurrent /Herwig/Decays/EtaGammaCurrent # e+e- -> eta phi create Herwig::MEee2Mesons MEee2EtaPhi HwMELeptonLowEnergy.so create Herwig::EtaPhiCurrent /Herwig/Decays/EtaPhiCurrent set MEee2EtaPhi:WeakCurrent /Herwig/Decays/EtaPhiCurrent # e+e- -> eta omega create Herwig::MEee2Mesons MEee2EtaOmega HwMELeptonLowEnergy.so create Herwig::EtaOmegaCurrent /Herwig/Decays/EtaOmegaCurrent set MEee2EtaOmega:WeakCurrent /Herwig/Decays/EtaOmegaCurrent # e+e- > p pbar -create Herwig::MEee2Mesons MEee2ppbar HwMELeptonLowEnergy.so +create Herwig::MEee2Mesons MEee2PPbar HwMELeptonLowEnergy.so create Herwig::WeakBaryonCurrent /Herwig/Decays/CzyzCurrent create Herwig::CzyzNucleonFormFactor /Herwig/Decays/CzyzFormFactor HwFormFactors.so set /Herwig/Decays/CzyzCurrent:FormFactor /Herwig/Decays/CzyzFormFactor -set MEee2ppbar:WeakCurrent /Herwig/Decays/CzyzCurrent +set MEee2PPbar:WeakCurrent /Herwig/Decays/CzyzCurrent # e+e- > hyperons create Herwig::MEee2Mesons MEee2LL HwMELeptonLowEnergy.so create Herwig::WeakBaryonCurrent /Herwig/Decays/KornerKurodaCurrent create Herwig::KornerKurodaFormFactor /Herwig/Decays/KornerKurodaFormFactor set /Herwig/Decays/KornerKurodaCurrent:FormFactor /Herwig/Decays/KornerKurodaFormFactor set MEee2LL:WeakCurrent /Herwig/Decays/KornerKurodaCurrent # e+e- -> KKpi create Herwig::MEee2Mesons MEee2KKPi HwMELeptonLowEnergy.so create Herwig::KKPiCurrent /Herwig/Decays/KKPiCurrent set MEee2KKPi:WeakCurrent /Herwig/Decays/KKPiCurrent # e+e- -> phi pi create Herwig::MEee2Mesons MEee2PhiPi HwMELeptonLowEnergy.so create Herwig::PhiPiCurrent /Herwig/Decays/PhiPiCurrent set MEee2PhiPi:WeakCurrent /Herwig/Decays/PhiPiCurrent # e+ e- -> omega pi pi create Herwig::MEee2Mesons MEee2OmegaPiPi HwMELeptonLowEnergy.so create Herwig::OmegaPiPiCurrent /Herwig/Decays/OmegaPiPiCurrent set MEee2OmegaPiPi:WeakCurrent /Herwig/Decays/OmegaPiPiCurrent ############################################################ # NLO (POWHEG e+e- matrix elements ############################################################ library HwPowhegMELepton.so create Herwig::MEee2gZ2qqPowheg PowhegMEee2gZ2qq newdef PowhegMEee2gZ2qq:MinimumFlavour 1 newdef PowhegMEee2gZ2qq:MaximumFlavour 5 newdef PowhegMEee2gZ2qq:AlphaQCD /Herwig/Shower/AlphaQCD newdef PowhegMEee2gZ2qq:AlphaQED /Herwig/Shower/AlphaQED create Herwig::MEee2gZ2llPowheg PowhegMEee2gZ2ll newdef PowhegMEee2gZ2ll:Allowed Charged set PowhegMEee2gZ2ll:AlphaQED /Herwig/Shower/AlphaQED ############################################################ # hadron-hadron matrix elements ############################################################ ################################### # Electroweak processes ################################### # q qbar -> gamma/Z -> l+l- create Herwig::MEqq2gZ2ff MEqq2gZ2ff newdef MEqq2gZ2ff:Process 3 newdef MEqq2gZ2ff:Coupling /Herwig/Shower/AlphaQCD # q qbar to W -> l nu create Herwig::MEqq2W2ff MEqq2W2ff newdef MEqq2W2ff:Process 2 newdef MEqq2W2ff:Coupling /Herwig/Shower/AlphaQCD # W+jet create Herwig::MEPP2WJet MEWJet newdef MEWJet:WDecay Leptons # Z+jet create Herwig::MEPP2ZJet MEZJet newdef MEZJet:ZDecay ChargedLeptons # PP->WW/WZ/ZZ create Herwig::MEPP2VV MEPP2VV # PP->WZ gamma create Herwig::MEPP2VGamma MEPP2VGamma ################################### # Photon and jet processes ################################### # qqbar/gg -> gamma gamma create Herwig::MEPP2GammaGamma MEGammaGamma # hadron-hadron to gamma+jet create Herwig::MEPP2GammaJet MEGammaJet # QCD 2-to-2 create Herwig::MEQCD2to2 MEQCD2to2 # MinBias create Herwig::MEMinBias MEMinBias newdef MEMinBias:csNorm 0.01 newdef MEMinBias:Scale 2.0 ################################### # Heavy Quark ################################### # qqbar/gg -> t tbar create Herwig::MEPP2QQ MEHeavyQuark create Herwig::MEPP2SingleTop MESingleTopTChannel set MESingleTopTChannel:Process tChannel create Herwig::MEPP2SingleTop MESingleTopSChannel set MESingleTopSChannel:Process sChannel create Herwig::MEPP2SingleTop MESingleTopTW set MESingleTopTW:Process tW ################################### # Higgs processes ################################### # hadron-hadron to higgs create Herwig::MEPP2Higgs MEHiggs newdef MEHiggs:ShapeScheme MassGenerator newdef MEHiggs:Process gg newdef MEHiggs:Coupling /Herwig/Shower/AlphaQCD # hadron-hadron to higgs+jet create Herwig::MEPP2HiggsJet MEHiggsJet # PP->ZH create Herwig::MEPP2ZH MEPP2ZH newdef MEPP2ZH:Coupling /Herwig/Shower/AlphaQCD # PP->WH create Herwig::MEPP2WH MEPP2WH newdef MEPP2WH:Coupling /Herwig/Shower/AlphaQCD # PP -> Higgs via VBF create Herwig::MEPP2HiggsVBF MEPP2HiggsVBF newdef MEPP2HiggsVBF:ShowerAlphaQCD /Herwig/Shower/AlphaQCD # PP -> t tbar Higgs create Herwig::MEPP2QQHiggs MEPP2ttbarH newdef MEPP2ttbarH:QuarkType Top # PP -> b bbar Higgs create Herwig::MEPP2QQHiggs MEPP2bbbarH newdef MEPP2bbbarH:QuarkType Bottom ########################################################## # Hadron-Hadron NLO matrix elements in the Powheg scheme ########################################################## library HwPowhegMEHadron.so # q qbar -> gamma/Z -> l+l- create Herwig::MEqq2gZ2ffPowheg PowhegMEqq2gZ2ff newdef PowhegMEqq2gZ2ff:Process 3 newdef PowhegMEqq2gZ2ff:Coupling /Herwig/Shower/AlphaQCD # q qbar to W -> l nu create Herwig::MEqq2W2ffPowheg PowhegMEqq2W2ff newdef PowhegMEqq2W2ff:Process 2 newdef PowhegMEqq2W2ff:Coupling /Herwig/Shower/AlphaQCD # PP->ZH create Herwig::MEPP2ZHPowheg PowhegMEPP2ZH newdef PowhegMEPP2ZH:Coupling /Herwig/Shower/AlphaQCD # PP->WH create Herwig::MEPP2WHPowheg PowhegMEPP2WH newdef PowhegMEPP2WH:Coupling /Herwig/Shower/AlphaQCD # hadron-hadron to higgs create Herwig::MEPP2HiggsPowheg PowhegMEHiggs newdef PowhegMEHiggs:ShapeScheme MassGenerator newdef PowhegMEHiggs:Process gg newdef PowhegMEHiggs:Coupling /Herwig/Shower/AlphaQCD # PP->VV create Herwig::MEPP2VVPowheg PowhegMEPP2VV newdef PowhegMEPP2VV:Coupling /Herwig/Shower/AlphaQCD # PP -> Higgs via VBF create Herwig::MEPP2HiggsVBFPowheg PowhegMEPP2HiggsVBF newdef PowhegMEPP2HiggsVBF:ShowerAlphaQCD /Herwig/Shower/AlphaQCD # PP -> diphoton NLO create Herwig::MEPP2GammaGammaPowheg MEGammaGammaPowheg set MEGammaGammaPowheg:Process 0 set MEGammaGammaPowheg:Contribution 1 set MEGammaGammaPowheg:ShowerAlphaQCD /Herwig/Shower/AlphaQCD set MEGammaGammaPowheg:ShowerAlphaQED /Herwig/Shower/AlphaQED ########################################################## # DIS matrix elements ########################################################## # neutral current create Herwig::MENeutralCurrentDIS MEDISNC newdef MEDISNC:Coupling /Herwig/Shower/AlphaQCD newdef MEDISNC:Contribution 0 # charged current create Herwig::MEChargedCurrentDIS MEDISCC newdef MEDISCC:Coupling /Herwig/Shower/AlphaQCD newdef MEDISCC:Contribution 0 # neutral current (POWHEG) create Herwig::MENeutralCurrentDIS PowhegMEDISNC newdef PowhegMEDISNC:Coupling /Herwig/Shower/AlphaQCD newdef PowhegMEDISNC:Contribution 1 # charged current (POWHEG) create Herwig::MEChargedCurrentDIS PowhegMEDISCC newdef PowhegMEDISCC:Coupling /Herwig/Shower/AlphaQCD newdef PowhegMEDISCC:Contribution 1 ########################################################## # Gamma-Gamma matrix elements ########################################################## # fermion-antiferimon create Herwig::MEGammaGamma2ff MEgg2ff HwMEGammaGamma.so # W+ W- create Herwig::MEGammaGamma2WW MEgg2WW HwMEGammaGamma.so ########################################################## # Gamma-Hadron matrix elements ########################################################## # gamma parton -> 2 jets create Herwig::MEGammaP2Jets MEGammaP2Jets HwMEGammaHadron.so ########################################################## # Set up the Subprocesses # # Generic for all colliders ########################################################## create ThePEG::SubProcessHandler SubProcess newdef SubProcess:PartonExtractor /Herwig/Partons/PPExtractor