diff --git a/Tests/Rivet/EE/EE-91.in b/Tests/Rivet/EE/EE-91.in --- a/Tests/Rivet/EE/EE-91.in +++ b/Tests/Rivet/EE/EE-91.in @@ -1,146 +1,152 @@ # -*- ThePEG-repository -*- ################################################## # LEP physics parameters (override defaults) ################################################## set /Herwig/Generators/EventGenerator:EventHandler:LuminosityFunction:Energy 91.2 ################################################## # select the analyses ################################################## # Validated ################################################## # ALEPH charged particle multiplicity insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1991_S2435284 # OPAL charged particle multiplicity insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1992_I321190 +# OPAL charged particle multiplicity (light flavour separated) +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2001_I536266 # DELPHI charged particle multiplicity insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1991_I301657 # ALEPH main LEP I QCD summary paper insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1996_S3486095 # ALEPH D* insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1999_S4193598 # OPAL D* insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1995_I382219 # OPAL charged hadron analysis insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1994_S2927284 # OPAL Delta++ analysis insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1995_S3198391 # OPAL J/Psi analysis analysis insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1996_S3257789 # ALEPH eta/omega analysis insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_2002_S4823664 # OPAL K*0 analysis insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1997_S3608263 # OPAL flavour specific charged multiplicities etc insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1998_S3780481 # OPAL f_0,f_2 and phi production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1998_S3702294 # OPAL gamma,pi0,eta,eta',rho+/-,a0+/- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1998_S3749908 # OPAL K0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2000_S4418603 # OPAL K* +/- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1993_I342766 # SLD flavour specific charged multiplicities etc insert /Herwig/Analysis/RivetAnalysis:Analyses 0 SLD_1996_S3398250 # SLD flavour specific charged multiplicities etc insert /Herwig/Analysis/RivetAnalysis:Analyses 0 SLD_1999_S3743934 # SLD flavour specific charged multiplicities etc insert /Herwig/Analysis/RivetAnalysis:Analyses 0 SLD_2004_S5693039 # OPAL event shapes and multiplicities at different energies insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2004_S6132243 # ALEPH jet and event shapes at many energies insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_2004_S5765862 # OPAL/JADE jet rates at many energies insert /Herwig/Analysis/RivetAnalysis:Analyses 0 JADE_OPAL_2000_S4300807 # DELPHI strange baryon production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_S3137023 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_2006_I719387 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_2000_I524694 # DELPHI f_0, rho_0 and f_2 production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1999_S3960137 # DELPHI K*0 and phi production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1996_I420528 # DELPHI f'_2 production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1996_I416741 # DELPHI pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1996_I401100 # DELPHI K- K*+/- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_I377487 # DELPHI K+/-, p pbar insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_I394052 # DELPHI Delta++ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_I399737 # DELPHI pi, K p, flavour seperated insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1998_I473409 # OPAL strange baryon production insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1997_S3396100 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1997_I421977 # DELPHI tuning paper insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1996_S3430090 # DELPHI b quark insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_2011_I890503 # ALEPH b quark insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_2001_S4656318 # SLD b quark insert /Herwig/Analysis/RivetAnalysis:Analyses 0 SLD_2002_S4869273 # OPAL b quark insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2003_I599181 # PDG hadron multiplicities and ratios insert /Herwig/Analysis/RivetAnalysis:Analyses 0 PDG_HADRON_MULTIPLICITIES insert /Herwig/Analysis/RivetAnalysis:Analyses 0 PDG_HADRON_MULTIPLICITIES_RATIOS # OPAL from gluon paper insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2004_I648738 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2004_I631361:PROCESS=QQ # L3 eta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 L3_1992_I336180 # L3 jet rates insert /Herwig/Analysis/RivetAnalysis:Analyses 0 L3_2004_I652683 # ALEPH pi+-, K+- and (p, anti-p) insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1995_I382179 # L3 eta' and omega specta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 L3_1997_I427107 # charged mult, different # of jets insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1992_I334948 # charged mult, different rapidity regions insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1991_I324035 # lambda, lambdabar correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1993_I360638 # delphi flavour sep charged # lambda, lambdabar correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1997_I428178 # ALEPH pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1997_I427131 # ALEPH Lambda polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1996_I415745 # OPAL Lambda polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1997_I447188 # DELPHI Lambda/K asymmetry insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_I382285 # OPAL b baryon polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1998_I474012 # DELPHI b baryon polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_2000_I513614 # ALEPH b baryon polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1996_I402895 # OPAL rho and comega polarization insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2000_I502750 # DELPHI charged particle distributions insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1999_I448370 # OPAL K0 spectrum insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1995_I393503 # OPAL phi, D*, B* spin insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1997_I440103 # DELPHI B* insert /Herwig/Analysis/RivetAnalysis:Analyses 0 DELPHI_1995_I395026 # ALEPH B* insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALEPH_1995_I398426 +# OPAL B* +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_1996_I428493 +# L3 B* +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 L3_1995_I381046 ################################################## # unvalidated ################################################## # OPAL 4 jet angles insert /Herwig/Analysis/RivetAnalysis:Analyses 0 OPAL_2001_S4553896 ################################################## # MC ################################################## insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_Bottom_Hadrons insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_Particle_Multiplicities diff --git a/Tests/python/LowEnergy.py b/Tests/python/LowEnergy.py --- a/Tests/python/LowEnergy.py +++ b/Tests/python/LowEnergy.py @@ -1,391 +1,392 @@ #! /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" :{} , "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"] # 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"] # 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"] # 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"] # 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"] # 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", "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"] # 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"] # 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", "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", "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"] # 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["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", "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", "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