00001 #ifndef MuonThresholds_h
00002 #define MuonThresholds_h
00003
00004 #include "TrigMuonEfficiency/MuonHypothesis.h"
00005
00006 namespace MuonHypoThresholds {
00007
00008 const float hypo_etaregion[] = {1.05, 1.5, 2.0, 9.9};
00009 size_t hypo_binsize = sizeof(hypo_etaregion)/sizeof(hypo_etaregion[0]);
00010
00011
00012 const double ef_cb_4[] = { 3.93, 3.91, 3.88, 3.88};
00013 const double ef_cb_6[] = { 5.88, 5.81, 5.78, 5.76};
00014 const double ef_cb_8[] = { 7.82, 7.74, 7.70, 7.72};
00015 const double ef_cb_10[] = { 9.77, 9.67, 9.62, 9.57};
00016 const double ef_cb_13[] = { 12.67, 12.55, 12.49, 12.46};
00017 const double ef_cb_15[] = { 14.63, 14.49, 14.42, 14.38};
00018 const double ef_cb_18[] = { 17.53, 17.39, 17.34, 17.28};
00019 const double ef_cb_24[] = { 23.34, 23.19, 23.14, 23.06};
00020 const double ef_cb_36[] = { 34.96, 34.78, 34.69, 34.63};
00021 const double ef_cb_40[] = { 38.76, 38.54, 38.38, 38.31};
00022
00023 const MuonHypothesis::EF_thresholdPair ef_cb_pair[] = {
00024 MuonHypothesis::EF_thresholdPair("4GeV", ef_cb_4),
00025 MuonHypothesis::EF_thresholdPair("6GeV", ef_cb_6),
00026 MuonHypothesis::EF_thresholdPair("8GeV", ef_cb_8),
00027 MuonHypothesis::EF_thresholdPair("10GeV", ef_cb_10),
00028 MuonHypothesis::EF_thresholdPair("13GeV", ef_cb_13),
00029 MuonHypothesis::EF_thresholdPair("15GeV", ef_cb_15),
00030 MuonHypothesis::EF_thresholdPair("18GeV", ef_cb_18),
00031 MuonHypothesis::EF_thresholdPair("24GeV", ef_cb_24),
00032 MuonHypothesis::EF_thresholdPair("36GeV", ef_cb_36),
00033 MuonHypothesis::EF_thresholdPair("40GeV", ef_cb_40),
00034 };
00035
00036 const MuonHypothesis::EF_thresholdMap ef_cb_map(ef_cb_pair, ef_cb_pair + sizeof(ef_cb_pair)/sizeof(ef_cb_pair[0]));
00037
00038
00039
00040
00041
00042 const double ef_sa_40[] = { 31.5, 30.0, 28.5, 32.5};
00043 const double ef_sa_40_barrelOnly[] = { 31.5,1000.0,1000.0,1000.0};
00044 const double ef_sa_40_uptoEC2[] = { 31.5, 30.0, 28.5, 1000.0};
00045 const double ef_sa_50_barrelOnly[] = { 45.0,1000.0,1000.0,1000.0};
00046
00047 const MuonHypothesis::EF_thresholdPair ef_sa_pair[] = {
00048 MuonHypothesis::EF_thresholdPair("40GeV", ef_sa_40),
00049 MuonHypothesis::EF_thresholdPair("40GeV_barrelOnly", ef_sa_40_barrelOnly),
00050 MuonHypothesis::EF_thresholdPair("40GeV_uptoEC2", ef_sa_40_uptoEC2),
00051 MuonHypothesis::EF_thresholdPair("50GeV_barrelOnly", ef_sa_50_barrelOnly),
00052 };
00053
00054 const MuonHypothesis::EF_thresholdMap ef_sa_map(ef_sa_pair, ef_sa_pair + sizeof(ef_sa_pair)/sizeof(ef_sa_pair[0]));
00055 }
00056 #endif
00057