diff --git a/macros/histCorr.C b/macros/histCorr.C index cffc5bf2b310989f47249e4e076a86cbefbeb7a1..debf8df01fa9e0eb36446899270492bab36b87d9 100644 --- a/macros/histCorr.C +++ b/macros/histCorr.C @@ -1,6 +1,6 @@ //#include "TCanvas.h" -void histCorr(const char *filename = "analysis_07_8.root", const char *foldername = "7_8", const char *ext = ".eps", const Bool_t can1 = 1, const Bool_t can2 = 1, const Bool_t can3 = 1, const Bool_t can4 = 1) +void histCorr(const char *filename = "analysis_07_8.root", const char *foldername = "7_8", const char *ext = ".eps", const Bool_t can1 = 0, const Bool_t can2 = 0, const Bool_t can3 = 0, const Bool_t can4 = 0, const Bool_t can5 = 1) { gSystem->Load("../libData.so"); TFile *f = new TFile(Form("../data/dataDSR4/%s", filename)); @@ -154,4 +154,21 @@ void histCorr(const char *filename = "analysis_07_8.root", const char *foldernam } + if (can5) { + TCanvas *c5 = new TCanvas("c5", "Sigma"); + c5->cd(); + TH1F *h11 = new TH1F("h11", "sigma", 1023, -10, 10); + + t->Draw("Ach1.fAmpPos >> h11","Ach1.fTime >0","col"); + //x + //h10->GetXaxis()->SetTitle("\\Delta \\tau_{1}^{(rise)} [ns]"); + //h10->GetXaxis()->CenterTitle(); + //y + //h10->GetYaxis()->SetTitle("fEdgeXi"); + //h10->GetYaxis()->CenterTitle(); + c5->Update(); + //c4->Print(Form("../macros/picsDRS4/file%s/XiSquares%s", foldername, ext)); + + } + } diff --git a/macros/testShowGraphs.cxx b/macros/testShowGraphs.cxx index af320168f569eb9a8f9ca6637da1056823582c48..6f9b399c41a413a687abcfd1030b5b0099a06bc8 100644 --- a/macros/testShowGraphs.cxx +++ b/macros/testShowGraphs.cxx @@ -7,7 +7,7 @@ void testShowGraphs() // TFile fr("../data/dataDSR4/analysis_07_1.root"); // TFile fr("../data/dataDSR4/analysis_07_8.root"); - TFile fr("../data/dataTektronix/exp18.root"); + TFile fr("../data/dataDSR4/Neurad_081216_2NEW.root"); // TFile fr("../data/dataTektronix/analysisExp7.root"); TTree *tr = (TTree*)fr.Get("atree"); @@ -45,19 +45,19 @@ void testShowGraphs() // cout << "warning 3" << endl; TCanvas *c1 = new TCanvas("c1","test",10,10,1000,600); - c1->Divide(3,2); + //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++) { + /*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"); - } + //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);