[yt-users] PHOX mock observation errors

[Barbara Ramirez]

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150915/f03986f3/attachment.htm>