diff --git a/PDF/Makefile.am b/PDF/Makefile.am --- a/PDF/Makefile.am +++ b/PDF/Makefile.am @@ -1,64 +1,65 @@ EXTRA_DIST = diffraction pkglib_LTLIBRARIES = HwPomeronPDF.la HwPomeronPDF_la_SOURCES = \ PomeronPDF.cc PomeronPDF.h ## add this to produce tests of the PDFs ## HwDIFFRACTIVEPDF_la_CPPFLAGS=$(AM_CPPFLAGS) -DDIFFRACTIVEPDF_TESTING HwPomeronPDF_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 4:0:0 pkglib_LTLIBRARIES += HwReggeonPDF.la HwReggeonPDF_la_SOURCES = \ ReggeonPDF.cc ReggeonPDF.h ## add this to produce tests of the PDFs ## HwDIFFRACTIVEPDF_la_CPPFLAGS=$(AM_CPPFLAGS) -DDIFFRACTIVEPDF_TESTING HwReggeonPDF_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 4:0:0 pkglib_LTLIBRARIES += HwPomeronFlux.la HwPomeronFlux_la_SOURCES = \ PomeronFlux.h PomeronFlux.cc HwPomeronFlux_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 5:0:0 ## bound into HwShower.la noinst_LTLIBRARIES = libHwRemDecayer.la libHwRemDecayer_la_SOURCES = \ HwRemDecayer.h HwRemDecayer.cc HwRemDecayer.fh ## bound into HwShower.la noinst_LTLIBRARIES += libHwMPIPDF.la libHwMPIPDF_la_SOURCES = \ MPIPDF.h MPIPDF.cc MPIPDF.fh \ +MPIPDF2.h MPIPDF2.cc MPIPDF2.fh \ MinBiasPDF.h MinBiasPDF.cc MinBiasPDF.fh pkglib_LTLIBRARIES += HwSatPDF.la HwSatPDF_la_SOURCES = \ SatPDF.h SatPDF.cc HwSatPDF_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 7:0:0 pkglib_LTLIBRARIES += HwPartonExtractor.la HwPartonExtractor_la_SOURCES = \ MultiPartonExtractor.h MultiPartonExtractor.cc HwPartonExtractor_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 7:0:0 pkglib_LTLIBRARIES += HwIncomingPhotonEvolver.la HwIncomingPhotonEvolver_la_SOURCES = \ IncomingPhotonEvolver.h IncomingPhotonEvolver.cc HwIncomingPhotonEvolver_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 7:0:0 pkglib_LTLIBRARIES += HwSaSPhotonPDF.la HwSaSPhotonPDF_la_SOURCES = \ SaSGamma.f SaSPhotonPDF.cc SaSPhotonPDF.h HwSaSPhotonPDF_la_LDFLAGS = $(AM_LDFLAGS) -module -version-info 2:0:0 install-data-local: for i in `find $(srcdir)/diffraction -name '*.data'`; \ do \ $(install_sh_DATA) $$i $(DESTDIR)$(pkgdatadir)/PDF/$${i#$(srcdir)/}; \ done uninstall-local: rm -rf $(DESTDIR)$(pkgdatadir)/PDF diff --git a/src/Matchbox/LO-DipoleShower.in.in b/src/Matchbox/LO-DipoleShower.in.in --- a/src/Matchbox/LO-DipoleShower.in.in +++ b/src/Matchbox/LO-DipoleShower.in.in @@ -1,24 +1,24 @@ # -*- ThePEG-repository -*- read Matchbox/LO.in cd /Herwig/EventHandlers set EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in cd /Herwig/Analysis @CREATE_RIVET@ ThePEG::RivetAnalysis Rivet RivetAnalysis.so @CREATE_HEPMC@ ThePEG::HepMCFile HepMC HepMCAnalysis.so cd /Herwig/Particles do c:UnsetHardProcessMass do cbar:UnsetHardProcessMass set c:NominalMass 0*GeV set cbar:NominalMass 0*GeV cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:StrictFlavourScheme Yes diff --git a/src/Matchbox/MCatLO-DipoleShower.in.in b/src/Matchbox/MCatLO-DipoleShower.in.in --- a/src/Matchbox/MCatLO-DipoleShower.in.in +++ b/src/Matchbox/MCatLO-DipoleShower.in.in @@ -1,28 +1,28 @@ # -*- ThePEG-repository -*- read Matchbox/LO.in cd /Herwig/EventHandlers set EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in set /Herwig/DipoleShower/DipoleShowerHandler:RestrictPhasespace Yes set /Herwig/DipoleShower/DipoleShowerHandler:MaxPtIsMuF Yes set /Herwig/DipoleShower/DipoleShowerHandler:HardScaleProfile /Herwig/MatrixElements/Matchbox/HardScaleProfile cd /Herwig/Analysis @CREATE_RIVET@ ThePEG::RivetAnalysis Rivet RivetAnalysis.so @CREATE_HEPMC@ ThePEG::HepMCFile HepMC HepMCAnalysis.so cd /Herwig/Particles do c:UnsetHardProcessMass do cbar:UnsetHardProcessMass set c:NominalMass 0*GeV set cbar:NominalMass 0*GeV cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:StrictFlavourScheme Yes diff --git a/src/Matchbox/MCatNLO-DipoleShower.in.in b/src/Matchbox/MCatNLO-DipoleShower.in.in --- a/src/Matchbox/MCatNLO-DipoleShower.in.in +++ b/src/Matchbox/MCatNLO-DipoleShower.in.in @@ -1,32 +1,32 @@ # -*- ThePEG-repository -*- cd /Herwig/EventHandlers set EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler set /Herwig/DipoleShower/DipoleShowerHandler:RestrictPhasespace Yes set /Herwig/DipoleShower/DipoleShowerHandler:MaxPtIsMuF Yes set /Herwig/DipoleShower/DipoleShowerHandler:HardScaleProfile /Herwig/MatrixElements/Matchbox/HardScaleProfile cd /Herwig/MatrixElements/Matchbox set Factory:ShowerApproximation DipoleMatching # tunes -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in cd /Herwig/Analysis @CREATE_RIVET@ ThePEG::RivetAnalysis Rivet RivetAnalysis.so @CREATE_HEPMC@ ThePEG::HepMCFile HepMC HepMCAnalysis.so cd /Herwig/Particles do c:UnsetHardProcessMass do cbar:UnsetHardProcessMass set c:NominalMass 0*GeV set cbar:NominalMass 0*GeV cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:StrictFlavourScheme Yes diff --git a/src/Matchbox/Powheg-DipoleShower.in.in b/src/Matchbox/Powheg-DipoleShower.in.in --- a/src/Matchbox/Powheg-DipoleShower.in.in +++ b/src/Matchbox/Powheg-DipoleShower.in.in @@ -1,30 +1,30 @@ # -*- ThePEG-repository -*- read Matchbox/Powheg.in set /Herwig/MatrixElements/Matchbox/MEMatching:TruncatedShower No cd /Herwig/EventHandlers set EventHandler:CascadeHandler /Herwig/DipoleShower/DipoleShowerHandler set /Herwig/DipoleShower/DipoleShowerHandler:RestrictPhasespace Yes set /Herwig/DipoleShower/DipoleShowerHandler:MaxPtIsMuF Yes set /Herwig/DipoleShower/DipoleShowerHandler:HardScaleProfile NULL # tunes -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in cd /Herwig/Analysis @CREATE_RIVET@ ThePEG::RivetAnalysis Rivet RivetAnalysis.so @CREATE_HEPMC@ ThePEG::HepMCFile HepMC HepMCAnalysis.so cd /Herwig/Particles do c:UnsetHardProcessMass do cbar:UnsetHardProcessMass set c:NominalMass 0*GeV set cbar:NominalMass 0*GeV cd /Herwig/UnderlyingEvent set MEQCD2to2Fast:StrictFlavourScheme Yes diff --git a/src/Merging/LEP-Merging.in b/src/Merging/LEP-Merging.in --- a/src/Merging/LEP-Merging.in +++ b/src/Merging/LEP-Merging.in @@ -1,155 +1,155 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/EECollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 91.2*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 0 set MergingFactory:OrderInAlphaEW 2 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process e- e+ -> j j [ j j ] set MergingFactory:NLOProcesses 2 set Merger:MergingScale 4.*GeV set Merger:MergingScaleSmearing 0.1 cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production # read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in # read Matchbox/OnShellHProduction.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package ## As massive b-quarks are currently not supported by the ## build in ME an external ME-Provider is needed (e.g. MG+OL). # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in read Matchbox/MadGraph-OpenLoops.in ################################################## ## Cut selection ## See the documentation for more options ################################################## ## cuts on additional jets cd /Herwig/Cuts # read Matchbox/DefaultEEJets.in # set NJetsCut:NJetsMin 3 ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice SHatScale ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for LEP ## at the Z Mass. (The collection might not be complete.) # read Merging/LEP91-Analysis.in ################################################## ## Do not apply profile scales for LEP as hard ## scale coincides with kinematic limit ################################################## set /Herwig/Shower/ShowerHandler:HardScaleProfile NULL set /Herwig/DipoleShower/DipoleShowerHandler:HardScaleProfile NULL ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode cd /Herwig/Generators saverun LEP-Merging EventGenerator diff --git a/src/Merging/LHC-H-Merging.in b/src/Merging/LHC-H-Merging.in --- a/src/Merging/LHC-H-Merging.in +++ b/src/Merging/LHC-H-Merging.in @@ -1,195 +1,195 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/PPCollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 8000*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 2 set MergingFactory:OrderInAlphaEW 1 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process p p -> h0 [ j j ] set MergingFactory:NLOProcesses 2 # Set the merging scale deviding the parton shower # from the matrix element region in phase space. set Merger:MergingScale 15.*GeV set Merger:MergingScaleSmearing 0.1 # The following line clear a preweighter # that is not working for Higgs as here # pt and HT are 0 before the shower. # See other merging files for more information. clear MergingFactory:Preweighters # The next line can switch of hadronization # and MPI modelling. Use with care!! # read Matchbox/PQCDLevel.in ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production # read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in read Matchbox/OnShellHProduction.in set /Herwig/Particles/h0:Stable Stable # Special settings for the VBF approximation # read Matchbox/VBFDiagramsOnly.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in read Matchbox/MadGraph-OpenLoops.in #read Matchbox/HJets.in # read Matchbox/VBFNLO.in ## Uncomment this to use ggh effective couplings ## currently only supported by MadGraph-GoSam and ## MadGraph-Openloops #read Matchbox/HiggsEffective.in cd /Herwig/MatrixElements/Matchbox/Amplitudes set OpenLoops:HiggsEff Yes set MadGraph:Model heft ################################################## ## Cut selection ## See the documentation for more options ################################################## cd /Herwig/Cuts/ set ChargedLeptonPairMassCut:MinMass 60*GeV set ChargedLeptonPairMassCut:MaxMass 120*GeV cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## cuts on additional jets cd /Herwig/Cuts/ # read Matchbox/DefaultPPJets.in # insert JetCuts:JetRegions 0 FirstJet # insert JetCuts:JetRegions 1 SecondJet # insert JetCuts:JetRegions 2 ThirdJet # insert JetCuts:JetRegions 3 FourthJet ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice FixedScale set FixedScale:FixedScale 125.*GeV ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## PDF choice ################################################## read Matchbox/FiveFlavourNoBMassScheme.in read Matchbox/CT14.in #read Matchbox/MMHT2014.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for Higgs at LHC ## at the 8 TeV. (The collection might not be complete.) # read Merging/LHC8-H-Analysis.in ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode set /Herwig/Generators/EventGenerator:IntermediateOutput Yes cd /Herwig/Generators saverun LHC-H-Merging EventGenerator diff --git a/src/Merging/LHC-J-Merging.in b/src/Merging/LHC-J-Merging.in --- a/src/Merging/LHC-J-Merging.in +++ b/src/Merging/LHC-J-Merging.in @@ -1,192 +1,192 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/PPCollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 7000*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 2 set MergingFactory:OrderInAlphaEW 0 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process p p -> j j [ j ] set MergingFactory:NLOProcesses 1 # Set the merging scale deviding the parton shower # from the matrix element region in phase space. set Merger:MergingScale 20.*GeV set Merger:MergingScaleSmearing 0.1 # The following lines control a preweighter, # that can be used to force more events in higher # HT or pt regions. The unweighted events are accepted # with a enhanced probability W, that is divided fron the # event weight once the event is accepted. # W = ((HT/scale)^HTPower + (pt_max/scale)^MaxPTPower) # with scale = MZ (can be changed) # Note that the weights will therefore differ from "1" # if the powers are not zero. set MPreWeight:HTPower 0 set MPreWeight:MaxPTPower 0 set MPreWeight:OnlyColoured No # The next line can switch of hadronization # and MPI modelling. Use with care!! # read Matchbox/PQCDLevel.in ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production # read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in # read Matchbox/OnShellHProduction.in # Special settings for the VBF approximation # read Matchbox/VBFDiagramsOnly.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in read Matchbox/MadGraph-OpenLoops.in # read Matchbox/HJets.in # read Matchbox/VBFNLO.in ## Uncomment this to use ggh effective couplings ## currently only supported by MadGraph-GoSam ## and MadGraph-OpenLoops # read Matchbox/HiggsEffective.in ################################################## ## Cut selection ## See the documentation for more options ################################################## cd /Herwig/Cuts/ set ChargedLeptonPairMassCut:MinMass 60*GeV set ChargedLeptonPairMassCut:MaxMass 120*GeV cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## cuts on additional jets cd /Herwig/Cuts/ read Matchbox/DefaultPPJets.in insert JetCuts:JetRegions 0 FirstJet # insert JetCuts:JetRegions 1 SecondJet # insert JetCuts:JetRegions 2 ThirdJet # insert JetCuts:JetRegions 3 FourthJet ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice MaxJetPtScale ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## PDF choice ################################################## read Matchbox/FiveFlavourNoBMassScheme.in read Matchbox/CT14.in # read Matchbox/MMHT2014.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for Jets at LHC ## at the 7 TeV. (The collection might not be complete.) # read Merging/LHC7-J-Analysis.in ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode set /Herwig/Generators/EventGenerator:IntermediateOutput Yes cd /Herwig/Generators saverun LHC-J-Merging EventGenerator diff --git a/src/Merging/LHC-T-Merging.in b/src/Merging/LHC-T-Merging.in --- a/src/Merging/LHC-T-Merging.in +++ b/src/Merging/LHC-T-Merging.in @@ -1,194 +1,194 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/PPCollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 7000*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 2 set MergingFactory:OrderInAlphaEW 0 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process p p -> t tbar [ j] set MergingFactory:NLOProcesses 1 # Set the merging scale deviding the parton shower # from the matrix element region in phase space. set Merger:MergingScale 20.*GeV set Merger:MergingScaleSmearing 0.1 # The following lines control a preweighter, # that can be used to force more events in higher # HT or pt regions. The unweighted events are accepted # with a enhanced probability W, that is divided fron the # event weight once the event is accepted. # W = ((HT/scale)^HTPower + (pt_max/scale)^MaxPTPower) # with scale = MZ (can be changed) # Note that the weights will therefore differ from "1" # if the powers are not zero. set MPreWeight:HTPower 0 set MPreWeight:MaxPTPower 0 set MPreWeight:OnlyColoured No # The next line can switch of hadronization # and MPI modelling. Use with care!! # read Matchbox/PQCDLevel.in ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in # read Matchbox/OnShellHProduction.in # Special settings for the VBF approximation # read Matchbox/VBFDiagramsOnly.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in read Matchbox/MadGraph-OpenLoops.in # read Matchbox/HJets.in # read Matchbox/VBFNLO.in ## Uncomment this to use ggh effective couplings ## currently only supported by MadGraph-GoSam and ## MadGraph-Openloops # read Matchbox/HiggsEffective.in ################################################## ## Cut selection ## See the documentation for more options ################################################## cd /Herwig/Cuts/ set ChargedLeptonPairMassCut:MinMass 60*GeV set ChargedLeptonPairMassCut:MaxMass 120*GeV cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## cuts on additional jets cd /Herwig/Cuts/ # read Matchbox/DefaultPPJets.in # insert JetCuts:JetRegions 0 FirstJet # insert JetCuts:JetRegions 1 SecondJet # insert JetCuts:JetRegions 2 ThirdJet # insert JetCuts:JetRegions 3 FourthJet ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice TopPairMassScale ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## PDF choice ################################################## read Matchbox/FiveFlavourNoBMassScheme.in read Matchbox/CT14.in # read Matchbox/MMHT2014.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for Tops at LHC ## at the 7 TeV. (The collection might not be complete.) # read Merging/LHC7-T-Analysis.in ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode set /Herwig/Generators/EventGenerator:IntermediateOutput Yes cd /Herwig/Generators saverun LHC-T-Merging EventGenerator diff --git a/src/Merging/LHC-W-Merging.in b/src/Merging/LHC-W-Merging.in --- a/src/Merging/LHC-W-Merging.in +++ b/src/Merging/LHC-W-Merging.in @@ -1,195 +1,195 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/PPCollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 7000*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 0 set MergingFactory:OrderInAlphaEW 2 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process p p -> l nu [ j j ] set MergingFactory:NLOProcesses 2 # Set the merging scale deviding the parton shower # from the matrix element region in phase space. set Merger:MergingScale 15.*GeV set Merger:MergingScaleSmearing 0.1 # The following lines control a preweighter, # that can be used to force more events in higher # HT or pt regions. The unweighted events are accepted # with a enhanced probability W, that is divided fron the # event weight once the event is accepted. # W = ((HT/scale)^HTPower + (pt_max/scale)^MaxPTPower) # with scale = MZ (can be changed) # Note that the weights will therefore differ from "1" # if the powers are not zero. set MPreWeight:HTPower 0 set MPreWeight:MaxPTPower 0 set MPreWeight:OnlyColoured No # The next line can switch of hadronization # and MPI modelling. Use with care!! # read Matchbox/PQCDLevel.in ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production # read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in # read Matchbox/OnShellHProduction.in # Special settings for the VBF approximation # read Matchbox/VBFDiagramsOnly.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in # read Matchbox/MadGraph-OpenLoops.in # read Matchbox/HJets.in # read Matchbox/VBFNLO.in ## Uncomment this to use ggh effective couplings ## currently only supported by MadGraph-GoSam and ## MadGraph-Openloops # read Matchbox/HiggsEffective.in ################################################## ## Cut selection ## See the documentation for more options ################################################## cd /Herwig/Cuts/ set /Herwig/Cuts/LeptonPairMassCut:MinMass 60*GeV set /Herwig/Cuts/LeptonPairMassCut:MaxMass 120*GeV cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## cuts on additional jets cd /Herwig/Cuts/ # read Matchbox/DefaultPPJets.in # insert JetCuts:JetRegions 0 FirstJet # insert JetCuts:JetRegions 1 SecondJet # insert JetCuts:JetRegions 2 ThirdJet # insert JetCuts:JetRegions 3 FourthJet ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice LeptonPairMassScale ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## PDF choice ################################################## read Matchbox/FiveFlavourNoBMassScheme.in read Matchbox/CT14.in # read Matchbox/MMHT2014.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for Ws at LHC ## at the 7 TeV. (The collection might not be complete.) # read Merging/LHC7-W-Analysis.in ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode set /Herwig/Generators/EventGenerator:IntermediateOutput Yes cd /Herwig/Generators saverun LHC-W-Merging EventGenerator diff --git a/src/Merging/LHC-Z-Merging.in b/src/Merging/LHC-Z-Merging.in --- a/src/Merging/LHC-Z-Merging.in +++ b/src/Merging/LHC-Z-Merging.in @@ -1,195 +1,195 @@ # -*- ThePEG-repository -*- ################################################## ## Herwig/Merging example input file ################################################## ################################################## ## Collider type ################################################## read snippets/DipoleMerging.in read snippets/PPCollider.in read snippets/MonacoSampler.in ################################################## ## Beam energy sqrt(s) ################################################## cd /Herwig/EventHandlers set EventHandler:LuminosityFunction:Energy 7000*GeV ################################################## ## Process selection ################################################## ## Note that event generation may fail if no matching matrix element has ## been found. Coupling orders are with respect to the Born process, ## i.e. NLO QCD does not require an additional power of alphas. ## Model assumptions read Matchbox/StandardModelLike.in read Matchbox/DiagonalCKM.in ## Set the order of the couplings cd /Herwig/Merging set MergingFactory:OrderInAlphaS 0 set MergingFactory:OrderInAlphaEW 2 ## Select the process ## You may use identifiers such as p, pbar, j, l, mu+, h0 etc. do MergingFactory:Process p p -> e+ e- [ j j ] set MergingFactory:NLOProcesses 2 # Set the merging scale deviding the parton shower # from the matrix element region in phase space. set Merger:MergingScale 15.*GeV set Merger:MergingScaleSmearing 0.1 # The following lines control a preweighter, # that can be used to force more events in higher # HT or pt regions. The unweighted events are accepted # with a enhanced probability W, that is divided fron the # event weight once the event is accepted. # W = ((HT/scale)^HTPower + (pt_max/scale)^MaxPTPower) # with scale = MZ (can be changed) # Note that the weights will therefore differ from "1" # if the powers are not zero. set MPreWeight:HTPower 0 set MPreWeight:MaxPTPower 0 set MPreWeight:OnlyColoured No # The next line can switch of hadronization # and MPI modelling. Use with care!! # read Matchbox/PQCDLevel.in ## Special settings required for on-shell production of unstable particles ## enable for on-shell top production # read Matchbox/OnShellTopProduction.in ## enable for on-shell W, Z or h production # read Matchbox/OnShellWProduction.in # read Matchbox/OnShellZProduction.in # read Matchbox/OnShellHProduction.in # Special settings for the VBF approximation # read Matchbox/VBFDiagramsOnly.in ################################################## ## Matrix element library selection ################################################## ## Select a generic tree/loop combination or a ## specialized NLO package # read Matchbox/MadGraph-GoSam.in # read Matchbox/MadGraph-MadGraph.in # read Matchbox/MadGraph-NJet.in # read Matchbox/MadGraph-OpenLoops.in # read Matchbox/HJets.in # read Matchbox/VBFNLO.in ## Uncomment this to use ggh effective couplings ## currently only supported by MadGraph-GoSam and ## MadGraph-Openloops # read Matchbox/HiggsEffective.in ################################################## ## Cut selection ## See the documentation for more options ################################################## cd /Herwig/Cuts/ set ChargedLeptonPairMassCut:MinMass 60*GeV set ChargedLeptonPairMassCut:MaxMass 120*GeV cd /Herwig/MatrixElements/Matchbox/Utility insert DiagramGenerator:ExcludeInternal 0 /Herwig/Particles/gamma ## cuts on additional jets cd /Herwig/Cuts/ # read Matchbox/DefaultPPJets.in # insert JetCuts:JetRegions 0 FirstJet # insert JetCuts:JetRegions 1 SecondJet # insert JetCuts:JetRegions 2 ThirdJet # insert JetCuts:JetRegions 3 FourthJet ################################################## ## Scale choice ## See the documentation for more options ################################################## cd /Herwig/MatrixElements/Matchbox/Scales/ set /Herwig/Merging/MergingFactory:ScaleChoice LeptonPairMassScale ################################################## ## Scale uncertainties ################################################## # read Matchbox/MuDown.in # read Matchbox/MuUp.in ################################################## ## Shower scale uncertainties ################################################## # read Matchbox/MuQDown.in # read Matchbox/MuQUp.in ################################################## ## CMW - Scheme ################################################## -read snippets/Dipole_AutoTune_prel.in +read snippets/Dipole_AutoTunes_gss.in ### Use factor in alpha_s argument: alpha_s(q) -> alpha_s(fac*q) ### with fac=exp(-(67-3pi^2-10/3*Nf)/(33-2Nf)) read Merging/FactorCMWScheme.in ### Linear CMW multiplication: ### alpha_s(q) -> alpha_s(q)(1+K_g*alpha_s(q)/2pi ) # read Merging/LinearCMWScheme.in ################################################## ## PDF choice ################################################## read Matchbox/FiveFlavourNoBMassScheme.in read Matchbox/CT14.in #read Matchbox/MMHT2014.in ################################################## ## Analyses ################################################## cd /Herwig/Analysis ## Write HepMC events. Modify the PrintEvent interface for your needs. # insert /Herwig/Generators/EventGenerator:AnalysisHandlers 0 HepMCFile ## Setup the Rivet analysis: #read snippets/Rivet.in #insert Rivet:Analyses 0 XXX_2017_ABC123 ## Here we collected a various Rivet analysis for Zs at LHC ## at the 8 TeV. (The collection might not be complete.) # read Merging/LHC7-Z-Analysis.in ################################################## ## Save the generator ################################################## do /Herwig/Merging/MergingFactory:ProductionMode set /Herwig/Generators/EventGenerator:IntermediateOutput Yes cd /Herwig/Generators saverun LHC-Z-Merging EventGenerator diff --git a/src/Merging/Makefile.am b/src/Merging/Makefile.am --- a/src/Merging/Makefile.am +++ b/src/Merging/Makefile.am @@ -1,34 +1,33 @@ BUILT_SOURCES = done-all-links Matchboxdir = ${pkgdatadir}/Merging INPUTFILES = \ \ LEP91-Analysis.in \ LHC7-Z-Analysis.in \ LHC7-W-Analysis.in \ LHC7-J-Analysis.in \ LHC7-T-Analysis.in \ LHC8-H-Analysis.in \ \ FactorCMWScheme.in \ LinearCMWScheme.in \ -Merging-Dipole-FactorCMWSchemeTune.in \ \ LEP-Merging.in \ LHC-W-Merging.in \ LHC-Z-Merging.in \ LHC-H-Merging.in \ LHC-T-Merging.in \ LHC-J-Merging.in dist_Matchbox_DATA = $(INPUTFILES) CLEANFILES = done-all-links done-all-links: $(INPUTFILES) @echo "Linking input files" @for i in $(INPUTFILES); do \ if test -f $(srcdir)/$$i -a ! -e $$i; then \ $(LN_S) -f $(srcdir)/$$i; fi; done @touch done-all-links diff --git a/src/Merging/Merging-Dipole-FactorCMWSchemeTune.in b/src/Merging/Merging-Dipole-FactorCMWSchemeTune.in deleted file mode 100644 --- a/src/Merging/Merging-Dipole-FactorCMWSchemeTune.in +++ /dev/null @@ -1,73 +0,0 @@ -# -*- ThePEG-repository -*- - -################################################## -## CMW Scheme with alpha_S=0.118 -################################################## - -read Merging/FactorCMWScheme.in - - -################################################## -## intrinsic pt -################################################## -set /Herwig/DipoleShower/IntrinsicPtGenerator:ValenceIntrinsicPtScale 1.67 -set /Herwig/DipoleShower/IntrinsicPtGenerator:SeaIntrinsicPtScale 1.67 - -################################################## -## (two-loop) alpha_s -################################################## -set /Herwig/DipoleShower/NLOAlphaS:input_alpha_s 0.118 -set /Herwig/DipoleShower/NLOAlphaS:input_scale 91.18*GeV -set /Herwig/Couplings/NLOAlphaS:input_alpha_s 0.118 -set /Herwig/Couplings/NLOAlphaS:input_scale 91.18*GeV - - -################################################## -## cutoffs -################################################## -set /Herwig/DipoleShower/Kinematics/FFMassiveKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/FFLightKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/FIMassiveKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/FILightKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/IFMassiveKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/IFLightKinematics:IRCutoff 0.604 -set /Herwig/DipoleShower/Kinematics/IILightKinematics:IRCutoff 0.604 - -################################################## -## hadronization parameters -################################################## - -set /Herwig/Particles/g:ConstituentMass 0.80 -set /Herwig/Hadronization/ClusterDecayer:ClSmrBottom 0.1 -set /Herwig/Hadronization/ClusterDecayer:ClSmrLight 0.26 -set /Herwig/Hadronization/ClusterDecayer:ClSmrCharm 0.16 - -set /Herwig/Hadronization/ClusterFissioner:ClPowBottom 0.7 -set /Herwig/Hadronization/ClusterFissioner:ClPowLight 0.9 -set /Herwig/Hadronization/ClusterFissioner:ClPowCharm 0.7 - -set /Herwig/Hadronization/ClusterFissioner:ClMaxBottom 5.8 -set /Herwig/Hadronization/ClusterFissioner:ClMaxCharm 2.9 -set /Herwig/Hadronization/ClusterFissioner:ClMaxLight 2.7 - -set /Herwig/Hadronization/ClusterFissioner:PSplitLight 1.0 -set /Herwig/Hadronization/ClusterFissioner:PSplitCharm 1.1 -set /Herwig/Hadronization/ClusterFissioner:PSplitBottom 0.9 - -set /Herwig/Hadronization/HadronSelector:SingleHadronLimitBottom 0.13 -set /Herwig/Hadronization/HadronSelector:SingleHadronLimitCharm 0.0 - -set /Herwig/Hadronization/HadronSelector:PwtUquark 0.65 -set /Herwig/Hadronization/HadronSelector:PwtDquark 0.65 -set /Herwig/Hadronization/HadronSelector:PwtSquark 0.74 -set /Herwig/Hadronization/HadronSelector:PwtCquark 0.30 -set /Herwig/Hadronization/HadronSelector:PwtBquark 0.40 -set /Herwig/Hadronization/HadronSelector:PwtDIquark 0.24 -set /Herwig/Hadronization/HadronSelector:SngWt 0.59 -set /Herwig/Hadronization/HadronSelector:DecWt 0.50 - -set /Herwig/Particles/b:NominalMass 5.3 -set /Herwig/Particles/b:ConstituentMass 4.8 - - -cd / diff --git a/src/snippets/Dipole_AutoTune_prel.in b/src/snippets/Dipole_AutoTune_prel.in deleted file mode 100644 --- a/src/snippets/Dipole_AutoTune_prel.in +++ /dev/null @@ -1,139 +0,0 @@ -# This is a preliminary Tune. -# With this version (7.1.3) various changes to the dipole shower are done. -# The dipole shower will be retuned for version 7.1.4. - -# Note: The gluon ConstituentMass is lower than the strange quark -# threshold assumed in the BaryonicReconnection. -# Using both snippets will produce no resonable results. - -set /Herwig/DipoleShower/NLOAlphaS:input_alpha_s 0.118 -set /Herwig/DipoleShower/NLOAlphaS:input_scale 91.18*GeV -set /Herwig/Couplings/NLOAlphaS:input_alpha_s 0.118 -set /Herwig/Couplings/NLOAlphaS:input_scale 91.18*GeV - -set /Herwig/DipoleShower/Kinematics/FFMassiveKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/FFLightKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/FIMassiveKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/FILightKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/IFMassiveKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/IFLightKinematics:IRCutoff 0.80399886545 -set /Herwig/DipoleShower/Kinematics/IILightKinematics:IRCutoff 0.80399886545 - -set /Herwig/Particles/g:ConstituentMass 0.762996296601 - -set /Herwig/Hadronization/ClusterDecayer:ClSmrLight 0.342544255421 -set /Herwig/Hadronization/ClusterDecayer:ClSmrCharm 0.237545842862 -set /Herwig/Hadronization/ClusterDecayer:ClSmrBottom 0.0620545338048 - -set /Herwig/Hadronization/ClusterFissioner:ClPowLight 0.950949552153 -set /Herwig/Hadronization/ClusterFissioner:ClPowCharm 0.725163909065 -set /Herwig/Hadronization/ClusterFissioner:ClPowBottom 0.809774254452 - -set /Herwig/Hadronization/ClusterFissioner:ClMaxLight 2.80478186219 -set /Herwig/Hadronization/ClusterFissioner:ClMaxCharm 3.13483791689 -set /Herwig/Hadronization/ClusterFissioner:ClMaxBottom 5.79181302131 - -set /Herwig/Hadronization/ClusterFissioner:PSplitLight 0.995391682928 -set /Herwig/Hadronization/ClusterFissioner:PSplitCharm 0.769749617516 -set /Herwig/Hadronization/ClusterFissioner:PSplitBottom 0.479199412481 - -set /Herwig/Hadronization/HadronSelector:SingleHadronLimitCharm 0.0485533655422 -set /Herwig/Hadronization/HadronSelector:SingleHadronLimitBottom 0.0968582906399 - -set /Herwig/Hadronization/HadronSelector:PwtUquark 0.718641969264 -set /Herwig/Hadronization/HadronSelector:PwtDquark 0.674967529289 -set /Herwig/Hadronization/HadronSelector:PwtSquark 0.542042947191 -set /Herwig/Hadronization/HadronSelector:PwtCquark 0.251544496683 -set /Herwig/Hadronization/HadronSelector:PwtBquark 0.426714635544 -set /Herwig/Hadronization/HadronSelector:PwtDIquark 0.296781084091 -set /Herwig/Hadronization/HadronSelector:SngWt 0.541287379336 -set /Herwig/Hadronization/HadronSelector:DecWt 0.569933528697 - -set /Herwig/Particles/b:NominalMass 5.38329996032 -set /Herwig/Particles/b:ConstituentMass 4.84231667563 - -####### -# CMW scheme: -# Tune was done with CMW for all kernels, as merging currently only is consistent if applied to all. -####### - - -cd /Herwig/DipoleShower/Kernels -set FFgx2ggxDipoleKernel:CMWScheme Factor -set FFqx2qgxDipoleKernel:CMWScheme Factor -set FFgx2ddxDipoleKernel:CMWScheme Factor -set FFgx2uuxDipoleKernel:CMWScheme Factor -set FFgx2ccxDipoleKernel:CMWScheme Factor -set FFgx2ssxDipoleKernel:CMWScheme Factor -set FFgx2bbxDipoleKernel:CMWScheme Factor -set FFMgx2ggxDipoleKernel:CMWScheme Factor -set FFMdx2dgxDipoleKernel:CMWScheme Factor -set FFMux2ugxDipoleKernel:CMWScheme Factor -set FFMcx2cgxDipoleKernel:CMWScheme Factor -set FFMsx2sgxDipoleKernel:CMWScheme Factor -set FFMbx2bgxDipoleKernel:CMWScheme Factor -set FFMtx2tgxDipoleKernel:CMWScheme Factor -set FFMgx2ddxDipoleKernel:CMWScheme Factor -set FFMgx2uuxDipoleKernel:CMWScheme Factor -set FFMgx2ccxDipoleKernel:CMWScheme Factor -set FFMgx2ssxDipoleKernel:CMWScheme Factor -set FFMgx2bbxDipoleKernel:CMWScheme Factor -set FIgx2ggxDipoleKernel:CMWScheme Factor -set FIqx2qgxDipoleKernel:CMWScheme Factor -set FIgx2ddxDipoleKernel:CMWScheme Factor -set FIgx2uuxDipoleKernel:CMWScheme Factor -set FIgx2ccxDipoleKernel:CMWScheme Factor -set FIgx2ssxDipoleKernel:CMWScheme Factor -set FIgx2bbxDipoleKernel:CMWScheme Factor -set FIMdx2dgxDipoleKernel:CMWScheme Factor -set FIMux2ugxDipoleKernel:CMWScheme Factor -set FIMcx2cgxDipoleKernel:CMWScheme Factor -set FIMsx2sgxDipoleKernel:CMWScheme Factor -set FIMbx2bgxDipoleKernel:CMWScheme Factor -set FIMtx2tgxDipoleKernel:CMWScheme Factor -set FIMgx2ddxDipoleKernel:CMWScheme Factor -set FIMgx2uuxDipoleKernel:CMWScheme Factor -set FIMgx2ccxDipoleKernel:CMWScheme Factor -set FIMgx2ssxDipoleKernel:CMWScheme Factor -set FIMgx2bbxDipoleKernel:CMWScheme Factor -#set FIMgx2ttxDipoleKernel:CMWScheme Factor -set IFgx2ggxDipoleKernel:CMWScheme Factor -set IFqx2qgxDipoleKernel:CMWScheme Factor -set IFqx2gqxDipoleKernel:CMWScheme Factor -set IFgx2ddbarxDipoleKernel:CMWScheme Factor -set IFgx2dbardxDipoleKernel:CMWScheme Factor -set IFgx2uubarxDipoleKernel:CMWScheme Factor -set IFgx2ubaruxDipoleKernel:CMWScheme Factor -set IFgx2ccbarxDipoleKernel:CMWScheme Factor -set IFgx2cbarcxDipoleKernel:CMWScheme Factor -set IFgx2ssbarxDipoleKernel:CMWScheme Factor -set IFgx2sbarsxDipoleKernel:CMWScheme Factor -set IFMgx2ggxDipoleKernel:CMWScheme Factor -set IFMqx2qgxDipoleKernel:CMWScheme Factor -set IFMqx2gqxDipoleKernel:CMWScheme Factor -set IFMgx2ddbarxDipoleKernel:CMWScheme Factor -set IFMgx2dbardxDipoleKernel:CMWScheme Factor -set IFMgx2uubarxDipoleKernel:CMWScheme Factor -set IFMgx2ubaruxDipoleKernel:CMWScheme Factor -set IFMgx2ccbarxDipoleKernel:CMWScheme Factor -set IFMgx2cbarcxDipoleKernel:CMWScheme Factor -set IFMgx2ssbarxDipoleKernel:CMWScheme Factor -set IFMgx2sbarsxDipoleKernel:CMWScheme Factor -set IIgx2ggxDipoleKernel:CMWScheme Factor -set IIqx2qgxDipoleKernel:CMWScheme Factor -set IIqx2gqxDipoleKernel:CMWScheme Factor -set IIgx2ddbarxDipoleKernel:CMWScheme Factor -set IIgx2dbardxDipoleKernel:CMWScheme Factor -set IIgx2uubarxDipoleKernel:CMWScheme Factor -set IIgx2ubaruxDipoleKernel:CMWScheme Factor -set IIgx2ccbarxDipoleKernel:CMWScheme Factor -set IIgx2cbarcxDipoleKernel:CMWScheme Factor -set IIgx2ssbarxDipoleKernel:CMWScheme Factor -set IIgx2sbarsxDipoleKernel:CMWScheme Factor -set IFgx2bbbarxDipoleKernel:CMWScheme Factor -set IFgx2bbarbxDipoleKernel:CMWScheme Factor -set IFMgx2bbbarxDipoleKernel:CMWScheme Factor -set IFMgx2bbarbxDipoleKernel:CMWScheme Factor -set IIgx2bbbarxDipoleKernel:CMWScheme Factor -set IIgx2bbarbxDipoleKernel:CMWScheme Factor - diff --git a/src/snippets/Dipole_AutoTunes_gss.in b/src/snippets/Dipole_AutoTunes_gss.in new file mode 100644 --- /dev/null +++ b/src/snippets/Dipole_AutoTunes_gss.in @@ -0,0 +1,119 @@ +set /Herwig/Shower/AlphaQCD:AlphaIn 0.118 +set /Herwig/Couplings/NLOAlphaS:input_alpha_s 0.118 +set /Herwig/DipoleShower/NLOAlphaS:input_alpha_s 0.118 +cd /Herwig/Particles/ +set g:ConstituentMass 0.95 +set b:NominalMass 4.700501 +set b:ConstituentMass 4.084889 +cd /Herwig/Hadronization/ +set ClusterDecayer:ClSmrBottom 0.085964 +set ClusterDecayer:ClSmrLight 0.698877 +set ClusterDecayer:ClSmrCharm 0.246296 +set ClusterFissioner:ClPowBottom 0.591646 +set ClusterFissioner:ClPowLight 0.99945 +set ClusterFissioner:ClPowCharm 3.386187 +set ClusterFissioner:ClMaxBottom 3.771649 +set ClusterFissioner:ClMaxCharm 4.780456 +set ClusterFissioner:ClMaxLight 3.055256 +set ClusterFissioner:PSplitLight 0.7779 +set ClusterFissioner:PSplitCharm 0.628766 +set ClusterFissioner:PSplitBottom 0.662911 +set HadronSelector:SingleHadronLimitBottom 0.000446 +set HadronSelector:SingleHadronLimitCharm 0.000508 +set HadronSelector:SngWt 0.927141 +set HadronSelector:DecWt 0.630787 + +set /Herwig/DipoleShower/Kinematics/FFMassiveKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/FFLightKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/FIMassiveKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/FILightKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/IFMassiveKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/IFLightKinematics:IRCutoff 1.000896 +set /Herwig/DipoleShower/Kinematics/IILightKinematics:IRCutoff 1.000896 + + +####### +# CMW scheme: +# Tune was done with CMW for all kernels, as merging currently only is consistent if applied to all. +####### + +cd /Herwig/DipoleShower/Kernels +set FFgx2ggxDipoleKernel:CMWScheme Factor +set FFqx2qgxDipoleKernel:CMWScheme Factor +set FFgx2ddxDipoleKernel:CMWScheme Factor +set FFgx2uuxDipoleKernel:CMWScheme Factor +set FFgx2ccxDipoleKernel:CMWScheme Factor +set FFgx2ssxDipoleKernel:CMWScheme Factor +set FFgx2bbxDipoleKernel:CMWScheme Factor +set FFMgx2ggxDipoleKernel:CMWScheme Factor +set FFMdx2dgxDipoleKernel:CMWScheme Factor +set FFMux2ugxDipoleKernel:CMWScheme Factor +set FFMcx2cgxDipoleKernel:CMWScheme Factor +set FFMsx2sgxDipoleKernel:CMWScheme Factor +set FFMbx2bgxDipoleKernel:CMWScheme Factor +set FFMtx2tgxDipoleKernel:CMWScheme Factor +set FFMgx2ddxDipoleKernel:CMWScheme Factor +set FFMgx2uuxDipoleKernel:CMWScheme Factor +set FFMgx2ccxDipoleKernel:CMWScheme Factor +set FFMgx2ssxDipoleKernel:CMWScheme Factor +set FFMgx2bbxDipoleKernel:CMWScheme Factor +set FIgx2ggxDipoleKernel:CMWScheme Factor +set FIqx2qgxDipoleKernel:CMWScheme Factor +set FIgx2ddxDipoleKernel:CMWScheme Factor +set FIgx2uuxDipoleKernel:CMWScheme Factor +set FIgx2ccxDipoleKernel:CMWScheme Factor +set FIgx2ssxDipoleKernel:CMWScheme Factor +set FIgx2bbxDipoleKernel:CMWScheme Factor +set FIMdx2dgxDipoleKernel:CMWScheme Factor +set FIMux2ugxDipoleKernel:CMWScheme Factor +set FIMcx2cgxDipoleKernel:CMWScheme Factor +set FIMsx2sgxDipoleKernel:CMWScheme Factor +set FIMbx2bgxDipoleKernel:CMWScheme Factor +set FIMtx2tgxDipoleKernel:CMWScheme Factor +set FIMgx2ddxDipoleKernel:CMWScheme Factor +set FIMgx2uuxDipoleKernel:CMWScheme Factor +set FIMgx2ccxDipoleKernel:CMWScheme Factor +set FIMgx2ssxDipoleKernel:CMWScheme Factor +set FIMgx2bbxDipoleKernel:CMWScheme Factor +#set FIMgx2ttxDipoleKernel:CMWScheme Factor +set IFgx2ggxDipoleKernel:CMWScheme Factor +set IFqx2qgxDipoleKernel:CMWScheme Factor +set IFqx2gqxDipoleKernel:CMWScheme Factor +set IFgx2ddbarxDipoleKernel:CMWScheme Factor +set IFgx2dbardxDipoleKernel:CMWScheme Factor +set IFgx2uubarxDipoleKernel:CMWScheme Factor +set IFgx2ubaruxDipoleKernel:CMWScheme Factor +set IFgx2ccbarxDipoleKernel:CMWScheme Factor +set IFgx2cbarcxDipoleKernel:CMWScheme Factor +set IFgx2ssbarxDipoleKernel:CMWScheme Factor +set IFgx2sbarsxDipoleKernel:CMWScheme Factor +set IFMgx2ggxDipoleKernel:CMWScheme Factor +set IFMqx2qgxDipoleKernel:CMWScheme Factor +set IFMqx2gqxDipoleKernel:CMWScheme Factor +set IFMgx2ddbarxDipoleKernel:CMWScheme Factor +set IFMgx2dbardxDipoleKernel:CMWScheme Factor +set IFMgx2uubarxDipoleKernel:CMWScheme Factor +set IFMgx2ubaruxDipoleKernel:CMWScheme Factor +set IFMgx2ccbarxDipoleKernel:CMWScheme Factor +set IFMgx2cbarcxDipoleKernel:CMWScheme Factor +set IFMgx2ssbarxDipoleKernel:CMWScheme Factor +set IFMgx2sbarsxDipoleKernel:CMWScheme Factor +set IIgx2ggxDipoleKernel:CMWScheme Factor +set IIqx2qgxDipoleKernel:CMWScheme Factor +set IIqx2gqxDipoleKernel:CMWScheme Factor +set IIgx2ddbarxDipoleKernel:CMWScheme Factor +set IIgx2dbardxDipoleKernel:CMWScheme Factor +set IIgx2uubarxDipoleKernel:CMWScheme Factor +set IIgx2ubaruxDipoleKernel:CMWScheme Factor +set IIgx2ccbarxDipoleKernel:CMWScheme Factor +set IIgx2cbarcxDipoleKernel:CMWScheme Factor +set IIgx2ssbarxDipoleKernel:CMWScheme Factor +set IIgx2sbarsxDipoleKernel:CMWScheme Factor +set IFgx2bbbarxDipoleKernel:CMWScheme Factor +set IFgx2bbarbxDipoleKernel:CMWScheme Factor +set IFMgx2bbbarxDipoleKernel:CMWScheme Factor +set IFMgx2bbarbxDipoleKernel:CMWScheme Factor +set IIgx2bbbarxDipoleKernel:CMWScheme Factor +set IIgx2bbarbxDipoleKernel:CMWScheme Factor + + diff --git a/src/snippets/Makefile.am b/src/snippets/Makefile.am --- a/src/snippets/Makefile.am +++ b/src/snippets/Makefile.am @@ -1,48 +1,48 @@ BUILT_SOURCES = done-all-links snippetsdir = ${pkgdatadir}/snippets INPUTFILES = \ CellGridSampler.in \ Diffraction.in \ DipoleMerging.in \ DipoleShowerFiveFlavours.in \ DipoleShowerFourFlavours.in \ EECollider.in \ EPCollider.in \ HepMCFixedOrder.in \ HepMC.in \ Matchbox.in \ MB-DipoleShower.in \ MB.in \ MonacoSampler.in \ Particles-SetLonglivedParticlesStable.in \ PDF-CT10.in \ PDF-NNPDF30NLO.in \ PPCollider.in \ RivetFixedOrder.in \ Rivet.in \ YFS.in \ CMWinQtiledShower.in \ -Dipole_AutoTune_prel.in \ +Dipole_AutoTunes_gss.in \ Tune-DotProduct.in \ Tune-DotProduct-Veto.in \ Tune-pT.in \ Tune-Q2.in \ EvolutionScheme-DotProduct.in \ EvolutionScheme-DotProduct-Veto.in \ EvolutionScheme-pT.in \ EvolutionScheme-Q2.in \ FixedTarget.in \ FixedTarget-PP.in dist_snippets_DATA = $(INPUTFILES) CLEANFILES = done-all-links done-all-links: $(INPUTFILES) @echo "Linking input files" @for i in $(INPUTFILES); do \ if test -f $(srcdir)/$$i -a ! -e $$i; then \ $(LN_S) -f $(srcdir)/$$i; fi; done @touch done-all-links