diff options
Diffstat (limited to 'utils/plot-muons')
| -rwxr-xr-x | utils/plot-muons | 7 |
1 files changed, 2 insertions, 5 deletions
diff --git a/utils/plot-muons b/utils/plot-muons index fc8427c..d14f8b6 100755 --- a/utils/plot-muons +++ b/utils/plot-muons @@ -153,11 +153,8 @@ if __name__ == '__main__': stopping_muons = pd.merge(ev[ev.muon & ev.stopping_muon],michel,left_on=['run','gtid'],right_on=['run','muon_gtid'],suffixes=('','_michel')) stopping_muons_mc = pd.merge(ev_mc[ev_mc.muon & ev_mc.stopping_muon],michel_mc,left_on=['run','gtid'],right_on=['run','muon_gtid'],suffixes=('','_michel')) - stopping_muons = stopping_muons[stopping_muons.nhit_cal < 4000] - stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.nhit_cal < 4000] - - stopping_muons = stopping_muons[stopping_muons.udotr < -0.5] - stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.udotr < -0.5] + stopping_muons = stopping_muons[stopping_muons.ke < 10e3] + stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.ke < 10e3] stopping_muons = stopping_muons[stopping_muons.cos_theta < -0.5] stopping_muons_mc = stopping_muons_mc[stopping_muons_mc.cos_theta < -0.5] |