Skip to content
Snippets Groups Projects
Commit 8e014d1a authored by Jihee Kim's avatar Jihee Kim
Browse files

fixed filter

parent 2137af76
Branches
Tags
1 merge request!36Update Analysis Script
This commit is part of merge request !36. Comments created here will be created in the context of that merge request.
Hide whitespace changes
Inline Side-by-side
ecal/scripts/emcal_electrons_analysis.cxx
+ 1
− 1
View file @ 8e014d1a
...@@ -138,7 +138,7 @@ void emcal_electrons_analysis(const char* input_fname = "rec_electron_10kEvt.roo ...@@ -138,7 +138,7 @@ void emcal_electrons_analysis(const char* input_fname = "rec_electron_10kEvt.roo
auto d2 = d1.Filter("ncluster==1"); 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 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 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"); .Histo1D({"hPthr_accepted", "Thrown momentum for reconstructed particle; p_{thr} [GeV]; Events", 100, -0.5, 30.5}, "p_thr");
// Draw Histograms // Draw Histograms
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment