diff --git a/Gen/LbPythia8/options/LHCbDefault.cmd b/Gen/LbPythia8/options/LHCbDefault.cmd index 17a3a9d9fb3c44388863865bdae3886e317514b3..723bcf713ac9a748f45125d01a514121547dcd21 100644 --- a/Gen/LbPythia8/options/LHCbDefault.cmd +++ b/Gen/LbPythia8/options/LHCbDefault.cmd @@ -49,27 +49,34 @@ SigmaProcess:alphaSorder = 2 # ------------------------------- PDF:pSet = LHAPDF6:CT09MCS/0 # ------------------------------- -# First tuning for pythia8 and ct09mcs: Tune 1 +# Tuning for pythia8 and ct09mcs: Sim10 Lambda0 tune # ------------------------------- SpaceShower:rapidityOrder = off # General MultipartonInteractions:bProfile = 1 # reset the default bProfile SpaceShower:alphaSvalue = 0.130 MultipartonInteractions:alphaSvalue = 0.130 MultipartonInteractions:ecmRef = 7000 -MultipartonInteractions:pT0Ref = 2.742289e+00' # Multiparton interactions +MultipartonInteractions:pT0Ref = 2.4188 # Multiparton interactions MultipartonInteractions:ecmPow = 0.238 # ------------------------------- -# flavour selection and excited states -# with seeting for light flavour as determined in LHCb Tune1: +# flavour selection and excited states with color reconnection +# with setting for light flavour as determined in LHCb Sim10 Lambda0 Tune # mesonSvector, probQQtoQ, probStoUD # ------------------------------- +BeamRemnants:remnantMode = 1 +ColourReconnection:mode = 1 +ColourReconnection:lambdaForm = 0 +ColourReconnection:allowDoubleJunRem = off +ColourReconnection:timeDilationPar = 1.9667 +ColourReconnection:m0 = 0.9021 +ColourReconnection:junctionCorrection = 0.8399 StringFlav:mesonUDvector = 0.6 -StringFlav:mesonSvector = 7.474387e-01 +StringFlav:mesonSvector = 0.3819 StringFlav:mesonCvector = 3.0 StringFlav:mesonBvector = 3.0 -StringFlav:probQQtoQ = 1.615701e-01 -StringFlav:probStoUD = 3.501613e-01 -StringFlav:probSQtoQQ = 0.4 +StringFlav:probQQtoQ = 0.1086 +StringFlav:probStoUD = 0.3590 +StringFlav:probSQtoQQ = 0.6451 StringFlav:probQQ1toQQ0 = 0.05 StringFlav:mesonUDL1S0J1 = 0.0989 StringFlav:mesonUDL1S1J0 = 0.0132 diff --git a/Gen/LbPythia8/options/Pythia8_Lambda2.py b/Gen/LbPythia8/options/Pythia8_Lambda2.py new file mode 100644 index 0000000000000000000000000000000000000000..083b8a6dcb444bddb144aa4c45957b49d2e57815 --- /dev/null +++ b/Gen/LbPythia8/options/Pythia8_Lambda2.py @@ -0,0 +1,72 @@ +from Configurables import Generation +from Configurables import MinimumBias +from Configurables import Inclusive +from Configurables import ( SignalPlain, SignalRepeatedHadronization, Special ) +from Configurables import Pythia8Production + +commandsTuning_Lambda2 = [ + 'MultipartonInteractions:pT0Ref = 2.3003', + 'ColourReconnection:lambdaForm = 2', + 'ColourReconnection:timeDilationPar = 2.3123', + 'ColourReconnection:m0 = 1.4563', + 'ColourReconnection:junctionCorrection = 1.0763' + 'StringFlav:probSQtoQQ = 0.5117', + 'StringFlav:mesonSvector = 0.3860', + 'StringFlav:probStoUD = 0.3162', + 'StringFlav:probQQtoQ = 0.0380', + ] + +Pythia8TurnOffFragmentation = [ "HadronLevel:all = off" ] + +gen = Generation("Generation") + +gen.addTool( MinimumBias , name = "MinimumBias" ) +gen.MinimumBias.ProductionTool = "Pythia8Production" +gen.MinimumBias.addTool( Pythia8Production , name = "Pythia8Production" ) +gen.MinimumBias.Pythia8Production.Tuning = "LHCbDefault.cmd" +gen.MinimumBias.Pythia8Production.Commands += commandsTuning_Lambda2 + +gen.addTool( Inclusive , name = "Inclusive" ) +gen.Inclusive.ProductionTool = "Pythia8Production" +gen.Inclusive.addTool( Pythia8Production , name = "Pythia8Production" ) +gen.Inclusive.Pythia8Production.Tuning = "LHCbDefault.cmd" +gen.MinimumBias.Pythia8Production.Commands += commandsTuning_Lambda2 + +gen.addTool( SignalPlain , name = "SignalPlain" ) +gen.SignalPlain.ProductionTool = "Pythia8Production" +gen.SignalPlain.addTool( Pythia8Production , name = "Pythia8Production" ) +gen.SignalPlain.Pythia8Production.Tuning = "LHCbDefault.cmd" +gen.SignalPlain.Pythia8Production.Commands += commandsTuning_Lambda2 + +gen.addTool( SignalRepeatedHadronization , name = "SignalRepeatedHadronization" ) +gen.SignalRepeatedHadronization.ProductionTool = "Pythia8Production" +gen.SignalRepeatedHadronization.addTool( Pythia8Production , name = "Pythia8Production" ) +gen.SignalRepeatedHadronization.Pythia8Production.Tuning = "LHCbDefault.cmd" +gen.SignalRepeatedHadronization.Pythia8Production.Commands += Pythia8TurnOffFragmentation +gen.SignalRepeatedHadronization.Pythia8Production.Commands += commandsTuning_Lambda2 + +gen.addTool( Special , name = "Special" ) +gen.Special.ProductionTool = "Pythia8Production" +gen.Special.addTool( Pythia8Production , name = "Pythia8Production" ) +gen.Special.Pythia8Production.Tuning = "LHCbDefault.cmd" +gen.Special.Pythia8Production.Commands += commandsTuning_Lambda2 +gen.Special.PileUpProductionTool = "Pythia8Production/Pythia8PileUp" +gen.Special.addTool(Pythia8Production, name = "Pythia8PileUp") +gen.Special.Pythia8PileUp.Tuning = "LHCbDefault.cmd" +gen.Special.Pythia8PileUp.Commands += commandsTuning_Lambda2 +gen.Special.ReinitializePileUpGenerator = False + +# Use same generator and configuration for spillover +from Configurables import Gauss +spillOverList = Gauss().getProp("SpilloverPaths") +for slot in spillOverList: + genSlot = Generation("Generation"+slot) + genSlot.addTool(MinimumBias, name = "MinimumBias") + genSlot.MinimumBias.ProductionTool = "Pythia8Production" + genSlot.MinimumBias.addTool(Pythia8Production, name = "Pythia8Production") + genSlot.MinimumBias.Pythia8Production.Commands += commandsTuning_Lambda2 + genSlot.MinimumBias.Pythia8Production.Tuning = "LHCbDefault.cmd" + + + +