SectorProcessorLUT Class Reference

#include <SectorProcessorLUT.h>

Public Member Functions
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 (unsigned pc_lut_version)
	SectorProcessorLUT ()
	~SectorProcessorLUT ()

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

Private Attributes
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_
unsigned	version_

Detailed Description

Definition at line 9 of file SectorProcessorLUT.h.

Constructor & Destructor Documentation

SectorProcessorLUT::SectorProcessorLUT ( )

explicit

Definition at line 11 of file SectorProcessorLUT.cc.

                                       :
    version_(0xFFFFFFFF)
{

}

SectorProcessorLUT::~SectorProcessorLUT

(

)

Definition at line 17 of file SectorProcessorLUT.cc.

                                        {

}

Member Function Documentation

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

Definition at line 140 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and ph_disp_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  size_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  return ph_disp_neighbor_.at(index);
}

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

Definition at line 135 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and ph_init_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  size_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  return ph_init_neighbor_.at(index);
}

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

Definition at line 207 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and ph_init_hard_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                {
  size_t index = fw_station * 16 + fw_cscid;
  return ph_init_hard_.at(index);
}

uint32_t SectorProcessorLUT::get_ph_patt_corr

(

int

pattern

)

const

Definition at line 194 of file SectorProcessorLUT.cc.

References ph_patt_corr_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                               {
  return ph_patt_corr_.at(pattern);
}

uint32_t SectorProcessorLUT::get_ph_patt_corr_sign

(

int

pattern

)

const

Definition at line 198 of file SectorProcessorLUT.cc.

References ph_patt_corr_sign_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                    {
  return ph_patt_corr_sign_.at(pattern);
}

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

Definition at line 202 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and ph_zone_offset_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                    {
  size_t index = pc_station * 9 + pc_chamber;
  return ph_zone_offset_.at(index);
}

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 169 of file SectorProcessorLUT.cc.

References Exception, diffTreeTool::index, 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 {
    throw cms::Exception("SectorProcessorLUT")
      << "get_th_corr_lut(): out of range pc_lut_id: " << pc_lut_id;
  }

  size_t index = ((fw_endcap * 6 + fw_sector) * 7 + pc_lut_id2) * 128 + pc_wire_strip_id;
  return th_corr_lut_neighbor_.at(index);
}

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

Definition at line 150 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and th_disp_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  size_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  return th_disp_neighbor_.at(index);
}

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

Definition at line 145 of file SectorProcessorLUT.cc.

References diffTreeTool::index, and th_init_neighbor_.

Referenced by PrimitiveConversion::convert_csc_details().

                                                                                          {
  size_t index = (fw_endcap * 6 + fw_sector) * 61 + pc_lut_id;
  return th_init_neighbor_.at(index);
}

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

Definition at line 155 of file SectorProcessorLUT.cc.

References diffTreeTool::index, 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;

  size_t index = ((fw_endcap * 6 + fw_sector) * 61 + pc_lut_id2) * 128 + pc_wire_id;
  return th_lut_neighbor_.at(index);
}

void SectorProcessorLUT::read

(

unsigned

pc_lut_version

)

Definition at line 21 of file SectorProcessorLUT.cc.

References Exception, ph_disp_neighbor_, ph_init_hard_, ph_init_neighbor_, ph_patt_corr_, ph_patt_corr_sign_, ph_zone_offset_, 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 TrackFinder::process().

                                                     {
  if (version_ == pc_lut_version)  return;

  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
  else 
    throw cms::Exception("SectorProcessorLUT")
      << "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 + "/";

  // std::cout << "coord_lut_path = " << coord_lut_path << std::endl;

  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_);

  if (ph_init_neighbor_.size() != 2*6*61) {  // [endcap_2][sector_6][chamber_61]
    throw cms::Exception("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "Expected th_corr_lut_neighbor_ to get " << 2*6*7*128 << " values, "
        << "got " << th_corr_lut_neighbor_.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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "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("SectorProcessorLUT")
        << "Expected ph_init_hard_ to get " << 5*16 << " values, "
        << "got " << ph_init_hard_.size() << " values.";
  }

  version_ = pc_lut_version;
  return;
}

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

private

Definition at line 212 of file SectorProcessorLUT.cc.

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::ph_disp_neighbor_

private

Definition at line 40 of file SectorProcessorLUT.h.

Referenced by get_ph_disp(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_init_hard_

private

Definition at line 49 of file SectorProcessorLUT.h.

Referenced by get_ph_init_hard(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_init_neighbor_

private

Definition at line 39 of file SectorProcessorLUT.h.

Referenced by get_ph_init(), and read().

std::vector<uint32_t> SectorProcessorLUT::ph_patt_corr_

private

Definition at line 46 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 47 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 48 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 44 of file SectorProcessorLUT.h.

Referenced by get_th_corr_lut(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_disp_neighbor_

private

Definition at line 42 of file SectorProcessorLUT.h.

Referenced by get_th_disp(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_init_neighbor_

private

Definition at line 41 of file SectorProcessorLUT.h.

Referenced by get_th_init(), and read().

std::vector<uint32_t> SectorProcessorLUT::th_lut_neighbor_

private

Definition at line 43 of file SectorProcessorLUT.h.

Referenced by get_th_lut(), and read().

unsigned SectorProcessorLUT::version_

private

Definition at line 51 of file SectorProcessorLUT.h.

Referenced by read().