Project

General

Profile

Bug #820 » field_configuration.py

Rogers, Chris, 14 December 2011 07:12

 
1
import os
2

    
3
geometry_root_dir = os.environ['ERIT_ROOT']
4

    
5
verbose_level = 0
6
physics_processes = "none"
7

    
8
simulation_geometry_filename = os.path.join(geometry_root_dir, "maus", "erit_geometry.dat")
9
geant4_visualisation = False
10
visualisation_viewer = "OGLIXm"
11
visualisation_theta = 90.
12
visualisation_phi = 90.
13

    
14
spill_generator_number_of_spills = 1
15
beam_defs = []
16
n_particles = 101
17
for particle_index in range(n_particles):
18
    if n_particles == 1:
19
        radius = 2350.
20
    else:
21
        radius = 2350.-260.+520.*float(particle_index)/float(n_particles-1)
22
    beam_defs.append({
23
        "reference":{
24
            "position":{"x":radius, "y":-0.0, "z":-500.0},
25
            "momentum":{"x":0.0, "y":0.0, "z":1.0},
26
            "particle_id":2112, "energy":10000., "time":0.0, "random_seed":1, # neutrons
27
         },
28
        "random_seed_algorithm":"incrementing_random", # algorithm for seeding MC
29
        "n_particles_per_spill":1, # probability of generating a particle
30
        "transverse":{
31
            "transverse_mode":"pencil", # transverse distribution matched to constant solenoid field
32
        },
33
        "longitudinal":{
34
            "longitudinal_mode":"pencil", # longitudinal distribution sawtooth in time
35
            "momentum_variable":"p", # Gaussian in total momentum (options energy, pz)
36
        }, # end time of sawtooth
37
        "coupling":{"coupling_mode":"none"} # no dispersion
38
    })
39

    
40
beam = {
41
    "particle_generator":"counter", # routine for generating empty primaries
42
    "random_seed":0, # random seed for beam generation; controls also how the MC
43
                     # seeds are generated
44
    "definitions":beam_defs
45
}
46

    
(4-4/5)