Skip to content
Snippets Groups Projects
Commit 61b5ff34 authored by Marshall Scott's avatar Marshall Scott
Browse files

Commented out benchmark.h, util.h, and mt.h; added TF1 for gaussian fit

parent 46b41f2d
No related branches found
No related tags found
No related merge requests found
This commit is part of merge request !24. Comments created here will be created in the context of that merge request.
...@@ -147,6 +147,8 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -147,6 +147,8 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -7.5, 7.5}, "dE_rel"); auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -7.5, 7.5}, "dE_rel");
auto f1 = hdE_rel->Fit("gaus", "S"); auto f1 = hdE_rel->Fit("gaus", "S");
const double* res = f1->GetParams(); const double* res = f1->GetParams();
auto tf1 = new TF1("", "TMath::Gaus(x, [0], [1], [2])", -7.5, 7.5);
tf1 -> SetParameters(res[0], res[1], res[2]);
//Pass/Fail //Pass/Fail
/* /*
...@@ -172,10 +174,10 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -172,10 +174,10 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
hdE_rel->GetYaxis()->SetTitleOffset(1.4); hdE_rel->GetYaxis()->SetTitleOffset(1.4);
hdE_rel->SetLineWidth(2); hdE_rel->SetLineWidth(2);
hdE_rel->SetLineColor(kBlue); hdE_rel->SetLineColor(kBlue);
f1->SetLineWidth(2); tf1->SetLineWidth(2);
f1->SetLineColor(kRed); tf1->SetLineColor(kRed);
hdE_rel->DrawClone(); hdE_rel->DrawClone();
f1->Draw("SAME"); tf1->Draw("SAME");
cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.png"); cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.png");
cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.pdf"); cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.pdf");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment