diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-10-09 12:55:51 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-10-09 12:55:51 -0500 |
| commit | 17e6899bcc57a5d24021cd4b3a9427df2e5628b9 (patch) | |
| tree | 7f18c73687772776a5c9d4192d0b57cbd29ba503 /utils/chi2 | |
| parent | fb7d6a0ccbeadfd9b05ac2d5bd82af5389e826d9 (diff) | |
| download | sddm-17e6899bcc57a5d24021cd4b3a9427df2e5628b9.tar.gz sddm-17e6899bcc57a5d24021cd4b3a9427df2e5628b9.tar.bz2 sddm-17e6899bcc57a5d24021cd4b3a9427df2e5628b9.zip | |
only include Monte Carlo files for the runs we are analyzing
Diffstat (limited to 'utils/chi2')
| -rwxr-xr-x | utils/chi2 | 9 |
1 files changed, 9 insertions, 0 deletions
@@ -482,9 +482,18 @@ if __name__ == '__main__': #muon = ev[ev.muon & ev.prompt & ~ev.atm] #muon_atm = ev[ev.muon & ev.prompt & ev.atm] + if (~rhdr.run.isin(ev_mc.run)).any(): + print_warning("Error! The following runs have no Monte Carlo: %s" % \ + rhdr.run[~rhdr.run.isin(ev_mc.run)].values) + sys.exit(1) + + ev_mc = ev_mc[ev_mc.run.isin(rhdr.run)] + data_mc = ev_mc[ev_mc.signal & ev_mc.prompt & ~ev_mc.atm] data_atm_mc = ev_mc[ev_mc.signal & ev_mc.prompt & ev_mc.atm] + data_mc = ev_mc[ev_mc.signal & ev_mc.prompt & ~ev_mc.atm] + bins = np.logspace(np.log10(20),np.log10(10e3),21) atmo_scale_factor = 100.0 |