diff --git a/Tests/Rivet/EE/EE-Etac.in b/Tests/Rivet/EE/EE-Etac.in --- a/Tests/Rivet/EE/EE-Etac.in +++ b/Tests/Rivet/EE/EE-Etac.in @@ -1,22 +1,24 @@ # -*- ThePEG-repository -*- create ThePEG::LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxA 5.2897*GeV set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxB 5.2897*GeV set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 10.5792 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 create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so set /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S) set /Herwig/MatrixElements/MEUpsilon:Coupling 96.72794 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon decaymode Upsilon(4S)->eta_c,eta_c; 1. 1 /Herwig/Decays/DecayME0 do /Herwig/Particles/Upsilon(4S):SelectDecayModes /Herwig/Particles/Upsilon(4S)/Upsilon(4S)->eta_c,eta_c; # BES multiplicty in eta_c decays insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1724880 # eta_c -> K+K-eta and eta_c -> K+K-pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2014_I1287632 # eta_c -> K0_SK+- pi-+ and eta_c -> K+K-pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2015_I1403544 +# eta_c -> K0_SK+- pi-+ +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2012_I1187787 # eta_c -> K+K-eta', pi+pi-eta and pi+pi-eta' insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2021_I1867843 diff --git a/Tests/Rivet/EE/EE-JPsi.in b/Tests/Rivet/EE/EE-JPsi.in --- a/Tests/Rivet/EE/EE-JPsi.in +++ b/Tests/Rivet/EE/EE-JPsi.in @@ -1,40 +1,52 @@ # -*- ThePEG-repository -*- # e+ e- -> J/Psi create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEJpsi HwMELepton.so set /Herwig/MatrixElements/MEJpsi:VectorMeson /Herwig/Particles/Jpsi set /Herwig/MatrixElements/MEJpsi:Coupling 11.43148 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEJpsi set EventGenerator:EventHandler:LuminosityFunction:Energy 3.096916 set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 0.2 set /Herwig/Particles/e-:PDF /Herwig/Partons/NoPDF set /Herwig/Particles/e+:PDF /Herwig/Partons/NoPDF # J/psi-> lambda anti-lambda and sigma anti-sigma insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1510563 # J/psi -> p pbar and n nbar insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2012_I1113599 # J/Psi -> xi- and Sigma*+/- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2016_I1422780 # J/Psi -> xi0 and Sigma*0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1506414 # J/psi -> xi- xi+bar correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I1864775 # J/psi-> lambda anti-lambda insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1691850 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2099126 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2071715 # J/psi-> lambda anti-Sigma0 +cc insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2012_I1121378 # J/Psi -> xi*- xibar+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1765606 # BES J/psi dalitz insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1697377 # J/Psi -> sigma+ sigmabar- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2020_I1791570 # J/psi -> pi+ pi- pi0, K+K-pi0 and KS0 K+-pi-+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2017_I1512302 # J/psi -> pi+ pi- pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2012_I1088606 # J/psi -> pi+ pi- eta' insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1621266 # J/psi -> K+K-pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1731057 +# J/psi -> p pbar eta' +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1709205 +# J/psi -> p pbar omega +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2013_I1223625 +# J/psi -> p pbar pi0 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESII_2009_I819937 +# J/psi -> p pbar eta +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BES_2001_I556330 +# J/psi -> p pbar phi +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2016_I1411448 +# J/psi -> p nbar pi- +cc +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESII_2006_I650381 diff --git a/Tests/Rivet/EE/EE-Lambdac.in b/Tests/Rivet/EE/EE-Lambdac.in --- a/Tests/Rivet/EE/EE-Lambdac.in +++ b/Tests/Rivet/EE/EE-Lambdac.in @@ -1,25 +1,33 @@ # -*- ThePEG-repository -*- create ThePEG::LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxA 5.2897*GeV set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxB 5.2897*GeV set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 10.5792 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 create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so set /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S) set /Herwig/MatrixElements/MEUpsilon:Coupling 96.72794 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon decaymode Upsilon(4S)->Lambda_c+,Lambdabar_c-; 1. 1 /Herwig/Decays/DecayME0 do /Herwig/Particles/Upsilon(4S):SelectDecayModes /Herwig/Particles/Upsilon(4S)/Upsilon(4S)->Lambda_c+,Lambdabar_c-; # Lambda_c decays (Lambda Pi+) insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2006_I693639 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ARGUS_1992_I319105 # Lambda_c decays (Lambda Pi+ and Sigma+ Pi0) insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_1995_I392704 # Lambda_c decays (e+ nu_e) +# correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ARGUS_1994_I371613 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOII_1994_I371611 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOII_2005_I668268 +# pK- e+ nu_e +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2122399 +# Lambda0 e+ nu_e +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2127373 +# other decays # Lambda_c -> Xi* K+ -> Xi- pi_ K+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2008_I781294 +# Lambda_c -> Lambda eta pi+ +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1711896 diff --git a/Tests/Rivet/EE/EE-Psi2S.in b/Tests/Rivet/EE/EE-Psi2S.in --- a/Tests/Rivet/EE/EE-Psi2S.in +++ b/Tests/Rivet/EE/EE-Psi2S.in @@ -1,35 +1,53 @@ # -*- ThePEG-repository -*- # e+ e- -> psi(2S) create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEPsi2S HwMELepton.so set /Herwig/MatrixElements/MEPsi2S:VectorMeson /Herwig/Particles/psi(2S) set /Herwig/MatrixElements/MEPsi2S:Coupling 19.25684 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEPsi2S set EventGenerator:EventHandler:LuminosityFunction:Energy 3.686097 set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 0.2 set /Herwig/Particles/e-:PDF /Herwig/Partons/NoPDF set /Herwig/Particles/e+:PDF /Herwig/Partons/NoPDF # psi(2S) -> lambda anti-lambda and sigma anti-sigma insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1510563 # psi(2S) -> p pbar and n nbar insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1658762 # psi(2S) -> xi- and Sigma+/- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2016_I1422780 # psi(2S) -> xi0 and Sigma*0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1506414 # psi(2S) -> xi*- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1747092 +# psi(2S) -> xi*0 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2021_I1921775 # MARKII psi(2s) -> J/Psi pi+pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MARKII_1979_I144382 # BES psi(2s) -> J/Psi pi+pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BES_1999_I507637 # BES psi(2s) -> sigma+ sigmabar- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2020_I1791570 # psi(2S) -> pi+ pi- pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2012_I1088606 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BES_2005_I689969 # psi(2S) -> pi+ pi- eta' insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1621266 # psi(2S) -> xi- xi+bar correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2099144 # psi(2S) -> gamma chi_c correlations insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1507887 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2009_I832707 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1624548 +# psi(2S) -> omega KS0 KS0 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2021_I1921801 +# psi(2S) -> K+ K- eta +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2020_I1771616 +# psi(2S) -> p pbar eta' +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1709205 +# psi(2S) -> p pbar eta +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2013_I1227512 +# psi(2S) -> p pbar pi0 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2013_I1120737 +# psi(2S) -> Lambda Lambdabar omega +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I2122392 +# J/psi -> p nbar pi- +cc +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESII_2006_I716589 diff --git a/Tests/Rivet/EE/EE-Psi3770.in b/Tests/Rivet/EE/EE-Psi3770.in --- a/Tests/Rivet/EE/EE-Psi3770.in +++ b/Tests/Rivet/EE/EE-Psi3770.in @@ -1,117 +1,118 @@ # -*- ThePEG-repository -*- # e+ e- -> psi(3770) create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEPsi3770 HwMELepton.so set /Herwig/MatrixElements/MEPsi3770:VectorMeson /Herwig/Particles/psi(3770) set /Herwig/MatrixElements/MEPsi3770:Coupling 58.12041 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEPsi3770 set EventGenerator:EventHandler:LuminosityFunction:Energy 3.7711 set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 0.2 set /Herwig/Particles/e-:PDF /Herwig/Partons/NoPDF set /Herwig/Particles/e+:PDF /Herwig/Partons/NoPDF ########## branching ratios ################ # D0 inclusive insert /Herwig/Analysis/RivetAnalysis:Analyses 0 PDG_D0 # D+ inclusive insert /Herwig/Analysis/RivetAnalysis:Analyses 0 PDG_DPLUS ########## semi-leptonic ################### # CLEO D lepton spectra insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2009_I823313 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2006_I715096 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2008_I769777 # D0 -> K- semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1697371 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2007_I1091435 # D+ -> eta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2011_I875526 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1662660 # D0/+ -> pi semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1655158 # D0 -> Kbar0 pi- e+ nu_e insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1705754 # D -> pi pi e+ nu_e insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1694530 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2013_I1081165 # D0 -> pi- semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2015_I1334693 # D0 -> pi-, K- semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2015_I1391138 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2004_I654843 # D+ -> K0 pi0 semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2017_I1519425 # D+ -> K- pi+ semi-leptonic insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2016_I1411645 # D+ -> omega l nu insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2015_I1386254 ######### Two body decays ################### # D0 -> omega phi insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2022_I1900094 ########## 3-body decays #################### # Dalitz plot analysis of D -> Kpipi decays insert /Herwig/Analysis/RivetAnalysis:Analyses 0 E691_1992_I342947 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MARKIII_1987_I247266 # Dalitz plot analysis of D0 -> KS0 pi0 pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2011_I913909 # Dalitz plot analysis of D+ -> K0S pi+ pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2014_I1277070 # Dalitz plot analysis of D0 -> K0S pi+ pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2003_I633196 # Dalitz plot analysis of D+ -> K- pi+ pi+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 E791_2002_I585322 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2007_I750701 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2008_I780363 # # Kinematic distributions in the decay D0 -> K-pi+pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOII_2001_I537154 # Dalitz decay of D0 -> K0S pi0 eta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2004_I649917 # Dalitz decay of D0 -> K-pi+eta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2020_I1785816 # D -> K pi omega insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2021_I1940222 # D -> K pi eta' insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1693610 # Dalitz plot analysis of D0 -> K0S K+ K- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2020_I1799437 # Dalitz plot analysis of D0 -> K0S pi+ pi- and D0 -> K0S K+ K- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2010_I853279 # Dalitz plot analysis of D+ -> K+K+K- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2019_I1720423 # Dalitz plot analysis of D+ -> K+ K- pi+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2008_I791716 # Dalitz plot analysis of D+ -> K+ K0S pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2021_I1859124 # Dalitz plot analysis of D0 -> K+ K- pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2007_I749390 # Dalitz plot analysis of D0 -> K0S K+- pi-+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2012_I1094160 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2016_I1394391 # Dalitz decay of D+ and D+_s -> K+pi+pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2004_I654030 # Dalitz plot analysis of D+ -> pi+ pi+ pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 E791_2001_I530320 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2007_I749602 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2003_I635446 # Dalitz plot analysis of D0 -> pi+ pi- pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2007_I747154 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2016_I1441203 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2005_I679349 # D0 -> eta pi+pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2008_I779705 # D0 -> pi0 eta eta insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2018_I1662665 ########## Four body decays #################### # D+ -> KS0 pi+pi+pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1714778 # D0 -> K- pi+ pi0 pi0 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BESIII_2019_I1725265 # D0 -> K-K-K+pi+ insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2003_I626320 # D0 -> K+ K- pi+ pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2012_I1086166 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2004_I663820 insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2018_I1704426 # D0 -> K+ K- pi+ pi- and 2pi+ 2pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2017_I1519168 # D0 -> 2pi+ 2pi- insert /Herwig/Analysis/RivetAnalysis:Analyses 0 FOCUS_2007_I741543 ########## Spectra ############################## # CLEO eta, eta' phi spectra in D decays insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEOC_2006_I728043 diff --git a/Tests/Rivet/EE/EE-Xic0.in b/Tests/Rivet/EE/EE-Xic0.in --- a/Tests/Rivet/EE/EE-Xic0.in +++ b/Tests/Rivet/EE/EE-Xic0.in @@ -1,20 +1,25 @@ # -*- ThePEG-repository -*- create ThePEG::LuminosityFunction /Herwig/EventHandlers/BFactoryLuminosity set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxA 5.2897*GeV set /Herwig/EventHandlers/BFactoryLuminosity:BeamEMaxB 5.2897*GeV set /Herwig/Generators/EventGenerator:EventHandler:Cuts:MHatMin 10.5792 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 create Herwig::MEee2VectorMeson /Herwig/MatrixElements/MEUpsilon HwMELepton.so set /Herwig/MatrixElements/MEUpsilon:VectorMeson /Herwig/Particles/Upsilon(4S) set /Herwig/MatrixElements/MEUpsilon:Coupling 96.72794 set /Herwig/MatrixElements/SubProcess:MatrixElements 0 /Herwig/MatrixElements/MEUpsilon decaymode Upsilon(4S)->Xi_c0,Xi_cbar0; 1. 1 /Herwig/Decays/DecayME0 do /Herwig/Particles/Upsilon(4S):SelectDecayModes /Herwig/Particles/Upsilon(4S)/Upsilon(4S)->Xi_c0,Xi_cbar0; -# Xi_c decays +# Xi_c -> Xi- pi+ decay asymmetry insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CLEO_2000_I537236 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2021_I1851126 +# Xi_c Lambda/Sigma K* asymmetries insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2021_I1859517 -insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2021_I1851126 -# Omega decays +# Xi_c0 -> K+ K- Xi0 +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2021_I1835729 +# Omega decay asymmetries insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BABAR_2006_I719581 +# Xi_c0 -> p K- K- pi+ +insert /Herwig/Analysis/RivetAnalysis:Analyses 0 BELLE_2005_I660759 diff --git a/Tests/python/LowEnergy-EE.py.in b/Tests/python/LowEnergy-EE.py.in --- a/Tests/python/LowEnergy-EE.py.in +++ b/Tests/python/LowEnergy-EE.py.in @@ -1,540 +1,543 @@ #! @PYTHON@ from __future__ import print_function 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().decode("UTF-8") #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" : {}, "EtaprimePiPi" : {}, "4Pi" : {}, "EtaPhi" : {}, "EtaOmega" : {}, "2K1Pi" : {}, "2K2Pi" : {}, "4K" : {}, "6m" : {}, "EtaPiPi" : {}, "OmegaPi" : {}, "PiGamma" : {}, "EtaGamma" : {}, "PhiPi" : {}, "OmegaPiPi" : {}, "DD" : {}, "BB" : {}, "5Pi" : {}, "LL" : {}, "Baryon" : {} } # 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']["SND_2020_I1789269" ] = ["d01-x01-y04"] # 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"] analyses['KK']["BABAR_2020_I1769654" ] = ["d01-x01-y01"] analyses['KK']["BESIII_2021_I1866051"] = ["d01-x01-y01"] # proton-antiproton analyses['PPbar']["BESIII_2021_I1966612"] = ["d01-x01-y01"] 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']["BABAR_2013_I1247058" ] = ["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_2021_I1853007"] = ["d01-x01-y01"] analyses['PPbar']["BESIII_2021_I1847766"] = ["d01-x01-y01"] # 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 pion analyses["3Pi"]["BABAR_2004_I656680" ] = ["d01-x01-y01"] analyses["3Pi"]["BABAR_2021_I1937349" ] = ["d01-x01-y01"] analyses["3Pi"]["BESIII_2019_I1773081" ] = ["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"] analyses["3Pi"]["SND_2020_I1809286" ] = ["d01-x01-y01","d02-x01-y01", "d03-x01-y01","d03-x01-y02","d03-x01-y03"] # eta pipi analyses["EtaPiPi"]["BABAR_2007_I758568" ] = ["d01-x01-y01","d02-x01-y01"] analyses["EtaPiPi"]["BABAR_2018_I1647139" ] = ["d01-x01-y01"] analyses["EtaPiPi"]["BABAR_2018_I1700745" ] = ["d01-x01-y01","d02-x01-y01"] analyses["EtaPiPi"]["BESIII_2022_I2039027"] = ["d01-x01-y01","d02-x01-y01"] analyses["EtaPiPi"]["CMD2_2000_I532970" ] = ["d02-x01-y01"] analyses['EtaPiPi']["CMD3_2019_I1744510" ] = ["d02-x01-y01"] analyses["EtaPiPi"]["DM2_1988_I264144" ] = ["d01-x01-y01"] analyses['EtaPiPi']["ND_1991_I321108" ] = ["d06-x01-y01","d14-x01-y01"] analyses["EtaPiPi"]["SND_2015_I1332929" ] = ["d01-x01-y01"] analyses["EtaPiPi"]["SND_2018_I1638368" ] = ["d01-x01-y01"] # eta' pipi analyses["EtaprimePiPi"]["BABAR_2007_I758568" ] = ["d05-x01-y01","d06-x01-y01"] analyses["EtaprimePiPi"]["BESIII_2020_I1836509"] = ["d01-x01-y01"] # Eta Phi analyses["EtaPhi"]["BABAR_2006_I709730" ] = ["d02-x01-y01"] analyses["EtaPhi"]["BABAR_2006_I731865" ] = ["d01-x01-y02"] analyses["EtaPhi"]["BABAR_2007_I758568" ] = ["d08-x01-y01","d09-x01-y01"] analyses["EtaPhi"]["BABAR_2008_I765258" ] = ["d04-x01-y01","d05-x01-y01"] analyses["EtaPhi"]["BABAR_2017_I1511276" ] = ["d03-x01-y01"] analyses["EtaPhi"]["BABAR_2022_I2120528" ] = ["d04-x01-y01","d05-x01-y01"] analyses["EtaPhi"]["BELLE_2009_I823878" ] = ["d01-x01-y01"] analyses["EtaPhi"]["BESII_2008_I801210" ] = ["d01-x01-y03"] analyses["EtaPhi"]["BESIII_2021_I1857930"] = ["d01-x01-y01"] analyses["EtaPhi"]["CMD3_2017_I1606078" ] = ["d01-x01-y01"] analyses["EtaPhi"]["CMD3_2019_I1740541" ] = ["d01-x01-y06","d02-x01-y06","d03-x01-y06"] analyses["EtaPhi"]["SND_2018_I1693737" ] = ["d01-x01-y01"] analyses["EtaPhi"]["SND_2019_I1726419" ] = ["d01-x01-y01","d01-x01-y03"] analyses["EtaPhi"]["SND_2021_I1942539" ] = ["d01-x01-y01"] # Eta Omega analyses["EtaOmega"]["BABAR_2006_I709730" ] = ["d02-x01-y01"] analyses["EtaOmega"]["BABAR_2021_I1938254" ] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01","d04-x01-y01"] analyses["EtaOmega"]["BESII_2008_I801210" ] = ["d01-x01-y03"] analyses["EtaOmega"]["BESIII_2020_I1817739"] = ["d01-x01-y01"] analyses["EtaOmega"]["BESIII_2022_I2047667"] = ["d01-x01-y02"] analyses["EtaOmega"]["CMD3_2017_I1606078" ] = ["d01-x01-y01","d01-x01-y02"] analyses["EtaOmega"]["SND_2016_I1473343" ] = ["d01-x01-y01"] analyses["EtaOmega"]["SND_2019_I1726419" ] = ["d01-x01-y01","d01-x01-y02"] analyses["EtaOmega"]["SND_2020_I1800477" ] = ["d01-x01-y01","d03-x01-y01"] # 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","d03-x01-y01"] analyses["4Pi"]["BESII_2008_I801210" ] = ["d01-x01-y01"] analyses["4Pi"]["BESIII_2022_I2047667" ] = ["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", "d01-x01-y01","d02-x01-y01","d03-x01-y01","d04-x01-y01","d10-x01-y03"] 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"] analyses["4Pi"]["BESIII_2020_I1817739" ] = ["d02-x01-y01"] analyses["4Pi"]["BESIII_2021_I1929314" ] = ["d01-x01-y03"] # (these are Omega(-> pi0 gamma) pi0) analyses["OmegaPi"]["CMD2_2003_I616446" ] = ["d01-x01-y01"] analyses["OmegaPi"]["SND_2000_I503946" ] = ["d01-x01-y01"] analyses["OmegaPi"]["SND_2000_I527752" ] = ["d01-x01-y01"] analyses["OmegaPi"]["SND_2016_I1489182" ] = ["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"] analyses["OmegaPi"]["BESIII_2020_I1817739"] = ["d02-x01-y01"] # omega 2 pi analyses["OmegaPiPi"]["BABAR_2007_I758568" ] = ["d01-x01-y01","d03-x01-y01","d04-x01-y01"] analyses["OmegaPiPi"]["BABAR_2018_I1700745" ] = ["d01-x01-y01","d03-x01-y01"] analyses['OmegaPiPi']["BESIII_2021_I1999208"] = ["d01-x01-y01"] analyses["OmegaPiPi"]["CMD2_2000_I532970" ] = ["d01-x01-y01"] analyses["OmegaPiPi"]["DM1_1981_I166964" ] = ["d01-x01-y01"] analyses["OmegaPiPi"]["DM2_1992_I339265" ] = ["d02-x01-y01"] analyses['OmegaPiPi']["ND_1991_I321108" ] = ["d14-x01-y01"] # 5 pion 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["5Pi"]["BESIII_2021_I1929314" ] = ["d01-x01-y07"] # 2K 1 pi analyses["2K1Pi"]["BABAR_2008_I765258" ] = ["d01-x01-y01","d02-x01-y01"] analyses["2K1Pi"]["BABAR_2017_I1511276" ] = ["d01-x01-y01"] analyses["2K1Pi"]["BESII_2008_I801208" ] = ["d01-x01-y01"] analyses["2K1Pi"]["BESIII_2018_I1691798"] = ["d01-x01-y01"] analyses["2K1Pi"]["BESIII_2022_I2033007"] = ["d01-x01-y01","d03-x01-y01","d04-x01-y01"] analyses["2K1Pi"]["DM1_1982_I176801" ] = ["d01-x01-y01"] analyses["2K1Pi"]["DM2_1991_I318558" ] = ["d01-x01-y01","d02-x01-y01"] analyses["2K1Pi"]["SND_2018_I1637194" ] = ["d01-x01-y01"] analyses["2K1Pi"]["SND_2020_I1806118" ] = ["d01-x01-y01"] # phi pi analyses["PhiPi"]["BABAR_2008_I765258" ] = ["d02-x01-y01","d03-x01-y01"] analyses["PhiPi"]["BABAR_2017_I1511276" ] = ["d01-x01-y01","d02-x01-y01"] analyses["PhiPi"]["BESIII_2022_I2033007"] = ["d01-x01-y01","d02-x01-y01"] analyses["PhiPi"]["SND_2020_I1806118" ] = ["d02-x01-y01"] # 2K 2 pi analyses["2K2Pi"]["BABAR_2005_I676691" ] = ["d02-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"]["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_2014_I1287920" ] = ["d09-x01-y01","d10-x01-y01","d11-x01-y01"] analyses["2K2Pi"]["BABAR_2017_I1511276" ] = ["d03-x01-y01","d04-x01-y01"] analyses["2K2Pi"]["BABAR_2017_I1591716" ] = ["d01-x01-y01","d02-x01-y01"] analyses['2K2Pi']["BESII_2007_I750713" ] = ["d01-x01-y04"] analyses["2K2Pi"]["BESII_2008_I801210" ] = ["d01-x01-y02"] analyses["2K2Pi"]["BESII_2008_I801208" ] = ["d01-x01-y02"] analyses['2K2Pi']["BESIII_2018_I1699641"] = ["d01-x01-y01","d02-x01-y01"] analyses['2K2Pi']["BESIII_2020_I1775344"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", "d04-x01-y01","d05-x01-y01","d06-x01-y01"] analyses['2K2Pi']["BESIII_2021_I1997451"] = ["d01-x01-y01"] analyses["2K2Pi"]["BELLE_2009_I809630" ] = ["d01-x01-y01"] analyses["2K2Pi"]["CMD3_2016_I1395968" ] = ["d01-x01-y01"] analyses['2K2Pi']["CMD3_2019_I1770428" ] = ["d01-x01-y06"] analyses["2K2Pi"]["DM1_1982_I169382" ] = ["d01-x01-y01"] analyses["2K2Pi"]["BESIII_2021_I1929314"] = ["d01-x01-y01"] # 4K analyses["4K"]["BESIII_2019_I1743841"] = ["d01-x01-y01","d02-x01-y01"] analyses["4K"]["BESIII_2021_I1929314"] = ["d01-x01-y02"] 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"]["BESIII_2021_I1929314"] = ["d01-x01-y05","d01-x01-y06"] 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["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"]["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"] analyses['6m']["BESIII_2020_I1788734"] = ["d01-x01-y01"] analyses['6m']["BABAR_2021_I1844422" ] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01", "d04-x01-y01","d05-x01-y01","d06-x01-y01"] analyses['6m']["BESIII_2020_I1837725" ] = ["d01-x01-y01"] analyses["6m"]["BABAR_2021_I1938254"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01","d04-x01-y01","d05-x01-y01"] analyses["6m"]["CMD3_2022_I2108984"] = ["d01-x01-y01","d02-x01-y01","d02-x01-y02"] analyses["6m"]["BABAR_2022_I2120528"] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01","d04-x01-y01","d05-x01-y01", "d06-x01-y01","d07-x01-y01","d08-x01-y01","d09-x01-y01","d10-x01-y01"] # other baryon processes analyses['Baryon']["BESIII_2017_I1509241" ] = ["d01-x01-y01"] analyses['Baryon']["BESIII_2021_I1845443" ] = ["d01-x01-y01","d02-x01-y01"] analyses['Baryon']["BESIII_2021_I1859248" ] = ["d01-x01-y01"] analyses["Baryon"]["BESIII_2021_I1929314" ] = ["d01-x01-y04","d01-x01-y08"] # 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"]["BELLE_2020_I1789775" ] = ["d01-x01-y01"] analyses["DD"]["BELLE_2019_I1762826" ] = ["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_2022_I1989527" ] = ["d01-x01-y03","d02-x01-y03"] analyses["DD"]["BESIII_2021_I1933191" ] = ["d01-x01-y03"] 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"] analyses["DD"]["BABAR_2007_I729388" ] = ["d01-x01-y01"] analyses["DD"]["BESIII_2015_I1329785" ] = ["d01-x01-y08","d02-x01-y08","d03-x01-y08"] analyses["DD"]["BESIII_2017_I1494065" ] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["BESIII_2017_I1596897" ] = ["d01-x01-y01"] analyses["DD"]["BESIII_2018_I1653121" ] = ["d01-x01-y01","d01-x01-y02"] analyses["DD"]["BESIII_2020_I1762922" ] = ["d01-x01-y01"] analyses["DD"]["BESIII_2018_I1633425" ] = ["d01-x01-y01"] analyses["DD"]["BESIII_2018_I1685535" ] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["BELLE_2011_I878228" ] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03"] analyses["DD"]["BABAR_2010_I864027" ] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03"] analyses["DD"]["BABAR_2009_I815035" ] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03","d02-x01-y01"] analyses["DD"]["BES_1999_I508349" ] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03","d01-x01-y04"] analyses["DD"]["BESIII_2020_I1795949" ] = ["d01-x01-y01","d02-x01-y01"] analyses["DD"]["BESIII_2021_I1867196" ] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] +analyses["DD"]["BESIII_2020_I1784442" ] = ["d01-x01-y01","d02-x01-y01"] +analyses["DD"]["BESIII_2021_I1849633" ] = ["d01-x01-y01"] +analyses["DD"]["BESIII_2021_I1880103" ] = ["d01-x01-y01","d02-x01-y01"] # BB analyses["BB"]["BELLE_2008_I764099" ] = ["d01-x01-y01","d02-x01-y01", "d03-x01-y01","d04-x01-y01"] analyses["BB"]["BELLE_2016_I1389855" ] = ["d01-x01-y02","d01-x01-y03"] analyses["BB"]["BELLE_2021_I1859137" ] = ["d01-x01-y01","d01-x01-y02","d01-x01-y03"] analyses["BB"]["CLEO_1991_I29927" ] = ["d01-x01-y01"] analyses["BB"]["CLEO_1999_I474676" ] = ["d01-x01-y01","d01-x01-y02"] analyses["BB"]["CUSB_1991_I325661" ] = ["d01-x01-y01"] # hyperons analyses["LL"]["BABAR_2007_I760730" ] = ["d01-x01-y01","d02-x01-y01","d03-x01-y01"] analyses["LL"]["BESIII_2018_I1627871"] = ["d01-x01-y01"] analyses["LL"]["BESIII_2019_I1726357"] = ["d01-x01-y01"] analyses["LL"]["BESIII_2019_I1758883"] = ["d01-x01-y05"] analyses["LL"]["BESIII_2020_I1814783"] = ["d01-x01-y01","d01-x01-y02", "d02-x01-y01","d02-x01-y02"] analyses["LL"]["BESIII_2020_I1823448"] = ["d01-x01-y04"] analyses["LL"]["BESIII_2021_I1866233"] = ["d01-x01-y01"] analyses["LL"]["BESIII_2021_I1940960"] = ["d01-x01-y01"] analyses["LL"]["BESIII_2021_I1900124"] = ["d01-x01-y01"] analyses["LL"]["DM2_1990_I297706" ] = ["d02-x01-y01"] # list the analysis if required and quit() allProcesses=False if "All" in opts.processes : allProcesses=True 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] : if(analysis=="CMD3_2019_I1770428") : for iy in range(1,3) : output+= " -m/%s/%s" % (analysis,"d02-x01-y0%s"%iy) elif(analysis=="BES_1999_I508349") : for ix in range(2,4) : for iy in range(1,3) : output+= " -m/%s/%s" % (analysis,"d0%s-x01-y0%s"%(ix,iy)) elif(analysis=="BESIII_2019_I1726357") : for ix in range(2,4) : output+= " -m/%s/%s" % (analysis,"d0%s-x01-y01"% ix) elif(analysis=="BESIII_2020_I1775344") : for ix in range(1,6) : output+= " -m/%s/%s" % (analysis,"d07-x01-y0%s"% ix) output+= " -m/%s/%s" % (analysis,"d08-x01-y0%s"% ix) elif(analysis=="BESIII_2020_I1814783") : for ix in range(1,3) : output+= " -m/%s/%s" % (analysis,"d03-x01-y0%s"% ix) elif(analysis=="SND_2020_I1809286") : for ix in range(1,5) : output+= " -m/%s/%s" % (analysis,"d04-x01-y0%s"% ix) 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", "LL" : "MEee2LL" , "2K1Pi" : "MEee2KKPi" } # get the particle masses from Herwig particles = { "pi+" : 0., "pi0" : 0. ,"eta" : 0. ,"eta'" : 0. ,"phi" : 0. ,"omega" : 0. ,"p+" : 0. ,"K+" : 0.} for val in particles : tempTxt = "get /Herwig/Particles/%s:NominalMass\nget /Herwig/Particles/%s:WidthLoCut\n" % (val,val) with open("temp.in",'w') as f: f.write(tempTxt) p = subprocess.Popen(["../src/Herwig", "read","temp.in"],stdout=subprocess.PIPE) vals = p.communicate()[0].split() mass = float(vals[0])-float(vals[1]) particles[val]=mass os.remove("temp.in") # minimum CMS energies for specific processes minMass = { "PiPi" : 2.*particles["pi+"], "KK" : 2.*particles["K+"], "3Pi" : 2.*particles["pi+"]+particles["pi0"], "4Pi" : 2.*particles["pi+"]+2.*particles["pi0"], "EtaPiPi" : particles["eta"]+2.*particles["pi+"], "EtaprimePiPi" : particles["eta'"]+2.*particles["pi+"], "EtaPhi" : particles["phi"]+particles["eta"], "EtaOmega" : particles["omega"]+particles["eta"], "OmegaPi" : particles["omega"]+particles["pi0"], "OmegaPiPi" : particles["omega"]+2.*particles["pi0"], "PhiPi" : particles["phi"]+particles["pi0"], "PiGamma" : particles["pi0"], "EtaGamma" : particles["eta"], "PPbar" : 2.*particles["p+"], "LL" : 0., "2K1Pi" : 2.*particles["K+"]+particles["pi0"] } # 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] and minMass[process]<=energy) : 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] and minMass[process]<=energy) : anals[0].append(matrix) else : if(matrix=="") : energies[energy]=[[],[analysis]] elif(minMass[process]<=energy) : 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="" matrices = energies[energy][0] if(allProcesses) : matrices = me.values() for matrix in matrices: 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, 'name' : "EE"}) 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, 'name' : "EE"}) 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)