Skip to content
Snippets Groups Projects
Commit 203ea0c6 authored by Dominic Jones's avatar Dominic Jones Committed by Dominic Jones
Browse files

add Cascadebc0 -> Lambda_b0+Ks and Sigmab+ -> Jpsi + Lambda + pi job options

parent 51892cdc
No related branches found
No related tags found
1 merge request!3154add Cascadebc0 -> Lambda_b0+Ks and Sigmab+ -> Jpsi + Lambda + pi job options
This commit is part of merge request !3154. Comments created here will be created in the context of that merge request.
- estimated CPU for CI job = 0.03 hrs
- using release = AthGeneration-23.6.29
- ecmEnergy = 13000.0
- randomSeed = 2113678
- EVNT to EVNT = False
- LHEonly = False
- platform = Linux-5.14.0-427.16.1.el9_4.x86_64-x86_64-with-glibc2.17
\ No newline at end of file
#--------------------------------------------------------------
#Pythia8 pp/gg -> Cascade_bc0 -> Lambda_b0 + K_S production
#--------------------------------------------------------------
evgenConfig.nEventsPerJob = 5000
evgenConfig.description = "Pythia8 pp/gg -> Cascade_bc0 -> Lambda_b0 + K_S"
evgenConfig.keywords = ["heavyFlavour","Jpsi","Muon"]
evgenConfig.contact = ["shaogang.peng@cern.ch"]
evgenConfig.process = "pp/gg -> X -> Cascade_bc0 -> Lambda_b0 + K_S"
evgenConfig.generators += ['Pythia8']
include('Pythia8_i/Pythia8_A14_NNPDF23LO_Common.py')
include("Pythia8_i/Pythia8_Photospp.py")
genSeq.Pythia8.Commands += [
'Higgs:useBSM = on',
'HiggsBSM:gg2H2 = on',
'HiggsH2:coup2d = 10.0',
'HiggsH2:coup2u = 10.0',
'HiggsH2:coup2Z = 0.0',
'HiggsH2:coup2W = 0.0',
'HiggsA3:coup2H2Z = 0.0',
'HiggsH2:coup2A3A3 = 0.0',
'HiggsH2:coup2H1H1 = 0.0',
'35:mMin = 0',
'35:mMax = 15.0',
############# For Fixed Mass Distribution#############
'35:m0 = 7.45', # Cascade_bc0,
'35:mWidth = 0.1',
'35:onMode = off',
'35:addChannel = 1 1.00 100 5122 310', #Lambda_b0 + K_S
'5122:onMode = off',
'5122:onIfMatch = 443 3122',
'443:onMode = off',
'443:onIfMatch = 13 13',
]
## Lambda and K_s have long lifetime and decay in Geant4
### Set lepton filters
if not hasattr(filtSeq, "MultiMuonFilter" ):
from GeneratorFilters.GeneratorFiltersConf import MultiMuonFilter
muonfilter1 = MultiMuonFilter("muonfilter1")
muonfilter2 = MultiMuonFilter("muonfilter2")
filtSeq += muonfilter1
filtSeq.muonfilter1.Ptcut = 3500.0 #MeV
filtSeq.muonfilter1.Etacut = 2.7
filtSeq.muonfilter1.NMuons = 2 #minimum
- estimated CPU for CI job = 0.16 hrs
- using release = AthGeneration-23.6.29
- ecmEnergy = 13000.0
- randomSeed = 12113
- EVNT to EVNT = False
- LHEonly = False
- platform = Linux-5.14.0-427.20.1.el9_4.x86_64-x86_64-with-glibc2.17
\ No newline at end of file
evgenConfig.description = "Pythia8B Sigma_b+ -> Jpsi + Lambda + pi"
evgenConfig.process = "Sigma_b+ -> Jpsi + Lambda + pi"
evgenConfig.keywords = ["muon","Jpsi","heavyFlavour"]
evgenConfig.generators = ["Pythia8B"]
evgenConfig.contact = ["shaogang.peng@cern.ch"]
evgenConfig.nEventsPerJob = 50
include("Pythia8B_i/Pythia8B_A14_NNPDF23LO_Common.py")
include("Pythia8B_i/Pythia8B_Photospp.py")
# Hard process
genSeq.Pythia8B.Commands += ['HardQCD:all = on'] # Equivalent of MSEL1
genSeq.Pythia8B.Commands += ['ParticleDecays:mixB = off']
genSeq.Pythia8B.Commands += ['HadronLevel:all = off']
# Event selection
genSeq.Pythia8B.SelectBQuarks = True
genSeq.Pythia8B.SelectCQuarks = False
genSeq.Pythia8B.VetoDoubleBEvents = True
#genSeq.Pythia8B.UserSelection = 'BJPSIINCLUSIVE'
genSeq.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 8.']
genSeq.Pythia8B.QuarkPtCut = 6.0
genSeq.Pythia8B.AntiQuarkPtCut = 6.0
genSeq.Pythia8B.QuarkEtaCut = 3.5
genSeq.Pythia8B.AntiQuarkEtaCut = 3.5
genSeq.Pythia8B.RequireBothQuarksPassCuts = False
#
genSeq.Pythia8B.NHadronizationLoops = 10 # 1 (old value)
# Bs
genSeq.Pythia8B.Commands += ['5222:onMode = off']
genSeq.Pythia8B.Commands += ['5222:addChannel = 1 1.00 0 443 3122 211']
genSeq.Pythia8B.Commands += ['443:onMode = off']
genSeq.Pythia8B.Commands += ['443:onIfMatch = 13 13']
genSeq.Pythia8B.SignalPDGCodes = [5222,443,13,-13,3122,211] ## need the full decay chain, select Sigma_b+ decay
genSeq.Pythia8B.SignalPtCuts = [0.0,0.0,3.5,3.5,0.0,0.48]
genSeq.Pythia8B.SignalEtaCuts = [102.5,102.5,2.7,2.7,102.5,2.6]
genSeq.Pythia8B.NumberOfSignalsRequiredPerEvent = 1
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment