void SumDiv1()
{
TFile *f=new TFile("SPEC_IONEXCElecDNA3L.root"); // opens the root file energy ion = ana = energy; count ion = ana1 = energy1; energy exc = ana3 = energy2; count exc = ana = energy3;
TH1D *trEI=(TH1D*)f->Get("energy"); //energy X
TH1D *trCI=(TH1D*)f->Get("energy1"); //count Y
TH1D *trCE=(TH1D*)f->Get("energy3"); // count Y
TH1D *trEE=(TH1D*)f->Get("energy2"); // energy X
//TH1F *tr; f->GetObject("energy", tr);

float energyE, energyI ,countE, countI, value_normalized;
double integral = trEI->Integral();
double integ = trCI->Integral();
double integ1 = trEE->Integral();
double integ2 = trCE->Integral();
   
ofstream myfile;
myfile.open ("SomaEI.txt");
//myfile.open ("values1.txt");
//myfile << "energy count \n"; 

for (int idx=0,idx3=0idx1=0,idx2=0; idx<trEI->GetNbinsX(),idx2<trEE->GetNbinsX() idx1<trCI->GetNbinsX(),  idx3<trCE->GetNbinsX();++idx,++idx3, ++idx1,++idx2) {
        energyI = trEI->GetBinCenter(idx);
        energyE = trEE->GetBinCenter(idx2);


                countI = (trCI->GetBinContent(idx1));
   
                countE = (trCE->GetBinContent(idx3));

   
//      
  
   //value_normalized = count * (integral);
   
   
   
   
   
   // value_normalized = count/integ;
   cout << energyE << " " << countE << " " << energyI << " " << countI << " " << endl; //print to the screen
   myfile << energyE << " " << countE << " " << energyI << " " << countI << " " << "\n"; //write to file
   
//   cout << energy << " " << value_normalized <<endl; print to the screen
//   myfile << energy << " " << value_normalized << "\n"; write to file

}

cout<<"integral = " << integral << endl;
cout<<"integ = " << integ << endl;
cout<<"integ = " << integ1 << endl;
cout<<"integ = " << integ2 << endl;

myfile.close();
}