Commit 3628b433 authored by Muzalevsky I.A's avatar Muzalevsky I.A

added script for converting data from txt files

parent e51e9398
void read() {
//Double_t N[5]; // массив для чисел ()
//ifstream myfile;
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");
//std::vector<Double_t> *A = new std::vector<Double_t>();
//tree->Branch("mybranch","vector<Double_t>",A);
ifstream myfile1;
myfile1.open("ch12016.12.07-02.11.54.dat");
for(i = 0; i<100000; i++){
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;
if(i-1000*n ==999) {
tree->Fill();
for(j=0;j<1000;j++){
T[j]=0;A[j]=0;
}
}
}
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();
}
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