payloadPlotting Namespace Reference

Functions
std::vector< std::tuple< int, int, int > >	maskedBarrelRocsToBins (int layer, int ladder, int module, std::bitset< 16 > bad_rocs, bool isFlipped)
std::vector< std::tuple< int, int, int > >	maskedForwardRocsToBins (int ring, int blade, int panel, int disk, std::bitset< 16 > bad_rocs, bool isFlipped)

Function Documentation

maskedBarrelRocsToBins()

std::vector<std::tuple<int, int, int> > payloadPlotting::maskedBarrelRocsToBins	(	int	layer,
		int	ladder,
		int	module,
		std::bitset< 16 >	bad_rocs,
		bool	isFlipped
	)

Definition at line 56 of file SiPixelQualityPlotter.cc.

References funct::abs(), heavyIonCSV_trainingSettings::idx, and PVValHelper::ladder.

Referenced by SiPixelQualityPlotter::analyze().

  {
    std::vector<std::tuple<int, int, int>> rocsToMask;

    int nlad_list[4] = {6, 14, 22, 32};
    int nlad = nlad_list[layer - 1];

    int start_x = module > 0 ? ((module + 4) * 8) + 1 : ((4 - (std::abs(module))) * 8) + 1;
    int start_y = ladder > 0 ? ((ladder + nlad) * 2) + 1 : ((nlad - (std::abs(ladder))) * 2) + 1;

    int roc0_x = ((layer == 1) || (layer > 1 && module > 0)) ? start_x + 7 : start_x;
    int roc0_y = start_y - 1;

    size_t idx = 0;
    while (idx < bad_rocs.size()) {
      if (bad_rocs.test(idx)) {
        int roc_x(0), roc_y(0);

        if ((layer == 1) || (layer > 1 && module > 0)) {
          if (!isFlipped) {
            roc_x = idx < 8 ? roc0_x - idx : (start_x - 8) + idx;
            roc_y = idx < 8 ? roc0_y + 1 : roc0_y + 2;
          } else {
            roc_x = idx < 8 ? roc0_x - idx : (start_x - 8) + idx;
            roc_y = idx < 8 ? roc0_y + 2 : roc0_y + 1;
          }
        } else {
          if (!isFlipped) {
            roc_x = idx < 8 ? roc0_x + idx : (roc0_x + 7) - (idx - 8);
            roc_y = idx < 8 ? roc0_y + 1 : roc0_y + 2;
          } else {
            roc_x = idx < 8 ? roc0_x + idx : (roc0_x + 7) - (idx - 8);
            roc_y = idx < 8 ? roc0_y + 2 : roc0_y + 1;
          }
        }

        edm::LogInfo("payloadPlotting") << bad_rocs << " : (idx)= " << idx << std::endl;
        edm::LogInfo("payloadPlotting") << " layer:  " << layer << std::endl;
        edm::LogInfo("payloadPlotting") << "module: " << module << " roc_x:" << roc_x << std::endl;
        edm::LogInfo("payloadPlotting") << "ladder: " << ladder << " roc_y:" << roc_y << std::endl;
        edm::LogInfo("payloadPlotting") << "=================================================================="
                                        << std::endl;

        rocsToMask.push_back(std::make_tuple(roc_x, roc_y, idx));
      }
      ++idx;
    }
    return rocsToMask;
  }

maskedForwardRocsToBins()

std::vector<std::tuple<int, int, int> > payloadPlotting::maskedForwardRocsToBins	(	int	ring,
		int	blade,
		int	panel,
		int	disk,
		std::bitset< 16 >	bad_rocs,
		bool	isFlipped
	)

Definition at line 109 of file SiPixelQualityPlotter.cc.

References funct::abs(), heavyIonCSV_trainingSettings::idx, and relativeConstraints::ring.

Referenced by SiPixelQualityPlotter::analyze().

  {
    std::vector<std::tuple<int, int, int>> rocsToMask;
    int nybins_list[2] = {92, 140};
    int nybins = nybins_list[ring - 1];

    int start_x = disk > 0 ? ((disk + 3) * 8) + 1 : ((3 - (std::abs(disk))) * 8) + 1;
    int start_y = blade > 0 ? (nybins / 2) + (blade * 4) - (panel * 2) + 3
                            : ((nybins / 2) - (std::abs(blade) * 4) - panel * 2) + 3;

    int roc0_x = disk > 0 ? start_x + 7 : start_x;
    int roc0_y = start_y - 1;

    size_t idx = 0;
    while (idx < bad_rocs.size()) {
      if (bad_rocs.test(idx)) {
        int roc_x(0), roc_y(0);

        if (disk > 0) {
          if (!isFlipped) {
            roc_x = idx < 8 ? roc0_x - idx : (start_x - 8) + idx;
            roc_y = idx < 8 ? roc0_y + 1 : roc0_y + 2;
          } else {
            roc_x = idx < 8 ? roc0_x - idx : (start_x - 8) + idx;
            roc_y = idx < 8 ? roc0_y + 2 : roc0_y + 1;
          }
        } else {
          if (!isFlipped) {
            roc_x = idx < 8 ? roc0_x + idx : (roc0_x + 7) - (idx - 8);
            roc_y = idx < 8 ? roc0_y + 1 : roc0_y + 2;
          } else {
            roc_x = idx < 8 ? roc0_x + idx : (roc0_x + 7) - (idx - 8);
            roc_y = idx < 8 ? roc0_y + 2 : roc0_y + 1;
          }
        }

        edm::LogInfo("payloadPlotting") << bad_rocs << " : (idx)= " << idx << std::endl;
        edm::LogInfo("payloadPlotting") << " panel: " << panel << " isFlipped: " << isFlipped << std::endl;
        edm::LogInfo("payloadPlotting") << " disk:  " << disk << " roc_x:" << roc_x << std::endl;
        edm::LogInfo("payloadPlotting") << " blade: " << blade << " roc_y:" << roc_y << std::endl;
        edm::LogInfo("payloadPlotting") << "===============================" << std::endl;

        rocsToMask.push_back(std::make_tuple(roc_x, roc_y, idx));
      }
      ++idx;
    }
    return rocsToMask;
  }