/*

FakePhysicist.com
admin@fakephysicist.com
14 April 2017

Compile with:
g++ genvidata.cpp -o genvidata `root-config --cflags --libs`

*/

#include <iostream>
#include <iomanip>
#include <cmath>
#include <fstream>
#include <vector>
#include <cstring>

/****** Include the required ROOT Classes *****/
#include "TRandom.h"

using namespace std;

int main() {

  int points = 11;		// No of voltage readings
  int sets = 100;		// Total no of measurement sets
  double maxV = 5.;		// Max Voltage
  double minV = 0.;		// Min Voltage
  double res = 1000.;		// Resistance Value
  double noise = 0.00005;	// Gaussian Noise of 0.05mA

  ofstream outfile("vidata.txt");
  
  for(int jk=0;jk<sets;jk++) {
    for(int ij=0;ij<points;ij++) {
      double volt = minV + ij*(maxV - minV)/(points - 1);
      // volt ==  0.,0.5,1.0,1.5,.....,5.0
      
      outfile << volt << "\t"
	      << gRandom->Gaus(volt/res,noise)
	// including Gaussian noise with mean at true position and sigma 0.05mA into the current measurement
	      << endl;
    }
  }

  return 0;
}