diff --git a/Tests/python/LowEnergy-EE-NonPerturbative.in b/Tests/python/LowEnergy-EE-NonPerturbative.in --- a/Tests/python/LowEnergy-EE-NonPerturbative.in +++ b/Tests/python/LowEnergy-EE-NonPerturbative.in @@ -1,109 +1,40 @@ # -*- ThePEG-repository -*- read snippets/EECollider.in ################################################## # Selected the hard process ################################################## # leading-order processes ################################################## cd /Herwig/MatrixElements -# e+ e- -> pi+pi- -create Herwig::MEee2Mesons MEee2Pions HwMELeptonLowEnergy.so -create Herwig::TwoPionCzyzCurrent /Herwig/Decays/TwoPionCzyzCurrent -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 -#create Herwig::TwoPionPhotonCurrent /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::WeakBaryonCurrent /Herwig/Decays/CzyzCurrent -create Herwig::CzyzNucleonFormFactor /Herwig/Decays/CzyzFormFactor -set /Herwig/Decays/CzyzCurrent:FormFactor /Herwig/Decays/CzyzFormFactor -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 -# default e+e- > q qbar (5 flavours d,u,s,c,b) +# set the process ${processes} ################################################## # LEP physics parameters (override defaults) ################################################## cd /Herwig/Generators set EventGenerator:EventHandler:LuminosityFunction:Energy ${ECMS} 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 set /Herwig/Particles/pi0:Stable Stable set /Herwig/Particles/K_S0:Stable Stable cd /Herwig/Generators ################################################## ## prepare for Rivet analysis or HepMC output ## when running with parton shower ################################################## create ThePEG::RivetAnalysis /Herwig/Analysis/Rivet RivetAnalysis.so insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/Rivet ${ANALYSES} ################################################### # Save run for later usage with 'Herwig run' ################################################## set EventGenerator:MaxErrors 10000 set EventGenerator:EventHandler:StatLevel Full set EventGenerator:EventHandler:CascadeHandler NULL saverun Rivet-LowEnergy-EE-NonPerturbative-${ECMS} EventGenerator diff --git a/src/defaults/MatrixElements.in b/src/defaults/MatrixElements.in --- a/src/defaults/MatrixElements.in +++ b/src/defaults/MatrixElements.in @@ -1,250 +1,324 @@ # -*- 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 +#create Herwig::TwoPionPhotonCurrent /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::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 +# 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