read.cpp

//#include <iostream>
//#include <fstream>
//
//#include "TFile.h"
//#include "TTree.h"
//
//#include "../dataClasses/RawEvent.h"

void read() {
	//Double_t N[5]; // массив для чисел ()
	//ifstream myfile;

	const Int_t par1 = 5;
	const Double_t par2 = 5.34;
	// ....
	// ....


	gSystem->Load("../libData.so");

	Double_t A[1000];
	Double_t T[1000];
	Int_t i,j,n;

	TFile *f1 = new TFile("exp2.root","RECREATE");
	TTree *tree = new TTree("tree","signal");
	tree->Branch("A",A,"A[1000]/D");
	tree->Branch("T",T,"T[1000]/D");

	RawEvent *event = new RawEvent();
	tree->Bronch("rawEvent", "RawEvent", &event);


	//std::vector<Double_t> *A = new std::vector<Double_t>();
	//tree->Branch("mybranch","vector<Double_t>",A);

	ifstream myfile1;
	myfile1.open("../data/rawDataTektronix/ch12016.12.07-02.11.54.dat");

	if (!myfile1.is_open()) {
		Error("read.c", "Some error when opening file");
		return;
	}

	Double_t amp = 0;
	Double_t time = 0;

	//rewrite using while-loop
	for(i = 0; i<100000; i++){
		n=i/1000;			// if ((i==1002)||(i==900)||(i==2002)||(i==8002)) {cout<<i<<" "<<n<<endl;}

		myfile1 >> amp;
		time = 5e-013+ (i-1000*n)*1e-010;

		A[i-1000*n] = amp;
		T[i-1000*n] = time;

		event->SetAmp(amp, i-1000*n);
		event->SetTime(time, i-1000*n);

		if(i-1000*n ==999) {
			tree->Fill();

			for(j=0;j<1000;j++){
				T[j]=0;
				A[j]=0;
			}

			event->Reset();
		}
	}
	myfile1.close();

	ifstream myfile2;
	myfile2.open("ch12016.12.07-08.48.05.dat");
	for(i = 0; i<100000; i++){	
		n=i/1000;			// if ((i==1002)||(i==900)||(i==2002)||(i==8002)) {cout<<i<<" "<<n<<endl;}
		myfile2 >> A[i-1000*n];	
		T[i-1000*n]=5e-013+ (i-1000*n)*1e-010;	

		if(i-1000*n ==999){
			tree->Fill();

			for(j=0;j<1000;j++){
				T[j]=0;A[j]=0;
			}
		}
	}
	myfile2.close();

	ifstream myfile3;
	myfile3.open("ch12016.12.07-10.03.01.dat");
	for(i = 0; i<29000; i++){	
		n=i/1000;// if ((i==1002)||(i==900)||(i==2002)||(i==8002)) {cout<<i<<" "<<n<<endl;}
		myfile3 >> A[i-1000*n];	
		T[i-1000*n]=5e-013+ (i-1000*n)*1e-010;	

		if(i-1000*n ==999) {
			tree->Fill();

			for(j=0;j<1000;j++){
				T[j]=0;A[j]=0;
			}
		}
	}
	myfile3.close();

	tree->Write();
	f1->Close();
}