#include #include #include #include #include #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_07_8.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); }