Analyze.
576 std::vector<std::vector<int>> probe_coll_GE11_region;
577 std::vector<std::vector<int>> probe_coll_GE11_lay;
578 std::vector<std::vector<int>> probe_coll_GE11_chamber;
579 std::vector<std::vector<float>> probe_coll_GE11_pt;
580 std::vector<std::vector<float>> probe_coll_GE11_eta;
581 std::vector<std::vector<float>> probe_coll_GE11_ieta;
582 std::vector<std::vector<float>> probe_coll_GE11_phi;
583 std::vector<std::vector<int>> probe_coll_GE11_sta;
584 std::vector<std::vector<float>> probe_coll_GE11_dx;
587 std::vector<std::vector<int>> probe_coll_GE21_region;
588 std::vector<std::vector<int>> probe_coll_GE21_lay;
589 std::vector<std::vector<int>> probe_coll_GE21_chamber;
590 std::vector<std::vector<float>> probe_coll_GE21_pt;
591 std::vector<std::vector<float>> probe_coll_GE21_eta;
592 std::vector<std::vector<float>> probe_coll_GE21_ieta;
593 std::vector<std::vector<float>> probe_coll_GE21_phi;
594 std::vector<std::vector<int>> probe_coll_GE21_sta;
595 std::vector<std::vector<float>> probe_coll_GE21_dx;
597 std::vector<uint8_t> probe_coll_GEM_staMatch;
600 std::vector<std::vector<int>> probe_coll_ME0_region;
601 std::vector<std::vector<int>> probe_coll_ME0_roll;
602 std::vector<std::vector<int>> probe_coll_ME0_lay;
603 std::vector<std::vector<int>> probe_coll_ME0_chamber;
604 std::vector<std::vector<float>> probe_coll_ME0_pt;
605 std::vector<std::vector<float>> probe_coll_ME0_eta;
606 std::vector<std::vector<float>> probe_coll_ME0_ieta;
607 std::vector<std::vector<float>> probe_coll_ME0_phi;
608 std::vector<std::vector<int>> probe_coll_ME0_sta;
609 std::vector<std::vector<float>> probe_coll_ME0_dx;
611 std::vector<unsigned> probe_indices;
616 for (
const auto i : probe_indices) {
618 std::vector<int> probe_GE11_region;
619 std::vector<int> probe_GE11_sta;
620 std::vector<int> probe_GE11_lay;
621 std::vector<int> probe_GE11_chamber;
622 std::vector<float> probe_GE11_pt;
623 std::vector<float> probe_GE11_eta;
624 std::vector<float> probe_GE11_ieta;
625 std::vector<float> probe_GE11_phi;
626 std::vector<float> probe_GE11_dx;
628 std::vector<int> probe_GE21_region;
629 std::vector<int> probe_GE21_sta;
630 std::vector<int> probe_GE21_lay;
631 std::vector<int> probe_GE21_chamber;
632 std::vector<float> probe_GE21_pt;
633 std::vector<float> probe_GE21_eta;
634 std::vector<float> probe_GE21_ieta;
635 std::vector<float> probe_GE21_phi;
636 std::vector<float> probe_GE21_dx;
638 uint8_t GEM_stationMatching = 0;
640 std::vector<int> probe_ME0_region;
641 std::vector<int> probe_ME0_roll;
642 std::vector<int> probe_ME0_sta;
643 std::vector<int> probe_ME0_lay;
644 std::vector<int> probe_ME0_chamber;
645 std::vector<float> probe_ME0_pt;
646 std::vector<float> probe_ME0_eta;
647 std::vector<float> probe_ME0_ieta;
648 std::vector<float> probe_ME0_phi;
649 std::vector<float> probe_ME0_dx;
651 bool gem_matched =
false;
653 for (
const auto& chambMatch : (*muons).at(
i).matches()) {
655 bool hit_matched =
false;
660 GEMDetId chId(chambMatch.id.rawId());
662 const int roll = chId.
roll();
663 const int region = chId.region();
664 const int station = chId.station();
665 const int layer = chId.layer();
666 const int chamber = chId.chamber();
667 const int ieta = chId.ieta();
668 const float pt = (*muons).at(
i).pt();
669 const float eta = (*muons).at(
i).eta();
670 const float phi = (*muons).at(
i).phi();
671 GEM_stationMatching = GEM_stationMatching | (1 <<
station);
675 double smallestDx = 99999.;
676 double matched_GEMHit_x = 99999.;
678 for (
auto& gemHit : chambMatch.gemHitMatches) {
680 if (
dx < smallestDx) {
682 closest_matchedHit = gemHit;
683 matched_GEMHit_x = gemHit.
x;
689 probe_GE11_region.push_back(
region);
690 probe_GE11_sta.push_back(
station);
691 probe_GE11_lay.push_back(
layer);
692 probe_GE11_chamber.push_back(
chamber);
693 probe_GE11_ieta.push_back(
ieta);
694 probe_GE11_pt.push_back(
pt);
695 probe_GE11_eta.push_back(
eta);
696 probe_GE11_phi.push_back(
phi);
697 probe_GE11_dx.push_back(smallestDx);
701 probe_GE21_region.push_back(
region);
702 probe_GE21_sta.push_back(
station);
703 probe_GE21_lay.push_back(
layer);
704 probe_GE21_chamber.push_back(
chamber);
705 probe_GE21_ieta.push_back(
ieta);
706 probe_GE21_pt.push_back(
pt);
707 probe_GE21_eta.push_back(
eta);
708 probe_GE21_phi.push_back(
phi);
709 probe_GE21_dx.push_back(smallestDx);
714 m_histos.find(
"GEMhit_dx_GE1")->second->Fill(smallestDx);
715 m_histos.find(
"GEMhit_x_GE1")->second->Fill(matched_GEMHit_x);
716 m_histos.find(
"Cham_x_GE1")->second->Fill(chambMatch.x);
719 m_histos.find(
"GEMhit_dx_GE2")->second->Fill(smallestDx);
720 m_histos.find(
"GEMhit_x_GE2")->second->Fill(matched_GEMHit_x);
721 m_histos.find(
"Cham_x_GE2")->second->Fill(chambMatch.x);
728 double smallestDx_seg = 99999.;
730 for (
auto& seg : chambMatch.gemMatches) {
731 float dx_seg =
std::abs(chambMatch.x - seg.x);
732 if (dx_seg < smallestDx_seg) {
733 smallestDx_seg = dx_seg;
734 closest_matchedSegment = seg;
739 probe_ME0_region.push_back(
region);
740 probe_ME0_roll.push_back(roll);
741 probe_ME0_sta.push_back(
station);
742 probe_ME0_lay.push_back(
layer);
743 probe_ME0_chamber.push_back(
chamber);
744 probe_ME0_ieta.push_back(
ieta);
745 probe_ME0_pt.push_back(
pt);
746 probe_ME0_eta.push_back(
eta);
747 probe_ME0_phi.push_back(
phi);
748 probe_ME0_dx.push_back(smallestDx_seg);
751 m_histos.find(
"GEMseg_dx_ME0")->second->Fill(smallestDx_seg);
752 m_histos.find(
"GEMseg_x_ME0")->second->Fill(closest_matchedSegment.
x);
753 m_histos.find(
"Cham_x_ME0")->second->Fill(chambMatch.x);
763 m_histos.find(
"probeEta")->second->Fill((*muons).at(
i).eta());
764 m_histos.find(
"probePhi")->second->Fill((*muons).at(
i).phi());
765 m_histos.find(
"probeNumberOfMatchedStations")->second->Fill((*muons).at(
i).numberOfMatchedStations());
766 m_histos.find(
"probePt")->second->Fill((*muons).at(
i).pt());
775 probe_coll_GE11_region.push_back(probe_GE11_region);
776 probe_coll_GE11_sta.push_back(probe_GE11_sta);
777 probe_coll_GE11_lay.push_back(probe_GE11_lay);
778 probe_coll_GE11_chamber.push_back(probe_GE11_chamber);
779 probe_coll_GE11_ieta.push_back(probe_GE11_ieta);
780 probe_coll_GE11_pt.push_back(probe_GE11_pt);
781 probe_coll_GE11_eta.push_back(probe_GE11_eta);
782 probe_coll_GE11_phi.push_back(probe_GE11_phi);
783 probe_coll_GE11_dx.push_back(probe_GE11_dx);
785 probe_coll_GEM_staMatch.push_back(GEM_stationMatching);
788 probe_coll_GE21_region.push_back(probe_GE21_region);
789 probe_coll_GE21_sta.push_back(probe_GE21_sta);
790 probe_coll_GE21_lay.push_back(probe_GE21_lay);
791 probe_coll_GE21_chamber.push_back(probe_GE21_chamber);
792 probe_coll_GE21_ieta.push_back(probe_GE21_ieta);
793 probe_coll_GE21_pt.push_back(probe_GE21_pt);
794 probe_coll_GE21_eta.push_back(probe_GE21_eta);
795 probe_coll_GE21_phi.push_back(probe_GE21_phi);
796 probe_coll_GE21_dx.push_back(probe_GE21_dx);
799 probe_coll_ME0_region.push_back(probe_ME0_region);
800 probe_coll_ME0_roll.push_back(probe_ME0_roll);
801 probe_coll_ME0_sta.push_back(probe_ME0_sta);
802 probe_coll_ME0_lay.push_back(probe_ME0_lay);
803 probe_coll_ME0_chamber.push_back(probe_ME0_chamber);
804 probe_coll_ME0_ieta.push_back(probe_ME0_ieta);
805 probe_coll_ME0_pt.push_back(probe_ME0_pt);
806 probe_coll_ME0_eta.push_back(probe_ME0_eta);
807 probe_coll_ME0_phi.push_back(probe_ME0_phi);
808 probe_coll_ME0_dx.push_back(probe_ME0_dx);
813 for (
unsigned i = 0;
i < probe_indices.size(); ++
i) {
817 unsigned nME0_matches = probe_coll_ME0_region.at(
i).size();
818 for (
unsigned j = 0;
j < nME0_matches; ++
j) {
820 int ME0_region = probe_coll_ME0_region.at(
i).at(
j);
824 int ME0_chamber = probe_coll_ME0_chamber.at(
i).at(
j);
826 float ME0_dx = probe_coll_ME0_dx.at(
i).at(
j);
831 m_histos.find(
"ME0_nPassingProbe_chamber_1D")->second->Fill(ME0_chamber);
833 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(4, ME0_chamber);
834 else if (ME0_region > 0)
835 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(5, ME0_chamber);
837 m_histos.find(
"ME0_nFailingProbe_chamber_1D")->second->Fill(ME0_chamber);
839 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(4, ME0_chamber);
840 else if (ME0_region > 0)
841 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(5, ME0_chamber);
847 unsigned nGE11_matches = probe_coll_GE11_region.at(
i).size();
848 for (
unsigned j = 0;
j < nGE11_matches; ++
j) {
850 int GEM_region = probe_coll_GE11_region.at(
i).at(
j);
851 int GEM_sta = probe_coll_GE11_sta.at(
i).at(
j);
852 int GEM_lay = probe_coll_GE11_lay.at(
i).at(
j);
853 int GEM_chamber = probe_coll_GE11_chamber.at(
i).at(
j);
854 int GEM_ieta = probe_coll_GE11_ieta.at(
i).at(
j);
855 float GEM_pt = probe_coll_GE11_pt.at(
i).at(
j);
856 float GEM_dx = probe_coll_GE11_dx.at(
i).at(
j);
857 float GEM_eta = probe_coll_GE11_eta.at(
i).at(
j);
858 float GEM_phi = probe_coll_GE11_phi.at(
i).at(
j);
861 m_histos.find(
"GE11_nPassingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
862 m_histos.find(
"GE11_nPassingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
863 m_histos.find(
"GE11_nPassingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
864 m_histos.find(
"GE11_nPassingProbe_allCh_1D")->second->Fill(GEM_region);
865 m_histos.find(
"GE11_nPassingProbe_chamber_1D")->second->Fill(GEM_chamber);
866 if (GEM_region < 0) {
868 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(2, GEM_chamber);
869 else if (GEM_lay == 1)
870 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(3, GEM_chamber);
872 if (GEM_region > 0) {
874 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(6, GEM_chamber);
875 else if (GEM_lay == 2)
876 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(7, GEM_chamber);
878 if (GEM_region == -1) {
879 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1_NoL")->second->Fill(0, GEM_chamber);
880 }
else if (GEM_region == 1) {
881 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1_NoL")->second->Fill(1, GEM_chamber);
884 if (GEM_region == 1 && GEM_lay == 1) {
885 m_histos.find(
"GEM_nPassingProbe_chamber_p1_1D")->second->Fill(GEM_chamber);
886 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1")->second->Fill(2, GEM_chamber);
887 m_histos.find(
"GEM_nPassingProbe_pt_p1_1D")->second->Fill(GEM_pt);
888 m_histos.find(
"GEM_nPassingProbe_eta_p1_1D")->second->Fill(
abs(GEM_eta));
889 m_histos.find(
"GEM_nPassingProbe_phi_p1_1D")->second->Fill(GEM_phi);
890 }
else if (GEM_region == 1 && GEM_lay == 2) {
891 m_histos.find(
"GEM_nPassingProbe_chamber_p2_1D")->second->Fill(GEM_chamber);
892 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1")->second->Fill(3, GEM_chamber);
893 m_histos.find(
"GEM_nPassingProbe_pt_p2_1D")->second->Fill(GEM_pt);
894 m_histos.find(
"GEM_nPassingProbe_eta_p2_1D")->second->Fill(
abs(GEM_eta));
895 m_histos.find(
"GEM_nPassingProbe_phi_p2_1D")->second->Fill(GEM_phi);
896 }
else if (GEM_region == -1 && GEM_lay == 1) {
897 m_histos.find(
"GEM_nPassingProbe_chamber_n1_1D")->second->Fill(GEM_chamber);
898 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1")->second->Fill(1, GEM_chamber);
899 m_histos.find(
"GEM_nPassingProbe_pt_n1_1D")->second->Fill(GEM_pt);
900 m_histos.find(
"GEM_nPassingProbe_eta_n1_1D")->second->Fill(
abs(GEM_eta));
901 m_histos.find(
"GEM_nPassingProbe_phi_n1_1D")->second->Fill(GEM_phi);
902 }
else if (GEM_region == -1 && GEM_lay == 2) {
903 m_histos.find(
"GEM_nPassingProbe_chamber_n2_1D")->second->Fill(GEM_chamber);
904 m_histos.find(
"GEM_nPassingProbe_Ch_region_GE1")->second->Fill(0, GEM_chamber);
905 m_histos.find(
"GEM_nPassingProbe_pt_n2_1D")->second->Fill(GEM_pt);
906 m_histos.find(
"GEM_nPassingProbe_eta_n2_1D")->second->Fill(
abs(GEM_eta));
907 m_histos.find(
"GEM_nPassingProbe_phi_n2_1D")->second->Fill(GEM_phi);
909 m_histos.find(
"GEM_nPassingProbe_pt_1D")->second->Fill(GEM_pt);
910 m_histos.find(
"GEM_nPassingProbe_eta_1D")->second->Fill(
abs(GEM_eta));
911 m_histos.find(
"GEM_nPassingProbe_phi_1D")->second->Fill(GEM_phi);
913 m_histos.find(
"GE11_nFailingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
914 m_histos.find(
"GE11_nFailingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
915 m_histos.find(
"GE11_nFailingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
916 m_histos.find(
"GE11_nFailingProbe_allCh_1D")->second->Fill(GEM_region);
917 m_histos.find(
"GE11_nFailingProbe_chamber_1D")->second->Fill(GEM_chamber);
918 if (GEM_region < 0) {
919 if (GEM_sta == 2 and GEM_lay == 2)
920 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
921 else if (GEM_sta == 2 and GEM_lay == 1)
922 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
923 else if (GEM_sta == 1 and GEM_lay == 2)
924 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(2, GEM_chamber);
925 else if (GEM_sta == 1 and GEM_lay == 1)
926 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(3, GEM_chamber);
928 if (GEM_region > 0) {
929 if (GEM_sta == 1 and GEM_lay == 1)
930 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(6, GEM_chamber);
931 else if (GEM_sta == 1 and GEM_lay == 2)
932 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(7, GEM_chamber);
933 else if (GEM_sta == 2 and GEM_lay == 1)
934 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
935 else if (GEM_sta == 2 and GEM_lay == 2)
936 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
938 if (GEM_region == -1) {
939 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1_NoL")->second->Fill(0, GEM_chamber);
940 }
else if (GEM_region == 1) {
941 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1_NoL")->second->Fill(1, GEM_chamber);
944 if (GEM_region == 1 && GEM_lay == 1) {
945 m_histos.find(
"GEM_nFailingProbe_chamber_p1_1D")->second->Fill(GEM_chamber);
946 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1")->second->Fill(2, GEM_chamber);
947 m_histos.find(
"GEM_nFailingProbe_pt_p1_1D")->second->Fill(GEM_pt);
948 m_histos.find(
"GEM_nFailingProbe_eta_p1_1D")->second->Fill(
abs(GEM_eta));
949 m_histos.find(
"GEM_nFailingProbe_phi_p1_1D")->second->Fill(GEM_phi);
950 }
else if (GEM_region == 1 && GEM_lay == 2) {
951 m_histos.find(
"GEM_nFailingProbe_chamber_p2_1D")->second->Fill(GEM_chamber);
952 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1")->second->Fill(3, GEM_chamber);
953 m_histos.find(
"GEM_nFailingProbe_pt_p2_1D")->second->Fill(GEM_pt);
954 m_histos.find(
"GEM_nFailingProbe_eta_p2_1D")->second->Fill(
abs(GEM_eta));
955 m_histos.find(
"GEM_nFailingProbe_phi_p2_1D")->second->Fill(GEM_phi);
956 }
else if (GEM_region == -1 && GEM_lay == 1) {
957 m_histos.find(
"GEM_nFailingProbe_chamber_n1_1D")->second->Fill(GEM_chamber);
958 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1")->second->Fill(1, GEM_chamber);
959 m_histos.find(
"GEM_nFailingProbe_pt_n1_1D")->second->Fill(GEM_pt);
960 m_histos.find(
"GEM_nFailingProbe_eta_n1_1D")->second->Fill(
abs(GEM_eta));
961 m_histos.find(
"GEM_nFailingProbe_phi_n1_1D")->second->Fill(GEM_phi);
962 }
else if (GEM_region == -1 && GEM_lay == 2) {
963 m_histos.find(
"GEM_nFailingProbe_chamber_n2_1D")->second->Fill(GEM_chamber);
964 m_histos.find(
"GEM_nFailingProbe_Ch_region_GE1")->second->Fill(0, GEM_chamber);
965 m_histos.find(
"GEM_nFailingProbe_pt_n2_1D")->second->Fill(GEM_pt);
966 m_histos.find(
"GEM_nFailingProbe_eta_n2_1D")->second->Fill(
abs(GEM_eta));
967 m_histos.find(
"GEM_nFailingProbe_phi_n2_1D")->second->Fill(GEM_phi);
969 m_histos.find(
"GEM_nFailingProbe_pt_1D")->second->Fill(GEM_pt);
970 m_histos.find(
"GEM_nFailingProbe_eta_1D")->second->Fill(
abs(GEM_eta));
971 m_histos.find(
"GEM_nFailingProbe_phi_1D")->second->Fill(GEM_phi);
976 unsigned nGE21_matches = probe_coll_GE21_region.at(
i).size();
977 for (
unsigned j = 0;
j < nGE21_matches; ++
j) {
979 int GEM_region = probe_coll_GE21_region.at(
i).at(
j);
980 int GEM_lay = probe_coll_GE21_lay.at(
i).at(
j);
981 int GEM_chamber = probe_coll_GE21_chamber.at(
i).at(
j);
982 float GEM_ieta = probe_coll_GE21_ieta.at(
i).at(
j);
983 float GEM_dx = probe_coll_GE21_dx.at(
i).at(
j);
984 float GEM_phi = probe_coll_GE21_phi.at(
i).at(
j);
988 m_histos.find(
"GE21_nPassingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
989 m_histos.find(
"GE21_nPassingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
990 m_histos.find(
"GE21_nPassingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
991 m_histos.find(
"GE21_nPassingProbe_allCh_1D")->second->Fill(GEM_region);
992 m_histos.find(
"GE21_nPassingProbe_chamber_1D")->second->Fill(GEM_chamber);
993 if (GEM_region < 0) {
995 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
996 else if (GEM_lay == 1)
997 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
999 if (GEM_region > 0) {
1001 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
1002 else if (GEM_lay == 2)
1003 m_histos.find(
"GEM_nPassingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
1006 m_histos.find(
"GE21_nFailingProbe_Ch_region")->second->Fill(GEM_region, GEM_chamber);
1007 m_histos.find(
"GE21_nFailingProbe_Ch_ieta")->second->Fill(GEM_ieta, GEM_chamber);
1008 m_histos.find(
"GE21_nFailingProbe_Ch_phi")->second->Fill(GEM_phi, GEM_chamber);
1009 m_histos.find(
"GE21_nFailingProbe_allCh_1D")->second->Fill(GEM_region);
1010 m_histos.find(
"GE21_nFailingProbe_chamber_1D")->second->Fill(GEM_chamber);
1011 if (GEM_region < 0) {
1013 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(0, GEM_chamber);
1014 else if (GEM_lay == 1)
1015 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(1, GEM_chamber);
1017 if (GEM_region > 0) {
1019 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(8, GEM_chamber);
1020 else if (GEM_lay == 2)
1021 m_histos.find(
"GEM_nFailingProbe_Ch_region_layer_phase2")->second->Fill(9, GEM_chamber);
const edm::EDGetTokenT< reco::MuonCollection > m_muToken
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::vector< std::vector< unsigned > > m_probeIndices
const bool m_detailedAnalysis
Abs< T >::type abs(const T &t)
std::map< std::string, MonitorElement * > m_histos
void analyze(const edm::Event &event, const edm::EventSetup &context) override
Analyze.
constexpr int roll() const