Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
bbToJpsimu4mu4.py
Go to the documentation of this file.
1 
4 
5 # General evgen config options
6 from Pythia8B_i.Pythia8B_iConf import Pythia8B_i
7 topAlg+=Pythia8B_i("Pythia8B")
8 evgenConfig.generators += ["Pythia8B"]
9 evgenConfig.description = "Inclusive bb->J/psi(mu4mu4) production"
10 evgenConfig.keywords = ["BPHYS"]
11 evgenConfig.minevents = 500
12 
13 # Hard process
14 topAlg.Pythia8B.Commands += ['HardQCD:all = on'] # Uncomment for MSEL1
15 #topAlg.Pythia8B.Commands += ['HardQCD:gg2bbbar = on'] # Uncomment for MSEL5
16 #topAlg.Pythia8B.Commands += ['HardQCD:qqbar2bbbar = on'] # Uncomment for MSEL5
17 topAlg.Pythia8B.Commands += ['PhaseSpace:pTHatMin = 6.'] # Equivalent of CKIN3
18 topAlg.Pythia8B.Commands += ['ParticleDecays:mixB = off']
19 topAlg.Pythia8B.Commands += ['HadronLevel:all = off']
20 
21 # Quark cuts
22 topAlg.Pythia8B.SelectBQuarks = True
23 topAlg.Pythia8B.SelectCQuarks = False
24 topAlg.Pythia8B.QuarkPtCut = 0.0
25 topAlg.Pythia8B.AntiQuarkPtCut = 4.0
26 topAlg.Pythia8B.QuarkEtaCut = 102.5
27 topAlg.Pythia8B.AntiQuarkEtaCut = 2.5
28 topAlg.Pythia8B.RequireBothQuarksPassCuts = True
29 topAlg.Pythia8B.VetoDoubleBEvents = True
30 topAlg.Pythia8B.UserSelection = 'BJPSIINCLUSIVE'
31 
32 # Job options for closing (anti)B decays but opening B(antiB) decays
33 include ("Pythia8B_i/CloseBDecays.py")
34 #include ("Pythia8B_i/CloseAntiBDecays.py")
35 
36 # Job options for opening inclusive B->J/psi decays
37 include ("Pythia8B_i/OpenBJpsiDecays.py")
38 
39 # Close all J/psi decays apart from J/psi->mumu
40 topAlg.Pythia8B.Commands += ['443:onMode = off']
41 topAlg.Pythia8B.Commands += ['443:2:onMode = on']
42 
43 # Signal topology - only events containing this sequence will be accepted
44 topAlg.Pythia8B.SignalPDGCodes = [443,-13,13]
45 
46 # Number of repeat-hadronization loops
47 topAlg.Pythia8B.NHadronizationLoops = 1
48 
49 # List of B-species - for counting purposes (no effect on generation)
50 include("Pythia8B_i/BPDGCodes.py")
51 
52 # Final state selections
53 topAlg.Pythia8B.TriggerPDGCode = 13
54 topAlg.Pythia8B.TriggerStatePtCut = [4.0]
55 topAlg.Pythia8B.TriggerStateEtaCut = [2.5]
56 topAlg.Pythia8B.MinimumCountPerCut = [2]
57 
58 
59 #Uncomment for DumpMC print-out
60 #from TruthExamples.TruthExamplesConf import DumpMC
61 #topAlg += DumpMC()
python.Include.include
include
Definition: Include.py:318