testShowGraphs.cxx

#include <TTree.h>
#include <TGraph.h>
#include <TFile.h>
#include <TSystem.h>
#include <TCanvas.h>
#include "TStyle.h"

void testShowGraphs()
{
	using std::cout;
	using std::endl;

	gStyle->SetCanvasDefH(900);
	gStyle->SetCanvasDefW(1500);

	const char *foldername = "7_8";
	const char *ext = ".gif";

	gSystem->Load("../libData.so");

	const Long64_t kFirstEvent = 128;

//	TFile fr("../data/dataDSR4/Neurad_7_08_NEW.root");
//	TFile fr("../data/dataDSR4/analysis_07_8.root");

//	TFile fr("../data/dataDSR4/Neurad_081216_2NEW.root");

//	TFile fr("../data/dataTektronix/exp18.root");

	TFile fr("../data/dataDSR4/analysis_08_2.root");

//	TFile fr("../data/dataDSR4/Neurad_081216_2NEW.root");
//	TFile fr("../data/dataTektronix/exp18.root");

//	TFile fr("../data/dataTektronix/analysisExp7.root");
	TTree *tr = (TTree*)fr.Get("atree");

	AEvent *revent = new AEvent();
	tr->SetBranchAddress("Ach0.",&revent);

	TGraph *gr[13];
	Double_t t10[13];
	Double_t t90[13];

	//loop over events
	for (Long64_t i = 0; i < 13; i++) {
		gr[i] = 0;
		t10[i] = 0;
		t90[i] = 0;

		tr->GetEntry(i+kFirstEvent);

		gr[i] = new TGraph(*revent->GetGraphSignal());
		t10[i] = revent->GetT_10();
		t90[i] = revent->GetT_90();

//		cout << t10[i] << "\t" << t90[i] << endl;
//		cout << revent->GetEdgeSlope() << endl;


	}//for over events


	TCanvas *c1 = new TCanvas("c1","test",10,10,1000,600);
	c1->Divide(2,2);


	TF1 *f1 = new TF1("f1name", "[0]+x*[1]");
//	TF1 *f1 = new TF1("pol1", );


	for (Int_t k = 0; k < 4; k++) {
		c1->cd(k+1);
		gr[k]->GetXaxis()->SetRangeUser(130, 175);
		gr[k]->GetXaxis()->SetTitle("Time [ns]");
		gr[k]->GetXaxis()->CenterTitle();
		gr[k]->GetYaxis()->SetTitle("Signal [V]");
		gr[k]->GetYaxis()->CenterTitle();
		gr[k]->Draw("AL*");
		//cout << t10[k] << "\t" << t90[k] << endl;
		//f1->SetRange(t10[k], t90[k]);
		//gr[k]->Fit(f1, "RQ");
	}
	c1->Print(Form("../macros/picsDRS4/file%s/signals%s", foldername, ext));
//	c1->cd(2);
//	tr->Draw("gAmp.Draw()","","goff",1,124);
}