-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathsettings.xml
More file actions
121 lines (111 loc) · 6.63 KB
/
Copy pathsettings.xml
File metadata and controls
121 lines (111 loc) · 6.63 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
<Settings>
<!-- TODO: add separation User-changable and "do not touch"? -->
<!--OPT means value is optional. If it is not specified the default value/behaviour will be used-->
<PhysicalConstants>
<!--See global_definitions.h for unit and meanings description-->
<electron_charge_SI>1.60217662e-19</electron_charge_SI>
<electron_mass_SI>9.10938356e-31</electron_mass_SI>
<electron_mass_eV>5.109989461e5</electron_mass_eV>
<h_bar_SI>1.054571800e-34</h_bar_SI>
<h_bar_eVs>6.582119514e-16</h_bar_eVs>
<a_bohr_SI>5.2917721092e-1</a_bohr_SI> <!--multiplied by e10 for XS to be in 1e-20 m2-->
<Ry_energy_eV>13.605693</Ry_energy_eV>
<boltzmann_SI>1.38064852e-23</boltzmann_SI>
<light_speed_SI>299792458</light_speed_SI>
<Argon>
<Ar_mass_eV>3.726e10</Ar_mass_eV>
<Ar_primal_line_nm>128</Ar_primal_line_nm>
<MERT5_Lmax>10</MERT5_Lmax>
<!--see Kurokawa Phys. Rev. A84 2011, MERT5+ fit http://dx.doi.org/10.1103/PhysRevA.84.062717 -->
<MERT5_A>-1.365</MERT5_A>
<MERT5_D>80.5</MERT5_D>
<MERT5_F>-153</MERT5_F>
<MERT5_G>31.0</MERT5_G>
<MERT5_A1>8.8</MERT5_A1>
<MERT5_H>29.7</MERT5_H>
<MERT5_alpha_d>11.08</MERT5_alpha_d>
<MERT5_alpha_q>0.0</MERT5_alpha_q>
<elastic_phases_En_max>12.0</elastic_phases_En_max> <!--Could not find phaseshift values for larger E-->
<elastic_phases_En_min>5e-3</elastic_phases_En_min> <!--MERT5 phaseshifts are questionable below this E-->
<elastic_phases_En_threshold>0.24</elastic_phases_En_threshold> <!--MERT 5 is applied only below this energy for differential cross-section-->
<elastic_XS_En_max>20.0</elastic_XS_En_max> <!--No data for larger energy. They do exist, simply not in "./data/"-->
<elastic_XS_En_min>1e-3</elastic_XS_En_min> <!--linear extrapolation to XS(E=0) is used below En_min-->
<elastic_XS_En_threshold>1.0</elastic_XS_En_threshold>
<elastic_XS_En_smooth_width>0.1</elastic_XS_En_smooth_width>
<elastic_XS_no_Ramsauer_minimum>false</elastic_XS_no_Ramsauer_minimum> <!--OPT-->
<elastic_XS_Ramsauer_minimum_En>0.3286</elastic_XS_Ramsauer_minimum_En> <!--OPT if no Ramsauer minimum is false or not specified-->
<elastic_XS_at_0_En>7.491</elastic_XS_at_0_En>
<scattering_angle_model>0</scattering_angle_model> <!--OPT, 0 - normal, 1 - uniform-->
<time_delay_model>None</time_delay_model> <!--OPT, Precise - default - using dt = d/dE Im arg(f), Rough - using d/dE phaseshift, None - no time delay -->
<Feshbach_resonance_3o2_En>11.103</Feshbach_resonance_3o2_En>
<Feshbach_resonance_1o2_En>11.270</Feshbach_resonance_1o2_En>
<Feshbach_resonance_3o2_Width>2.3e-3</Feshbach_resonance_3o2_Width>
<Feshbach_resonance_1o2_Width>2.2e-3</Feshbach_resonance_1o2_Width>
<Dissoc_attachment_En_thresh>11.548</Dissoc_attachment_En_thresh>
<Dissoc_attachment_XS>1</Dissoc_attachment_XS>
<Argon_ion_decay_time>360e-9</Argon_ion_decay_time> <!--In seconds-->
<Feshbach_resonance_NBrS_XS>0.0</Feshbach_resonance_NBrS_XS><!--in e-20 m^2-->
<Feshbach_resonance_NBrS_En_loss>0.0</Feshbach_resonance_NBrS_En_loss> <!--OPT-->
</Argon>
</PhysicalConstants>
<ProgramConstants>
<is_test_version>false</is_test_version> <!--OPT for running tests before simulation. TODO: test modules and mb even submodules-->
<temperature>87</temperature> <!--Kelvins-->
<pressure>1.05e5</pressure> <!--Pascals-->
<thread_number>1</thread_number>
<maximal_energy>16.0</maximal_energy>
<angle_discretization>1001</angle_discretization>
<drift_distance_ignore_history>0</drift_distance_ignore_history> <!--OPT-->
<skip_history_rate>0</skip_history_rate> <!--OPT-->
<skip_history_time>2e-8</skip_history_time> <!--in s, OPT. Write history every ... seconds. Has higher priority than skip_history_rate, which will be ignored if skip_history_time exists and is not zero -->
<RecordedValues> <!-- if field is not specified, its value is taken as false by default (do not record) -->
<!-- TODO: add profiles: for E specrtras, for V drift, Custom -->
<process_type>false</process_type>
<time_initial>true</time_initial> <!-- Used for separation of dirrerent electrons! TODO: change TTree structure so that electrons are separate -->
<time_delta>true</time_delta>
<time_delta_full>true</time_delta_full>
<energy_initial>true</energy_initial>
<energy_collision>true</energy_collision>
<energy_final>false</energy_final>
<energy_average>false</energy_average> <!-- useless -->
<position_initial>false</position_initial> <!-- position - coordinate against electric field (along e drift direction)-->
<position_delta>false</position_delta>
<position_final>true</position_final> <!-- there is redundancy here -->
<path_delta>true</path_delta> <!-- full covered path in 3D -->
<theta_initial>false</theta_initial>
<theta_collision>false</theta_collision>
<theta_delta>false</theta_delta>
<theta_final>false</theta_final>
<photon_energy>false</photon_energy> <!-- required for NBrS study. -->
<deb_log_rand>false</deb_log_rand> <!-- TODO: 3D table solver does not support uncertainty tracking -->
<deb_solver_y_left>false</deb_solver_y_left>
<deb_solver_y_right>false</deb_solver_y_right>
<deb_solver_E_left>false</deb_solver_E_left>
<deb_solver_E_right>false</deb_solver_E_right>
</RecordedValues>
<data_location>data</data_location><!--OPT-->
<DataFiles>
<elastic_XS_data>data/ArScatteringCross.dat</elastic_XS_data> <!--OPT-->
<elastic_XS_phaseshifts_data>data/McEachranArPhaseShifts.dat</elastic_XS_phaseshifts_data> <!--OPT-->
<excitation_XS_data>data/ArExcitations_Magboltz.dat</excitation_XS_data> <!--OPT-->
<ionization_XS_data>data/ArIonizations_Magboltz.dat</ionization_XS_data> <!--OPT-->
</DataFiles>
<cache_data_folder>data_derived</cache_data_folder>
<test_folder>tests</test_folder> <!-- OPT -->
<output_file>Output/v01.t/eData_($Td)Td.root</output_file>
<!--Ec_spectrum_data>Output/v12.1/Ec_forms/($Td)Td.dat</Ec_spectrum_data> <!--OPT. If set, the simulation is run in fixed Ec spectrum mode.-->
<drift_distance>1e-4</drift_distance> <!--default value for each run-->
<n_electrons>10</n_electrons> <!--default value for each run-->
<random_seed>42</random_seed> <!--default value for each run-->
<random_generator>TRandom1</random_generator> <!-- TRandom1 (best, slowest), TRandom2 (BAD), TRandom3 (faster, good) -->
<Mixture>
<Argon>0.99</Argon>
<ArgonVanDerWaals>0.01</ArgonVanDerWaals>
</Mixture>
<Runs>
<Run>
<Td>7.0</Td>
</Run>
</Runs>
</ProgramConstants>
</Settings>