gROOT->Reset();

#include "Riostream.h"
#include <stdio.h>
#include <fstream>

void U62() {
   ifstream in;
   in.open("OutPut10.dat");

   Float_t n,Et,ECsI1,ECsI2,ECsI3,ECsI4,ECsI5,ECsI6,ECsI7,ECsI8,ECsI9,ECsI10,ECsI11,ECsI12,ECsIF,ECsIB,ECMO;

   Int_t nlines = 0;
   TFile *f = new TFile("OutPut10.root","RECREATE");
//   TH1F *h1 = new TH1F("h1","dN/dE",500,0,4999);
//   TH1F *h2 = new TH1F("h2","dN/dE",500,0,4999);
 
   TNtuple *ntuple = new TNtuple("ntuple","data from ascii file","n:Et:ECsI1:ECsI2:ECsI3:ECsI4:ECsI5:ECsI6:ECsI7:ECsI8:ECsI9:ECsI10:ECsI11:ECsI12:ECsIF:ECsIB:ECMO");

while (1) {
 
in >>n>>Et>>ECsI1>>ECsI2>>ECsI3>>ECsI4>>ECsI5>>ECsI6>>ECsI7>>ECsI8>>ECsI9>>ECsI10>>ECsI11>>ECsI12>>ECsIF>>ECsIB>>ECMO ;        
   
    if (!in.good()) break;
    if (nlines < 10) printf("n=%8f, Et=%8f, ECsI1=%8f, ECsI2=%8f, ECsI3=%8f, ECsI4=%8f, ECsI5=%8f, ECsI6=%8f, ECsI7=%8f, ECsI8=%8f, ECsI9=%8f, ECsI10=%8f, ECsI11=%8f, ECsI12=%8f, ECsIF=%8f, ECsIB=%8f, ECMO=%8f\n",n,Et,ECsI1,ECsI2,ECsI3,ECsI4,ECsI5,ECsI6,ECsI7,ECsI8,ECsI9,ECsI10,ECsI11,ECsI12,ECsIF,ECsIB,ECMO);

 //   if (ECsI1 + ECsI2 + ECsI3 + ECsI4 + ECsI5 + ECsI6 + ECsI7 + ECsI8 + ECsI9 + ECsI10 + ECsI11 + ECsI12 + ECsIF + ECsIB = 0)
  //    {h1->Fill(Et,ECMO);} //anticoincidences with CsI

      ntuple->Fill(n,Et,ECsI1,ECsI2,ECsI3,ECsI4,ECsI5,ECsI6,ECsI7,ECsI8,ECsI9,ECsI10,ECsI11,ECsI12,ECsIF,ECsIB,ECMO);

//      ntuple->Fill(n:Et:ECsI1:ECsI2:ECsI3:ECsI4:ECsI5:ECsI6:ECsI7:ECsI8:ECsI9:ECsI10:ECsI11:ECsI12:ECsIF:ECsIB:ECMO);

      nlines++;
          }
   printf(" found %d points\n",nlines);

   in.close();

   TCanvas *c=new TCanvas("c","",800,800); 
   ntuple->Draw("n:Et:ECsI1:ECsI2:ECsI3:ECsI4:ECsI5:ECsI6:ECsI7:ECsI8:ECsI9:ECsI10:ECsI11:ECsI12:ECsIF:ECsIB:ECMO");
//	 h1->Draw();
	// h2->Draw("same");
   f->Write();
}