← Back to branch summary

~ajdobbs/maus/rootfitter

« back to all changes in this revision

Viewing changes to bin/Global/datacard_mu_plus_through_tracks.py

  • Committer: Adam Dobbs
  • Date: 2017-02-28 17:58:27 UTC
  • mfrom: (1149.2.21 maus-primarychain)
  • Revision ID: phuccj@gmail.com-20170228175827-vhunrks5hzqsayb5
Merge: global recon updates

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