From 2f558d9e1c3ccc2c2b5ba086a4e0b352f7e1c99e Mon Sep 17 00:00:00 2001
From: Maria Zurek <zurek@anl.gov>
Date: Thu, 15 Jul 2021 00:59:55 +0000
Subject: [PATCH] Add plots with ScFi and Img energy deposits
---
.../emcal_barrel_particles_analysis.cxx | 47 +++++++++++++++++--
1 file changed, 44 insertions(+), 3 deletions(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
index 0e35ac33..536cb195 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_particles_analysis.cxx
@@ -109,9 +109,21 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
auto hEsim = d1.Histo1D(
{"hEsim", "Energy Deposit; Energy Deposit [GeV]; Events", 500, 0.0, 0.5},
"Esim");
+ auto hEsimImg = d1.Histo1D(
+ {"hEsimImg", "Energy Deposit; Energy Deposit [GeV]; Events", 500, 0.0, 0.5},
+ "EsimImg");
+ auto hEsimScFi = d1.Histo1D(
+ {"hEsimScFi", "Energy Deposit; Energy Deposit [GeV]; Events", 500, 0.0, 0.5},
+ "EsimScFi");
auto hfsam = d1.Histo1D(
- {"hfsam", "Sampling Fraction; Sampling Fraction; Events", 800, 0.0, 0.2},
+ {"hfsam", "Sampling Fraction; Sampling Fraction; Events", 400, 0.0, 0.2},
"fsam");
+ auto hfsamImg = d1.Histo1D(
+ {"hfsamImg", "Sampling Fraction; Sampling Fraction; Events", 400, 0.0, 0.2},
+ "fsamImg");
+ auto hfsam = d1.Histo1D(
+ {"hfsamScFi", "Sampling Fraction; Sampling Fraction; Events", 400, 0.0, 0.2},
+ "fsamScFi");
addDetectorName(detector_name, hEthr.GetPtr());
addDetectorName(detector_name, hEsim.GetPtr());
@@ -155,14 +167,13 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
{
TCanvas* c4 = new TCanvas("c4", "c4", 700, 500);
- c4->SetLogy(1);
auto h = hfsam->DrawCopy();
h->SetLineWidth(2);
h->SetLineColor(kBlue);
double up_fit = h->GetMean() + 5*h->GetStdDev();
double down_fit = h->GetMean() - 5*h->GetStdDev();
h->Fit("gaus", "", "", down_fit, up_fit);
- h->GetXaxis()->SetRangeUser(0.,up_fit);
+ h->GetXaxis()->SetRangeUser(down_fit,up_fit);
TF1 *gaus = h->GetFunction("gaus");
fSam_mean = gaus->GetParameter(1);
gaus->SetLineWidth(2);
@@ -170,6 +181,36 @@ void emcal_barrel_particles_analysis(std::string particle_name = "electron", boo
save_canvas(c4,"fsam",particle_name);
}
+ {
+ TCanvas* c5 = new TCanvas("c5", "c5", 700, 500);
+ auto h = hfsamImg->DrawCopy();
+ h->SetLineWidth(2);
+ h->SetLineColor(kBlue);
+ double up_fit = h->GetMean() + 5*h->GetStdDev();
+ double down_fit = h->GetMean() - 5*h->GetStdDev();
+ h->Fit("gaus", "", "", down_fit, up_fit);
+ h->GetXaxis()->SetRangeUser(down_fit,up_fit);
+ TF1 *gaus = h->GetFunction("gaus");
+ gaus->SetLineWidth(2);
+ gaus->SetLineColor(kRed);
+ save_canvas(c4,"fsamImg",particle_name);
+ }
+
+ {
+ TCanvas* c6 = new TCanvas("c6", "c6", 700, 500);
+ auto h = hfsamScFi->DrawCopy();
+ h->SetLineWidth(2);
+ h->SetLineColor(kBlue);
+ double up_fit = h->GetMean() + 5*h->GetStdDev();
+ double down_fit = h->GetMean() - 5*h->GetStdDev();
+ h->Fit("gaus", "", "", down_fit, up_fit);
+ h->GetXaxis()->SetRangeUser(down_fit,up_fit);
+ TF1 *gaus = h->GetFunction("gaus");
+ gaus->SetLineWidth(2);
+ gaus->SetLineColor(kRed);
+ save_canvas(c4,"fsamScFi",particle_name);
+ }
+
j[particle_name] = {
{"particle_name", particle_name},
{"thrown_energy", Ethr_mean},
--
GitLab