#include<iostream>
#include<cmath>
#include<fstream>
#include<string>
#include<TCanvas.h>
#include <TFile.h>
#include<TApplication.h>
#include<TGraphErrors.h>
#include<TH1D.h>
#include<TF1.h>
#include<TStyle.h>
#include<TAxis.h>
#include<TLatex.h>
#include<TMath.h>
#include<TTree.h>
#include<cstdlib>
#include<stdlib.h>
#include<TRandom.h>
#include "TH2F.h"
#include "TCanvas.h"
#include "TPad.h"
#include<cstring>
#include<sstream>  
#include<iomanip>

using namespace std;


void edepsplit() {

  
  		std::vector<double> *Edep  = new vector <double>();
  		double Edep0[1];
  		double Edep1[1];
  
  	TFile *RootIn = TFile::Open("Sim.root","READ");
  	TTree *hIn = (TTree*) (RootIn->Get("Tree1"));
  	TFile *RootOut = new TFile("SimSplit.root","RECREATE");
  	TTree *hOut = new TTree("Tree1", "EdepSplit");
    		hIn->SetBranchAddress("Edep", &Edep);
    		hOut->Branch("EdepGePD", &Edep0[0]);
    		hOut->Branch("EdepGeDD", &Edep1[0]);
  	

  	long int NEvents = hIn->GetEntries();
  	cout << "Number of Events: " << NEvents << endl;
  	long int Step = (long int)(NEvents/10);
  	for(long int i=0; i<NEvents; i++) {
    		if (i%Step==0) cout << (int)(i*1e2/NEvents) << "% events processed.." << endl;
    		hIn->GetEntry(i);
    		if (Edep->size()>0 && (Edep->at(0)>0 || Edep->at(1)>0)){
    			Edep0[0]=Edep->at(0);
				Edep1[0]=Edep->at(1); 
    		}
    	
    	hOut->Fill();
  	}
  	hOut->Write();
  	RootIn->Close();
  	RootOut->Close();
}