void testShowGraphs() { gSystem->Load("../libData.so"); const Long64_t kFirstEvent = 50; TFile fr("../data/dataDSR4/analysis_08_2.root"); // TFile fr("../data/dataDSR4/analysis_07_8.root"); // TFile fr("../data/dataTektronix/exp18.root"); // TFile fr("../data/dataTektronix/analysisExp7.root"); TTree *tr = (TTree*)fr.Get("atree"); AEvent *revent = new AEvent(); cout << tr->SetBranchAddress("Ach0.",&revent) << endl; TGraph *gr[10]; Double_t t10[10]; Double_t t90[10]; // cout << "warning" << endl; cout << tr->GetEntries() << endl; //loop over events for (Long64_t i = 0; i < 10; i++) { gr[i] = 0; t10[i] = 0; t90[i] = 0; // cout << "warning 1" << endl; // cout << i+kFirstEvent << endl; tr->GetEntry(i+kFirstEvent); // tr->GetEntry(i); // cout << "warning 2" << endl; 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 // cout << "warning 3" << endl; TCanvas *c1 = new TCanvas("c1","test",10,10,1000,600); c1->Divide(3,2); TF1 *f1 = new TF1("f1name", "[0]+x*[1]"); // TF1 *f1 = new TF1("pol1", ); for (Int_t k = 0; k < 6; k++) { c1->cd(k+1); gr[k]->Draw("AL*"); cout << t10[k] << "\t" << t90[k] << endl; f1->SetRange(t10[k], t90[k]); gr[k]->Fit(f1, "RQ"); } // c1->cd(2); // tr->Draw("gAmp.Draw()","","goff",1,124); }