From 442ef59d9332e70f23235601436e5a465075fc91 Mon Sep 17 00:00:00 2001
From: Maria <zurek@anl.gov>
Date: Fri, 17 Sep 2021 13:58:10 -0500
Subject: [PATCH] Add momentum smearing

---
 .../scripts/emcal_barrel_energy_scan_analysis.cxx           | 2 +-
 .../scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx      | 6 +++---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis.cxx b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis.cxx
index c47e4f71..1febeccd 100644
--- a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis.cxx
+++ b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis.cxx
@@ -170,7 +170,7 @@ std::tuple <double, double, double, double, double, double, double, double> extr
   hErecImg->SetLineWidth(2);
   hErecImg->SetLineColor(kBlue);
   hErecImg->DrawClone();
-  set_histo_range(hErecImg);
+  //set_histo_range(hErecImg);
   c2->cd(4);
   hErecScFi->GetYaxis()->SetTitleOffset(1.4);
   hErecScFi->SetLineWidth(2);
diff --git a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
index 977215bd..e8e7a093 100644
--- a/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
+++ b/benchmarks/imaging_ecal/scripts/emcal_barrel_energy_scan_analysis_pi_e.cxx
@@ -83,15 +83,15 @@ std::tuple <double, double> extract_sampling_fraction_parameters(std::string E_l
     auto p = input[2];
     auto momentum = TMath::Sqrt(p.ps.x * p.ps.x + p.ps.y * p.ps.y + p.ps.z * p.ps.z);
     return momentum;
-  }
+  };
 
   // Smeared Momentum [GeV/c]
   auto momentum_smeared = [](const double momentum) {
     TRandom *smear = new TRandom();
     auto sigma_relative = 0.005;
-    auto sigma = smear->Gauss(0, sigma_relative*momentum);
+    auto sigma = smear->Gaus(0, sigma_relative*momentum);
     return momentum + sigma;
-  }
+  };
 
   // Reconstructed Energy [GeV]
   auto Erec = [](const std::vector<eic::CalorimeterHitData>& evt) {
-- 
GitLab