#!/bin/env python

import h5py
import numpy as np
import healpy as hp

import read_lightcone as rl

plt.subplot(2, 2, 1)
totalmass_map = rl.read_map(
    "./lightcones/", "lightcone0", shell_nr=0, map_name="TotalMass"
)
hp.mollview(totalmass_map + 1, norm="log", title="Projected mass", hold=True)

plt.subplot(2, 2, 2)
gasmass_map = rl.read_map(
    "./lightcones/", "lightcone0", shell_nr=0, map_name="SmoothedGasMass"
)
hp.mollview(gasmass_map + 1, norm="log", title="Gas mass", hold=True)

plt.subplot(2, 2, 3)
stellarmass_map = rl.read_map(
    "./lightcones/", "lightcone0", shell_nr=0, map_name="StellarMass"
)
hp.mollview(stellarmass_map + 1, norm="log", title="Stellar mass", hold=True)

plt.subplot(2, 2, 4)
xray_map = rl.read_map(
    "./lightcones/", "lightcone0", shell_nr=0, map_name="XrayROSATIntrinsicPhotons"
)
hp.mollview(xray_map + 1e50, norm="log", title="ROSAT photons", hold=True)

plt.suptitle("SmallCosmoVolume lightcone")