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, 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"
# 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")
#-----------------------------------------------------------------------------
# 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"
# 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
# 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
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" ]
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