diff --git a/convertTektronix/exp2.root b/convertTektronix/exp2.root index 89c1b247325b02523a365e7641a7ca488e23610e..862ca8338b506ba016b1f4d7bba94a7b118cdf54 100644 Binary files a/convertTektronix/exp2.root and b/convertTektronix/exp2.root differ diff --git a/convertTektronix/infiles.dat b/convertTektronix/infiles.dat new file mode 100644 index 0000000000000000000000000000000000000000..7e7f5493c9a5dfe5baad37657ff23b78487d4831 --- /dev/null +++ b/convertTektronix/infiles.dat @@ -0,0 +1,12 @@ +ch1 +../data/rawDataTektronix/ch12016.12.07-02.11.54.dat +../data/rawDataTektronix/ch12016.12.07-08.48.05.dat +../data/rawDataTektronix/ch12016.12.07-10.03.01.dat +ch2 +... +... +... +ch3 +... +... +... \ No newline at end of file diff --git a/convertTektronix/read.cpp b/convertTektronix/read.cpp index 8387096c081b1ebb9913364c42281705127579c6..551bde153ca67e0915c0cde6e25ba48c928cc959 100644 --- a/convertTektronix/read.cpp +++ b/convertTektronix/read.cpp @@ -1,15 +1,23 @@ -#include -#include - -#include "TFile.h" -#include "TTree.h" - -#include "../dataClasses/RawEvent.h" +//#include +//#include +// +//#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; @@ -19,7 +27,7 @@ void read() { tree->Branch("A",A,"A[1000]/D"); tree->Branch("T",T,"T[1000]/D"); - RawEvent *event = new RawEvent; + RawEvent *event = new RawEvent(); tree->Bronch("rawEvent", "RawEvent", &event); @@ -34,17 +42,31 @@ void read() { return; } - for(i = 0; i<100000; i++){ + 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<> A[i-1000*n]; - T[i-1000*n]=5e-013+ (i-1000*n)*1e-010; + + 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; + T[j]=0; + A[j]=0; } + + event->Reset(); } } myfile1.close(); diff --git a/dataClasses/RawEvent.cpp b/dataClasses/RawEvent.cpp index 37435800c08772c7d25b7f6636b86555e532ddac..6e9be6be8ec270cf6a951228f071128f104f7810 100644 --- a/dataClasses/RawEvent.cpp +++ b/dataClasses/RawEvent.cpp @@ -45,8 +45,9 @@ void RawEvent::PrintAmp(Int_t i) { } void RawEvent::SetAmp(Double_t a, Int_t i) { +// cout << fNPoints << endl; if (i >= fNPoints) { - cout << "Error: array with raw amplitudes is overloaded!" << endl; + Error("RawEvent::SetAmp", "Array with raw amplitudes is overloaded!"); return; } fAmp[i] = a; @@ -55,7 +56,7 @@ void RawEvent::SetAmp(Double_t a, Int_t i) { void RawEvent::SetTime(Double_t t, Int_t i) { if (i >=fNPoints) { - cout << "Error: array with raw times is overloaded!" << endl; + Error("RawEvent::SetTime", "Array with raw times is overloaded!"); return; } fTime[i] = t;