10 theApp.setup( MONTECARLO )
11 include(
"PartPropSvc/PartPropSvc.py" )
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" ]
25 theApp.TopAlg = [
"PythiaB" ,
"EvtLbPolDecay",
"LambdabFilter"]
32 EventSelector.RunNumber = 1
33 EventSelector.FirstEvent = 1
36 include(
"EvtGen_i/StopPytWeakBdecays.py" )
43 PythiaB = Algorithm(
"PythiaB" )
44 PythiaB.PythiaCommand += [
"pysubs ckin 3 9.",
47 "pysubs ckin 11 -3.5",
51 PythiaB.cutbq = [
"0. 100. or 10. 3.5"]
55 PythiaB.ForceDecayChannel =
"LambdabJpsimumuLambda"
56 PythiaB.DecayChannelParameters = [7.,3.]
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
70 EvtLbPolDecay = Algorithm(
"EvtLbPolDecay" )
71 EvtLbPolDecay.userDecayTableName =
"LAMBDAB.DEC"
73 EvtLbPolDecay.EvtGenPol=0.4
75 EvtLbPolDecay.MakeHisto=
False
80 theApp.Dlls += [
"RootHistCnv" ]
81 theApp.HistogramPersistency =
"ROOT"
82 NtupleSvc = Service(
"NtupleSvc" )
83 NtupleSvc.Output = [
"FILE1 DATAFILE='ntuple_LbPol.root' TYP='ROOT' OPT='NEW'"]
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"