[yt-users] PHOX mock observation errors

[John ZuHone]

Hi Barbara,

Sorry to hear you’re having trouble at the moment. I have a few other questions. 

1) Can you post the entire script? Use http://paste.yt-project.org and send a link. 
2) Can you give me a sense of what the mass and temperature of your cluster is? It sounds like it is fairly large. And the radius? 
3) What code did the simulation come from?

After knowing these, I think we can get somewhere. I have some ideas. 

Best,

John 

> On Sep 15, 2015, at 10:59 AM, Barbara Ramirez <Barbara.Ramirez-Mosquera at uibk.ac.at> wrote:
> 
> 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
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

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