cosmoglobe / c3pp

A Commander3 post processing tool.

Commander3 postprocessing tool (c3pp)

A python code for processing and plotting commander3 files.

Temperature	Polarization
Thermal dust	Synchrotron
Spinning dust	Free-free

Installation

Install tools by running

pip install git+https://github.com/Cosmoglobe/c3pp.git -U

or

pip install git+https://github.com/Cosmoglobe/c3pp.git --user

and make sure that

PATH=$PATH:~/.local/bin

is in your path.

Usage

The tool uses command line arguments and has many different tools built in (and more coming). In order to get an overview of the available tools and how to access them, simply type

c3pp --help

in your terminal. This will display a list of all available tools. For more information on each specific tool, simply type

c3pp [toolname] --help

and a description will be printed:

[command prompt]$ c3pp
Usage: c3pp [OPTIONS] COMMAND [ARGS]...

Options:
  --help  Show this message and exit.

Commands:
  alm2fits          Converts c3 alms in .h5 file to fits.
  crosspec          Calculates a powerspectrum from polspice.
  dlbin2dat         Outputs binned powerspectra averaged over a range of output samples.
  fits-mean         Calculates the mean over sample range from fits-files.
  fits-stddev       Calculates the standard deviation over sample range from fits-files.
  generate-sky      Generate sky maps from separate input maps.
  gnomplot          Gnomonic view plotting.
  h52fits           Outputs a .h5 map to fits on the form 000001_cmb_amp_n1024.fits
  make-diff-plots   Produces difference maps between output directories.
  mean              Calculates the mean over sample range from .h5 file.
  output-sky-model  Outputs spectrum plots.
  pixreg2trace      Outputs the values of the pixel regions for each sample to a dat file.
  plot              Plots map from .fits or h5 file.
  plotrelease       Plots all release files.
  printdata         Prints the data of a fits file
  printheader       Prints the header of a fits file.
  qu2ang            Calculates polarization angle map from QU signals.
  release           Creates a release file-set on the BeyondPlanck format.
  rmcolumn          Removes columns in fits file
  sigma-l2fits      Converts c3-h5 dataset to fits for c1 BR and GBR estimator analysis.
  specplot          Plots the file output by the Crosspec function.
  stddev            Calculates the stddev over sample range from .h5 file.
  traceplot         Traceplot of samples from .dat.

Contributing

Pull requests are welcome! For feature requests open an issue.

---

Funding

This work has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreements No 776282 (COMPET-4; BeyondPlanck), 772253 (ERC; bits2cosmology) and 819478 (ERC; Cosmoglobe).

---

License

MIT

c3pp's Issues

Add option for WMAP colormap

Just because sometimes it's nice to compare to older figures

Colorbar labels

The support for specifying which values the labels on the color bar takes have been removed. This leads to none or badly places labels on color bar.

Suggestion: Add the feature back to the code if possible.

Polarization plots

First, can we make a plot of Q and U as an option? A la this:

Second, would it be possible to plot headless vectors as an option?

Userdefined color bar units

Would be nice if there was a possibility to use a user defined label for the units of a color bar in a plot. E.g. when plotting the synchrotron amplitude at reference frequency 408 MHz, the intuitive unit would be K_RJ (not uK) as the values typically range between 10 and 300 K_RJ.

Minor ticks so it is more obvious when things are log/linear scale

Remake colorbar labels if using logscale

If using log-scale on user defined limits (e.g. "-min 1 -max 4") the color bar returns the log-limits and not 10 to the power of the limits in the labels of the color bar. It would be nice if the labels of the limits looks like 10^1 and 10^4 (not just 1 and 4), especially when the units of the color bar is not the logarithm.

Add possibility to make GIFs

Some kind of functionality where you could give a deterministic filename-scheme and make GIFs out of them would be cool!

Possibility for map cutouts

Using the gnomview functionality to allow the creation of cut-outs of HEALPix maps would be cool!

make remove_dipole argument less opaque

I need to type -remove_dipole TEXT and I don't know what the TEXT argument is supposed to be. I feel like this should be an optional argument

"plot" not working properly

Running c3pp plot while specifying -min and/or -max gives an error message:

File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/click/core.py", line 764, in call
return self.main(*args, **kwargs)
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/click/core.py", line 717, in main
rv = self.invoke(ctx)
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/click/core.py", line 1137, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/click/core.py", line 956, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/click/core.py", line 555, in invoke
return callback(*args, **kwargs)
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/c3postproc/commands.py", line 214, in plot
Plotter(input, dataset, nside, auto, min, max, minmax, range, colorbar, lmax, fwhm, mask, mfill, sig, remove_dipole, logscale, size, white_background, darkmode, pdf, cmap, title, ltitle, unit, scale, verbose,)
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/c3postproc/plotter.py", line 371, in Plotter
ticklabels = [fmt(i, 1) for i in ticklabels]
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/c3postproc/plotter.py", line 371, in
ticklabels = [fmt(i, 1) for i in ticklabels]
File "/mn/stornext/u3/krisjand/.local/lib/python3.6/site-packages/c3postproc/plotter.py", line 786, in fmt
TypeError: bad operand type for abs(): 'str'

Apparently, the min/max values are stored as string (or something).

Furthermore, running "plot" with the -range option does nothing (the range is not applied).