comparison_plots.py
add systematics