Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
LbJpsiL.py
Go to the documentation of this file.
1 
10 theApp.setup( MONTECARLO )
11 include( "PartPropSvc/PartPropSvc.py" )
12 #--------------------------------------------------------------
13 # Private Application Configuration options
14 #--------------------------------------------------------------
15 theApp.Dlls += [ "GaudiAlg" ]
16 theApp.Dlls += [ "TruthExamples", "PythiaB" ]
17 theApp.Dlls += [ "GaudiAud" ]
18 theApp.Dlls += [ "HbookCnv" ]
19 theApp.Dlls += [ "GeneratorFilters" ,"EvtGen_i" ]
20 theAuditorSvc = AuditorSvc()
21 theAuditorSvc.Auditors = [ "ChronoAuditor" ]
22 #--------------------------------------------------------------
23 # Algorithms
24 #--------------------------------------------------------------
25 theApp.TopAlg = ["PythiaB" , "EvtLbPolDecay", "LambdabFilter"]
26 #--------------------------------------------------------------
27 # Number of events to be accepted !! (default is 10)
28 # re-written if use B job submition script
29 # RunNumber, FirstEvent re-written if use B job submition script
30 #--------------------------------------------------------------
31 theApp.EvtMax = 1000
32 EventSelector.RunNumber = 1
33 EventSelector.FirstEvent = 1
34 #-------------------------------------------------------------- # PARAMETERS SPECIFIC TO PYTHIAB
35 #--------------------------------------------------------------
36 include( "EvtGen_i/StopPytWeakBdecays.py" )
37 #--------------------------------------------------------------
38 # -------- PYTHIA PARAMETERS OPTIMAL FOR BEAUTY PRODUCTION --
39 #--------------------------------------------------------------
40 # 'msel 5' is only for fast tests!
41 # for correct b-producion you should use 'msel 1'
42 include( "PythiaB/Btune.py" )
43 PythiaB = Algorithm( "PythiaB" )
44 PythiaB.PythiaCommand += ["pysubs ckin 3 9.",
45  "pysubs ckin 9 -3.5",
46  "pysubs ckin 10 3.5",
47  "pysubs ckin 11 -3.5",
48  "pysubs ckin 12 3.5",
49  "pysubs msel 1"]
50 # ------------- Selections on b quarks -------------
51 PythiaB.cutbq = ["0. 100. or 10. 3.5"]
52 # ------------- Number of repeated hadronization mhadr -------------
53 PythiaB.mhadr = 15
54 # ------------- Control flags on the userFinsel file -------------
55 PythiaB.ForceDecayChannel = "LambdabJpsimumuLambda"
56 PythiaB.DecayChannelParameters = [7.,3.] # Minimum pt of Lambdab = 7, and maximum abs(etha)=3.
57 # ------------- For how many events store B-chain in NTUPLE -------------
58 LambdabFilter = Algorithm( "LambdabFilter" )
59 LambdabFilter.SignaltoNtup = 20000
60 LambdabFilter.EtaCut = 2.7;
61 LambdabFilter.PtCut = 500.;
62 LambdabFilter.LowPtMuCut = 2500.;
63 LambdabFilter.HighPtMuCut = 4000.;
64 LambdabFilter.WriteNtuple = False
65 #--------------------------------------------------------------
66 # YOUR CONTROL OF EVTGEN
67 #--------------------------------------------------------------
68 # If you do not define a name to your Decay-table file here
69 # then EvtGen will work with a default decay table DECAY.DEC
70 EvtLbPolDecay = Algorithm( "EvtLbPolDecay" )
71 EvtLbPolDecay.userDecayTableName = "LAMBDAB.DEC"
72 # Here you set the polarization value
73 EvtLbPolDecay.EvtGenPol=0.4
74 # This control the production of histograms for certification
75 EvtLbPolDecay.MakeHisto=False
76 #--------------------------------------------------------------
77 # Histogram & Ntuple Persistency
78 #--------------------------------------------------------------
79 #Change the following line to "RootHistCnv" for ROOT persistency
80 theApp.Dlls += [ "RootHistCnv" ]
81 theApp.HistogramPersistency = "ROOT"
82 NtupleSvc = Service( "NtupleSvc" )
83 NtupleSvc.Output = ["FILE1 DATAFILE='ntuple_LbPol.root' TYP='ROOT' OPT='NEW'"]
84 #----------------------------------------------------------------
85 # Pool persistency
86 #----------------------------------------------------------------
87 include( "AthenaPoolCnvSvc/WriteAthenaPool_jobOptions.py" )
88 theApp.Dlls += [ "GeneratorObjectsAthenaPoolPoolCnv" ]
89 Stream1 = Algorithm( "Stream1" )
90 Stream1.ItemList += ["2101#*", "133273#*" ]
91 Stream1.AcceptAlgs += ["LambdabFilter"]
92 PoolSvc = Service( "PoolSvc" )
93 Stream1.OutputFile = "LbPol.root.pool"
94 #==============================================================
95 #
96 # End of job options file
97 #
98 
99 
python.Include.include
include
Definition: Include.py:318