#include <TStyle.h>
#include <TCanvas.h>
#include <TMath.h>
#include <TGraph2D.h>
#include <TH2Poly.h>
#include <TClonesArray.h>
#include <iostream>
#include <sstream>
#include <fstream>
#include <vector>
#include <TVector3.h>
#include <TString.h>
#include <TROOT.h>
#include <TRint.h>
#include <TTree.h>
#include <TFile.h>
#include <TCutG.h>
#include <TChain.h>
#include <TH1F.h>
#include <TH2.h>
#include <TH2F.h>
#include <TF1.h>
#include <TH3.h>
#include <TH3F.h>
#include <TF3.h>
#include <TMarker3DBox.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <TMath.h>
#include <TGraph.h>
#include <TGraph2D.h>
#include <TRandom2.h>
#include <TStyle.h>
#include <TF2.h>
#include <TH1.h>
#include <TMinuit.h>
#include <TVirtualFitter.h>
#include <TFitter.h>
#include <TPolyLine3D.h>
#include <TH2Poly.h>
#include <TLine.h>
#include <TArc.h>
#include <sstream>
#include <TPolyMarker.h>
#include <algorithm>
using namespace std;
void analysis_SRC()
{
/*	TFile*f1=new TFile("../S018_rootfiles/physics/analys_v3_physics_run460.root");
	TFile*f2=new TFile("../S018_rootfiles/physics/analys_v3_p2p_physics_run460.root");*/
	//TTree *t=new TTree("t","t");
	TChain*t1=new TChain("tree","tree");
	//TChain*t2=new TChain("tree","tree");
	for(int i=19; i<20; i++)
	{
		//t2->Add(Form("rootfiles/psp15%04d_minos.root",i));
		t1->Add(Form("rootfiles/psp15%04d_brips.root",i));
	} 



//parfit b,d, number tracks; z, aoq
	vector<double> *parFit2=0;
	vector<double> *parFit4=0;
	int NumberTracks;
	double BigRIPSBeam_aoq[6];
	double BigRIPSBeam_zet[6];

    //t1->Print();
	TBranch* bparFit2;
	TBranch* bparFit4;
	TBranch* bNumberTracks;
	TBranch* bBigRIPSBeam_zet;
	TBranch* bBigRIPSBeam_aoq;
	t1->GetBranch("BigRIPSBeam")->SetMakeClass(true);
	t1->SetBranchAddress("BigRIPSBeam.zet",&(BigRIPSBeam_zet[0]),&bBigRIPSBeam_zet);
	//t1->SetBranchAddress("BigRIPSBeam.zet",BigRIPSBeam_zet,&(bBigRIPSBeam_zet[0]));
	//t1->SetBranchAddress("BigRIPSBeam.aoq",BigRIPSBeam_aoq);//,&bBigRIPSBeam_aoq);
	/*t2->SetBranchAddress("parFit2",&parFit2,&bparFit2);
	t2->SetBranchAddress("parFit4",&parFit4,&bparFit4);
	t2->SetBranchAddress("NumberTracks",&NumberTracks,&bNumberTracks);
	TH1F*h_phi=new TH1F("h_phi","h_phi",180,0,180);
	double phi[3];
	TFile*fout=new TFile("out.root","recreate");

	TTree*tout=new TTree("tout","tout");
	tout->Branch("phi",phi,"phi[4]/D");
	tout->Branch("NumberTracks",&NumberTracks,"NumberTracks/I");


	TFile*ff=new TFile("rootfiles/cut.root");
	TCutG*cutbr=(TCutG*)ff->Get("brGe85");
	TCutG*cutzd=(TCutG*)ff->Get("zdZn83");
	cout<<"last printout "<<t1->GetEntries()<<endl;*/
	for(int i=0; i<t1->GetEntries(); i++)
	{
		t1->GetEntry(i);
		cout<<i<<" "<<BigRIPSBeam_zet[4]<<endl;
		//t2->GetEntry(i);
		/*if(cutbr->IsInside(BigRIPSBeam_aoq[2],BigRIPSBeam_zet[2])&&cutzd->IsInside(BigRIPSBeam_aoq[4],BigRIPSBeam_zet[4])&&NumberTracks==3) 
		{
			//TVector3 v(0,0,1);
			TVector3 v0(parFit2->at(0),parFit4->at(0),1);
			TVector3 v1(parFit2->at(1),parFit4->at(1),1);
			TVector3 v2(parFit2->at(2),parFit4->at(2),1);
			phi[0]=v0.Phi();phi[1]=v1.Phi();phi[2]=v2.Phi();
			for(int j=0; j<NumberTracks; j++)
			{
				for(int k=j+1; k<NumberTracks; k++)
				{
					h_phi->Fill(abs(phi[j]-phi[k]));
				}
			}				
		}*/
		//tout->Fill();
	}
	//fout->cd();
	//h_phi->Write();tout->Write();

	

}