← Back to branch summary

~christopher-hunt08/maus/maus_integrated_kalman

« back to all changes in this revision

Viewing changes to bin/Global/datacard_mu_plus_through_tracks.py

merging in changes in merge branch

Show diffs side-by-side

added added

removed removed

Lines of Context:
bin/Global/datacard_mu_plus_through_tracks.py
 
1
"""
 
2
Datacard to run with simulation of 200MeV/c muons
 
3
"""
 
4
# pylint: disable=C0103, W0611
 
5
 
 
6
import os
 
7
job_number = int(os.getenv("LSB_JOBINDEX", 10))
 
8
 
 
9
output_json_file_name = "200MeV_mu_plus_through" + str(job_number) + ".json"
 
10
output_json_file_type = "text"
 
11
efficiency_filename = "efficiency" + str(job_number) + ".csv"
 
12
# Used by MapPyRemoveTracks.
 
13
keep_only_muon_tracks = False
 
14
 
 
15
# Used by MapCppSimulation
 
16
keep_tracks = True # set to true to keep start and end point of every track
 
17
keep_steps = False # set to true to keep start and end point of every track and
 
18
                   # every step point
 
19
simulation_geometry_filename = "Stage4_NoField.dat"
 
20
maximum_number_of_steps = 10000
 
21
simulation_reference_particle = {
 
22
    "position":{"x":0.0, "y":-0.0, "z":10500.0},
 
23
    "momentum":{"x":0.0, "y":0.0, "z":1.0},
 
24
    "particle_id":-13, "energy":266.0, "time":0.0, "random_seed":job_number
 
25
}
 
26
 
 
27
muon_half_life = 1.e24
 
28
# used by InputPySpillGenerator to determine the number of empty spills that
 
29
# will be generated by the simulation
 
30
spill_generator_number_of_spills = 100
 
31
# used by MapPyBeamMaker to generate input particle data
 
32
# This is a sample beam distribution based on guesses by Chris Rogers of what
 
33
# an optimised beam might look like
 
34
beam = {
 
35
    "particle_generator":"counter", # routine for generating empty primaries
 
36
    "random_seed":job_number, # random seed for beam generation; controls also
 
37
                              # how the MC seeds are generated
 
38
    "definitions":[
 
39
    {
 
40
       "reference":simulation_reference_particle, # reference particle
 
41
       "random_seed_algorithm":"incrementing_random", # algorithm for seeding MC
 
42
       "n_particles_per_spill":1, # probability of generating a particle
 
43
       "transverse":{
 
44
          "transverse_mode":"twiss",
 
45
          "beta_x":8000.,
 
46
          "alpha_x":1.,
 
47
          "emittance_x":3.0, # 4d emittance
 
48
          "beta_y":8000.,
 
49
          "alpha_y":1.0,
 
50
          "emittance_y":3.0,
 
51
          },
 
52
       "longitudinal":{"longitudinal_mode":"sawtooth_time",
 
53
                   "momentum_variable":"p",
 
54
                   "sigma_p":25., # RMS total momentum
 
55
                   "t_start":-1.e6, # start time of sawtooth
 
56
                   "t_end":+1.e6}, # end time of sawtooth
 
57
       "coupling":{"coupling_mode":"none"} # no dispersion
 
58
    }]
 
59
}
 
60
 
 
61
particle_decay = True
 
62
physics_processes = "mean_energy_loss"
 
63
TOF_calib_by = "date"
 
64
TOF_calib_date_from = "2013-10-01 00:00:00"
 
65
SciFiPRHelicalOn = False