def CkvPVTimeSmearActivate(): from Configurables import Generation from Configurables import BeamSpotMarkovChainSampleVertex Generation("Generation").addTool(BeamSpotMarkovChainSampleVertex, name="BeamSpotMarkovChainSampleVertex") Generation( "Generation").VertexSmearingTool = "BeamSpotMarkovChainSampleVertex"
# Import the necessary modules. from Configurables import Generation, MinimumBias, Special, Pythia8Production from Configurables import BcVegPyProduction # Add Pythia as minimum bias production tool. Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "Pythia8Production" Generation().MinimumBias.addTool(Pythia8Production) # Add BcVegPy as special production tool. Generation().addTool(Special) Generation().Special.ProductionTool = "BcVegPyProduction" Generation().Special.addTool(BcVegPyProduction) Generation().Special.BcVegPyProduction.ShowerToolName = "Pythia8Production" Generation().Special.PileUpProductionTool = "Pythia8Production/Pythia8PileUp" Generation().Special.ReinitializePileUpGenerator = False
# Import the necessary modules. from Configurables import Generation, HerwigppProduction from Configurables import MinimumBias, Inclusive, Special, SignalPlain from Configurables import SignalRepeatedHadronization # Add Herwig++ as minimum bias production tool. Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "HerwigppProduction" Generation().MinimumBias.addTool(HerwigppProduction) # Add Herwig++ as inclusive production tool. Generation().addTool(Inclusive) Generation().Inclusive.ProductionTool = "HerwigppProduction" Generation().Inclusive.addTool(HerwigppProduction) # Add Herwig++ as special production tool. Generation().addTool(Special) Generation().Special.ProductionTool = "HerwigppProduction" Generation().Special.addTool(HerwigppProduction) Generation().Special.HerwigppProduction.Nlo = True # Add Herwig++ as plain signal production tool. Generation().addTool(SignalPlain) Generation().SignalPlain.ProductionTool = "HerwigppProduction" Generation().SignalPlain.addTool(HerwigppProduction) # Add Herwig++ as repeated signal hadronization production tool. Generation().addTool(SignalRepeatedHadronization) Generation().SignalRepeatedHadronization.ProductionTool = "HerwigppProduction" Generation().SignalRepeatedHadronization.addTool(HerwigppProduction)
# Import the necessary modules. from Configurables import Generation, MinimumBias, Special, PythiaProduction from Configurables import BcVegPyProduction # Add Pythia as minimum bias production tool. Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "PythiaProduction" Generation().MinimumBias.addTool(PythiaProduction) # Add BcVegPy as special production tool. Generation().addTool(Special) Generation().Special.ProductionTool = "BcVegPyProduction" Generation().Special.addTool(BcVegPyProduction) Generation().Special.BcVegPyProduction.ShowerToolName = "PythiaProduction"
#----------------------------------------------------------------------------- # Generator phase #----------------------------------------------------------------------------- GaussGen = GenInit("GaussGen") # # Set the random numbers - these fix the random seed. # GaussGen.FirstEventNumber = 1 GaussGen.RunNumber = 1082 # The output is managed below, so we disable the standard Gauss output. Gauss().OutputType = 'NONE' Gauss().Histograms = 'NONE' # Switch off pileup (recommended!). Generation().addTool(FixedNInteractions, name="FixedNInteractions") Generation().FixedNInteractions.NInteractions = 0 Generation().PileUpTool = "FixedNInteractions" # Define "special" production. Generation().addTool(Special, "Special") Generation().SampleGenerationTool = "Special" # No cuts. Generation().Special.CutTool = "" # Define the production tool Generation().Special.addTool(ReadLHEfileProduction, name="ReadLHEfileProduction") Generation().Special.ProductionTool = "ReadLHEfileProduction"
"pypars parp 91 1.0" , "pypars parp 149 0.02" , "pypars parp 150 0.085" , "pydat1 parj 11 0.4" , "pydat1 parj 12 0.4" , "pydat1 parj 13 0.769" , "pydat1 parj 14 0.09" , "pydat1 parj 15 0.018" , "pydat1 parj 16 0.0815" , "pydat1 parj 17 0.0815" , "pydat1 mstj 26 0" , "pydat1 parj 33 0.4" ] from Configurables import Generation, Special, PythiaProduction, MinimumBias, Inclusive, SignalPlain, SignalRepeatedHadronization gen = Generation() gen.addTool( Special ) gen.addTool( MinimumBias ) gen.addTool( Inclusive ) gen.addTool( SignalPlain ) gen.addTool( SignalRepeatedHadronization ) gen.Special.ProductionTool = "PythiaProduction" gen.MinimumBias.ProductionTool = "PythiaProduction" gen.Inclusive.ProductionTool = "PythiaProduction" gen.SignalPlain.ProductionTool = "PythiaProduction" gen.SignalRepeatedHadronization.ProductionTool = "PythiaProduction" # Use same generator and configuration for spillover from Configurables import Gauss spillOverList = Gauss().getProp("SpilloverPaths")
# Generator Algos generator = GaudiSequencer("Generator") generator.OutputLevel = 2 genSlotMainSeq = GaudiSequencer("GeneratorSlotMainSeq") genSlotMainSeq.OutputLevel = 2 from Configurables import GenInit gaussGen = GenInit("GaussGen") gaussGen.OutputLevel = 2 gaussGen.FirstEventNumber = 1 gaussGen.RunNumber = 1082 from Configurables import Generation generation = Generation("Generation") generation.OutputLevel = 2 genMonitor = GaudiSequencer("GenMonitor") genMonitor.OutputLevel = 2 from Configurables import GenMonitorAlg genMonitorAlg = GenMonitorAlg("GenMonitorAlg") genMonitorAlg.OutputLevel = 2 # Simulation Algos simulation = GaudiSequencer("Simulation") simulation.OutputLevel = 2 from Configurables import GiGaDataStoreAlgorithm gigaStore = GiGaDataStoreAlgorithm("GiGaStore")
from Gaudi.Configuration import * from Configurables import Generation, MinimumBias, HijingProduction from Configurables import ToolSvc, EvtGenDecay Generation().addTool( MinimumBias ) Generation().MinimumBias.addTool( HijingProduction ) # Configuration for pA collisions Generation().MinimumBias.HijingProduction.Commands += [ "hijinginit frame CMS", "hijinginit proj A", "hijinginit targ P", "hijinginit beam1", "hijinginit iap 208", "hijinginit izp 82", "hijinginit iat 1", "hijinginit izt 1" ]
# file /build/jenkins-tests/workspace/nightly-builds/checkout/tmp/checkout/DBASE/Gen/DecFiles/v30r41/options/11102432.py generated: Mon, 16 Dec 2019 12:50:06 # # Event Type: 11102432 # # ASCII decay Descriptor: {[[B0]nos -> K+ pi- (pi0 -> gamma gamma)]cc, [[B0]os -> K- pi+ (pi0 -> gamma gamma)]cc} # from Configurables import Generation Generation().EventType = 11102432 Generation().SampleGenerationTool = "SignalRepeatedHadronization" from Configurables import SignalRepeatedHadronization Generation().addTool(SignalRepeatedHadronization) Generation().SignalRepeatedHadronization.ProductionTool = "PythiaProduction" from Configurables import ToolSvc from Configurables import EvtGenDecay ToolSvc().addTool(EvtGenDecay) ToolSvc( ).EvtGenDecay.UserDecayFile = "$DECFILESROOT/dkfiles/Bd_K+pi-pi0=DecProdCut,sqDalitz.dec" Generation().SignalRepeatedHadronization.CutTool = "DaughtersInLHCb" Generation().SignalRepeatedHadronization.SignalPIDList = [511, -511] Generation().PileUpTool = "FixedNInteractions"
GaussGen = GenInit("GaussGen") GaussGen.RunNumber = 1082 LHCbApp().EvtMax = 10000 #LHCbApp().setProp("OutputLevel", 3) #LHCbApp().setProp("OutputLevel", 4) LHCbApp().OutputLevel = WARNING Gauss().Histograms = "NONE" Gauss().OutputType = "NONE" Gauss().DatasetName = "GaussDevWHepMC" #Special settings for beam conditions: #** ignore multiple pp collisions, force nPV===1 (algorithm will complain if not) from Configurables import Generation gaussGen = Generation("Generation") gaussGen.PileUpTool = "FixedNInteractions" #** force head-on pp collisions: # - theoreticians always assume pp collision is in CM frame # - you may experience some problems in generating event with settings below (issue under investigation) # - see CorrectCrossingAngles option below to control CM boost at runtime #Gauss().BeamCrossingAngle = 0.0 #Gauss().BeamBetaStar = 0.0 GenMonitor = GaudiSequencer("GenMonitor") GenMonitor.Members += ["RivetAnalysisHandler"] rivet = RivetAnalysisHandler() rivet.BaseFileName = "testBird" rivet.RunName = "LHCbMBPy8" rivet.Analyses = ["MC_GENERIC", "MC_IDENTIFIED"]
from Configurables import Generation, Gauss from GaudiKernel import SystemOfUnits Generation().Special.OniaPairsProduction.Ecm = 2 * Gauss( ).BeamMomentum / SystemOfUnits.GeV Generation().Special.OniaPairsProduction.Psi1S1S = 1 Generation().Special.OniaPairsProduction.Psi1S2S = 1 Generation().Special.OniaPairsProduction.Psi2S2S = 1
# file /build/jenkins-tests/workspace/nightly-builds/checkout/tmp/checkout/DBASE/Gen/DecFiles/v30r43/options/15102430.py generated: Wed, 26 Feb 2020 18:42:35 # # Event Type: 15102430 # # ASCII decay Descriptor: [ Lambda_b0 -> p+ K- pi0 ]cc # from Configurables import Generation Generation().EventType = 15102430 Generation().SampleGenerationTool = "SignalPlain" from Configurables import SignalPlain Generation().addTool(SignalPlain) Generation().SignalPlain.ProductionTool = "PythiaProduction" from Configurables import ToolSvc from Configurables import EvtGenDecay ToolSvc().addTool(EvtGenDecay) ToolSvc( ).EvtGenDecay.UserDecayFile = "$DECFILESROOT/dkfiles/Lb_pKpi0=DecProdCut.dec" Generation().SignalPlain.CutTool = "DaughtersInLHCb" Generation().SignalPlain.SignalPIDList = [5122, -5122] Generation().PileUpTool = "FixedNInteractions"
# Import the necessary modules. from Configurables import Generation, Special, SuperChicProduction # Add SuperChic as special production tool. Generation().addTool(Special) Generation().Special.ProductionTool = "SuperChicProduction" Generation().Special.addTool(SuperChicProduction)
# Generate Pile Up from Gaudi.Configuration import * from Configurables import Generation, Pythia8Production, Gauss gen = Generation() gen.Special.PileUpProductionTool = "Pythia8Production/Pythia8PileUp" gen.PileUpTool = "FixedLuminosityForRareProcess" gen.Special.addTool(Pythia8Production, name="Pythia8PileUp") gen.Special.Pythia8PileUp.Tuning = "LHCbDefault.cmd" gen.Special.ReinitializePileUpGenerator = False # AlpGen specific options from Configurables import AlpGenProductionwqq, AlpGenProductionzbb, AlpGenProductiontop, AlpGenProduction4q # AlpGen conf dict try: AlpGenDict = Generation.AlpGenDict except: raise NameError("AlpGenDict not defined") # File type try: alpgen_FileLabel = AlpGenDict["alpgen_FileLabel"] except: raise NameError('value for alpgen_FileLabel property not set') else: if alpgen_FileLabel == "wbb" or alpgen_FileLabel == "wcc": gen.Special.ProductionTool = "AlpGenProductionwqq" gen.Special.addTool(AlpGenProductionwqq) alpgen = gen.Special.AlpGenProductionwqq
# Import the necessary modules. from Configurables import Generation, MinimumBias, Special, PythiaProduction from Configurables import GenXiccProduction # Add Pythia as minimum bias production tool. Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "PythiaProduction" Generation().MinimumBias.addTool(PythiaProduction) # Add GenXicc as special production tool. Generation().addTool(Special) Generation().Special.ProductionTool = "GenXiccProduction" Generation().Special.addTool(GenXiccProduction) Generation().Special.GenXiccProduction.ShowerToolName = "PythiaProduction"
# Beam options as in Beam6500GeV-md100-2016-nu1.6.py importOptions("$APPCONFIGOPTS/Gauss/Beam7000GeV-mu100-nu7.6-HorExtAngle.py") from Configurables import Gauss, CondDB CondDB().Upgrade = True # Add also the DDDB tag for now from Configurables import LHCbApp LHCbApp().DDDBtag = "dddb-20190223" LHCbApp().CondDBtag = "sim-20180530-vc-mu100" LHCbApp().DataType = "Upgrade" # Configure the generation phase as usual Generation().EventType = 30000000 GenPhase().SampleGenerationTool = "MinimumBias" GenPhase().ProductionTool = "Pythia8ProductionMT" GenPhase().DecayTool = "EvtGenDecay" from Configurables import ToolSvc from Configurables import EvtGenDecay ToolSvc().addTool(EvtGenDecay) ToolSvc( ).EvtGenDecay.UserDecayFile = "$DECFILESROOT/dkfiles/Dsst_ppbarpi=TightCut.dec" ToolSvc().EvtGenDecay.DecayFile = "$DECFILESROOT/dkfiles/DECAY.DEC" from Configurables import Gauss, Gaussino Gaussino().ConvertEDM = True nthreads = 4 from Configurables import GiGaMT
#----------------------------------------------------------------------------- # Generator phase #----------------------------------------------------------------------------- GaussGen = GenInit("GaussGen") # # Set the random numbers - these fix the random seed. # GaussGen.FirstEventNumber = FIRST_EVENT_NUMBER GaussGen.RunNumber = RUN_NUMBER # The output is managed below, so we disable the standard Gauss output. Gauss().OutputType = 'NONE' Gauss().Histograms = 'NONE' # Switch off pileup (recommended!). Generation().addTool(FixedNInteractions, name="FixedNInteractions") Generation().FixedNInteractions.NInteractions = 0 Generation().PileUpTool = "FixedNInteractions" # Define "special" production. Generation().addTool(Special, "Special") Generation().SampleGenerationTool = "Special" # No cuts. Generation().Special.CutTool = "" # Define the production tool Generation().Special.addTool(ReadLHEfileProduction, name="ReadLHEfileProduction") Generation().Special.ProductionTool = "ReadLHEfileProduction"
from Configurables import Generation, MinimumBias, Pythia8Production from Configurables import Inclusive, SignalPlain, SignalRepeatedHadronization from Configurables import Special Pythia8TurnOffFragmentation = ["HadronLevel:all = off"] # this is only a preliminary tuning that should probably need to be updated Pythia8NewFrag = [ "StringZ:useNonstandardB = on", "StringZ:aNonstandardB = 1.01", "StringZ:bNonstandardB = 0.46", "StringZ:rFactB = 0", "StringZ:useNonstandardC = on", "StringZ:aNonstandardC = 1.075", "StringZ:bNonstandardC = 0.694", "StringZ:rFactC = 0" ] 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 += Pythia8NewFrag gen.addTool(Inclusive, name="Inclusive") gen.Inclusive.ProductionTool = "Pythia8Production" gen.Inclusive.addTool(Pythia8Production, name="Pythia8Production") gen.Inclusive.Pythia8Production.Tuning = "LHCbDefault.cmd" gen.Inclusive.Pythia8Production.Commands += Pythia8NewFrag gen.addTool(SignalPlain, name="SignalPlain") gen.SignalPlain.ProductionTool = "Pythia8Production" gen.SignalPlain.addTool(Pythia8Production, name="Pythia8Production")
from Configurables import Generation, MinimumBias Generation().SampleGenerationTool = "MinimumBias" Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "HijingProduction"
else: print 'using flat momentum spectrum!' from Configurables import MomentumRange ParticleGun().addTool(MomentumRange) from GaudiKernel import SystemOfUnits ParticleGun().MomentumRange.MomentumMin = 1.0 * SystemOfUnits.GeV from GaudiKernel import SystemOfUnits ParticleGun().MomentumRange.MomentumMax = 100. * SystemOfUnits.GeV ParticleGun().EventType = eventType ParticleGun().ParticleGunTool = "MomentumRange" ParticleGun().NumberOfParticlesTool = "FlatNParticles" #figure this out from Configurables import Generation pid_list = [] if hasattr(Generation(), 'SignalRepeatedHadronization'): pid_list = Generation().SignalRepeatedHadronization.SignalPIDList elif hasattr(Generation(), "SignalPlain"): pid_list = Generation().SignalPlain.SignalPIDList else: print 'major configuration problem, please fix!' import sys sys.exit() print 'got signal PID list', pid_list ParticleGun().MomentumRange.PdgCodes = pid_list ParticleGun().SignalPdgCode = abs(pid_list[0]) ParticleGun().DecayTool = "EvtGenDecay" else: importOptions('$DECFILESOPTS/%s.py' % eventType) importOptions("$LBPYTHIA8ROOT/options/Pythia8.py")
from Configurables import Gauss from Configurables import Generation from Configurables import Special from GaudiKernel import SystemOfUnits from Configurables import OniaPairsProduction Generation().SampleGenerationTool = "Special" Generation().addTool(Special) Generation().Special.ProductionTool = "OniaPairsProduction" Generation().Special.addTool(OniaPairsProduction) Generation().Special.OniaPairsProduction.Ecm = 2 * Gauss( ).BeamMomentum / SystemOfUnits.GeV
# file /nightlies/jenkins/workspace/nightly-slot-checkout/tmp/checkout/DBASE/Gen/DecFiles/v27r39/options/13244202.py generated: Mon, 09 Feb 2015 18:32:21 # # Event Type: 13244??? # # ASCII decay Descriptor: [Bs0 -> ( [chi_c1(1P) , chi_c2(1P)] -> (J/psi(1S) -> mu+ mu-) gamma ) pi+ pi- ]cc # from Configurables import Generation Generation().EventType = 13244??? Generation().SampleGenerationTool = "SignalRepeatedHadronization" from Configurables import SignalRepeatedHadronization Generation().addTool( SignalRepeatedHadronization ) Generation().SignalRepeatedHadronization.ProductionTool = "PythiaProduction" from Configurables import ToolSvc from Configurables import EvtGenDecay ToolSvc().addTool( EvtGenDecay ) ToolSvc().EvtGenDecay.UserDecayFile = "$DECFILESROOT/dkfiles/Bs_chicpipi,Jpsig,mm=TightCuts.dec" Generation().SignalRepeatedHadronization.CutTool = "LoKi::GenCutTool/TightCut" Generation().SignalRepeatedHadronization.SignalPIDList = [ 531,-531 ] # from Configurables import LoKi__GenCutTool gen = Generation() gen.SignalRepeatedHadronization.addTool ( LoKi__GenCutTool , 'TightCut' ) # tightCut = gen.SignalRepeatedHadronization.TightCut tightCut.Decay = '[ B_s0 -> ( Meson -> ^( J/psi(1S) => ^mu+ ^mu-) ^gamma ) pi+ pi-) ]CC' tightCut.Cuts = { 'gamma' : ' goodGamma ' , '[mu+]cc' : ' goodMuon ' , '[pi+]cc' : ' goodPion ' , 'J/psi(1S)' : ' goodPsi ' }
GaussGen.FirstEventNumber = 1 #--Number of events nEvts = 10 LHCbApp().EvtMax = nEvts # LHCbApp().OutputLevel = DEBUG #========================================================== from Gauss.Configuration import * from Configurables import Generation from Configurables import MinimumBias, Special from Configurables import Pythia8Gun from Configurables import FixedNInteractions Generation().PileUpTool = 'FixedNInteractions' Generation().addTool(FixedNInteractions) Generation().FixedNInteractions.NInteractions = 1 Generation().SampleGenerationTool = 'MinimumBias' Generation().addTool(MinimumBias) Generation().MinimumBias.ProductionTool = "Pythia8Gun" Generation().MinimumBias.addTool(Pythia8Gun) #========================================================== from GaudiKernel.SystemOfUnits import GeV # List of particle to be used in the gun. PDG-code. Generation().MinimumBias.Pythia8Gun.PDGList = [23] # Z-boson