Commit 66cd88ed authored by Vratislav Chudoba's avatar Vratislav Chudoba

convertTektronix fill branch with RawEvent class.

parent 11a99d97
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
#include <iostream> //#include <iostream>
#include <fstream> //#include <fstream>
//
#include "TFile.h" //#include "TFile.h"
#include "TTree.h" //#include "TTree.h"
//
#include "../dataClasses/RawEvent.h" //#include "../dataClasses/RawEvent.h"
void read() { void read() {
//Double_t N[5]; // массив для чисел () //Double_t N[5]; // массив для чисел ()
//ifstream myfile; //ifstream myfile;
const Int_t par1 = 5;
const Double_t par2 = 5.34;
// ....
// ....
gSystem->Load("../libData.so");
Double_t A[1000]; Double_t A[1000];
Double_t T[1000]; Double_t T[1000];
Int_t i,j,n; Int_t i,j,n;
...@@ -19,7 +27,7 @@ void read() { ...@@ -19,7 +27,7 @@ void read() {
tree->Branch("A",A,"A[1000]/D"); tree->Branch("A",A,"A[1000]/D");
tree->Branch("T",T,"T[1000]/D"); tree->Branch("T",T,"T[1000]/D");
RawEvent *event = new RawEvent; RawEvent *event = new RawEvent();
tree->Bronch("rawEvent", "RawEvent", &event); tree->Bronch("rawEvent", "RawEvent", &event);
...@@ -34,17 +42,31 @@ void read() { ...@@ -34,17 +42,31 @@ void read() {
return; return;
} }
Double_t amp = 0;
Double_t time = 0;
//rewrite using while-loop
for(i = 0; i<100000; i++){ for(i = 0; i<100000; i++){
n=i/1000; // if ((i==1002)||(i==900)||(i==2002)||(i==8002)) {cout<<i<<" "<<n<<endl;} n=i/1000; // if ((i==1002)||(i==900)||(i==2002)||(i==8002)) {cout<<i<<" "<<n<<endl;}
myfile1 >> 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) { if(i-1000*n ==999) {
tree->Fill(); tree->Fill();
for(j=0;j<1000;j++){ for(j=0;j<1000;j++){
T[j]=0;A[j]=0; T[j]=0;
A[j]=0;
} }
event->Reset();
} }
} }
myfile1.close(); myfile1.close();
......
...@@ -45,8 +45,9 @@ void RawEvent::PrintAmp(Int_t i) { ...@@ -45,8 +45,9 @@ void RawEvent::PrintAmp(Int_t i) {
} }
void RawEvent::SetAmp(Double_t a, Int_t i) { void RawEvent::SetAmp(Double_t a, Int_t i) {
// cout << fNPoints << endl;
if (i >= fNPoints) { if (i >= fNPoints) {
cout << "Error: array with raw amplitudes is overloaded!" << endl; Error("RawEvent::SetAmp", "Array with raw amplitudes is overloaded!");
return; return;
} }
fAmp[i] = a; fAmp[i] = a;
...@@ -55,7 +56,7 @@ void RawEvent::SetAmp(Double_t a, Int_t i) { ...@@ -55,7 +56,7 @@ void RawEvent::SetAmp(Double_t a, Int_t i) {
void RawEvent::SetTime(Double_t t, Int_t i) { void RawEvent::SetTime(Double_t t, Int_t i) {
if (i >=fNPoints) { if (i >=fNPoints) {
cout << "Error: array with raw times is overloaded!" << endl; Error("RawEvent::SetTime", "Array with raw times is overloaded!");
return; return;
} }
fTime[i] = t; fTime[i] = t;
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment