code: def round_trip(
      simu_decays=[0.1, 0.5, 1.0],
      simu_baseline=[0.12, 0.07],
      simu_adjacency=[
          [[0, 0.1, 0.4], [0.2, 0.0, 0.2]],
          [[0, 0, 0], [0.6, 0.3, 0]],
      ],
      gfit_decays=[0.2, 0.6, 1.0],
      gfit_baseline=[0.12, 0.07],
      gfit_adjacency=[
          [[0, 0.1, 0.8], [0.2, 0.8, 0.2]],
          [[0.3, 0.3, 0], [0.6, 0.3, 0]],
      ],
      end_time=1000,
      max_jumps=5000,
      verbose=True,
      force_simulation=False,
  ):
      """
      If any of `gfit_` parameters do not match with the `simu_` parameters,
      the qqplots of `gfit_model` are expected to show poor fit.
      """
      simu_model = simulate_hawkes_sum_exp_kern(
          adjacency=simu_adjacency,
          decays=simu_decays,
          baseline=simu_baseline,
          end_time=end_time,
          max_jumps=max_jumps,
          verbose=verbose,
          force_simulation=force_simulation,
      )
      gfit_model = model_from_timestamps(
          simu_model.timestamps,
          adjacency=gfit_adjacency,
          decays=gfit_decays,
          baseline=gfit_baseline,
          end_time=end_time,
      )
      return simu_model, gfit_model
  
imports: from tick.hawkes import SimuHawkesSumExpKernels
  from tick.plot import qq_plots
label_desc: Use image node_id 4 for calling a global function with example usage: round_trip(simu_decays, simu_baseline, simu_adjacency, gfit_decays, gfit_baseline, gfit_adjacency, end_time, max_jumps, verbose, force_simulation) and returns: simu_model,  gfit_model