/*
 * gui_peaksearch.c
 *
 * Peak search parts of GUI
 *
 * Copyright © 2020-2021 Deutsches Elektronen-Synchrotron DESY,
 *                       a research centre of the Helmholtz Association.
 *
 * Authors:
 *   2020-2021 Thomas White <taw@physics.org>
 *
 * This file is part of CrystFEL.
 *
 * CrystFEL is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * CrystFEL is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with CrystFEL.  If not, see <http://www.gnu.org/licenses/>.
 *
 */


#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include <stdarg.h>
#include <stdio.h>
#include <getopt.h>
#include <string.h>
#include <gtk/gtk.h>
#include <gdk/gdkkeysyms-compat.h>
#include <assert.h>

#include <datatemplate.h>
#include <peaks.h>
#include <peakfinder8.h>

#include "crystfelimageview.h"
#include "gui_project.h"
#include "crystfel_gui.h"


void update_peaks(struct crystfelproject *proj)
{
	if ( proj->n_frames == 0 ) return;
	if ( proj->cur_image == NULL ) return;

	crystfel_image_view_set_peak_box_size(CRYSTFEL_IMAGE_VIEW(proj->imageview),
	                                      proj->peak_search_params.pk_inn);

	image_feature_list_free(proj->cur_image->features);
	proj->cur_image->features = NULL;

	switch ( proj->peak_search_params.method ) {

		ImageFeatureList *peaks;

		case PEAK_ZAEF:
		proj->cur_image->features = search_peaks(proj->cur_image,
		                                         proj->peak_search_params.threshold,
		                                         proj->peak_search_params.min_sq_gradient,
		                                         proj->peak_search_params.min_snr,
		                                         proj->peak_search_params.pk_inn,
		                                         proj->peak_search_params.pk_mid,
		                                         proj->peak_search_params.pk_out,
		                                         1);
		break;

		case PEAK_PEAKFINDER8:
		proj->cur_image->features = peakfinder8(proj->cur_image, 2048,
		                                        proj->peak_search_params.threshold,
		                                        proj->peak_search_params.min_snr,
		                                        proj->peak_search_params.min_pix_count,
		                                        proj->peak_search_params.max_pix_count,
		                                        proj->peak_search_params.local_bg_radius,
		                                        proj->peak_search_params.min_res,
		                                        proj->peak_search_params.max_res,
		                                        1, 0, NULL);
		break;

		case PEAK_PEAKFINDER9:
		proj->cur_image->features = search_peaks_peakfinder9(proj->cur_image,
		                                                     proj->peak_search_params.min_snr_biggest_pix,
		                                                     proj->peak_search_params.min_snr_peak_pix,
		                                                     proj->peak_search_params.min_snr,
		                                                     proj->peak_search_params.min_sig,
		                                                     proj->peak_search_params.min_peak_over_neighbour,
		                                                     proj->peak_search_params.local_bg_radius);
		break;

		case PEAK_HDF5:
		case PEAK_CXI:
		peaks = image_read_peaks(proj->dtempl,
		                         proj->cur_image->filename,
		                         proj->cur_image->ev,
		                         proj->peak_search_params.half_pixel_shift);
		if ( proj->peak_search_params.revalidate ) {
			proj->cur_image->features = validate_peaks(proj->cur_image, peaks,
			                                           proj->peak_search_params.min_snr,
			                                           proj->peak_search_params.pk_inn,
			                                           proj->peak_search_params.pk_mid,
			                                           proj->peak_search_params.pk_out,
			                                           1, 0);
			image_feature_list_free(peaks);
		}
		break;

		default:
		ERROR("This peak detection method not implemented!\n");
		break;

	}
}


struct param_callback_vals
{
	float *pfval;
	int *pival;
	char *str_val;
	struct crystfelproject *proj;
};


static void check_param_callback(GtkWidget *checkbox,
                                 struct param_callback_vals *vals)
{
	int val = gtk_toggle_button_get_active(GTK_TOGGLE_BUTTON(checkbox));
	*(vals->pival) = val;
	update_imageview(vals->proj);
}


static void int_param_callback(GtkWidget *entry,
                               struct param_callback_vals *cbvals)
{
	const char *text;
	int val;

	text = gtk_entry_get_text(GTK_ENTRY(entry));
	if (sscanf(text, "%i", &val) != 1) {
		ERROR("Invalid value\n");
		return;
	}

	/* Only update peaks if value has changed */
	if ( strcmp(cbvals->str_val, text) != 0 ) {
		*(cbvals->pival) = val;
		cbvals->proj->unsaved = 1;
		update_imageview(cbvals->proj);
		free(cbvals->str_val);
		cbvals->str_val = strdup(text);
	}
}


static gboolean int_param_focus_callback(GtkWidget *entry,
                                         GdkEvent *event,
                                         struct param_callback_vals *cbvals)
{
	int_param_callback(entry, cbvals);
	return FALSE;
}


static void float_param_callback(GtkWidget *entry,
                                 struct param_callback_vals *cbvals)
{
	const char *text;
	float val;

	text = gtk_entry_get_text(GTK_ENTRY(entry));
	if (sscanf(text, "%f", &val) != 1) {
		ERROR("Invalid value\n");
		return;
	}

	/* Only update peaks if value has changed */
	if ( strcmp(cbvals->str_val, text) != 0 ) {
		*(cbvals->pfval) = val;
		cbvals->proj->unsaved = 1;
		update_imageview(cbvals->proj);
		free(cbvals->str_val);
		cbvals->str_val = strdup(text);
	}
}


static gboolean float_param_focus_callback(GtkWidget *entry,
                                           GdkEvent *event,
                                           struct param_callback_vals *cbvals)
{
	float_param_callback(entry, cbvals);
	return FALSE;
}


static void free_callback_params(gpointer cbvals,
                                 GClosure *closure)
{
	free(cbvals);
}


static void add_int_param(GtkWidget *params_box, const char *labeltext,
                          int *pval, struct crystfelproject *proj,
                          const char *tooltip)
{
	GtkWidget *hbox;
	GtkWidget *label;
	GtkWidget *entry;
	struct param_callback_vals *cbvals;
	char tmp[64];

	hbox = gtk_hbox_new(FALSE, 0.0);
	gtk_box_pack_start(GTK_BOX(params_box),
	                   GTK_WIDGET(hbox), FALSE, FALSE, 8.0);
	label = gtk_label_new(labeltext);
	gtk_misc_set_alignment(GTK_MISC(label), 1.0, 0.5);
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(label),
	                   FALSE, FALSE, 2.0);
	entry = gtk_entry_new();
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(entry),
	                   TRUE, TRUE, 2.0);
	snprintf(tmp, 63, "%i", *pval);
	gtk_entry_set_text(GTK_ENTRY(entry), tmp);

	cbvals = malloc(sizeof(struct param_callback_vals));
	if ( cbvals != NULL ) {
		cbvals->proj = proj;
		cbvals->pival = pval;
		cbvals->str_val = strdup(tmp);
		g_signal_connect_data(G_OBJECT(entry),
		                      "activate",
		                      G_CALLBACK(int_param_callback),
		                      cbvals, free_callback_params, 0);
		g_signal_connect(G_OBJECT(entry),
		                 "focus-out-event",
		                 G_CALLBACK(int_param_focus_callback),
		                 cbvals);
	} else {
		ERROR("Failed to connect parameter callback\n");
	}

	gtk_widget_set_tooltip_text(hbox, tooltip);
}


static void add_float_param(GtkWidget *params_box, const char *labeltext,
                            float *pval, struct crystfelproject *proj,
                            const char *tooltip)
{
	GtkWidget *hbox;
	GtkWidget *label;
	GtkWidget *entry;
	struct param_callback_vals *cbvals;
	char tmp[64];

	hbox = gtk_hbox_new(FALSE, 0.0);
	gtk_box_pack_start(GTK_BOX(params_box),
	                   GTK_WIDGET(hbox), FALSE, FALSE, 8.0);
	label = gtk_label_new(labeltext);
	gtk_misc_set_alignment(GTK_MISC(label), 1.0, 0.5);
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(label),
	                   FALSE, FALSE, 2.0);
	entry = gtk_entry_new();
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(entry),
	                   TRUE, TRUE, 2.0);
	snprintf(tmp, 63, "%.2f", *pval);
	gtk_entry_set_text(GTK_ENTRY(entry), tmp);

	cbvals = malloc(sizeof(struct param_callback_vals));
	if ( cbvals != NULL ) {
		cbvals->proj = proj;
		cbvals->pfval = pval;
		cbvals->str_val = strdup(tmp);
		g_signal_connect_data(G_OBJECT(entry),
		                      "activate",
		                      G_CALLBACK(float_param_callback),
		                      cbvals, free_callback_params, 0);
		g_signal_connect(G_OBJECT(entry),
		                 "focus-out-event",
		                 G_CALLBACK(float_param_focus_callback),
		                 cbvals);

	} else {
		ERROR("Failed to connect parameter callback\n");
	}

	gtk_widget_set_tooltip_text(hbox, tooltip);
}


static void add_check_param(GtkWidget *params_box, const char *labeltext,
                            int *pval, struct crystfelproject *proj,
                            const char *tooltip)
{
	GtkWidget *checkbox;
	struct param_callback_vals *cbvals;

	checkbox = gtk_check_button_new_with_label(labeltext);
	gtk_box_pack_start(GTK_BOX(params_box),
	                   GTK_WIDGET(checkbox), FALSE, FALSE, 8.0);
	gtk_toggle_button_set_active(GTK_TOGGLE_BUTTON(checkbox),
	                             *pval);

	cbvals = malloc(sizeof(struct param_callback_vals));
	if ( cbvals != NULL ) {
		cbvals->proj = proj;
		cbvals->pival = pval;
		g_signal_connect_data(G_OBJECT(checkbox), "toggled",
		                      G_CALLBACK(check_param_callback),
		                      cbvals, free_callback_params, 0);
	} else {
		ERROR("Failed to connect parameter callback\n");
	}

	gtk_widget_set_tooltip_text(checkbox, tooltip);
}


static void add_radii(GtkWidget *params_box,
                      struct crystfelproject *proj)
{
	add_float_param(params_box, "Peak radius (inner):",
	                &proj->peak_search_params.pk_inn, proj,
	                "--peak-radius");
	add_float_param(params_box, "Peak radius (middle):",
	                &proj->peak_search_params.pk_mid, proj,
	                "--peak-radius");
	add_float_param(params_box, "Peak radius (outer):",
	                &proj->peak_search_params.pk_out, proj,
	                "--peak-radius");
}


static void peaksearch_algo_changed(GtkWidget *combo,
                                    struct crystfelproject *proj)
{
	const char *algo_id;

	if ( proj->peak_params != NULL ) {
		gtk_container_remove(GTK_CONTAINER(proj->peak_vbox),
		                     proj->peak_params);
		proj->peak_params = NULL;
	}

	proj->peak_params = gtk_vbox_new(FALSE, 0.0);
	gtk_box_pack_start(GTK_BOX(proj->peak_vbox),
	                   GTK_WIDGET(proj->peak_params),
	                   FALSE, FALSE, 8.0);

	algo_id = gtk_combo_box_get_active_id(GTK_COMBO_BOX(combo));

	if ( strcmp(algo_id, "zaef") == 0 ) {

		proj->peak_search_params.method = PEAK_ZAEF;

		add_float_param(proj->peak_params, "Threshold:",
		                &proj->peak_search_params.threshold, proj,
		                "--threshold");
		add_float_param(proj->peak_params, "Minimum squared gradient:",
		                &proj->peak_search_params.min_sq_gradient, proj,
		                "--min-squared-gradient (--min-gradient)");
		add_float_param(proj->peak_params, "Minimum signal/noise ratio:",
		                &proj->peak_search_params.min_snr, proj,
		                "--min-snr");
		add_radii(proj->peak_params, proj);

	} else if ( strcmp(algo_id, "peakfinder8") == 0 ) {

		proj->peak_search_params.method = PEAK_PEAKFINDER8;

		add_float_param(proj->peak_params, "Threshold:",
		                &proj->peak_search_params.threshold, proj,
		                "--threshold");
		add_float_param(proj->peak_params, "Minimum signal/noise ratio:",
		                &proj->peak_search_params.min_snr, proj,
		                "--min-snr");
		add_int_param(proj->peak_params, "Minimum number of pixels:",
		              &proj->peak_search_params.min_pix_count, proj,
		              "--min-pix-count");
		add_int_param(proj->peak_params, "Maximum number of pixels:",
		              &proj->peak_search_params.max_pix_count, proj,
		              "--max-pix-count");
		add_int_param(proj->peak_params, "Local background radius:",
		              &proj->peak_search_params.local_bg_radius, proj,
		              "--local-bg-radius");
		add_int_param(proj->peak_params, "Minimum resolution (pixels):",
		              &proj->peak_search_params.min_res, proj,
		              "--min-res");
		add_int_param(proj->peak_params, "Maximum resolution (pixels):",
		              &proj->peak_search_params.max_res, proj,
		              "--max-res");

	} else if ( strcmp(algo_id, "peakfinder9") == 0 ) {

		proj->peak_search_params.method = PEAK_PEAKFINDER9;

		add_float_param(proj->peak_params, "Minimum SNR of brightest pixel in peak:",
		                &proj->peak_search_params.min_snr_biggest_pix, proj,
		                "--min-snr-biggest-pix");
		add_float_param(proj->peak_params, "Minimum SNR of peak pixel:",
		                &proj->peak_search_params.min_snr_peak_pix, proj,
		                "--min-snr-peak-pix");
		add_float_param(proj->peak_params, "Minimum signal/noise ratio:",
		                &proj->peak_search_params.min_snr, proj,
		                "-min-snr");
		add_float_param(proj->peak_params, "Minimum background standard deviation:",
		                &proj->peak_search_params.min_sig, proj,
		                "--min-sig");
		add_float_param(proj->peak_params, "Brightest pixel cutoff (just for speed):",
		                &proj->peak_search_params.min_peak_over_neighbour, proj,
		                "--min-peak-over-neighbour");
		add_int_param(proj->peak_params, "Local background radius:",
		              &proj->peak_search_params.local_bg_radius, proj,
		              "--local-bg-radius");

	} else if ( strcmp(algo_id, "hdf5") == 0 ) {

		proj->peak_search_params.method = PEAK_HDF5;

		add_check_param(proj->peak_params, "Half pixel shift",
		                &proj->peak_search_params.half_pixel_shift,
		                proj, "--no-half-pixel-shift");
		add_check_param(proj->peak_params, "Check peaks first",
		                &proj->peak_search_params.revalidate,
		                proj, "--no-revalidate");
		add_radii(proj->peak_params, proj);

	} else if ( strcmp(algo_id, "cxi") == 0 ) {
		ERROR("algo_id should be hdf5, not cxi\n");

	} else {
		ERROR("Unrecognised peak search '%s'\n", algo_id);
	}

	gtk_widget_show_all(proj->peak_vbox);
	update_peaks(proj);
	update_imageview(proj);
	proj->unsaved = 1;
}


static void peaksearch_response_sig(GtkWidget *dialog, gint resp,
                                    struct crystfelproject *proj)
{
	if ( (resp == GTK_RESPONSE_DELETE_EVENT)
	  || (resp == GTK_RESPONSE_CANCEL) )
	{
		proj->peak_search_params = proj->original_params;
	}

	update_imageview(proj);
	gtk_widget_destroy(dialog);
	proj->peak_vbox = NULL;
	proj->peak_params = NULL;
}


gint peaksearch_sig(GtkWidget *widget, struct crystfelproject *proj)
{
	GtkWidget *dialog;
	GtkWidget *content_area;
	GtkWidget *vbox;
	GtkWidget *hbox;
	GtkWidget *label;
	GtkWidget *combo;

	if ( proj->peak_params != NULL ) return FALSE;

	force_peaks_on(proj);
	gtk_combo_box_set_active(GTK_COMBO_BOX(proj->results_combo), 0);

	/* Take a copy of the original parameters, for reverting */
	proj->original_params = proj->peak_search_params;

	dialog = gtk_dialog_new_with_buttons("Peak search",
	                                     GTK_WINDOW(proj->window),
	                                     GTK_DIALOG_DESTROY_WITH_PARENT,
	                                     "Discard changes", GTK_RESPONSE_CANCEL,
	                                     "Confirm", GTK_RESPONSE_ACCEPT,
	                                     NULL);

	g_signal_connect(G_OBJECT(dialog), "response",
	                 G_CALLBACK(peaksearch_response_sig), proj);

	vbox = gtk_vbox_new(FALSE, 0.0);
	content_area = gtk_dialog_get_content_area(GTK_DIALOG(dialog));
	gtk_container_add(GTK_CONTAINER(content_area), vbox);
	gtk_container_set_border_width(GTK_CONTAINER(content_area), 8);
	proj->peak_vbox = vbox;

	hbox = gtk_hbox_new(FALSE, 0.0);
	gtk_box_pack_start(GTK_BOX(vbox), GTK_WIDGET(hbox), FALSE, FALSE, 8.0);
	label = gtk_label_new("Peak search algorithm");
	gtk_misc_set_alignment(GTK_MISC(label), 1.0, 0.5);
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(label), FALSE, FALSE, 2.0);
	combo = gtk_combo_box_text_new();
	gtk_box_pack_start(GTK_BOX(hbox), GTK_WIDGET(combo), TRUE, TRUE, 2.0);
	gtk_combo_box_text_append(GTK_COMBO_BOX_TEXT(combo), "zaef",
	                "Zaefferer gradient search (zaef)");
	gtk_combo_box_text_append(GTK_COMBO_BOX_TEXT(combo), "hdf5",
	                "Use the peak lists in the data files (hdf5/cxi)");
	gtk_combo_box_text_append(GTK_COMBO_BOX_TEXT(combo), "peakfinder8",
	                "Radial background estimation (peakfinder8)");

	if ( crystfel_has_peakfinder9() ) {
		gtk_combo_box_text_append(GTK_COMBO_BOX_TEXT(combo), "peakfinder9",
		                          "Local background estimation (peakfinder9)");
	}

	g_signal_connect(G_OBJECT(combo), "changed",
	                 G_CALLBACK(peaksearch_algo_changed), proj);
	gtk_combo_box_set_active_id(GTK_COMBO_BOX(combo),
	                            str_peaksearch(proj->peak_search_params.method));
	proj->type_combo = combo;

	gtk_widget_show_all(dialog);

	return FALSE;
}