Garfield++ 3.0
A toolkit for the detailed simulation of particle detectors based on ionisation measurement in gases and semiconductors
Loading...
Searching...
No Matches
PhotoAbsCSLib.cpp
Go to the documentation of this file.
3#ifdef NOT_INCLUDE_GASLIB_IN_PACSLIB
4#else
6#endif
7
8#include <iostream>
9
10// 2004, I. Smirnov
11
12/*
13Remark 14.10.2005:
14
15The text in this file does not depend on wcpplib/matter/GasLib.c
16or wcpplib/matter/GasLib.h.
17The first was included earlier when it was thought to depend.
18But there is anyway a sense in such inclusion.
19If the user needs to provide sertain order of initialization,
20he will need to include both PhotoAbsCSLib.c and GasLib.c
21in certain his source file.
22If GasLib.c is already included here, there is no need to include
23it again there and there will be only one line of initializations.
24*/
25
26namespace {
27
28std::string getDataBasePath() {
29
30 std::string path;
31 // First try if the environment variable HEED_DATABASE is defined.
32 char* heed_database = std::getenv("HEED_DATABASE");
33 if (heed_database == NULL) {
34 // If HEED_DATABASE is not defined, try GARFIELD_HOME.
35 char* garfield_home = std::getenv("GARFIELD_HOME");
36 if (garfield_home == NULL) {
37 std::cerr << "Heed:\n Could not retrieve database path.\n";
38 } else {
39 path = std::string(garfield_home) + "/Heed/heed++/database";
40 }
41 } else {
42 path = std::string(heed_database);
43 }
44 if (!path.empty()) {
45 std::cout << "Heed:\n Database path: " << path << "\n";
46 }
47 return path;
48}
49
50Heed::ExAtomPhotoAbsCS generate_Ar_PACS(const std::string& shelllist_dir,
51 const std::string& pacs_table_dir) {
52
53 Heed::ExAtomPhotoAbsCS Argon_PACS_mod_esc(18, shelllist_dir + "shelllist.dat",
54 pacs_table_dir + "Ar.dat");
55
56 // ExAtomPhotoAbsCS Argon_PACS_mod_esc(18,
57 // shelllist_dir + "shelllist.dat",
58 // shelllist_dir + "mw3.dat");
59
60 // ExAtomPhotoAbsCS Argon_PACS_mod_esc(18, "argon",
61 // shelllist_dir + "ftbf18.dat", 2);
62
63 Heed::AtomicSecondaryProducts* asp = Argon_PACS_mod_esc.get_asp(1);
64 std::vector<double> electron_energy;
65 std::vector<double> photon_energy;
66 electron_energy.emplace_back(0.000200);
67 // electron_energy.emplace_back(0.002670);
68 asp->add_channel(0.65, electron_energy, photon_energy);
69 electron_energy.resize(2);
70 electron_energy[0] = 0.000050;
71 electron_energy[1] = 0.000200;
72 asp->add_channel(0.35, electron_energy, photon_energy, 1);
73 // mcout<<"L1:\n";
74 // asp->print(mcout, 2);
75
76 asp = Argon_PACS_mod_esc.get_asp(2);
77 electron_energy.resize(1);
78 electron_energy[0] = 0.000200;
79 asp->add_channel(1.0, electron_energy, photon_energy, 1);
80 // mcout<<"L2:\n";
81 // asp->print(mcout, 2);
82
83 asp = Argon_PACS_mod_esc.get_asp(3);
84 electron_energy.resize(1);
85 electron_energy[0] = 0.000200;
86 asp->add_channel(1.0, electron_energy, photon_energy, 1);
87 // mcout<<"L3:\n";
88 // asp->print(mcout, 2);
89
90 return Argon_PACS_mod_esc;
91}
92}
93
94namespace Heed {
95
96using CLHEP::gram;
97using CLHEP::mole;
98
99const std::string shelllist_dir_name = getDataBasePath() + "/";
100const std::string pacs_table_dir_name = shelllist_dir_name + "henke/";
101
102// Hydrogen
103SimpleAtomPhotoAbsCS Hydrogen_PACS(1, std::make_shared<HydrogenPhotoAbsCS>());
104
105// SimpleTablePhotoAbsCS Hydrogen_for_CH4_shell_PACS("Hydrogen_for_CH4",
106// 1, 12.65e-6,
107// shelllist_dir_name + "H_for_CH4.dat");
109 std::make_shared<PhenoPhotoAbsCS>("Hydrogen_for_H2", 1, 15.43e-6, 3.228));
111 std::make_shared<PhenoPhotoAbsCS>("Hydrogen_for_CH4", 1, 12.65e-06, 3.228));
113 std::make_shared<PhenoPhotoAbsCS>("Hydrogen_for_NH4", 1, 10.0e-06, 3.228));
114
115// ExAtomPhotoAbsCS Hydrogen_for_H2_PACS(1,
116// shelllist_dir_name + "shelllist.dat",
117// pacs_table_dir_name + "H.dat");
118// ExAtomPhotoAbsCS Hydrogen_for_CH4_PACS(1,
119// shelllist_dir_name + "shelllist.dat",
120// shelllist_dir_name + "H_for_CH4.dat",
121// "H_for_CH4",
122// 12.65e-6);
123// ExAtomPhotoAbsCS Hydrogen_for_NH4_PACS(1,
124// shelllist_dir_name + "shelllist.dat",
125// pacs_table_dir_name + "H.dat",
126// 10.0e-06);
127
129 pacs_table_dir_name + "He.dat");
131 pacs_table_dir_name + "Li.dat");
133 pacs_table_dir_name + "Be.dat");
135 pacs_table_dir_name + "B.dat");
137 pacs_table_dir_name + "C.dat");
138// for debug, FitBT
139// ExAtomPhotoAbsCS Carbon_PACS(6, "carbon",
140// shelllist_dir_name + "shelltscf.dat",
141// 2, 0, 0.0);
142
143// ExAtomPhotoAbsCS Carbon_for_CH4_PACS(6,
144// shelllist_dir_name + "shelllist.dat",
145// pacs_table_dir_name + "C.dat",
146// "C_for_CH4",
147// 12.65e-06);
149 shelllist_dir_name + "C_for_CH4.dat",
150 "C_for_CH4", 12.65e-6);
152 pacs_table_dir_name + "C.dat",
153 "C_for_C2H4", 10.51e-06);
155 pacs_table_dir_name + "C.dat",
156 "C_for_C2H6", 11.52e-06);
158 pacs_table_dir_name + "C.dat",
159 "C_for_C4H10", 10.55e-06);
161 shelllist_dir_name + "shelllist.dat",
162 pacs_table_dir_name + "C.dat",
163 "C_for_Methylal", 10.0e-06);
165 pacs_table_dir_name + "C.dat", "C_for_CF4",
166 16.23e-06);
168 pacs_table_dir_name + "C.dat", "C_for_CO2",
169 13.79e-06);
171 pacs_table_dir_name + "N.dat", "N_for_N2",
172 15.581e-6);
174 pacs_table_dir_name + "O.dat");
176 pacs_table_dir_name + "O.dat", "O_for_CO2",
177 13.79e-6);
179 pacs_table_dir_name + "F.dat");
180ExAtomPhotoAbsCS Neon_PACS(10, shelllist_dir_name + "shelllist.dat",
181 pacs_table_dir_name + "Ne.dat");
183 pacs_table_dir_name + "Na.dat");
185 pacs_table_dir_name + "Mg.dat");
187 pacs_table_dir_name + "Al.dat");
189 pacs_table_dir_name + "Si.dat");
191 "shelllist_solid.dat",
192 pacs_table_dir_name + "Si.dat",
193 "Si_crystal");
194ExAtomPhotoAbsCS Silicon_G4_PACS(14, shelllist_dir_name + "shelllist_solid.dat",
195 shelllist_dir_name + "Si_G4.dat", "Si_G4");
197 pacs_table_dir_name + "P.dat");
199 pacs_table_dir_name + "S.dat");
201 pacs_table_dir_name + "Cl.dat");
202// "Standard" Argon:
203// ExAtomPhotoAbsCS Argon_PACS(18,
204// shelllist_dir_name + "shelllist.dat",
205// pacs_table_dir_name + "Ar.dat");
206// Optional variants:
207// ExAtomPhotoAbsCS Argon_PACS(18,
208// shelllist_dir_name + "shelllist.dat",
209// shelllist_dir_name + "mw3.dat");
210// Variant for debug, pointwise cross section
211// ExAtomPhotoAbsCS Argon_PACS(18, "argon",
212// shelllist_dir_name + "ftbf18.dat", 2);
213// Variant for debug, fitted cross section
214// ExAtomPhotoAbsCS Argon_PACS(18, "argon",
215// shelllist_dir_name + "shelltscf.dat",
216// 2, 0, 0.0);
217// Variant for debug, fitted cross section with replacement from Henke
218// ExAtomPhotoAbsCS Argon_PACS(18, "argon",
219// shelllist_dir_name + "shelltscf.dat",
220// pacs_table_dir_name + "Ar.dat",
221// 40.0e-6, 2, 0.0);
222// Another variant for debug, fitted cross section with replacement from
223// Marr and West, should be similar to old Fortran verion
224// ExAtomPhotoAbsCS Argon_PACS(18, "argon",
225// shelllist_dir_name + "shelltscf.dat",
226// shelllist_dir_name + "mw3.dat",
227// 40.0e-6, 2, 0.0);
228
230 generate_Ar_PACS(shelllist_dir_name, pacs_table_dir_name);
231
233 pacs_table_dir_name + "Ga.dat");
235 pacs_table_dir_name + "Ga.dat",
236 "Ga_for_GaAs");
238 pacs_table_dir_name + "Ge.dat");
240 shelllist_dir_name + "shelllist.dat",
241 pacs_table_dir_name + "Ge.dat",
242 "Ge_crystal", 0.67e-06);
244 pacs_table_dir_name + "As.dat");
246 pacs_table_dir_name + "As.dat",
247 "As_for_GaAs");
249 pacs_table_dir_name + "Br.dat");
251 pacs_table_dir_name + "Kr.dat");
253 pacs_table_dir_name + "Cd.dat");
255 pacs_table_dir_name + "Cd.dat",
256 "Cd_for_CdTe");
258 pacs_table_dir_name + "Te.dat");
260 shelllist_dir_name + "shelllist_solid.dat",
261 pacs_table_dir_name + "Te.dat",
262 "Te_for_CdTe");
264 pacs_table_dir_name + "Xe.dat");
266 pacs_table_dir_name + "Cs.dat");
268 pacs_table_dir_name + "Hg.dat");
270 pacs_table_dir_name + "U.dat");
271
273// MolecPhotoAbsCS H2_MPACS(Hydrogen_PACS, 2);
284 33.0e-6);
285
286// MolecPhotoAbsCS CH4_MPACS(&Carbon_for_CH4_PACS, 1,
287// &Hydrogen_for_H2_PACS, 4, 27.3e-6);
289 27.3e-6);
291
292// !!! The following line may need to be refined
293// (to adjust outer shell energies).
295// !!! The following line may need to be refined
296// (to adjust outer shell energies).
298
300 25.8e-6);
302 25.8e-6);
304 25.0e-6);
306 24.0e-6);
309 23.4e-6);
310
311// !!! The following line may need to be refined
312// (to adjust outer shell energies).
315
316// MolecPhotoAbsCS C2H2_MPACS(&Carbon_for_CH4_PACS, 2, &Hydrogen_for_CH4_PACS, 2);
317// MolecPhotoAbsCS C2H4_MPACS(&Carbon_for_CH4_PACS, 2, &Hydrogen_for_CH4_PACS, 4);
318// MolecPhotoAbsCS C2H6_MPACS(&Carbon_for_CH4_PACS, 2, &Hydrogen_for_CH4_PACS, 6);
319// MolecPhotoAbsCS C3H8_MPACS(&Carbon_for_CH4_PACS, 3, &Hydrogen_for_CH4_PACS, 8);
320// MolecPhotoAbsCS C4H10_MPACS(&Carbon_for_CH4_PACS, 4, &Hydrogen_for_CH4_PACS,
321// 10);
324 10.0e-6 * 23.4 / 10.55); // similar to C4H10
325/*
326The value of W for noble gases is slightly less than
327twice the ionization potential.
328For organic gases it is very close to mean ionization potential
329averaged with taking into account of atomic charges of carbon and hydrogen.
330and assuming that the ionization potential of the hydrogen is the same
331as in pure molecular hydrogen H2.
332*/
333
334// Additional molecular photoabsorption-cross sections
335// for consistency with Magboltz
336// Where available, the W values are taken from ICRU report 31
338 23.2e-6);
343 &Oxygen_PACS, 1);
347 &Oxygen_PACS, 1, 24.7e-6);
349 &Oxygen_PACS, 1, 24.8e-6);
351 &Oxygen_PACS, 1);
356 1);
359 &Fluorine_PACS, 5);
363 &Fluorine_PACS, 3);
365 &Bromine_PACS, 1);
372}
void add_channel(double fchannel_prob_dens, const std::vector< double > &felectron_energy, const std::vector< double > &fphoton_energy, int s_all_rest=0)
Definition: PhotoAbsCS.cpp:469
Atomic photo-absorption with excitation.
Definition: PhotoAbsCS.h:430
Definition: BGMesh.cpp:6
ExAtomPhotoAbsCS Gallium_for_GaAs_PACS(31, shelllist_dir_name+"shelllist_solid.dat", pacs_table_dir_name+"Ga.dat", "Ga_for_GaAs")
Definition: PhotoAbsCSLib.h:51
ExAtomPhotoAbsCS Tellurium_for_CdTe_PACS(52, shelllist_dir_name+"shelllist_solid.dat", pacs_table_dir_name+"Te.dat", "Te_for_CdTe")
Definition: PhotoAbsCSLib.h:61
SimpleAtomPhotoAbsCS Hydrogen_for_CH4_PACS(1, std::make_shared< PhenoPhotoAbsCS >("Hydrogen_for_CH4", 1, 12.65e-06, 3.228))
Definition: PhotoAbsCSLib.h:18
ExAtomPhotoAbsCS Phosphorus_PACS(15, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"P.dat")
Definition: PhotoAbsCSLib.h:46
SimpleAtomPhotoAbsCS Hydrogen_for_H2_PACS(1, std::make_shared< PhenoPhotoAbsCS >("Hydrogen_for_H2", 1, 15.43e-6, 3.228))
MolecPhotoAbsCS C2H2_MPACS
const std::string pacs_table_dir_name
ExAtomPhotoAbsCS Mercury_PACS(80, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Hg.dat")
MolecPhotoAbsCS Methylal_MPACS
MolecPhotoAbsCS CH3OH_MPACS
ExAtomPhotoAbsCS Magnesium_PACS(12, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Mg.dat")
Definition: PhotoAbsCSLib.h:41
MolecPhotoAbsCS C2F4H2_MPACS
MolecPhotoAbsCS CF4_MPACS
ExAtomPhotoAbsCS Carbon_for_C2H4_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_C2H4", 10.51e-06)
ExAtomPhotoAbsCS Silicon_crystal_PACS(14, shelllist_dir_name+"shelllist_solid.dat", pacs_table_dir_name+"Si.dat", "Si_crystal")
Definition: PhotoAbsCSLib.h:44
MolecPhotoAbsCS CS2_MPACS
MolecPhotoAbsCS NH3_MPACS
ExAtomPhotoAbsCS Nitrogen_PACS(7, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"N.dat", "N_for_N2", 15.581e-6)
MolecPhotoAbsCS CF3Br_MPACS
MolecPhotoAbsCS Cs_MPACS
MolecPhotoAbsCS C2H6_MPACS
MolecPhotoAbsCS O2_MPACS
ExAtomPhotoAbsCS Beryllium_PACS(4, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Be.dat")
Definition: PhotoAbsCSLib.h:25
MolecPhotoAbsCS SiH4_MPACS
MolecPhotoAbsCS SF6_MPACS
MolecPhotoAbsCS F2_MPACS
MolecPhotoAbsCS C2H2F4_MPACS
ExAtomPhotoAbsCS Cadmium_PACS(48, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Cd.dat")
Definition: PhotoAbsCSLib.h:58
ExAtomPhotoAbsCS Bromine_PACS(35, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Br.dat")
Definition: PhotoAbsCSLib.h:56
ExAtomPhotoAbsCS Krypton_PACS(36, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Kr.dat")
MolecPhotoAbsCS H2O_MPACS
MolecPhotoAbsCS SF4_MPACS
MolecPhotoAbsCS C2H4_MPACS
ExAtomPhotoAbsCS Fluorine_PACS(9, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"F.dat")
MolecPhotoAbsCS C2F6_MPACS
MolecPhotoAbsCS CO2_MPACS
MolecPhotoAbsCS O3_MPACS
MolecPhotoAbsCS N2_MPACS
ExAtomPhotoAbsCS Boron_PACS(5, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"B.dat")
ExAtomPhotoAbsCS Oxygen_PACS(8, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"O.dat")
MolecPhotoAbsCS Hg_MPACS
ExAtomPhotoAbsCS Arsenic_PACS(33, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"As.dat")
Definition: PhotoAbsCSLib.h:54
MolecPhotoAbsCS COS_MPACS
ExAtomPhotoAbsCS Carbon_for_C2H6_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_C2H6", 11.52e-06)
ExAtomPhotoAbsCS Silicon_PACS(14, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Si.dat")
MolecPhotoAbsCS N2O_MPACS
ExAtomPhotoAbsCS Silicon_G4_PACS(14, shelllist_dir_name+"shelllist_solid.dat", shelllist_dir_name+"Si_G4.dat", "Si_G4")
Definition: PhotoAbsCSLib.h:45
ExAtomPhotoAbsCS Sodium_PACS(11, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Na.dat")
Definition: PhotoAbsCSLib.h:40
MolecPhotoAbsCS Kr_MPACS
ExAtomPhotoAbsCS Carbon_for_C4H10_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_C4H10", 10.55e-06)
ExAtomPhotoAbsCS Arsenic_for_GaAs_PACS(33, shelllist_dir_name+"shelllist_solid.dat", pacs_table_dir_name+"As.dat", "As_for_GaAs")
Definition: PhotoAbsCSLib.h:55
MolecPhotoAbsCS He_MPACS
ExAtomPhotoAbsCS Argon_PACS
ExAtomPhotoAbsCS Tellurium_PACS(52, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Te.dat")
Definition: PhotoAbsCSLib.h:60
ExAtomPhotoAbsCS Helium_PACS(2, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"He.dat")
ExAtomPhotoAbsCS Sulfur_PACS(16, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"S.dat")
ExAtomPhotoAbsCS Oxygen_for_CO2_PACS(8, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"O.dat", "O_for_CO2", 13.79e-6)
Definition: PhotoAbsCSLib.h:37
MolecPhotoAbsCS C2HF5_MPACS
ExAtomPhotoAbsCS Neon_PACS(10, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Ne.dat")
ExAtomPhotoAbsCS Carbon_for_CO2_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_CO2", 13.79e-06)
ExAtomPhotoAbsCS Aluminium_PACS(13, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Al.dat")
Definition: PhotoAbsCSLib.h:42
ExAtomPhotoAbsCS Chlorine_PACS(17, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Cl.dat")
Definition: PhotoAbsCSLib.h:48
ExAtomPhotoAbsCS Germanium_PACS(32, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Ge.dat")
ExAtomPhotoAbsCS Cadmium_for_CdTe_PACS(48, shelllist_dir_name+"shelllist_solid.dat", pacs_table_dir_name+"Cd.dat", "Cd_for_CdTe")
Definition: PhotoAbsCSLib.h:59
const std::string shelllist_dir_name
MolecPhotoAbsCS BF3_MPACS
ExAtomPhotoAbsCS Germanium_crystal_PACS(32, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Ge.dat", "Ge_crystal", 0.67e-06)
Definition: PhotoAbsCSLib.h:53
MolecPhotoAbsCS H2_MPACS
MolecPhotoAbsCS H2S_MPACS
ExAtomPhotoAbsCS Lithium_PACS(3, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Li.dat")
Definition: PhotoAbsCSLib.h:24
ExAtomPhotoAbsCS Xenon_PACS(54, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Xe.dat")
MolecPhotoAbsCS DME_MPACS
ExAtomPhotoAbsCS Uranium_PACS(92, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"U.dat")
Definition: PhotoAbsCSLib.h:65
MolecPhotoAbsCS C3H7OH_MPACS
MolecPhotoAbsCS C3F8_MPACS
MolecPhotoAbsCS NO_MPACS
MolecPhotoAbsCS Ne_MPACS
MolecPhotoAbsCS C3H6_MPACS
MolecPhotoAbsCS C4H10_MPACS
MolecPhotoAbsCS Xe_MPACS
ExAtomPhotoAbsCS Gallium_PACS(31, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Ga.dat")
Definition: PhotoAbsCSLib.h:50
ExAtomPhotoAbsCS Carbon_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat")
Definition: PhotoAbsCSLib.h:27
MolecPhotoAbsCS Ar_MPACS
MolecPhotoAbsCS C3H8_MPACS
MolecPhotoAbsCS C5H12_MPACS
SimpleAtomPhotoAbsCS Hydrogen_for_NH4_PACS(1, std::make_shared< PhenoPhotoAbsCS >("Hydrogen_for_NH4", 1, 10.0e-06, 3.228))
Definition: PhotoAbsCSLib.h:19
ExAtomPhotoAbsCS Caesium_PACS(55, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"Cs.dat")
SimpleAtomPhotoAbsCS Hydrogen_PACS(1, std::make_shared< HydrogenPhotoAbsCS >())
Definition: PhotoAbsCSLib.h:16
MolecPhotoAbsCS C2H5OH_MPACS
MolecPhotoAbsCS CHF3_MPACS
MolecPhotoAbsCS CH4_MPACS
MolecPhotoAbsCS CO_MPACS
ExAtomPhotoAbsCS Carbon_for_CF4_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_CF4", 16.23e-06)
ExAtomPhotoAbsCS Carbon_for_Methylal_PACS(6, shelllist_dir_name+"shelllist.dat", pacs_table_dir_name+"C.dat", "C_for_Methylal", 10.0e-06)
MolecPhotoAbsCS GeH4_MPACS
ExAtomPhotoAbsCS Carbon_for_CH4_PACS(6, shelllist_dir_name+"shelllist.dat", shelllist_dir_name+"C_for_CH4.dat", "C_for_CH4", 12.65e-6)