Update Analysis Script

Adding analysis script

• Put a cut on the energy difference, | E_{\mathrm{cluster}} - E_{\mathrm{thrown}} | < \delta_{E} or something similar.
• acceptance as a function of P_{T} , \eta, etc., for bins in different bins in E_{\mathrm{thrown}}
• resolution
• mass reconstruction

ecal/scripts/emcal_electrons_analysis.cxx

@@ -138,7 +138,7 @@ void emcal_electrons_analysis(const char* input_fname = "rec_electron_10kEvt.roo
   auto d2 = d1.Filter("ncluster==1");
   auto hClusterE1       = d2.Histo1D({"hClusterE1", "One Cluster Energy; Cluster Energy [GeV]; Events", 100, -0.5, 30.5},                      "clusterE");
   auto hEres            = d2.Histo1D({"hEres",      "Energy Resolution; #DeltaE/E; Events", 100,-0.3,0.3},                                     "E_res");
-  auto hPthr_accepted   = d2.Filter([] (float x, float y) { return x/y > 0.90; }, {"clusterE", "E_thr"})
+  auto hPthr_accepted   = d2.Filter("clusterE/E_thr > 0.9")
 	  		    .Histo1D({"hPthr_accepted", "Thrown momentum for reconstructed particle; p_{thr} [GeV]; Events", 100, -0.5, 30.5}, "p_thr");
 
   // Draw Histograms