Hello!
I have written before about trying to implement the xray photon
generator.
This time I've updated all the libraries, YT, Heasoft... etc.
I'm just trying to generate a photon list for the moment but I'm
getting two sorts of errors:
1)
---------------------------------------------------------------------------
RuntimeError Traceback (most
recent call last)
/home/barbara/hires_final/photon_trial.py in <module>()
90 cosmo = Cosmology()
91
---> 92 photons = PhotonList.from_scratch(Roi, redshift,
A, exp_time,thermal_model,center="center",cosmology=cosmo)
93 photons.write_h5_file("my_photons.h5")
94
/home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_simulator.pyc
in from_scratch(cls, data_source, redshift, area, exp_time,
photon_model, parameters, center, dist, cosmology)
329 parameters["Width"] =
2.*width.in_units("kpc")
330
--> 331 photons = photon_model(data_source,
parameters)
332
333 mylog.info("Finished generating photons.")
/home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_models.pyc
in __call__(self, data_source, parameters)
188 raise RuntimeError("Number of
photons generated for this chunk "+
189 "exceeds
photons_per_chunk (%d)! " % self.photons_per_chunk +
--> 190 "Increase
photons_per_chunk!")
191
192 if self.method == "invert_cdf":
RuntimeError: Number of photons generated for this chunk
exceeds photons_per_chunk (1000000000)! Increase
photons_per_chunk!
Ok, so I increase the number of photons per chunk and I get:
2)
---------------------------------------------------------------------------
MemoryError Traceback (most
recent call last)
/home/barbara/hires_final/photon_trial.py in <module>()
90 cosmo = Cosmology()
91
---> 92 photons = PhotonList.from_scratch(Roi, redshift,
A, exp_time,thermal_model,center="center",cosmology=cosmo)
93 photons.write_h5_file("my_photons.h5")
94
/home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_simulator.pyc
in from_scratch(cls, data_source, redshift, area, exp_time,
photon_model, parameters, center, dist, cosmology)
329 parameters["Width"] =
2.*width.in_units("kpc")
330
--> 331 photons = photon_model(data_source,
parameters)
332
333 mylog.info("Finished generating photons.")
/home/barbara/yt-x86_64/src/yt-hg/yt/analysis_modules/photon_simulator/photon_models.pyc
in __call__(self, data_source, parameters)
156
157 number_of_photons = np.zeros(num_cells,
dtype="uint64")
--> 158 energies =
np.zeros(self.photons_per_chunk)
159
160 start_e = 0
MemoryError:
This is just a trial, so I'm using almost a copy of the code that's
on the documentation:
Roi = ds.sphere("center",radius=width_half)
apec_model =
TableApecModel("/home/barbara/yt-x86_64/atomdb_v2.0.2",0.5,10.0,100,thermal_broad=False,apec_vers="2.0.2")
thermal_model = ThermalPhotonModel(apec_model,
X_H=0.75,photons_per_chunk=1000000000, Zmet=0.3,method="invert_cdf")
A = 6000.
exp_time = 100.0e3
redshift = 0.05
cosmo = Cosmology()
photons = PhotonList.from_scratch(Roi, redshift, A,
exp_time,thermal_model,center="center",cosmology=cosmo)
photons.write_h5_file("my_photons.h5")
My snapshot information:
length unit 0.999752515846 cm
time unit 1.0 s
mass unit 1.0 g
velocity unit 1.0 cm/s
Enter total Width of plot in Mpccm: 0.5
smallest cell: 14.6469091426 kpccm/h
minima and maxima 1.04339659244e-29 g/cm**3 1.5187912088e-24 g/cm**3
The simulation is [ 89.40546647 89.40546647 89.40546647] Mpccm
Centering region in:, [ 44.58268195 44.53902694 44.51719943] Mpccm
yt : [INFO ] 2015-09-15 16:49:43,219 Cosmology: h = 0.71,
omega_matter = 0.27, omega_lambda = 0.73
Am I just using wrong parameters to try to simulate the Xray
emission of this simulated cluster?
Thanks beforehand for the help!
--
Barbara L. Ramirez MSc.
Institut für Astro- und Teilchenphysik
University of Innsbruck