Index: trunk/share/debug/Makefile_full
===================================================================
--- trunk/share/debug/Makefile_full	(revision 8303)
+++ trunk/share/debug/Makefile_full	(revision 8304)
@@ -1,504 +1,545 @@
-FC=ifort2018beta
-FCFLAGS= 
+FC=pgfortran_2019
+FCFLAGS=-Mbackslash
 CC=gcc
 CCFLAGS=
 
 MODELS = \
 	SM.mdl \
 	SM_hadrons.mdl \
 	Test.mdl
 
 CC_SRC = \
 	sprintf_interface.c \
 	signal_interface.c
 
 F77_SRC = \
 	pythia.f \
 	pythia_pdf.f \
 	pythia6_up.f \
 	toppik.f \
 	toppik_axial.f
 
 FC0_SRC = 
 
 FC_SRC = \
 	format_defs.f90 \
 	io_units.f90 \
 	kinds.f90 \
 	constants.f90 \
 	iso_varying_string.f90 \
 	unit_tests.f90 \
 	system_dependencies.f90 \
 	string_utils.f90 \
 	system_defs.f90 \
+	debug_master.f90 \
 	diagnostics.f90 \
 	sorting.f90 \
 	physics_defs.f90 \
 	pdg_arrays.f90 \
 	bytes.f90 \
 	hashes.f90 \
 	md5.f90 \
 	model_data.f90 \
 	numeric_utils.f90 \
 	auto_components.f90 \
 	var_base.f90 \
 	model_testbed.f90 \
 	auto_components_uti.f90 \
 	auto_components_ut.f90 \
 	os_interface.f90 \
 	c_particles.f90 \
 	format_utils.f90 \
 	lorentz.f90 \
 	colors.f90 \
 	flavors.f90 \
 	helicities.f90 \
 	quantum_numbers.f90 \
 	state_matrices.f90 \
 	interactions.f90 \
 	CppStringsWrap_dummy.f90 \
 	FastjetWrap_dummy.f90 \
 	cpp_strings.f90 \
 	fastjet.f90 \
 	jets.f90 \
 	subevents.f90 \
 	su_algebra.f90 \
 	bloch_vectors.f90 \
 	polarizations.f90 \
 	particles.f90 \
 	event_base.f90 \
 	eio_data.f90 \
 	eio_base.f90 \
 	eio_base_uti.f90 \
 	eio_base_ut.f90 \
 	variables.f90 \
 	rng_base.f90 \
 	tao_random_numbers.f90 \
 	rng_tao.f90 \
 	rng_stream.f90 \
 	rng_base_uti.f90 \
 	rng_base_ut.f90 \
 	dispatch_rng.f90 \
 	dispatch_rng_uti.f90 \
 	dispatch_rng_ut.f90 \
 	beam_structures.f90 \
 	evaluators.f90 \
 	beams.f90 \
 	sm_physics.f90 \
 	file_registries.f90 \
 	sf_aux.f90 \
 	sf_mappings.f90 \
 	sf_base.f90 \
+	electron_pdfs.f90 \
 	sf_isr.f90 \
 	sf_epa.f90 \
 	sf_ewa.f90 \
 	sf_escan.f90 \
 	sf_gaussian.f90 \
 	sf_beam_events.f90 \
 	circe1.f90 \
 	sf_circe1.f90 \
 	circe2.f90 \
 	selectors.f90 \
 	sf_circe2.f90 \
 	sm_qcd.f90 \
 	mrst2004qed.f90 \
 	cteq6pdf.f90 \
 	mstwpdf.f90 \
 	ct10pdf.f90 \
 	CJpdf.f90 \
 	ct14pdf.f90 \
 	pdf_builtin.f90 \
 	LHAPDFWrap_dummy.f90 \
 	lhapdf5_full_dummy.f90 \
 	lhapdf5_has_photon_dummy.f90 \
 	lhapdf.f90 \
 	hoppet_dummy.f90 \
 	hoppet_interface.f90 \
 	sf_pdf_builtin.f90 \
 	sf_lhapdf.f90 \
 	dispatch_beams.f90 \
 	process_constants.f90 \
 	prclib_interfaces.f90 \
 	prc_core_def.f90 \
 	particle_specifiers.f90 \
 	process_libraries.f90 \
 	prc_test.f90 \
 	prc_core.f90 \
 	prc_test_core.f90 \
 	prc_omega.f90 \
 	phs_base.f90 \
 	ifiles.f90 \
 	lexers.f90 \
 	syntax_rules.f90 \
 	parser.f90 \
 	expr_base.f90 \
 	formats.f90 \
 	analysis.f90 \
 	user_code_interface.f90 \
 	observables.f90 \
 	eval_trees.f90 \
 	interpolation.f90 \
 	nr_tools.f90 \
 	ttv_formfactors.f90 \
 	ttv_formfactors_uti.f90 \
 	ttv_formfactors_ut.f90 \
 	models.f90 \
 	prclib_stacks.f90 \
 	user_files.f90 \
 	cputime.f90 \
 	mci_base.f90 \
 	integration_results.f90 \
-	pcm_base.f90 \
-	nlo_data.f90 \
+	integration_results_uti.f90 \
+	integration_results_ut.f90 \
 	mappings.f90 \
 	permutations.f90 \
 	resonances.f90 \
 	phs_trees.f90 \
 	phs_forests.f90 \
+	prc_external.f90 \
+	blha_config.f90 \
+	blha_olp_interfaces.f90 \
+	prc_openloops.f90 \
+	prc_threshold.f90 \
+	process_config.f90 \
+	process_counter.f90 \
+	process_mci.f90 \
+	pcm_base.f90 \
+	nlo_data.f90 \
 	cascades.f90 \
+	cascades2_lexer.f90 \
+	cascades2_lexer_uti.f90 \
+	cascades2_lexer_ut.f90 \
+	cascades2.f90 \
+	cascades2_uti.f90 \
+	cascades2_ut.f90 \
+	phs_none.f90 \
+	phs_rambo.f90 \
 	phs_wood.f90 \
 	phs_fks.f90 \
 	phs_single.f90 \
-	prc_user_defined.f90 \
 	fks_regions.f90 \
-	blha_config.f90 \
-	blha_olp_interfaces.f90 \
 	virtual.f90 \
 	pdf.f90 \
 	real_subtraction.f90 \
 	dglap_remnant.f90 \
-	nlo_color_data.f90 \
+	dispatch_fks.f90 \
+	dispatch_phase_space.f90 \
 	pcm.f90 \
-	process_counter.f90 \
-	core_manager.f90 \
-	prc_openloops.f90 \
 	recola_wrapper_dummy.f90 \
 	prc_recola.f90 \
-	prc_threshold.f90 \
-	process_config.f90 \
-	process_mci.f90 \
 	subevt_expr.f90 \
 	parton_states.f90 \
+	prc_template_me.f90 \
 	process.f90 \
 	process_stacks.f90 \
 	iterations.f90 \
 	rt_data.f90 \
-	dispatch_phase_space.f90 \
-	prc_template_me.f90 \
 	file_utils.f90 \
 	prc_gosam.f90 \
 	dispatch_me_methods.f90 \
 	sf_base_uti.f90 \
 	sf_base_ut.f90 \
 	dispatch_uti.f90 \
 	dispatch_ut.f90 \
 	formats_uti.f90 \
 	formats_ut.f90 \
 	md5_uti.f90 \
 	md5_ut.f90 \
 	os_interface_uti.f90 \
 	os_interface_ut.f90 \
 	sorting_uti.f90 \
 	sorting_ut.f90 \
 	grids.f90 \
 	grids_uti.f90 \
 	grids_ut.f90 \
 	solver.f90 \
 	solver_uti.f90 \
 	solver_ut.f90 \
 	cputime_uti.f90 \
 	cputime_ut.f90 \
 	sm_qcd_uti.f90 \
 	sm_qcd_ut.f90 \
 	sm_physics_uti.f90 \
 	sm_physics_ut.f90 \
 	lexers_uti.f90 \
 	lexers_ut.f90 \
 	parser_uti.f90 \
 	parser_ut.f90 \
 	xml.f90 \
 	xml_uti.f90 \
 	xml_ut.f90 \
 	colors_uti.f90 \
 	colors_ut.f90 \
 	state_matrices_uti.f90 \
 	state_matrices_ut.f90 \
 	analysis_uti.f90 \
 	analysis_ut.f90 \
 	particles_uti.f90 \
 	particles_ut.f90 \
 	radiation_generator.f90 \
 	radiation_generator_uti.f90 \
 	radiation_generator_ut.f90 \
 	blha_uti.f90 \
 	blha_ut.f90 \
 	evaluators_uti.f90 \
 	evaluators_ut.f90 \
 	models_uti.f90 \
 	models_ut.f90 \
 	eval_trees_uti.f90 \
 	eval_trees_ut.f90 \
 	resonances_uti.f90 \
 	resonances_ut.f90 \
 	phs_trees_uti.f90 \
 	phs_trees_ut.f90 \
 	phs_forests_uti.f90 \
 	phs_forests_ut.f90 \
 	beams_uti.f90 \
 	beams_ut.f90 \
 	su_algebra_uti.f90 \
 	su_algebra_ut.f90 \
 	bloch_vectors_uti.f90 \
 	bloch_vectors_ut.f90 \
 	polarizations_uti.f90 \
 	polarizations_ut.f90 \
 	sf_aux_uti.f90 \
 	sf_aux_ut.f90 \
 	sf_mappings_uti.f90 \
 	sf_mappings_ut.f90 \
 	sf_pdf_builtin_uti.f90 \
 	sf_pdf_builtin_ut.f90 \
 	sf_lhapdf_uti.f90 \
 	sf_lhapdf_ut.f90 \
 	sf_isr_uti.f90 \
 	sf_isr_ut.f90 \
 	sf_epa_uti.f90 \
 	sf_epa_ut.f90 \
 	sf_ewa_uti.f90 \
 	sf_ewa_ut.f90 \
 	sf_circe1_uti.f90 \
 	sf_circe1_ut.f90 \
 	sf_circe2_uti.f90 \
 	sf_circe2_ut.f90 \
 	sf_gaussian_uti.f90 \
 	sf_gaussian_ut.f90 \
 	sf_beam_events_uti.f90 \
 	sf_beam_events_ut.f90 \
 	sf_escan_uti.f90 \
 	sf_escan_ut.f90 \
 	phs_base_uti.f90 \
 	phs_base_ut.f90 \
+	phs_none_uti.f90 \
+	phs_none_ut.f90 \
 	phs_single_uti.f90 \
 	phs_single_ut.f90 \
+	phs_rambo_uti.f90 \
+	phs_rambo_ut.f90 \
 	phs_wood_uti.f90 \
 	phs_wood_ut.f90 \
 	phs_fks_uti.f90 \
 	phs_fks_ut.f90 \
 	fks_regions_uti.f90 \
 	fks_regions_ut.f90 \
 	mci_midpoint.f90 \
 	mci_base_uti.f90 \
 	mci_base_ut.f90 \
 	mci_midpoint_uti.f90 \
 	mci_midpoint_ut.f90 \
 	kinematics.f90 \
 	instances.f90 \
+	mci_none.f90 \
+	mci_none_uti.f90 \
+	mci_none_ut.f90 \
 	processes_uti.f90 \
 	processes_ut.f90 \
 	process_stacks_uti.f90 \
 	process_stacks_ut.f90 \
 	prc_recola_uti.f90 \
 	prc_recola_ut.f90 \
-	nlo_color_data_uti.f90 \
-	nlo_color_data_ut.f90 \
 	rng_tao_uti.f90 \
 	rng_tao_ut.f90 \
 	rng_stream_uti.f90 \
 	rng_stream_ut.f90 \
 	selectors_uti.f90 \
 	selectors_ut.f90 \
 	vegas.f90 \
 	vegas_uti.f90 \
 	vegas_ut.f90 \
 	vamp2.f90 \
 	vamp2_uti.f90 \
 	vamp2_ut.f90 \
 	exceptions.f90 \
 	vamp_stat.f90 \
 	utils.f90 \
 	divisions.f90 \
 	linalg.f90 \
 	vamp.f90 \
 	mci_vamp.f90 \
 	mci_vamp_uti.f90 \
 	mci_vamp_ut.f90 \
 	mci_vamp2.f90 \
 	mci_vamp2_uti.f90 \
 	mci_vamp2_ut.f90 \
 	prclib_interfaces_uti.f90 \
 	prclib_interfaces_ut.f90 \
 	particle_specifiers_uti.f90 \
 	particle_specifiers_ut.f90 \
 	process_libraries_uti.f90 \
 	process_libraries_ut.f90 \
 	prclib_stacks_uti.f90 \
 	prclib_stacks_ut.f90 \
 	slha_interface.f90 \
 	slha_interface_uti.f90 \
 	slha_interface_ut.f90 \
 	cascades_uti.f90 \
 	cascades_ut.f90 \
 	prc_test_uti.f90 \
 	prc_test_ut.f90 \
 	prc_template_me_uti.f90 \
 	prc_template_me_ut.f90 \
 	prc_omega_uti.f90 \
 	prc_omega_ut.f90 \
 	event_transforms.f90 \
 	event_transforms_uti.f90 \
 	event_transforms_ut.f90 \
 	hep_common.f90 \
 	hepev4_aux.f90 \
 	tauola_dummy.f90 \
 	tauola_interface.f90 \
 	shower_base.f90 \
 	shower_partons.f90 \
 	muli_base.f90 \
 	muli_momentum.f90 \
 	muli_trapezium.f90 \
 	muli_interactions.f90 \
 	muli_cuba.f90 \
 	muli_fibonacci_tree.f90 \
 	muli_aq.f90 \
 	muli_dsigma.f90 \
 	muli_mcint.f90 \
 	muli_remnant.f90 \
 	muli.f90 \
 	matching_base.f90 \
 	powheg_matching.f90 \
 	shower_core.f90 \
 	shower_base_uti.f90 \
 	shower_base_ut.f90 \
 	shower.f90 \
 	shower_uti.f90 \
 	shower_ut.f90 \
 	shower_pythia6.f90 \
+	whizard_lha.f90 \
+	whizard_lha_uti.f90 \
+	whizard_lha_ut.f90 \
+	LHAWhizard_dummy.f90 \
+	Pythia8Wrap_dummy.f90 \
+	pythia8.f90 \
+	pythia8_uti.f90 \
+	pythia8_ut.f90 \
+	shower_pythia8.f90 \
 	hadrons.f90 \
 	ktclus.f90 \
 	mlm_matching.f90 \
 	ckkw_matching.f90 \
 	jets_uti.f90 \
 	jets_ut.f90 \
 	pdg_arrays_uti.f90 \
 	pdg_arrays_ut.f90 \
 	interactions_uti.f90 \
 	interactions_ut.f90 \
 	decays.f90 \
 	decays_uti.f90 \
 	decays_ut.f90 \
 	evt_nlo.f90 \
 	events.f90 \
 	events_uti.f90 \
 	events_ut.f90 \
 	HepMCWrap_dummy.f90 \
 	hepmc_interface.f90 \
 	hepmc_interface_uti.f90 \
 	hepmc_interface_ut.f90 \
 	LCIOWrap_dummy.f90 \
 	lcio_interface.f90 \
 	lcio_interface_uti.f90 \
 	lcio_interface_ut.f90 \
 	hep_events.f90 \
 	hep_events_uti.f90 \
 	hep_events_ut.f90 \
 	expr_tests_uti.f90 \
 	expr_tests_ut.f90 \
 	parton_states_uti.f90 \
 	parton_states_ut.f90 \
 	eio_data_uti.f90 \
 	eio_data_ut.f90 \
 	eio_raw.f90 \
 	eio_raw_uti.f90 \
 	eio_raw_ut.f90 \
 	eio_checkpoints.f90 \
 	eio_checkpoints_uti.f90 \
 	eio_checkpoints_ut.f90 \
 	eio_lhef.f90 \
 	eio_lhef_uti.f90 \
 	eio_lhef_ut.f90 \
 	eio_hepmc.f90 \
 	eio_hepmc_uti.f90 \
 	eio_hepmc_ut.f90 \
 	eio_lcio.f90 \
 	eio_lcio_uti.f90 \
 	eio_lcio_ut.f90 \
 	stdhep_dummy.f90 \
 	xdr_wo_stdhep.f90 \
 	eio_stdhep.f90 \
 	eio_stdhep_uti.f90 \
 	eio_stdhep_ut.f90 \
 	eio_ascii.f90 \
 	eio_ascii_uti.f90 \
 	eio_ascii_ut.f90 \
 	eio_weights.f90 \
 	eio_weights_uti.f90 \
 	eio_weights_ut.f90 \
 	eio_dump.f90 \
 	eio_dump_uti.f90 \
 	eio_dump_ut.f90 \
 	eio_callback.f90 \
+	real_subtraction_uti.f90 \
+	real_subtraction_ut.f90 \
 	iterations_uti.f90 \
 	iterations_ut.f90 \
 	rt_data_uti.f90 \
 	rt_data_ut.f90 \
 	dispatch_mci.f90 \
 	dispatch_mci_uti.f90 \
 	dispatch_mci_ut.f90 \
 	dispatch_phs_uti.f90 \
 	dispatch_phs_ut.f90 \
+	resonance_insertion.f90 \
+	resonance_insertion_uti.f90 \
+	resonance_insertion_ut.f90 \
+	recoil_kinematics.f90 \
+	recoil_kinematics_uti.f90 \
+	recoil_kinematics_ut.f90 \
+	isr_epa_handler.f90 \
+	isr_epa_handler_uti.f90 \
+	isr_epa_handler_ut.f90 \
 	dispatch_transforms.f90 \
 	dispatch_transforms_uti.f90 \
 	dispatch_transforms_ut.f90 \
 	beam_structures_uti.f90 \
 	beam_structures_ut.f90 \
 	process_configurations.f90 \
 	process_configurations_uti.f90 \
 	process_configurations_ut.f90 \
 	compilations.f90 \
 	compilations_uti.f90 \
 	compilations_ut.f90 \
-	dispatch_fks.f90 \
 	integrations.f90 \
 	integrations_uti.f90 \
 	integrations_ut.f90 \
 	event_streams.f90 \
 	event_streams_uti.f90 \
 	event_streams_ut.f90 \
+	restricted_subprocesses.f90 \
+	eio_direct.f90 \
+	eio_direct_uti.f90 \
+	eio_direct_ut.f90 \
 	simulations.f90 \
+	restricted_subprocesses_uti.f90 \
+	restricted_subprocesses_ut.f90 \
 	simulations_uti.f90 \
 	simulations_ut.f90 \
 	commands.f90 \
 	commands_uti.f90 \
 	commands_ut.f90 \
 	cmdline_options.f90 \
 	libmanager.f90 \
 	features.f90 \
 	whizard.f90
 
 FC_OBJ = $(FC0_SRC:.f90=.o) $(F77_SRC:.f=.o) $(FC_SRC:.f90=.o) 
 CC_OBJ = $(CC_SRC:.c=.o)
 
-all: whizard_test whizard
+all: whizard_test 
 check: whizard_test
 	./whizard_test --check resonances
 
 whizard_test: $(FC_OBJ) $(CC_OBJ) main_ut.f90
 	$(FC) $(FC_OBJ) $(CC_OBJ) -ldl -o $@ main_ut.f90
 
 whizard: $(FC_OBJ) $(CC_OBJ) main.f90
 	$(FC) $(FC_OBJ) $(CC_OBJ) -ldl -o $@ main.f90
 
 %.o: %.f90
 	$(FC) $(FCFLAGS) -c $<
 
 %.o: %.f
 	$(FC) $(FCFLAGS) -c $<
 
 %.o: %.c
 	$(CC) $(CCFLAGS) -c $<
 
 tar: $(FC_SRC) $(F77_SRC) $(FC0_SRC) $(CC_SRC) $(MODELS)
 	tar cvvzf whizard-`date +%y%m%d`-`date +%H%M`.tar.gz $(FC_SRC) $(FC0_SRC) \
 	$(F77_SRC) $(CC_SRC) main_ut.f90 Makefile $(MODELS)
 
 clean:
 	rm -f *.mod *.o whizard_test