###################################################
theApp.setup( MONTECARLO )
include( "PartPropSvc/PartPropSvc.py" )
EventSelector = Service( "EventSelector" )
EventSelector.RunNumber   = 11

EventSelector.FirstEvent  = 1
# load relevant libraries
theApp.Dlls += [ "GaudiAlg" ]
# ApplicationMgr.DLLs += { "Isajet_i"};
theApp.Dlls += [ "Pythia_i"]
theApp.Dlls += [ "GeneratorFilters" ]
theApp.Dlls += [ "GaudiAud" ]
#//ApplicationMgr.DLLs += { "AtlfastAlgs"};
theApp.Dlls += [ "HbookCnv" ]
theAuditorSvc = AuditorSvc()
theAuditorSvc.Auditors  = [ "ChronoAuditor" ]
# No histogramming at the Generator stage
#--------------------------------------------------------------
# Event related parameters
#--------------------------------------------------------------
# Number of events to be processed (default is 0)
theApp.EvtMax = 100
#--------------------------------------------------------------
# Private Application Configuration options
#--------------------------------------------------------------
# OUTPUT PRINTOUT LEVEL
# Set output level threshold (2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
# you can override this for individual modules if necessary
MessageSvc = Service( "MessageSvc" )
MessageSvc.OutputLevel               = 2
# Execution of algorithms given in sequencer
#
theApp.TopAlg = ["Sequencer/TopSequence"]
#
#GENERATORS SETUP (select by uncommenting/commenting)
#----------------
TopSequence = Algorithm( "TopSequence" )
TopSequence.Members  = ["Sequencer/Generator"]
#//TopSequence.Members  += {"Sequencer/Atlfast"};
TopSequence.Members  += ["Sequencer/Filter"]
Generator = Algorithm( "Generator" )
Generator.Members    = ["Pythia"]
Filter = Algorithm( "Filter" )
Filter.Members = ["MultiLeptonFilter"]
AtRndmGenSvc = Service( "AtRndmGenSvc" )
AtRndmGenSvc.Seeds = [" PYTHIA 4789849 989220611 ", " PYTHIA_INIT 824221
2327511 "]
###############################################################


Pythia = Algorithm( "Pythia" )

Pythia.PythiaCommand = [
    "pysubs msel 0",
    "pysubs msub 81 1",#qqbar->qqbar
    "pysubs msub 82 1",#gg->qqbar

    "pydat3 mdme 153 1 0",
    "pydat3 mdme 154 1 0",
    "pydat3 mdme 155 1 0",
    "pydat3 mdme 156 1 1",#charm
    "pydat3 mdme 157 1 0",#bottom
    "pydat3 mdme 158 1 0",
    "pydat3 mdme 159 1 0",
    "pydat3 mdme 160 1 0",
    "pydat3 mdme 161 1 0",
     #output!
    "pyinit pylisti 12",
    "pyinit pylistf 1",
    "pyinit dumpr 1 5",
    "pystat 1 2 3 4 5"]
#
# Filter to get 3 leptons
# chiedo anche un Ds!!!!!
#MultiLeptonFilter = Algorithm( "MultiLeptonFilter" )
#MultiLeptonFilter.NLeptons = 3
#MultiLeptonFilter.Etacut = 2.5
#MultiLeptonFilter.Ptcut = 3000.;  # Note this is 4 GeV
Generator.Members += ["EventCounter"]

#----------------------------------------------------------------
# Ntuple service output
#----------------------------------------------------------------

theApp.Dlls += [ 'RootHistCnv' ]                               #1
theApp.HistogramPersistency = 'ROOT'
#HistogramPersistencySvc = Service( "HistogramPersistencySvc" )
#HistogramPersistencySvc.OutputFile  = "atlas.root"
NTupleSvc = Service( 'NTupleSvc' )
NTupleSvc.Output = [ "FILE1 DATAFILE='ccbar.cbnt.root' OPT='NEW'" ]

include( 'CBNT_Athena/CBNT_Athena_jobOptions.py' )             #2
include( 'CBNT_Athena/CBNT_EventInfo_jobOptions.py' )          #3
include( 'RecExCommon/CBNT_Truth_jobOptions.py' )

CBNT_Athena.NtupleLocID = '/FILE1/CBNT/t3333'

###################################################