rm(list=ls())


# Sample
# -----------------------
sample.freq  <- 1     # 1 Sek.
sample.T     <- 1 / sample.freq

sample.dauer <- 10    # Sample-Time
sample_num   <- sample.dauer * sample.freq  # Anzahl der Samples
# ------------------------------------------------------------------
measurement.freq  <- 6.528
measurement.T     <- 1 / measurement.freq

impulse_num <- measurement.freq * sample.dauer
####################################################################

measurement.impulse_times <- (0:impulse_num) * measurement.T

sample_times <- 0:(sample_num-1) * sample.T

impuls.messreihe <- numeric(0)

make_one_sample_set <- function()
{
cat("*** make_one_sample_set()\n")
  # für alle (start_intervall ... end_intervall) mache:
  # ---------------------------------------------------
  for( start_int in 1 : length(sample_times) - 1 )
  {
    end_int <- start_int + 1   # End-Intervall
    print( start_int)

    selected_impulses <-
      measurement.impulse_times[ (measurement.impulse_times >= start_int) & (measurement.impulse_times < end_int) ]

    impuls.anzahl <- length(selected_impulses)

    print( paste("Anzahl der Impulse im Samplebereich von", start_int, "bis", end_int, ":", length(selected_impulses) ) )
    impuls.messreihe [end_int] <- impuls.anzahl

    # print(selected_impulses)

    #print( paste(start_int, end_int) )
  }
cat("** after loop\n")

  cat("--------------------------------------------------------\n")
  cat(paste(" Sample-Frequenz:           ", sample.freq, "\n"))
  cat(paste(" Sample-Periodendauer:      ", sample.T), "\n")
  cat(paste(" Sample-Dauer               ", sample.dauer), "\n")
  cat(paste(" Impuls-Frequenz Meas.Equip.", measurement.freq, "\n"))

cat("** creating res.txt\n")

  res.txt <- paste(#"\n\n Anzahl der Messungen:", sample_times, "\n", # !!!!!HERE!!!!!
                   "Summe der Impulse: ", sum(impuls.messreihe), "\n",
                   "Impulse pro Messdauer (+/- 0):", sum(impuls.messreihe) / sample.dauer, "\n", "\n",
                   "Mittelwert der Impulse:", mean(impuls.messreihe), "\n",
                   "Standardabweichung:", sd(impuls.messreihe),"\n"
                  )

cat("** printing res.txt\n")
  cat(res.txt)
  cat("================\n")
  cat("Attention, output of \"sample_times\" only:\n")
  cat(sample_times)
  cat("\nOK, this was \"sample_times\"\n")
  cat("================\n")

cat("*** make_one_sample_set()\n")

}



make_one_sample_set()