Testing results

benchmarks/barrel_ecal/scripts/emcal_barrel_pion_rejection_analysis.cxx

@@ -10,9 +10,9 @@
 #include "dd4pod/Geant4ParticleCollection.h"
 #include "dd4pod/CalorimeterHitCollection.h"
 
-#include "benchmark.h"
-#include "mt.h"
-#include "util.h"
+//#include "benchmark.h"
+//#include "mt.h"
+//#include "util.h"
 
 #include "TCanvas.h"
 #include "TStyle.h"
@@ -44,7 +44,7 @@ void emcal_barrel_pion_rejection_analysis(
                                           const char* input_fname4 = "../sim_output/sim_emcal_barrel_uniform_piminus_E5.root", 
                                           const char* input_fname5 = "../sim_output/sim_emcal_barrel_uniform_piminus_E10.root",
                                           const char* input_fname6 = "../sim_output/sim_emcal_barrel_uniform_piminus_E18.root"
-                                          */
+                                          /**/
                                           )
 {
   // Setting for graphs
@@ -233,8 +233,8 @@ void emcal_barrel_pion_rejection_analysis(
 
       // Gather ratio of pi/e for each Energy and eta bin
       // Then plot the distributions
-      rejRatios[i][j] = hp->Integral() / he->Integral();
-      std::cout << rejRatios[i][j] << std::endl;
+      rejRatios[i][j] = (double)hp->Integral() / (double)he->Integral();
+      std::cout << hp->Integral() << " " << he->Integral()  << " " << rejRatios[i][j] << std::endl;
 
       hp->Divide(he.GetPtr());
       title = "#pi^{-}/e^{-} (E = " + std::to_string((int)E[i]) + " GeV) : " + etaTitle[j];