diff --git a/src/Matchbox/FiveFlavourNoBMassScheme.in b/src/Matchbox/FiveFlavourNoBMassScheme.in --- a/src/Matchbox/FiveFlavourNoBMassScheme.in +++ b/src/Matchbox/FiveFlavourNoBMassScheme.in @@ -1,126 +1,126 @@ # -*- ThePEG-repository -*- cd /Herwig/MatrixElements/Matchbox do Factory:StartParticleGroup p insert Factory:ParticleGroup 0 /Herwig/Particles/b insert Factory:ParticleGroup 0 /Herwig/Particles/bbar insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup do Factory:StartParticleGroup pbar insert Factory:ParticleGroup 0 /Herwig/Particles/b insert Factory:ParticleGroup 0 /Herwig/Particles/bbar insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup cd /Herwig/Merging/ do MergingFactory:StartParticleGroup p insert MergingFactory:ParticleGroup 0 /Herwig/Particles/b insert MergingFactory:ParticleGroup 0 /Herwig/Particles/bbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup do MergingFactory:StartParticleGroup pbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/b insert MergingFactory:ParticleGroup 0 /Herwig/Particles/bbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup cd /Herwig/Particles do b:UnsetHardProcessMass do bbar:UnsetHardProcessMass do c:UnsetHardProcessMass do cbar:UnsetHardProcessMass set b:NominalMass 0*GeV set bbar:NominalMass 0*GeV set c:NominalMass 0*GeV set cbar:NominalMass 0*GeV cd /Herwig/DipoleShower/Kernels set IFgx2bbbarxDipoleKernel:UseKernel Yes set IFgx2bbarbxDipoleKernel:UseKernel Yes ## Even though the b-quarks are massless in this scheme, ## the spectator can be massive. set IFMgx2bbbarxDipoleKernel:UseKernel Yes set IFMgx2bbarbxDipoleKernel:UseKernel Yes set IIgx2bbbarxDipoleKernel:UseKernel Yes set IIgx2bbarbxDipoleKernel:UseKernel Yes ## First remove the bbar splittings. cd /Herwig/Shower do SplittingGenerator:DeleteInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:DeleteInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting g->b,bbar; GtobbbarSudakov -do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov +#do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:DeleteInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->g,bbar; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov ## Now insert the bbar splittings. do SplittingGenerator:AddInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting g->b,bbar; GtobbbarSudakov -do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov +#do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:AddInitialSplitting bbar->g,bbar; QtoGQSudakov do SplittingGenerator:AddInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:MaximumFlavour 5 diff --git a/src/Matchbox/FiveFlavourScheme.in b/src/Matchbox/FiveFlavourScheme.in --- a/src/Matchbox/FiveFlavourScheme.in +++ b/src/Matchbox/FiveFlavourScheme.in @@ -1,110 +1,110 @@ # -*- ThePEG-repository -*- cd /Herwig/MatrixElements/Matchbox do Factory:StartParticleGroup p insert Factory:ParticleGroup 0 /Herwig/Particles/b insert Factory:ParticleGroup 0 /Herwig/Particles/bbar insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup do Factory:StartParticleGroup pbar insert Factory:ParticleGroup 0 /Herwig/Particles/b insert Factory:ParticleGroup 0 /Herwig/Particles/bbar insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup cd /Herwig/Merging/ do MergingFactory:StartParticleGroup p insert MergingFactory:ParticleGroup 0 /Herwig/Particles/b insert MergingFactory:ParticleGroup 0 /Herwig/Particles/bbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup do MergingFactory:StartParticleGroup pbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/b insert MergingFactory:ParticleGroup 0 /Herwig/Particles/bbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup cd /Herwig/Particles set b:HardProcessMass 0*GeV set bbar:HardProcessMass 0*GeV cd /Herwig/DipoleShower/Kernels set IFgx2bbbarxDipoleKernel:UseKernel Yes set IFgx2bbarbxDipoleKernel:UseKernel Yes set IFMgx2bbbarxDipoleKernel:UseKernel Yes set IFMgx2bbarbxDipoleKernel:UseKernel Yes set IIgx2bbbarxDipoleKernel:UseKernel Yes set IIgx2bbarbxDipoleKernel:UseKernel Yes ## First remove the bbar splittings. In case someone read a FourFlavour snippet. cd /Herwig/Shower do SplittingGenerator:DeleteInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:DeleteInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting g->b,bbar; GtobbbarSudakov -do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov +#do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:DeleteInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->g,bbar; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov ## Now insert the bbar splittings. do SplittingGenerator:AddInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting g->b,bbar; GtobbbarSudakov -do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov +#do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:AddInitialSplitting bbar->g,bbar; QtoGQSudakov do SplittingGenerator:AddInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:MaximumFlavour 5 diff --git a/src/Matchbox/FourFlavourScheme.in b/src/Matchbox/FourFlavourScheme.in --- a/src/Matchbox/FourFlavourScheme.in +++ b/src/Matchbox/FourFlavourScheme.in @@ -1,105 +1,105 @@ # -*- ThePEG-repository -*- cd /Herwig/MatrixElements/Matchbox do Factory:StartParticleGroup p insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup do Factory:StartParticleGroup pbar insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup do Factory:StartParticleGroup j insert Factory:ParticleGroup 0 /Herwig/Particles/c insert Factory:ParticleGroup 0 /Herwig/Particles/cbar insert Factory:ParticleGroup 0 /Herwig/Particles/s insert Factory:ParticleGroup 0 /Herwig/Particles/sbar insert Factory:ParticleGroup 0 /Herwig/Particles/d insert Factory:ParticleGroup 0 /Herwig/Particles/dbar insert Factory:ParticleGroup 0 /Herwig/Particles/u insert Factory:ParticleGroup 0 /Herwig/Particles/ubar insert Factory:ParticleGroup 0 /Herwig/Particles/g do Factory:EndParticleGroup cd /Herwig/Merging/ do MergingFactory:StartParticleGroup p insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup do MergingFactory:StartParticleGroup pbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup do MergingFactory:StartParticleGroup j insert MergingFactory:ParticleGroup 0 /Herwig/Particles/c insert MergingFactory:ParticleGroup 0 /Herwig/Particles/cbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/s insert MergingFactory:ParticleGroup 0 /Herwig/Particles/sbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/d insert MergingFactory:ParticleGroup 0 /Herwig/Particles/dbar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/u insert MergingFactory:ParticleGroup 0 /Herwig/Particles/ubar insert MergingFactory:ParticleGroup 0 /Herwig/Particles/g do MergingFactory:EndParticleGroup cd /Herwig/DipoleShower/Kernels set IFgx2bbbarxDipoleKernel:UseKernel No set IFgx2bbarbxDipoleKernel:UseKernel No set IFMgx2bbbarxDipoleKernel:UseKernel No set IFMgx2bbarbxDipoleKernel:UseKernel No set IIgx2bbbarxDipoleKernel:UseKernel No set IIgx2bbarbxDipoleKernel:UseKernel No cd /Herwig/Shower do SplittingGenerator:DeleteInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:DeleteInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:DeleteInitialSplitting g->b,bbar; GtobbbarSudakov -do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov +#do SplittingGenerator:DeleteInitialSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:DeleteInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->g,bbar; QtoGQSudakov do SplittingGenerator:DeleteInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:DeleteInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:MaximumFlavour 4