###############################################################
#
# Job options file for testing ParticleGenerator.
#
#==============================================================
theApp.setup( MONTECARLO )

#--------------------------------------------------------------
# Private Application Configuration options
#--------------------------------------------------------------
#load relevant libraries
include( "PartPropSvc/PartPropSvc.py" )

theApp.Dlls  += [ "ParticleGenerator" ]      
theApp.TopAlg = [ "ParticleGenerator" ]

# The following is needed to load the Athena Random
# Number Generation Service.
theApp.ExtSvc += ["AtRndmGenSvc"]
AtRndmGenSvc = Service( "AtRndmGenSvc" )
AtRndmGenSvc.Seeds = ["SINGLE 2040160768 443921183"]
# AtRndmGenSvc.ReadFromFile = TRUE

#--------------------------------------------------------------
# Set output level threshold (1=VERBOSE, 2=DEBUG, 3=INFO, 4=WARNING, 5=ERROR, 6=FATAL )
#--------------------------------------------------------------
MessageSvc.OutputLevel      = 3

#--------------------------------------------------------------
# Event related parameters
#--------------------------------------------------------------
StoreGateSvc = Service( "StoreGateSvc" )
StoreGateSvc.Dump = FALSE

#--------------------------------------------------------------
# ParticleGenerator parameters
#--------------------------------------------------------------

# e- = 11, e+ = -11, mu- = 13, mu+ = -13, pi+ = 211, pi- = -211, pi0 = 111, geantino = 999
# See PDGTABLE or Rev.Part.Properties for the complete list.

# The orders are in the form of an array of strings.  Each string
# must contain exactly one command.  In most cases, if ParticleGenerator
# can't figure out what a command means, it will terminate during its
# initialization.  Energies are in GeV, distances and time (ct) are in mm.

# If there is more than one definition of 'ParticleGenerator.orders',
# only the last one will be used.

# For an example of how to use user-defined histograms to generate a
# distribution, see jobOptions_ParticleGeneratorExampleHistogram.py.

ParticleGenerator = Algorithm( "ParticleGenerator" )

ParticleGenerator.orders = [
 "PDGcode: constant 11",
 "energy: constant 629000",
 "vertX: gauss 0.0 10",
 "vertY: gauss 0.0 10",
 "vertZ: gauss 0.0 50",
 "targetX: constant 420",
 "targetY: constant -4.2e+02",
 "targetZ: fixed 52000E-01",
 "Time=flat -1 1"
]
ParticleGenerator.orders = [
 "pdgcode: constant -11",
 "energy: constant 543000",
 "momX: fixed 1",
 "momY: Fixed 1",
 "momZ: fiXed 1"
]
ParticleGenerator.orders = [
 "id: constant 13",
 "Et: fixed 50000.0",
 "eta: sequence 0.1 0.3 0.6 1.2 2.8 3.0 5",
 "phi: Loop after Eta 0 6.283 .785398"
]
ParticleGenerator.orders = [
 "id: sequence -13 13",
 "Et: fixed 50000.0",
 "eta: fixed 0.5",
 "phi: Loop 5 -6 -1"
]
ParticleGenerator.orders = [
 "id: fixed 22 gamma",
 "Energy=g00fy This is ignored 10000 By ParticleGenerator 1 Most of the Time",
 "Et: constant 10000",
 "VertZ: fixed 0",
 "Pz: sequence -1 1",
 "Phi: Loop after Pz 0 6.283 0.785398"
]
#====Run this.==== 13 myu-minus 
ParticleGenerator.orders = [
 "id: fixed 13",
 "Et: fixed 50000.0",
 "eta: uniform 1.0 2.7",
 "phi: uniform 0.0 6.2832"
]
#==============================================================
#
# End of job options file
#
###############################################################