diff --git a/Hadronization/ColourReconnector.h b/Hadronization/ColourReconnector.h --- a/Hadronization/ColourReconnector.h +++ b/Hadronization/ColourReconnector.h @@ -1,265 +1,271 @@ // -*- C++ -*- // // ColourReconnector.h is a part of Herwig - A multi-purpose Monte Carlo event generator // Copyright (C) 2002-2017 The Herwig Collaboration // // Herwig is licenced under version 3 of the GPL, see COPYING for details. // Please respect the MCnet academic guidelines, see GUIDELINES for details. // #ifndef HERWIG_ColourReconnector_H #define HERWIG_ColourReconnector_H #include #include "CluHadConfig.h" #include "ColourReconnector.fh" namespace Herwig { using namespace ThePEG; /** \ingroup Hadronization * \class ColourReconnector * \brief Class for changing colour reconnections of partons. * \author Alberto Ribon, Christian Roehr * * This class does the nonperturbative colour rearrangement, after the * nonperturbative gluon splitting and the "normal" cluster formation. * It uses the list of particles in the event record, and the collections of * "usual" clusters which is passed to the main method. If the colour * reconnection is actually accepted, then the previous collections of "usual" * clusters is first deleted and then the new one is created. * * * @see \ref ColourReconnectorInterfaces "The interfaces" * defined for ColourReconnector. */ class ColourReconnector: public Interfaced { public: /** * Does the colour rearrangement, starting out from the list of particles in * the event record and the collection of "usual" clusters passed as * arguments. If the actual rearrangement is accepted, the initial collection of * clusters is overridden by the old ones. */ void rearrange(ClusterVector & clusters); using CluVecIt = ClusterVector::iterator; private: /** PRIVATE MEMBER FUNCTIONS */ /** * @brief Calculates the sum of the squared cluster masses. * @arguments q, aq vectors containing the quarks and antiquarks respectively * @return Sum of cluster squared masses M^2_{q[i],aq[i]}. */ Energy2 _clusterMassSum(const PVector & q, const PVector & aq) const; /** * @brief Examines whether the cluster vector (under the given permutation of * the antiquarks) contains colour-octet clusters * @param cv Cluster vector * @param P Permutation, a vector of permutated indices from 0 to * cv.size()-1 */ bool _containsColour8(const ClusterVector & cv, const vector & P) const; /** * @brief A Metropolis-type algorithm which finds a local minimum in the * total sum of cluster masses * @arguments cv cluster vector */ void _doRecoStatistical(ClusterVector & cv) const; /** * @brief Plain colour reconnection as used in Herwig 2.5.0 * @arguments cv cluster vector */ void _doRecoPlain(ClusterVector & cv) const; /** * Baryonic Colour Reconnection model */ void _doRecoBaryonic(ClusterVector & cv) const; void _makeBaryonicClusters(ClusterPtr &c1, ClusterPtr &c2, ClusterPtr &c3, ClusterPtr &newcluster1, ClusterPtr &newcluster2) const; /** * @brief Finds the cluster in cv which, if reconnected with the given * cluster cl, would result in the smallest sum of cluster masses. * If no reconnection partner can be found, a pointer to the * original Cluster cl is returned. * @arguments cv cluster vector * cl cluster iterator (must be from cv) which wants to have a reconnection partner * @return iterator to the found cluster, or the original cluster pointer if * no mass-reducing combination can be found */ CluVecIt _findRecoPartner(CluVecIt cl, ClusterVector & cv) const; CluVecIt _findPartnerRapidity(CluVecIt cl, ClusterVector & cv) const; CluVecIt _findPartnerBaryonic(CluVecIt cl, ClusterVector & cv, bool & tetraCand, const ClusterVector& a, CluVecIt &baryonic1, CluVecIt &baryonic2 ) const; /** * @brief Reconnects the constituents of the given clusters to the (only) * other possible cluster combination. * @return pair of pointers to the two new clusters */ pair _reconnect(ClusterPtr &c1, ClusterPtr &c2) const; /** * Reconnection method for baryonic reconenction model */ pair _reconnectBaryonic(ClusterPtr &c1, ClusterPtr &c2) const; /** * @brief At random, swap two antiquarks, if not excluded by the * constraint that there must not be any colour-octet clusters. * @arguments q, aq vectors containing the quarks and antiquarks respectively * maxtries maximal number of tries to find a non-colour-octet * reconfiguration * @return Pair of ints indicating the indices of the antiquarks to be * swapped. Returns (-1,-1) if no valid reconfiguration could be * found after maxtries trials */ pair _shuffle(const PVector & q, const PVector & aq, unsigned maxtries = 10) const; /** DATA MEMBERS */ /** * Specifies the colour reconnection algorithm to be used. */ int _algorithm = 0; /** - * The annealing factor is the ratio of two successive temperature steps: - * T_n = _annealingFactor * T_(n-1) - */ - double _annealingFactor = 0.9; - - /** - * Number of temperature steps in the statistical annealing algorithm - */ - unsigned int _annealingSteps = 50; - - /** * Do we do colour reconnections? */ int _clreco = 0; + /** + * Statistical Reco: * Factor used to determine the initial temperature according to * InitialTemperature = _initTemp * median {energy changes in a few random * rearrangements} */ - double _initTemp = 0.1; + double _initTemp = 0.01; + + /** + * Statistical Reco: + * The annealing factor is the ratio of two successive temperature steps: + * T_n = _annealingFactor * T_(n-1) + */ + double _annealingFactor = 0.21; + + /** + * Statistical Reco: + * Number of temperature steps in the statistical annealing algorithm + */ + unsigned int _annealingSteps = 10; + + /** + * Statistical Reco: + * The number of tries per temperature steps is the number of clusters times + * this factor. + */ + double _triesPerStepFactor = 0.66; /** * Probability that a found reconnection possibility is actually accepted. */ double _preco = 0.5; double _precoBaryonic = 0.5; - - /** - * The number of tries per temperature steps is the number of clusters times - * this factor. - */ - double _triesPerStepFactor = 5.0; + /** * maximum allowed distance in the eta phi space for reconnection to occur */ /** * Maximium distance for reconnections */ Length _maxDistance = picometer; /** * @return true, if the two partons are splitting products of the same * gluon */ bool _isColour8(tcPPtr p, tcPPtr q) const; /** * Option for handling octets */ unsigned int _octetOption = 0; public: /** @name Functions used by the persistent I/O system. */ //@{ /** * Function used to write out object persistently. * @param os the persistent output stream written to. */ void persistentOutput(PersistentOStream & os) const; /** * Function used to read in object persistently. * @param is the persistent input stream read from. * @param version the version number of the object when written. */ void persistentInput(PersistentIStream & is, int version); //@} /** * Standard Init function used to initialize the interfaces. */ static void Init(); protected: /** @name Clone Methods. */ //@{ /** * Make a simple clone of this object. * @return a pointer to the new object. */ virtual IBPtr clone() const; /** Make a clone of this object, possibly modifying the cloned object * to make it sane. * @return a pointer to the new object. */ virtual IBPtr fullclone() const; //@} private: /** * Private and non-existent assignment operator. */ ColourReconnector & operator=(const ColourReconnector &) = delete; }; } #endif /* HERWIG_ColourReconnector_H */ diff --git a/src/LHC-MB.in b/src/LHC-MB.in --- a/src/LHC-MB.in +++ b/src/LHC-MB.in @@ -1,98 +1,62 @@ # -*- ThePEG-repository -*- ################################################################################ # This file contains our best tune to UE data from ATLAS at 7 TeV. More recent # tunes and tunes for other centre-of-mass energies as well as more usage # instructions can be obtained from this Herwig wiki page: # http://projects.hepforge.org/herwig/trac/wiki/MB_UE_tunes # The model for soft interactions and diffractions is explained in # [S. Gieseke, P. Kirchgaesser, F. Loshaj, arXiv:1612.04701] ################################################################################ read snippets/PPCollider.in ################################################## # Technical parameters for this run ################################################## cd /Herwig/Generators ################################################## # LHC physics parameters (override defaults here) ################################################## set EventGenerator:EventHandler:LuminosityFunction:Energy 7000.0 # Intrinsic pT tune extrapolated to LHC energy set /Herwig/Shower/ShowerHandler:IntrinsicPtGaussian 2.2*GeV # Minimum Bias read snippets/MB.in -# Read in parameters of the soft model recommended for MB/UE simulations -read snippets/SoftTune.in - - #Diffraction model read snippets/Diffraction.in -#Turn on Baryonic reconnection -read snippets/BaryonicReconnection.in - -# Normalization of the Min bias cross section for correct diffractive cross section -set /Herwig/MatrixElements/MEMinBias:csNorm 0.01 - - -# Use LHC parametrization of the cross section -set /Herwig/UnderlyingEvent/MPIHandler:DLmode 2 - - -#some preliminary parameters for the MPI model which need to be tuned -#TODO -set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.02 -set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.4 -set /Herwig/UnderlyingEvent/MPIHandler:Power 0.308 -set /Herwig/Partons/RemnantDecayer:ladderMult 0.45 -set /Herwig/Partons/RemnantDecayer:ladderbFactor 1.0 -set /Herwig/UnderlyingEvent/MPIHandler:DiffractiveRatio 0.17 - -#ordered pTs -set /Herwig/Partons/RemnantDecayer:PtDistribution 4 - -# set the correct PDFs -read Matchbox/CT14.in - -#set /Herwig/Shower/ShowerHandler:CascadeHandler NULL # Switch off parton shower -#set /Herwig/Hadronization/HadronizationHandler:HadronizationHandler NULL -# do Soft interactions - -#set /Herwig/Hadronization/ColourReconnector:ColourReconnection No ################################################# # Analyses ################################################## -#Comment these lines out in order to use rivet analyses -cd /Herwig/Analysis -create ThePEG::RivetAnalysis RivetAnalysis RivetAnalysis.so - -cd /Herwig/Generators -insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/RivetAnalysis - -#Some example analyses -insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1084540 -insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562 +## Comment these lines out in order to use rivet analyses +# cd /Herwig/Analysis +# create ThePEG::RivetAnalysis RivetAnalysis RivetAnalysis.so + +# cd /Herwig/Generators +# insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/RivetAnalysis +## Some example analyses +# insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1084540 +# insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562 #set /Herwig/Analysis/Plot:EventNumber 54 #cd /Herwig/Generators #insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/Plot #insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile #set /Herwig/Analysis/HepMCFile:PrintEvent 1000000 #set /Herwig/Analysis/HepMCFile:Format GenEvent #set /Herwig/Analysis/HepMCFile:Units GeV_mm #set /Herwig/Analysis/HepMCFile:Filename events.fifo ################################################## # Save run for later usage with 'Herwig run' ################################################## cd /Herwig/Generators saverun LHC-MB EventGenerator diff --git a/src/defaults/Cuts.in b/src/defaults/Cuts.in --- a/src/defaults/Cuts.in +++ b/src/defaults/Cuts.in @@ -1,118 +1,118 @@ # -*- ThePEG-repository -*- ########################################################### # Default cuts (applied to the hard subprocess) # # Don't change values here, re-set them in your own input # files using these as examples. ########################################################### mkdir /Herwig/Matchers cd /Herwig/Matchers create ThePEG::Matcher Lepton create ThePEG::Matcher Neutrino create ThePEG::Matcher ChargedLepton create ThePEG::Matcher LightQuark create ThePEG::Matcher LightParticle create ThePEG::Matcher LightAntiQuark create ThePEG::Matcher StandardQCDParton create ThePEG::Matcher Photon create ThePEG::Matcher Top create ThePEG::Matcher Bottom create ThePEG::Matcher WBoson create ThePEG::Matcher ZBoson create ThePEG::Matcher HiggsBoson mkdir /Herwig/Cuts cd /Herwig/Cuts # cut on jet pt create ThePEG::SimpleKTCut JetKtCut SimpleKTCut.so newdef JetKtCut:Matcher /Herwig/Matchers/StandardQCDParton newdef JetKtCut:MinKT 20.0*GeV # cut on photon create ThePEG::SimpleKTCut PhotonKtCut SimpleKTCut.so newdef PhotonKtCut:Matcher /Herwig/Matchers/Photon newdef PhotonKtCut:MinKT 20.0*GeV newdef PhotonKtCut:MinEta -3. newdef PhotonKtCut:MaxEta 3. # cut on leptons create ThePEG::SimpleKTCut LeptonKtCut SimpleKTCut.so newdef LeptonKtCut:Matcher /Herwig/Matchers/Lepton newdef LeptonKtCut:MinKT 0.0*GeV # cut on charged leptons create ThePEG::SimpleKTCut ChargedLeptonKtCut SimpleKTCut.so newdef ChargedLeptonKtCut:Matcher /Herwig/Matchers/ChargedLepton newdef ChargedLeptonKtCut:MinKT 0.0*GeV # cut on top quarks create ThePEG::KTRapidityCut TopKtCut KTRapidityCut.so newdef TopKtCut:Matcher /Herwig/Matchers/Top newdef TopKtCut:MinKT 0.0*GeV # cut on W bosons create ThePEG::KTRapidityCut WBosonKtCut KTRapidityCut.so newdef WBosonKtCut:Matcher /Herwig/Matchers/WBoson newdef WBosonKtCut:MinKT 0.0*GeV # cut on Z bosons create ThePEG::KTRapidityCut ZBosonKtCut KTRapidityCut.so newdef ZBosonKtCut:Matcher /Herwig/Matchers/ZBoson newdef ZBosonKtCut:MinKT 0.0*GeV # cut on Higgs bosons create ThePEG::KTRapidityCut HiggsBosonKtCut KTRapidityCut.so newdef HiggsBosonKtCut:Matcher /Herwig/Matchers/HiggsBoson newdef HiggsBosonKtCut:MinKT 0.0*GeV # cut on Higgs bosons create ThePEG::KTRapidityCut LightParticleKtCut KTRapidityCut.so newdef LightParticleKtCut:Matcher /Herwig/Matchers/LightParticle newdef LightParticleKtCut:MinKT 0.0*GeV # create a cut on the invariant mass of lepton pairs create ThePEG::V2LeptonsCut MassCut V2LeptonsCut.so newdef MassCut:Families All newdef MassCut:CComb All newdef MassCut:MinM 20.*GeV newdef MassCut:MaxM 14000.*GeV # create a cut on Q^2 for neutral current DIS create ThePEG::SimpleDISCut NeutralCurrentCut SimpleDISCut.so newdef NeutralCurrentCut:MinQ2 20. newdef NeutralCurrentCut:Current Neutral # create a cut on Q^2 for charged current DIS create ThePEG::SimpleDISCut ChargedCurrentCut SimpleDISCut.so newdef ChargedCurrentCut:MinQ2 20. newdef ChargedCurrentCut:Current Charged # create a cut of Q^2 for charged current DIS # Default Cuts object create ThePEG::Cuts Cuts newdef /Herwig/Generators/EventGenerator:EventHandler:Cuts /Herwig/Cuts/Cuts newdef Cuts:MHatMin 20.*GeV # insert into hadron cuts (by defaults) insert Cuts:OneCuts[0] JetKtCut insert Cuts:OneCuts[1] PhotonKtCut insert Cuts:OneCuts[2] LeptonKtCut insert Cuts:OneCuts[3] TopKtCut insert Cuts:OneCuts[4] WBosonKtCut insert Cuts:OneCuts[5] ZBosonKtCut insert Cuts:OneCuts[6] HiggsBosonKtCut insert Cuts:OneCuts[7] ChargedLeptonKtCut insert Cuts:MultiCuts[0] MassCut # create diffrent cuts object for MinBias to avoid numerical problems create ThePEG::Cuts MinBiasCuts newdef MinBiasCuts:ScaleMin 2.0*GeV2 -newdef MinBiasCuts:X1Min 0.055 -newdef MinBiasCuts:X2Min 0.055 +newdef MinBiasCuts:X1Min 0.011 +newdef MinBiasCuts:X2Min 0.011 newdef MinBiasCuts:MHatMin 0.0*GeV diff --git a/src/defaults/Hadronization.in b/src/defaults/Hadronization.in --- a/src/defaults/Hadronization.in +++ b/src/defaults/Hadronization.in @@ -1,95 +1,96 @@ # -*- ThePEG-repository -*- ############################################################ # Setup of default hadronization # # There are no user servicable parts inside. # # Anything that follows below should only be touched if you # know what you're doing. ############################################################# cd /Herwig/Particles create ThePEG::ParticleData Cluster setup Cluster 81 Cluster 0.00990 0.0 0.0 0.0 0 0 0 1 create ThePEG::ParticleData Remnant setup Remnant 82 Remnant 0.00990 0.0 0.0 0.0 0 0 0 1 mkdir /Herwig/Hadronization cd /Herwig/Hadronization create Herwig::ClusterHadronizationHandler ClusterHadHandler create Herwig::PartonSplitter PartonSplitter create Herwig::ClusterFinder ClusterFinder create Herwig::ColourReconnector ColourReconnector create Herwig::ClusterFissioner ClusterFissioner create Herwig::LightClusterDecayer LightClusterDecayer create Herwig::ClusterDecayer ClusterDecayer create Herwig::HwppSelector HadronSelector newdef ClusterHadHandler:PartonSplitter PartonSplitter newdef ClusterHadHandler:ClusterFinder ClusterFinder newdef ClusterHadHandler:ColourReconnector ColourReconnector newdef ClusterHadHandler:ClusterFissioner ClusterFissioner newdef ClusterHadHandler:LightClusterDecayer LightClusterDecayer newdef ClusterHadHandler:ClusterDecayer ClusterDecayer newdef ClusterHadHandler:MinVirtuality2 0.1*GeV2 newdef ClusterHadHandler:MaxDisplacement 1.0e-10*millimeter newdef ClusterHadHandler:UnderlyingEventHandler NULL newdef ClusterFissioner:HadronSelector HadronSelector newdef LightClusterDecayer:HadronSelector HadronSelector newdef ClusterDecayer:HadronSelector HadronSelector newdef ColourReconnector:ColourReconnection Yes -newdef ColourReconnector:ReconnectionProbability 0.652710 -newdef ColourReconnector:Algorithm Plain -newdef ColourReconnector:InitialTemperature 0.01 -newdef ColourReconnector:AnnealingFactor 0.21 -newdef ColourReconnector:AnnealingSteps 10 -newdef ColourReconnector:TriesPerStepFactor 0.66 +newdef ColourReconnector:ReconnectionProbabilityBaryonic 0.57 +newdef ColourReconnector:ReconnectionProbability 0.89 +newdef ColourReconnector:Algorithm Baryonic newdef ColourReconnector:OctetTreatment All # Clustering parameters for light quarks newdef ClusterFissioner:ClMaxLight 3.649 newdef ClusterFissioner:ClPowLight 2.780 newdef ClusterFissioner:PSplitLight 0.899 newdef ClusterDecayer:ClDirLight 1 newdef ClusterDecayer:ClSmrLight 0.78 # Clustering parameters for b-quarks newdef ClusterFissioner:ClMaxBottom 3.757 newdef ClusterFissioner:ClPowBottom 0.547 newdef ClusterFissioner:PSplitBottom 0.625 newdef ClusterDecayer:ClDirBottom 1 newdef ClusterDecayer:ClSmrBottom 0.078 newdef HadronSelector:SingleHadronLimitBottom 0.000 # Clustering parameters for c-quarks newdef ClusterFissioner:ClMaxCharm 3.950 newdef ClusterFissioner:ClPowCharm 2.559 newdef ClusterFissioner:PSplitCharm 0.994 newdef ClusterDecayer:ClDirCharm 1 newdef ClusterDecayer:ClSmrCharm 0.163 newdef HadronSelector:SingleHadronLimitCharm 0.000 # Clustering parameters for exotic quarks # (e.g. hadronizing Susy particles) newdef ClusterFissioner:ClMaxExotic 2.7*GeV newdef ClusterFissioner:ClPowExotic 1.46 newdef ClusterFissioner:PSplitExotic 1.00 newdef ClusterDecayer:ClDirExotic 1 newdef ClusterDecayer:ClSmrExotic 0. newdef HadronSelector:SingleHadronLimitExotic 0. # +newdef PartonSplitter:SplitPwtSquark 0.824135 +newdef PartonSplitter:Split uds + +# newdef HadronSelector:PwtDquark 1.0 newdef HadronSelector:PwtUquark 1.0 -newdef HadronSelector:PwtSquark 0.700 +newdef HadronSelector:PwtSquark 0.291717 newdef HadronSelector:PwtCquark 1.0 newdef HadronSelector:PwtBquark 1.0 newdef HadronSelector:PwtDIquark 0.298 newdef HadronSelector:SngWt 0.74 newdef HadronSelector:DecWt 0.62 newdef HadronSelector:Mode 1 newdef HadronSelector:BelowThreshold All diff --git a/src/defaults/MatrixElements.in b/src/defaults/MatrixElements.in --- a/src/defaults/MatrixElements.in +++ b/src/defaults/MatrixElements.in @@ -1,249 +1,250 @@ # -*- 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 ############################################################ # 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 diff --git a/src/defaults/Partons.in b/src/defaults/Partons.in --- a/src/defaults/Partons.in +++ b/src/defaults/Partons.in @@ -1,36 +1,37 @@ # -*- ThePEG-repository -*- mkdir /Herwig/Partons cd /Herwig/Partons create ThePEG::NoRemnants NoRemnants create ThePEG::NoPDF NoPDF newdef NoPDF:RemnantHandler NoRemnants create ThePEG::UnResolvedRemnant LeptonRemnants UnResolvedRemnant.so create Herwig::HwRemDecayer RemnantDecayer HwShower.so newdef RemnantDecayer:PomeronStructure Gluon -newdef RemnantDecayer:colourDisrupt 0.7983 +newdef RemnantDecayer:colourDisrupt 0.0 newdef RemnantDecayer:ladderPower -0.088983 newdef RemnantDecayer:ladderNorm 1.086029 +newdef RemnantDecayer:PtDistribution 4 create ThePEG::SoftRemnantHandler HadronRemnants newdef HadronRemnants:RemnantDecayer RemnantDecayer read PDF.in # Make particle extractors create ThePEG::PartonExtractor EEExtractor newdef EEExtractor:NoPDF NoPDF cp EEExtractor PPExtractor newdef PPExtractor:FlatSHatY 1 newdef PPExtractor:FirstPDF HardLOPDF newdef PPExtractor:SecondPDF HardLOPDF cp EEExtractor EPExtractor newdef EPExtractor:FlatSHatY 1 newdef PPExtractor:FirstPDF HardLOPDF newdef PPExtractor:SecondPDF HardLOPDF newdef EPExtractor:SecondPDF HardLOPDF diff --git a/src/defaults/Shower.in b/src/defaults/Shower.in --- a/src/defaults/Shower.in +++ b/src/defaults/Shower.in @@ -1,321 +1,321 @@ # -*- ThePEG-repository -*- ############################################################ # Setup of default parton shower # # Useful switches for users are marked near the top of # this file. # # Don't edit this file directly, but reset the switches # in your own input files! ############################################################ library HwMPI.so library HwShower.so library HwMatching.so mkdir /Herwig/Shower cd /Herwig/Shower create Herwig::QTildeShowerHandler ShowerHandler newdef ShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler newdef ShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer # use LO PDFs for Shower, can be changed later newdef ShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF newdef ShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF newdef ShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF newdef ShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF ##################################### # initial setup, don't change these! ##################################### create Herwig::SplittingGenerator SplittingGenerator create Herwig::ShowerAlphaQCD AlphaQCD create Herwig::ShowerAlphaQED AlphaQED create Herwig::PartnerFinder PartnerFinder newdef PartnerFinder:PartnerMethod 0 newdef PartnerFinder:ScaleChoice 0 create Herwig::KinematicsReconstructor KinematicsReconstructor newdef KinematicsReconstructor:ReconstructionOption Colour3 newdef KinematicsReconstructor:InitialStateReconOption SofterFraction newdef KinematicsReconstructor:InitialInitialBoostOption LongTransBoost newdef KinematicsReconstructor:FinalFinalWeight Yes newdef /Herwig/Partons/RemnantDecayer:AlphaS AlphaQCD newdef /Herwig/Partons/RemnantDecayer:AlphaEM AlphaQED newdef ShowerHandler:PartnerFinder PartnerFinder newdef ShowerHandler:KinematicsReconstructor KinematicsReconstructor newdef ShowerHandler:SplittingGenerator SplittingGenerator newdef ShowerHandler:SpinCorrelations Yes newdef ShowerHandler:SoftCorrelations Singular ################################################################## # Intrinsic pT # # Recommended: # 1.9 GeV for Tevatron W/Z production. # 2.1 GeV for LHC W/Z production at 10 TeV # 2.2 GeV for LHC W/Z production at 14 TeV # # Set all parameters to 0 to disable ################################################################## -newdef ShowerHandler:IntrinsicPtGaussian 1.3*GeV +newdef ShowerHandler:IntrinsicPtGaussian 2.2*GeV newdef ShowerHandler:IntrinsicPtBeta 0 newdef ShowerHandler:IntrinsicPtGamma 0*GeV newdef ShowerHandler:IntrinsicPtIptmax 0*GeV ############################################################# # Set up truncated shower handler. ############################################################# create Herwig::PowhegShowerHandler PowhegShowerHandler set PowhegShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set PowhegShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer newdef PowhegShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF newdef PowhegShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF newdef PowhegShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF newdef PowhegShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF newdef PowhegShowerHandler:MPIHandler /Herwig/UnderlyingEvent/MPIHandler newdef PowhegShowerHandler:RemDecayer /Herwig/Partons/RemnantDecayer newdef PowhegShowerHandler:PDFA /Herwig/Partons/ShowerLOPDF newdef PowhegShowerHandler:PDFB /Herwig/Partons/ShowerLOPDF newdef PowhegShowerHandler:PDFARemnant /Herwig/Partons/RemnantPDF newdef PowhegShowerHandler:PDFBRemnant /Herwig/Partons/RemnantPDF newdef PowhegShowerHandler:PartnerFinder PartnerFinder newdef PowhegShowerHandler:KinematicsReconstructor KinematicsReconstructor newdef PowhegShowerHandler:SplittingGenerator SplittingGenerator newdef PowhegShowerHandler:Interactions QCDandQED newdef PowhegShowerHandler:SpinCorrelations Yes newdef PowhegShowerHandler:SoftCorrelations Singular -newdef PowhegShowerHandler:IntrinsicPtGaussian 1.3*GeV +newdef PowhegShowerHandler:IntrinsicPtGaussian 2.2*GeV newdef PowhegShowerHandler:IntrinsicPtBeta 0 newdef PowhegShowerHandler:IntrinsicPtGamma 0*GeV newdef PowhegShowerHandler:IntrinsicPtIptmax 0*GeV newdef PowhegShowerHandler:EvolutionScheme DotProduct ############################################################# # End of interesting user servicable section. # # Anything that follows below should only be touched if you # know what you're doing. # # Really. ############################################################# # # a few default values newdef ShowerHandler:MECorrMode 1 newdef ShowerHandler:EvolutionScheme DotProduct newdef AlphaQCD:ScaleFactor 1.0 newdef AlphaQCD:NPAlphaS 2 newdef AlphaQCD:Qmin 0.935 newdef AlphaQCD:NumberOfLoops 2 newdef AlphaQCD:AlphaIn 0.1186 # # # Lets set up all the splittings create Herwig::HalfHalfOneSplitFn QtoQGammaSplitFn set QtoQGammaSplitFn:InteractionType QED set QtoQGammaSplitFn:ColourStructure ChargedChargedNeutral set QtoQGammaSplitFn:AngularOrdered Yes set QtoQGammaSplitFn:StrictAO Yes create Herwig::HalfHalfOneSplitFn QtoQGSplitFn newdef QtoQGSplitFn:InteractionType QCD newdef QtoQGSplitFn:ColourStructure TripletTripletOctet set QtoQGSplitFn:AngularOrdered Yes set QtoQGSplitFn:StrictAO Yes create Herwig::OneOneOneSplitFn GtoGGSplitFn newdef GtoGGSplitFn:InteractionType QCD newdef GtoGGSplitFn:ColourStructure OctetOctetOctet set GtoGGSplitFn:AngularOrdered Yes set GtoGGSplitFn:StrictAO Yes create Herwig::OneOneOneMassiveSplitFn WtoWGammaSplitFn newdef WtoWGammaSplitFn:InteractionType QED newdef WtoWGammaSplitFn:ColourStructure ChargedChargedNeutral set WtoWGammaSplitFn:AngularOrdered Yes set WtoWGammaSplitFn:StrictAO Yes create Herwig::OneHalfHalfSplitFn GtoQQbarSplitFn newdef GtoQQbarSplitFn:InteractionType QCD newdef GtoQQbarSplitFn:ColourStructure OctetTripletTriplet set GtoQQbarSplitFn:AngularOrdered Yes set GtoQQbarSplitFn:StrictAO Yes create Herwig::OneHalfHalfSplitFn GammatoQQbarSplitFn newdef GammatoQQbarSplitFn:InteractionType QED newdef GammatoQQbarSplitFn:ColourStructure NeutralChargedCharged set GammatoQQbarSplitFn:AngularOrdered Yes set GammatoQQbarSplitFn:StrictAO Yes create Herwig::HalfOneHalfSplitFn QtoGQSplitFn newdef QtoGQSplitFn:InteractionType QCD newdef QtoGQSplitFn:ColourStructure TripletOctetTriplet set QtoGQSplitFn:AngularOrdered Yes set QtoGQSplitFn:StrictAO Yes create Herwig::HalfOneHalfSplitFn QtoGammaQSplitFn newdef QtoGammaQSplitFn:InteractionType QED newdef QtoGammaQSplitFn:ColourStructure ChargedNeutralCharged set QtoGammaQSplitFn:AngularOrdered Yes set QtoGammaQSplitFn:StrictAO Yes # # Now the Sudakovs create Herwig::PTCutOff PTCutOff newdef PTCutOff:pTmin 0.958*GeV create Herwig::SudakovFormFactor SudakovCommon newdef SudakovCommon:Alpha AlphaQCD newdef SudakovCommon:Cutoff PTCutOff newdef SudakovCommon:PDFmax 1.0 cp SudakovCommon QtoQGSudakov newdef QtoQGSudakov:SplittingFunction QtoQGSplitFn newdef QtoQGSudakov:PDFmax 1.9 cp SudakovCommon QtoQGammaSudakov set QtoQGammaSudakov:SplittingFunction QtoQGammaSplitFn set QtoQGammaSudakov:Alpha AlphaQED set QtoQGammaSudakov:PDFmax 1.9 cp QtoQGammaSudakov LtoLGammaSudakov cp PTCutOff LtoLGammaPTCutOff # Technical parameter to stop evolution. set LtoLGammaPTCutOff:pTmin 0.000001 set LtoLGammaSudakov:Cutoff LtoLGammaPTCutOff cp SudakovCommon GtoGGSudakov newdef GtoGGSudakov:SplittingFunction GtoGGSplitFn newdef GtoGGSudakov:PDFmax 2.0 cp SudakovCommon WtoWGammaSudakov newdef WtoWGammaSudakov:SplittingFunction WtoWGammaSplitFn set WtoWGammaSudakov:Alpha AlphaQED cp SudakovCommon GtoQQbarSudakov newdef GtoQQbarSudakov:SplittingFunction GtoQQbarSplitFn newdef GtoQQbarSudakov:PDFmax 120.0 cp SudakovCommon GammatoQQbarSudakov newdef GammatoQQbarSudakov:SplittingFunction GammatoQQbarSplitFn set GammatoQQbarSudakov:Alpha AlphaQED newdef GammatoQQbarSudakov:PDFmax 120.0 cp SudakovCommon GtobbbarSudakov newdef GtobbbarSudakov:SplittingFunction GtoQQbarSplitFn newdef GtobbbarSudakov:PDFmax 40000.0 cp SudakovCommon GtoccbarSudakov newdef GtoccbarSudakov:SplittingFunction GtoQQbarSplitFn newdef GtoccbarSudakov:PDFmax 2000.0 cp SudakovCommon QtoGQSudakov newdef QtoGQSudakov:SplittingFunction QtoGQSplitFn cp SudakovCommon QtoGammaQSudakov newdef QtoGammaQSudakov:SplittingFunction QtoGammaQSplitFn set QtoGammaQSudakov:Alpha AlphaQED cp SudakovCommon utoGuSudakov newdef utoGuSudakov:SplittingFunction QtoGQSplitFn newdef utoGuSudakov:PDFFactor OverOneMinusZ newdef utoGuSudakov:PDFmax 5.0 cp SudakovCommon dtoGdSudakov newdef dtoGdSudakov:SplittingFunction QtoGQSplitFn newdef dtoGdSudakov:PDFFactor OverOneMinusZ # # Now add the final splittings # do SplittingGenerator:AddFinalSplitting u->u,g; QtoQGSudakov do SplittingGenerator:AddFinalSplitting d->d,g; QtoQGSudakov do SplittingGenerator:AddFinalSplitting s->s,g; QtoQGSudakov do SplittingGenerator:AddFinalSplitting c->c,g; QtoQGSudakov do SplittingGenerator:AddFinalSplitting b->b,g; QtoQGSudakov do SplittingGenerator:AddFinalSplitting t->t,g; QtoQGSudakov # do SplittingGenerator:AddFinalSplitting g->g,g; GtoGGSudakov # do SplittingGenerator:AddFinalSplitting g->u,ubar; GtoQQbarSudakov do SplittingGenerator:AddFinalSplitting g->d,dbar; GtoQQbarSudakov do SplittingGenerator:AddFinalSplitting g->s,sbar; GtoQQbarSudakov do SplittingGenerator:AddFinalSplitting g->c,cbar; GtoccbarSudakov do SplittingGenerator:AddFinalSplitting g->b,bbar; GtobbbarSudakov do SplittingGenerator:AddFinalSplitting g->t,tbar; GtoQQbarSudakov # do SplittingGenerator:AddFinalSplitting gamma->u,ubar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->d,dbar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->s,sbar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->c,cbar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->b,bbar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->t,tbar; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->e-,e+; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->mu-,mu+; GammatoQQbarSudakov do SplittingGenerator:AddFinalSplitting gamma->tau-,tau+; GammatoQQbarSudakov # do SplittingGenerator:AddFinalSplitting u->u,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting d->d,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting s->s,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting c->c,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:AddFinalSplitting e-->e-,gamma; LtoLGammaSudakov do SplittingGenerator:AddFinalSplitting mu-->mu-,gamma; LtoLGammaSudakov do SplittingGenerator:AddFinalSplitting tau-->tau-,gamma; LtoLGammaSudakov do SplittingGenerator:AddFinalSplitting W+->W+,gamma; WtoWGammaSudakov # # Now lets add the initial splittings. Remember the form a->b,c; means # that the current particle b is given and we backward branch to new # particle a which is initial state and new particle c which is final state # do SplittingGenerator:AddInitialSplitting u->u,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting d->d,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting s->s,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting c->c,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting b->b,g; QtoQGSudakov do SplittingGenerator:AddInitialSplitting u->u,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting d->d,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting s->s,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting c->c,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting b->b,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting t->t,gamma; QtoQGammaSudakov do SplittingGenerator:AddInitialSplitting g->g,g; GtoGGSudakov # do SplittingGenerator:AddInitialSplitting g->d,dbar; GtoQQbarSudakov do SplittingGenerator:AddInitialSplitting g->u,ubar; GtoQQbarSudakov do SplittingGenerator:AddInitialSplitting g->s,sbar; GtoQQbarSudakov do SplittingGenerator:AddInitialSplitting g->c,cbar; GtoccbarSudakov do SplittingGenerator:AddInitialSplitting g->b,bbar; GtobbbarSudakov # do SplittingGenerator:AddInitialSplitting gamma->d,dbar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting gamma->u,ubar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting gamma->s,sbar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting gamma->c,cbar; GammatoQQbarSudakov do SplittingGenerator:AddInitialSplitting gamma->b,bbar; GammatoQQbarSudakov # do SplittingGenerator:AddInitialSplitting d->g,d; dtoGdSudakov do SplittingGenerator:AddInitialSplitting u->g,u; utoGuSudakov do SplittingGenerator:AddInitialSplitting s->g,s; QtoGQSudakov do SplittingGenerator:AddInitialSplitting c->g,c; QtoGQSudakov do SplittingGenerator:AddInitialSplitting b->g,b; QtoGQSudakov do SplittingGenerator:AddInitialSplitting dbar->g,dbar; dtoGdSudakov do SplittingGenerator:AddInitialSplitting ubar->g,ubar; utoGuSudakov do SplittingGenerator:AddInitialSplitting sbar->g,sbar; QtoGQSudakov do SplittingGenerator:AddInitialSplitting cbar->g,cbar; QtoGQSudakov do SplittingGenerator:AddInitialSplitting bbar->g,bbar; QtoGQSudakov # do SplittingGenerator:AddInitialSplitting d->gamma,d; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting u->gamma,u; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting s->gamma,s; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting c->gamma,c; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting b->gamma,b; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting dbar->gamma,dbar; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting ubar->gamma,ubar; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting sbar->gamma,sbar; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting cbar->gamma,cbar; QtoGammaQSudakov do SplittingGenerator:AddInitialSplitting bbar->gamma,bbar; QtoGammaQSudakov diff --git a/src/defaults/UnderlyingEvent.in b/src/defaults/UnderlyingEvent.in --- a/src/defaults/UnderlyingEvent.in +++ b/src/defaults/UnderlyingEvent.in @@ -1,75 +1,82 @@ # -*- ThePEG-repository -*- ################################################ # Set up the handler for the underlying event. ################################################ mkdir /Herwig/UnderlyingEvent cd /Herwig/UnderlyingEvent ######## cuts ################################################# library SimpleKTCut.so # cut on pt. Without a specific matcher object, it works on all # particles create ThePEG::SimpleKTCut KtCut newdef KtCut:MinKT 3.06 # create the cuts object for the Underlying Event create ThePEG::Cuts UECuts # This should always be 2*MinKT!! newdef UECuts:MHatMin 6.12 insert UECuts:OneCuts 0 KtCut ######## subprocess ########################################### create ThePEG::SubProcessHandler FastQCD create Herwig::MEQCD2to2Fast MEQCD2to2Fast HwMEHadronFast.so insert FastQCD:MatrixElements 0 MEQCD2to2Fast cp /Herwig/Partons/PPExtractor /Herwig/Partons/MPIExtractor newdef /Herwig/Partons/MPIExtractor:FlatSHatY 1 newdef FastQCD:PartonExtractor /Herwig/Partons/MPIExtractor # always use LO PDF for MPI and remnants newdef /Herwig/Partons/MPIExtractor:FirstPDF /Herwig/Partons/MPIPDF newdef /Herwig/Partons/MPIExtractor:SecondPDF /Herwig/Partons/MPIPDF ######## MPI Handler ########################################## library JetCuts.so library HwMPI.so create Herwig::MPIHandler MPIHandler # set the subprocesses and corresponding cuts insert MPIHandler:SubProcessHandlers 0 FastQCD insert MPIHandler:Cuts 0 UECuts # energy extrapolation for the MinKT cut newdef MPIHandler:EnergyExtrapolation Power -newdef MPIHandler:Power 0.420445 -newdef MPIHandler:pTmin0 3.568157 +newdef MPIHandler:Power 0.322 +newdef MPIHandler:pTmin0 2.55 newdef MPIHandler:ReferenceScale 7000.*GeV newdef MPIHandler:EEparamA 0.585*GeV newdef MPIHandler:EEparamB 37.45*GeV +# Paramters for the ladders for soft MPIs +newdef /Herwig/Partons/RemnantDecayer:ladderMult 0.466 +newdef /Herwig/Partons/RemnantDecayer:ladderbFactor 0.46 + +# Ratio of diffractive events in MinBias runs +newdef MPIHandler:DiffractiveRatio 0.133 + # The inverse hadron radius -newdef MPIHandler:InvRadius 1.489997 +newdef MPIHandler:InvRadius 1.58 # Set the details of the soft model # Flag to decide whether additional soft interactions (i.e. pt < ptmin) # should be simulated at all newdef MPIHandler:softInt Yes # Flag to decide whether to use the 2 component model. In this model, # the hadron radius (or better: mean parton-parton separation) can be # different for hard and soft scatters. The radius for the soft scatters # is computed during the run startup so that the elastic t-slope, B_el, # coincides with current measurements (for CMenergy < 1.8 TeV) or # Donnachie-Landshoff pomeron fit (CMenergy > 1.8 TeV). If this model is # disabled one has to take care that this observable is correctly # described by adjusting InvRadius properly. newdef MPIHandler:twoComp Yes newdef MPIHandler:DLmode 2 diff --git a/src/defaults/quarks.in b/src/defaults/quarks.in --- a/src/defaults/quarks.in +++ b/src/defaults/quarks.in @@ -1,35 +1,35 @@ # -*- ThePEG-repository -*- # # file containing the particle data for the quarks # create ThePEG::ConstituentParticleData d setup d 1 d 0.005 0. 0 0 -1 3 2 1 0.325 create ThePEG::ConstituentParticleData dbar setup dbar -1 dbar 0.005 0. 0 0 1 -3 2 1 0.325 makeanti dbar d create ThePEG::ConstituentParticleData ubar setup ubar -2 ubar 0.0023 0. 0 0 -2 -3 2 1 0.325 create ThePEG::ConstituentParticleData u setup u 2 u 0.0023 0. 0 0 2 3 2 1 0.325 makeanti u ubar create ThePEG::ConstituentParticleData c setup c 4 c 1.25 0. 0 0 2 3 2 1 1.6 create ThePEG::ConstituentParticleData cbar setup cbar -4 cbar 1.25 0. 0 0 -2 -3 2 1 1.6 makeanti cbar c create ThePEG::ConstituentParticleData sbar -setup sbar -3 sbar 0.095 0. 0 0 1 -3 2 1 0.5 +setup sbar -3 sbar 0.095 0. 0 0 1 -3 2 1 0.45 create ThePEG::ConstituentParticleData s -setup s 3 s 0.095 0. 0 0 -1 3 2 1 0.5 +setup s 3 s 0.095 0. 0 0 -1 3 2 1 0.45 makeanti s sbar create ThePEG::ConstituentParticleData b setup b 5 b 4.2 0. 0 0 -1 3 2 1 5 create ThePEG::ConstituentParticleData bbar setup bbar -5 bbar 4.2 0. 0 0 1 -3 2 1 5 makeanti bbar b create ThePEG::ParticleData t setup t 6 t 174.2 1.4 14 0 2 3 2 0 create ThePEG::ParticleData tbar setup tbar -6 tbar 174.2 1.4 14 0 -2 -3 2 0 makeanti tbar t diff --git a/src/snippets/BaryonicReconnection.in b/src/snippets/BaryonicReconnection.in --- a/src/snippets/BaryonicReconnection.in +++ b/src/snippets/BaryonicReconnection.in @@ -1,17 +1,10 @@ # Set strange quark mass to 0.45 in order to allow alternative gluon splitting set /Herwig/Particles/s:ConstituentMass 0.45*GeV set /Herwig/Particles/sbar:ConstituentMass 0.45*GeV # Use Baryonic Colour Reconnection Model set /Herwig/Hadronization/ColourReconnector:Algorithm Baryonic # Allow alternative gluon splitting set /Herwig/Hadronization/PartonSplitter:Split uds -# Parameters for the Baryonic Reconnection Model -set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.772606 -set /Herwig/Hadronization/ColourReconnector:ReconnectionProbabilityBaryonic 0.477612 -set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.053252 -set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.282032 -set /Herwig/Hadronization/HadronSelector:PwtSquark 0.291717 -set /Herwig/Hadronization/PartonSplitter:SplitPwtSquark 0.824135 diff --git a/src/snippets/SoftTune.in b/src/snippets/SoftTune.in --- a/src/snippets/SoftTune.in +++ b/src/snippets/SoftTune.in @@ -1,8 +1,16 @@ -# Recommended set of parameters for MB/UE runs +# Tund for UE and MB (these are the default parameters) -set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.5 -set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 3.502683 -set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.402055 -set /Herwig/UnderlyingEvent/MPIHandler:Power 0.416852 -set /Herwig/Partons/RemnantDecayer:ladderPower -0.08 -set /Herwig/Partons/RemnantDecayer:ladderNorm 0.95 +set /Herwig/MatrixElements/MEMinBias:csNorm 0.01 +set /Herwig/UnderlyingEvent/MPIHandler:DLmode 2 +set /Herwig/UnderlyingEvent/MPIHandler:pTmin0 2.55 +set /Herwig/UnderlyingEvent/MPIHandler:InvRadius 1.58 +set /Herwig/UnderlyingEvent/MPIHandler:Power 0.322 +set /Herwig/Partons/RemnantDecayer:ladderMult 0.466 +set /Herwig/Partons/RemnantDecayer:ladderbFactor 0.46 +set /Herwig/UnderlyingEvent/MPIHandler:DiffractiveRatio 0.133 +set /Herwig/Partons/RemnantDecayer:PtDistribution 4 +set /Herwig/Hadronization/ColourReconnector:ReconnectionProbability 0.89 +set /Herwig/Hadronization/ColourReconnector:ReconnectionProbabilityBaryonic 0.57 +set /Herwig/Hadronization/HadronSelector:PwtSquark 0.291717 +set /Herwig/Hadronization/PartonSplitter:SplitPwtSquark 0.824135 +