SectorProcessorLUT Class Reference

#include <SectorProcessorLUT.h>

Public Member Functions
uint32_t	get_cppf_lut_id (int rpc_region, int rpc_sector, int rpc_station, int rpc_ring, int rpc_subsector, int rpc_roll) const
uint32_t	get_cppf_ph_lut (int rpc_region, int rpc_sector, int rpc_station, int rpc_ring, int rpc_subsector, int rpc_roll, int halfstrip, bool is_neighbor) const
uint32_t	get_cppf_th_lut (int rpc_region, int rpc_sector, int rpc_station, int rpc_ring, int rpc_subsector, int rpc_roll) const
uint32_t	get_ph_disp (int fw_endcap, int fw_sector, int pc_lut_id) const
uint32_t	get_ph_init (int fw_endcap, int fw_sector, int pc_lut_id) const
uint32_t	get_ph_init_hard (int fw_station, int fw_cscid) const
uint32_t	get_ph_patt_corr (int pattern) const
uint32_t	get_ph_patt_corr_sign (int pattern) const
uint32_t	get_ph_zone_offset (int pc_station, int pc_chamber) const
uint32_t	get_th_corr_lut (int fw_endcap, int fw_sector, int pc_lut_id, int pc_wire_strip_id) const
uint32_t	get_th_disp (int fw_endcap, int fw_sector, int pc_lut_id) const
uint32_t	get_th_init (int fw_endcap, int fw_sector, int pc_lut_id) const
uint32_t	get_th_lut (int fw_endcap, int fw_sector, int pc_lut_id, int pc_wire_id) const
void	read (bool pc_lut_data, int pc_lut_version)
	SectorProcessorLUT ()
	~SectorProcessorLUT ()

Private Member Functions
void	read_cppf_file (const std::string &filename, std::vector< uint32_t > &vec1, std::vector< uint32_t > &vec2, bool local)
void	read_file (const std::string &filename, std::vector< uint32_t > &vec)

Private Attributes
std::vector< uint32_t >	cppf_ph_lut_
std::vector< uint32_t >	cppf_th_lut_
std::vector< uint32_t >	ph_disp_neighbor_
std::vector< uint32_t >	ph_init_hard_
std::vector< uint32_t >	ph_init_neighbor_
std::vector< uint32_t >	ph_patt_corr_
std::vector< uint32_t >	ph_patt_corr_sign_
std::vector< uint32_t >	ph_zone_offset_
std::vector< uint32_t >	th_corr_lut_neighbor_
std::vector< uint32_t >	th_disp_neighbor_
std::vector< uint32_t >	th_init_neighbor_
std::vector< uint32_t >	th_lut_neighbor_
int	version_

Detailed Description

Definition at line 8 of file SectorProcessorLUT.h.

Constructor & Destructor Documentation

SectorProcessorLUT::SectorProcessorLUT ( )

explicit

Definition at line 10 of file SectorProcessorLUT.cc.

: version_(0xFFFFFFFF) {}

SectorProcessorLUT::~SectorProcessorLUT

(

)

Definition at line 12 of file SectorProcessorLUT.cc.

{}

Member Function Documentation

uint32_t SectorProcessorLUT::get_cppf_lut_id	(	int	rpc_region,
		int	rpc_sector,
		int	rpc_station,
		int	rpc_ring,
		int	rpc_subsector,
		int	rpc_roll
	)		const

Definition at line 298 of file SectorProcessorLUT.cc.

Referenced by get_cppf_ph_lut(), get_cppf_th_lut(), and read_cppf_file().

                                                                                                          {
  uint32_t iendcap = (rpc_region == -1) ? 1 : 0;
  uint32_t isector = (rpc_sector - 1);
  uint32_t istationring = (rpc_station >= 3) ? ((rpc_station - 3) * 2 + (rpc_ring - 2) + 2) : (rpc_station - 1);
  uint32_t isubsector = (rpc_subsector - 1);
  uint32_t iroll = (rpc_roll - 1);
  return ((((iendcap * 6 + isector) * 6 + istationring) * 6 + isubsector) * 3 + iroll);
}

uint32_t SectorProcessorLUT::get_cppf_ph_lut	(	int	rpc_region,
		int	rpc_sector,
		int	rpc_station,
		int	rpc_ring,
		int	rpc_subsector,
		int	rpc_roll,
		int	halfstrip,
		bool	is_neighbor
	)		const

Definition at line 308 of file SectorProcessorLUT.cc.

References cppf_ph_lut_, and get_cppf_lut_id().

Referenced by PrimitiveConversion::convert_rpc_details().

                                                                     {
  const uint32_t th_index = get_cppf_lut_id(rpc_region, rpc_sector, rpc_station, rpc_ring, rpc_subsector, rpc_roll);
  const uint32_t ph_index = (th_index * 64) + (halfstrip - 1);
  uint32_t ph = 0;

  if (ph_index < cppf_ph_lut_.size()) {
    ph = cppf_ph_lut_.at(ph_index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. rpc_region: " << rpc_region
                         << ", rpc_sector: " << rpc_sector << ", rpc_station: " << rpc_station
                         << ", rpc_ring: " << rpc_ring << ", rpc_subsector: " << rpc_subsector
                         << ", rpc_roll: " << rpc_roll << ", halfstrip: " << halfstrip
                         << ", is_neighbor: " << is_neighbor;
  }

  if (!is_neighbor && rpc_subsector == 2)
    ph += 900;
  return ph;
}

uint32_t SectorProcessorLUT::get_cppf_th_lut	(	int	rpc_region,
		int	rpc_sector,
		int	rpc_station,
		int	rpc_ring,
		int	rpc_subsector,
		int	rpc_roll
	)		const

Definition at line 335 of file SectorProcessorLUT.cc.

References cppf_th_lut_, and get_cppf_lut_id().

Referenced by PrimitiveConversion::convert_rpc_details().

                                                                                                          {
  const uint32_t th_index = get_cppf_lut_id(rpc_region, rpc_sector, rpc_station, rpc_ring, rpc_subsector, rpc_roll);
  uint32_t th = 0;

  if (th_index < cppf_th_lut_.size()) {
    th = cppf_th_lut_.at(th_index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. rpc_region: " << rpc_region
                         << ", rpc_sector: " << rpc_sector << ", rpc_station: " << rpc_station
                         << ", rpc_ring: " << rpc_ring << ", rpc_subsector: " << rpc_subsector
                         << ", rpc_roll: " << rpc_roll;
  }
  return th;
}

uint32_t SectorProcessorLUT::get_ph_disp	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id
	)		const

Definition at line 153 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_disp_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  const uint32_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  uint32_t entry = 0;

  if (index < ph_disp_neighbor_.size()) {
    entry = ph_disp_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_ph_init	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id
	)		const

Definition at line 140 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_init_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  const uint32_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  uint32_t entry = 0;

  if (index < ph_init_neighbor_.size()) {
    entry = ph_init_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_ph_init_hard	(	int	fw_station,
		int	fw_cscid
	)		const

Definition at line 285 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_init_hard_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                {
  const uint32_t index = fw_station * 16 + fw_cscid;
  uint32_t entry = 0;

  if (index < ph_init_hard_.size()) {
    entry = ph_init_hard_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_station: " << fw_station
                         << ", fw_cscid: " << fw_cscid;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_ph_patt_corr

(

int

pattern

)

const

Definition at line 248 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_patt_corr_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                               {
  const uint32_t index = pattern;
  uint32_t entry = 0;

  if (index < ph_patt_corr_.size()) {
    entry = ph_patt_corr_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. pattern: " << pattern;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_ph_patt_corr_sign

(

int

pattern

)

const

Definition at line 260 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_patt_corr_sign_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                    {
  const uint32_t index = pattern;
  uint32_t entry = 0;

  if (index < ph_patt_corr_sign_.size()) {
    entry = ph_patt_corr_sign_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. pattern: " << pattern;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_ph_zone_offset	(	int	pc_station,
		int	pc_chamber
	)		const

Definition at line 272 of file SectorProcessorLUT.cc.

References mps_splice::entry, and ph_zone_offset_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                    {
  const uint32_t index = pc_station * 9 + pc_chamber;
  uint32_t entry = 0;

  if (index < ph_zone_offset_.size()) {
    entry = ph_zone_offset_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. pc_station: " << pc_station
                         << ", pc_chamber: " << pc_chamber;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_th_corr_lut	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id,
		int	pc_wire_strip_id
	)		const

Definition at line 215 of file SectorProcessorLUT.cc.

References mps_splice::entry, and th_corr_lut_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                                                    {
  int pc_lut_id2 = pc_lut_id;

  // Make ME1/1a the same as ME1/1b
  if ((9 <= pc_lut_id2 && pc_lut_id2 < 12) || (25 <= pc_lut_id2 && pc_lut_id2 < 28))
    pc_lut_id2 -= 9;
  // Make ME1/1a neighbor the same as ME1/1b
  if (pc_lut_id2 == 15)
    pc_lut_id2 -= 3;

  if (pc_lut_id2 <= 3) {
    pc_lut_id2 -= 0;
  } else if (pc_lut_id2 == 12) {
    pc_lut_id2 -= 9;
  } else if (16 <= pc_lut_id2 && pc_lut_id2 < 19) {
    pc_lut_id2 -= 12;
  } else {
    edm::LogError("L1T") << "get_th_corr_lut(): out of range pc_lut_id: " << pc_lut_id;
  }

  const uint32_t index = ((fw_endcap * 6 + fw_sector) * 7 + pc_lut_id2) * 128 + pc_wire_strip_id;
  uint32_t entry = 0;

  if (index < th_corr_lut_neighbor_.size()) {
    entry = th_corr_lut_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id
                         << ", pc_wire_strip_id: " << pc_wire_strip_id;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_th_disp	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id
	)		const

Definition at line 179 of file SectorProcessorLUT.cc.

References mps_splice::entry, and th_disp_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  const uint32_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  uint32_t entry = 0;

  if (index < th_disp_neighbor_.size()) {
    entry = th_disp_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_th_init	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id
	)		const

Definition at line 166 of file SectorProcessorLUT.cc.

References mps_splice::entry, and th_init_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  const uint32_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  uint32_t entry = 0;

  if (index < th_init_neighbor_.size()) {
    entry = th_init_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id;
  }
  return entry;
}

uint32_t SectorProcessorLUT::get_th_lut	(	int	fw_endcap,
		int	fw_sector,
		int	pc_lut_id,
		int	pc_wire_id
	)		const

Definition at line 192 of file SectorProcessorLUT.cc.

References mps_splice::entry, and th_lut_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                                         {
  int pc_lut_id2 = pc_lut_id;

  // Make ME1/1a the same as ME1/1b
  if ((9 <= pc_lut_id2 && pc_lut_id2 < 12) || (25 <= pc_lut_id2 && pc_lut_id2 < 28))
    pc_lut_id2 -= 9;
  // Make ME1/1a neighbor the same as ME1/1b
  if (pc_lut_id2 == 15)
    pc_lut_id2 -= 3;

  const uint32_t index = ((fw_endcap * 6 + fw_sector) * 61 + pc_lut_id2) * 128 + pc_wire_id;
  uint32_t entry = 0;

  if (index < th_lut_neighbor_.size()) {
    entry = th_lut_neighbor_.at(index);
  } else {
    edm::LogError("L1T") << "Could not retrieve entry from LUT. fw_endcap: " << fw_endcap
                         << ", fw_sector: " << fw_sector << ", pc_lut_id: " << pc_lut_id
                         << ", pc_wire_id: " << pc_wire_id;
  }
  return entry;
}

void SectorProcessorLUT::read	(	bool	pc_lut_data,
		int	pc_lut_version
	)

Definition at line 14 of file SectorProcessorLUT.cc.

References cppf_ph_lut_, cppf_th_lut_, Exception, ph_disp_neighbor_, ph_init_hard_, ph_init_neighbor_, ph_patt_corr_, ph_patt_corr_sign_, ph_zone_offset_, read_cppf_file(), read_file(), AlCaHLTBitMon_QueryRunRegistry::string, th_corr_lut_neighbor_, th_disp_neighbor_, th_init_neighbor_, th_lut_neighbor_, and version_.

Referenced by edmIntegrityCheck.PublishToFileSystem::get(), and EMTFSetup::reload().

                                                                  {
  if (version_ == pc_lut_version)
    return;

  edm::LogInfo("L1T") << "EMTF using pc_lut_ver: " << pc_lut_version << ", configured for "
                      << (pc_lut_data ? "data" : "MC");

  std::string coord_lut_dir = "";
  if (pc_lut_version == 0)
    coord_lut_dir = "ph_lut_v1";  // All year 2016
  else if (pc_lut_version == 1)
    coord_lut_dir = "ph_lut_v2";  // Beginning of 2017, improved alignment from ideal CMS geometry (MC)
  else if (pc_lut_version == 2 && pc_lut_data)
    coord_lut_dir = "ph_lut_v3_data";  // Update in September 2017 from ReReco alignment, data only
  else if (pc_lut_version == 2)
    coord_lut_dir = "ph_lut_v2";  // MC still uses ideal CMS aligment
  else if (pc_lut_version == -1 && pc_lut_data)
    coord_lut_dir = "ph_lut_v3_data";  // September 2017 data LCT alignment, but use local CPPF LUTs for RPC
  else if (pc_lut_version == -1)
    coord_lut_dir = "ph_lut_v2";  // MC using ideal CMS LCT alignment, but use local CPPF LUTs for RPC
  else
    throw cms::Exception("L1TMuonEndCap")
        << "Trying to use EMTF pc_lut_version = " << pc_lut_version << ", does not exist!";
  // Will catch user trying to run with Global Tag settings on 2016 data, rather than fakeEmtfParams. - AWB 08.06.17

  std::string coord_lut_path = "L1Trigger/L1TMuon/data/emtf_luts/" + coord_lut_dir + "/";

  read_file(coord_lut_path + "ph_init_neighbor.txt", ph_init_neighbor_);
  read_file(coord_lut_path + "ph_disp_neighbor.txt", ph_disp_neighbor_);
  read_file(coord_lut_path + "th_init_neighbor.txt", th_init_neighbor_);
  read_file(coord_lut_path + "th_disp_neighbor.txt", th_disp_neighbor_);
  read_file(coord_lut_path + "th_lut_neighbor.txt", th_lut_neighbor_);
  read_file(coord_lut_path + "th_corr_lut_neighbor.txt", th_corr_lut_neighbor_);

  std::string cppf_coord_lut_path = "L1Trigger/L1TMuon/data/cppf/";  // Coordinate LUTs actually used by CPPF
  bool use_local_cppf_files = (pc_lut_version == -1);
  if (use_local_cppf_files) {  // More accurate coordinate transformation LUTs from Jia Fu
    cppf_coord_lut_path = "L1Trigger/L1TMuon/data/cppf_luts/angleScale_v1/";
  }

  read_cppf_file(cppf_coord_lut_path,
                 cppf_ph_lut_,
                 cppf_th_lut_,
                 use_local_cppf_files);  // cppf filenames are hardcoded in the function

  if (ph_init_neighbor_.size() != 2 * 6 * 61) {  // [endcap_2][sector_6][chamber_61]
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_init_neighbor_ to get " << 2 * 6 * 61 << " values, "
                                          << "got " << ph_init_neighbor_.size() << " values.";
  }

  if (ph_disp_neighbor_.size() != 2 * 6 * 61) {  // [endcap_2][sector_6][chamber_61]
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_disp_neighbor_ to get " << 2 * 6 * 61 << " values, "
                                          << "got " << ph_disp_neighbor_.size() << " values.";
  }

  if (th_init_neighbor_.size() != 2 * 6 * 61) {  // [endcap_2][sector_6][chamber_61]
    throw cms::Exception("L1TMuonEndCap") << "Expected th_init_neighbor_ to get " << 2 * 6 * 61 << " values, "
                                          << "got " << th_init_neighbor_.size() << " values.";
  }

  if (th_disp_neighbor_.size() != 2 * 6 * 61) {  // [endcap_2][sector_6][chamber_61]
    throw cms::Exception("L1TMuonEndCap") << "Expected th_disp_neighbor_ to get " << 2 * 6 * 61 << " values, "
                                          << "got " << th_disp_neighbor_.size() << " values.";
  }

  if (th_lut_neighbor_.size() != 2 * 6 * 61 * 128) {  // [endcap_2][sector_6][chamber_61][wire_128]
    throw cms::Exception("L1TMuonEndCap") << "Expected th_lut_neighbor_ to get " << 2 * 6 * 61 * 128 << " values, "
                                          << "got " << th_lut_neighbor_.size() << " values.";
  }

  if (th_corr_lut_neighbor_.size() != 2 * 6 * 7 * 128) {  // [endcap_2][sector_6][chamber_61][strip_wire_128]
    throw cms::Exception("L1TMuonEndCap") << "Expected th_corr_lut_neighbor_ to get " << 2 * 6 * 7 * 128 << " values, "
                                          << "got " << th_corr_lut_neighbor_.size() << " values.";
  }

  if (cppf_ph_lut_.size() !=
      2 * 6 * 6 * 6 * 3 *
          64) {  // [endcap_2][rpc_sector_6][rpc_station_ring_6][rpc_subsector_6][rpc_roll_3][rpc_halfstrip_64]
    throw cms::Exception("L1TMuonEndCap") << "Expected cppf_ph_lut_ to get " << 2 * 6 * 6 * 6 * 3 * 64 << " values, "
                                          << "got " << cppf_ph_lut_.size() << " values.";
  }

  if (cppf_th_lut_.size() !=
      2 * 6 * 6 * 6 * 3) {  // [endcap_2][rpc_sector_6][rpc_station_ring_6][rpc_subsector_6][rpc_roll_3]
    throw cms::Exception("L1TMuonEndCap") << "Expected cppf_th_lut_ to get " << 2 * 6 * 6 * 6 * 3 << " values, "
                                          << "got " << cppf_th_lut_.size() << " values.";
  }

  // clct pattern convertion array from CMSSW
  //{0.0, 0.0, -0.60,  0.60, -0.64,  0.64, -0.23,  0.23, -0.21,  0.21, 0.0}
  // 0    0    -5      +5    -5      +5    -2      +2    -2      +2    0
  ph_patt_corr_ = {0, 0, 5, 5, 5, 5, 2, 2, 2, 2, 0};
  if (ph_patt_corr_.size() != 11) {
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_patt_corr_ to get " << 11 << " values, "
                                          << "got " << ph_patt_corr_.size() << " values.";
  }

  ph_patt_corr_sign_ = {0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0};
  if (ph_patt_corr_sign_.size() != 11) {
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_patt_corr_sign_ to get " << 11 << " values, "
                                          << "got " << ph_patt_corr_sign_.size() << " values.";
  }

  ph_zone_offset_ = {39, 57, 76,  39, 58, 76, 41, 60,  79,  95, 114, 132, 95, 114, 133, 98, 116, 135,
                     38, 76, 113, 39, 58, 76, 95, 114, 132, 38, 76,  113, 39, 58,  76,  95, 114, 132,
                     38, 76, 113, 38, 57, 76, 95, 113, 132, 21, 21,  23,  1,  21,  1,   21, 1,   20};
  if (ph_zone_offset_.size() != 6 * 9) {
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_zone_offset_ to get " << 6 * 9 << " values, "
                                          << "got " << ph_zone_offset_.size() << " values.";
  }

  // start phi of each chamber in reduced precision, for zone hits,
  // with negative offset to allow for possible chamber movement
  ph_init_hard_ = {39,  57,  76, 39,  58,  76, 41,  60,  79, 39, 57,  76,  21,  21, 23,  21,  95,  114, 132, 95,
                   114, 133, 98, 116, 135, 95, 114, 132, 0,  0,  0,   0,   38,  76, 113, 39,  58,  76,  95,  114,
                   132, 1,   21, 0,   0,   0,  0,   0,   38, 76, 113, 39,  58,  76, 95,  114, 132, 1,   21,  0,
                   0,   0,   0,  0,   38,  76, 113, 38,  57, 76, 95,  113, 132, 1,  20,  0,   0,   0,   0,   0};
  if (ph_init_hard_.size() != 5 * 16) {
    throw cms::Exception("L1TMuonEndCap") << "Expected ph_init_hard_ to get " << 5 * 16 << " values, "
                                          << "got " << ph_init_hard_.size() << " values.";
  }

  version_ = pc_lut_version;
  return;
}

void SectorProcessorLUT::read_cppf_file ( const std::string &  filename, std::vector< uint32_t > &  vec1, std::vector< uint32_t > &  vec2, bool  local  )

private

Definition at line 365 of file SectorProcessorLUT.cc.

References Exception, contentValuesFiles::fullPath, get_cppf_lut_id(), mps_fire::i, gpuClustering::id, and timingPdfMaker::infile.

Referenced by read().

                                                    {
  auto get_rpc_region = [](uint32_t id) { return (static_cast<int>((id >> 0) & 0X3) + (-1)); };
  auto get_rpc_sector = [](uint32_t id) { return (static_cast<int>((id >> 7) & 0XF) + (1)); };
  auto get_rpc_ring = [](uint32_t id) { return (static_cast<int>((id >> 2) & 0X7) + (1)); };
  auto get_rpc_station = [](uint32_t id) { return (static_cast<int>((id >> 5) & 0X3) + (1)); };
  auto get_rpc_subsector = [](uint32_t id) { return (static_cast<int>((id >> 12) & 0X7) + (1)); };
  auto get_rpc_roll = [](uint32_t id) { return (static_cast<int>((id >> 15) & 0X7) + (0)); };

  std::vector<std::string> cppf_filenames = {
      "angleScale_RPC_CPPFp1.txt",
      "angleScale_RPC_CPPFp2.txt",
      "angleScale_RPC_CPPFp3.txt",
      "angleScale_RPC_CPPFp4.txt",
      "angleScale_RPC_CPPFn1.txt",
      "angleScale_RPC_CPPFn2.txt",
      "angleScale_RPC_CPPFn3.txt",
      "angleScale_RPC_CPPFn4.txt",
  };

  vec1.clear();
  vec2.clear();
  vec1.resize(2 * 6 * 6 * 6 * 3 * 64, 0);
  vec2.resize(2 * 6 * 6 * 6 * 3, 0);

  for (size_t i = 0; i < cppf_filenames.size(); ++i) {
    std::ifstream infile;
    infile.open(edm::FileInPath(filename + cppf_filenames.at(i)).fullPath().c_str());

    // std::cout << "\n\nOpening CPPF LUT file " << cppf_filenames.at(i) << std::endl;

    int buf1, buf2, buf3, buf4, buf5, buf6;
    // Special variables for transforming centrally-provided CPPF LUTs
    int buf1_prev = 0, buf2_prev = 0, halfstrip_prev = 0;  // Values from previous line in file
    int line_num = 0;                                      // Line number in file
    int count_dir = -1;                                    // Direction of half-strip counting: +1 is up, -1 is down
    int dStrip = 0;                                        // Offset for half-strip from full strip
    while ((infile >> buf1) && (infile >> buf2) && (infile >> buf3) && (infile >> buf4) && (infile >> buf5) &&
           (infile >> buf6)) {
      if ((line_num % 192) == 191)
        line_num += 1;  // Gap in central files vs. Jia Fu's files
      line_num += 1;
      // On roughly every-other line, files in L1Trigger/L1TMuon/data/cppf have 0 in the first three columns
      // Skips a "0 0 0" line once every 192 lines
      if ((line_num % 2) == 1) {
        buf1_prev = buf1;
        buf2_prev = buf2;
      }

      if (local && (buf1 == 0 || buf2 == 0)) {
        throw cms::Exception("L1TMuonEndCap") << "Expected non-0 values, got buf1 = " << buf1 << ", buf2 = " << buf2;
      }
      if (!local && (buf1_prev == 0 || buf2_prev == 0)) {
        throw cms::Exception("L1TMuonEndCap")
            << "Expected non-0 values, got buf1_prev = " << buf1_prev << ", buf2_prev = " << buf2_prev;
      }

      uint32_t id = (local ? buf1 : buf1_prev);
      int32_t rpc_region = get_rpc_region(id);
      int32_t rpc_sector = get_rpc_sector(id);
      int32_t rpc_station = get_rpc_station(id);
      int32_t rpc_ring = get_rpc_ring(id);
      int32_t rpc_subsector = get_rpc_subsector(id);
      int32_t rpc_roll = get_rpc_roll(id);

      // Offset into halfstrips from centrally-provided LUTs
      if (buf2_prev * 2 > halfstrip_prev + 8 ||
          buf2_prev * 2 < halfstrip_prev - 8) {  // Starting a new series of strips
        if (buf2_prev == 1)
          count_dir = +1;  // Starting from a low number, counting up
        else
          count_dir = -1;  // Starting from a high number, counting down
      }
      if (count_dir == -1)
        dStrip = (buf2_prev * 2 == halfstrip_prev ? 1 : 0);
      if (count_dir == +1)
        dStrip = (buf2_prev * 2 == halfstrip_prev + 2 ? 1 : 0);
      if (buf2_prev * 2 < halfstrip_prev - 8 && buf2_prev == 1)
        dStrip = 1;

      //uint32_t strip        = buf2;
      uint32_t halfstrip =
          (local ? buf2
                 : buf2_prev * 2 -
                       dStrip);  // I modified the local text files to use 'halfstrip' instead of 'strip' in column 2
      halfstrip_prev = halfstrip;

      uint32_t ph = buf5;
      uint32_t th = buf6;

      const uint32_t th_index = get_cppf_lut_id(rpc_region, rpc_sector, rpc_station, rpc_ring, rpc_subsector, rpc_roll);
      const uint32_t ph_index = (th_index * 64) + (halfstrip - 1);

      // std::cout << id << " " << rpc_region << " " << rpc_sector << " " << rpc_station << " " << rpc_ring << " "
      //    << rpc_subsector << " " << rpc_roll << " " << halfstrip << " " << th_index << " " << ph_index << std::endl;

      vec1.at(ph_index) = ph;
      if (halfstrip == 1)
        vec2.at(th_index) = th;

      // Fill gap in centrally-provided LUTs once every 192 lines
      if (!local && (line_num % 192) == 191)
        vec1.at(ph_index + 1) = ph;

    }  // End while ((infile >> buf1) && ... && (infile >> buf6))
    infile.close();
  }  // End loop over CPPF LUT files
}

void SectorProcessorLUT::read_file ( const std::string &  filename, std::vector< uint32_t > &  vec  )

private

Definition at line 351 of file SectorProcessorLUT.cc.

References visDQMUpload::buf, contentValuesFiles::fullPath, and timingPdfMaker::infile.

Referenced by read().

                                                                                      {
  vec.clear();

  std::ifstream infile;
  infile.open(edm::FileInPath(filename).fullPath().c_str());

  int buf;
  while (infile >> buf) {
    buf = (buf == -999) ? 0 : buf;
    vec.push_back(buf);
  }
  infile.close();
}

Member Data Documentation

std::vector<uint32_t> SectorProcessorLUT::cppf_ph_lut_

private

Definition at line 67 of file SectorProcessorLUT.h.

Referenced by get_cppf_ph_lut(), and read().

std::vector<uint32_t> SectorProcessorLUT::cppf_th_lut_

private

Definition at line 68 of file SectorProcessorLUT.h.

Referenced by get_cppf_th_lut(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_disp_neighbor_

private

Definition at line 56 of file SectorProcessorLUT.h.

Referenced by get_ph_disp(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_init_hard_

private

Definition at line 65 of file SectorProcessorLUT.h.

Referenced by get_ph_init_hard(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_init_neighbor_

private

Definition at line 55 of file SectorProcessorLUT.h.

Referenced by get_ph_init(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_patt_corr_

private

Definition at line 62 of file SectorProcessorLUT.h.

Referenced by get_ph_patt_corr(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_patt_corr_sign_

private

Definition at line 63 of file SectorProcessorLUT.h.

Referenced by get_ph_patt_corr_sign(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_zone_offset_

private

Definition at line 64 of file SectorProcessorLUT.h.

Referenced by get_ph_zone_offset(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_corr_lut_neighbor_

private

Definition at line 60 of file SectorProcessorLUT.h.

Referenced by get_th_corr_lut(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_disp_neighbor_

private

Definition at line 58 of file SectorProcessorLUT.h.

Referenced by get_th_disp(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_init_neighbor_

private

Definition at line 57 of file SectorProcessorLUT.h.

Referenced by get_th_init(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_lut_neighbor_

private

Definition at line 59 of file SectorProcessorLUT.h.

Referenced by get_th_lut(), and read().

int SectorProcessorLUT::version_

private

Definition at line 70 of file SectorProcessorLUT.h.

Referenced by read().