trklet::VMRouterTable Class Reference

#include <VMRouterTable.h>

Public Member Functions
int	getLookup (unsigned int layerdisk, double z, double r, int iseed=-1)
void	init (unsigned int layerdisk, std::string const &name)
int	lookup (int zbin, int rbin)
int	lookupdisk (int zbin, int rbin)
int	lookupinner (int zbin, int rbin)
int	lookupinneroverlap (int zbin, int rbin)
int	lookupinnerThird (int zbin, int rbin)
	VMRouterTable (Settings const &settings)
	VMRouterTable (Settings const &settings, unsigned int layerdisk, std::string const &name)
void	writeVMTable (std::string const &name, std::vector< int > const &table)
	~VMRouterTable ()=default

Private Attributes
double	dr_
double	dz_
int	rbins_
int	rbits_
double	rmax_
double	rmin_
Settings const &	settings_
std::vector< int >	vmrtable_
std::vector< int >	vmrtabletedisk_
std::vector< int >	vmrtableteinner_
std::vector< int >	vmrtableteinneroverlap_
std::vector< int >	vmrtableteinnerThird_
int	zbins_
int	zbits_
double	zmax_
double	zmin_

Detailed Description

Definition at line 16 of file VMRouterTable.h.

Constructor & Destructor Documentation

VMRouterTable() [1/2]

VMRouterTable::VMRouterTable

(

Settings const &

settings

)

Definition at line 10 of file VMRouterTable.cc.

: settings_(settings) {}

VMRouterTable() [2/2]

VMRouterTable::VMRouterTable	(	Settings const &	settings,
		unsigned int	layerdisk,
		std::string const &	name
	)

Definition at line 12 of file VMRouterTable.cc.

    : settings_(settings) {
  init(layerdisk, name);
}

References init(), and Skims_PA_cff::name.

~VMRouterTable()

trklet::VMRouterTable::~VMRouterTable ( )

default

Member Function Documentation

getLookup()

int VMRouterTable::getLookup	(	unsigned int	layerdisk,
		double	z,
		double	r,
		int	iseed = -1
	)

Definition at line 141 of file VMRouterTable.cc.

                                                                                  {
  double z0cut = settings_.z0cut();
 
  if (layerdisk < N_LAYER) {
    if (iseed == 1 && std::abs(z) < 52.0)
      return -1;
 
    double rmean = settings_.rmean(layerdisk);
 
    double rratio1 = rmean / (r + 0.5 * dr_);
    double rratio2 = rmean / (r - 0.5 * dr_);
 
    double z1 = (z - 0.5 * dz_) * rratio1 + z0cut * (rratio1 - 1.0);
    double z2 = (z + 0.5 * dz_) * rratio1 + z0cut * (rratio1 - 1.0);
    double z3 = (z - 0.5 * dz_) * rratio2 + z0cut * (rratio2 - 1.0);
    double z4 = (z + 0.5 * dz_) * rratio2 + z0cut * (rratio2 - 1.0);
    double z5 = (z - 0.5 * dz_) * rratio1 - z0cut * (rratio1 - 1.0);
    double z6 = (z + 0.5 * dz_) * rratio1 - z0cut * (rratio1 - 1.0);
    double z7 = (z - 0.5 * dz_) * rratio2 - z0cut * (rratio2 - 1.0);
    double z8 = (z + 0.5 * dz_) * rratio2 - z0cut * (rratio2 - 1.0);
 
    double zmin = std::min({z1, z2, z3, z4, z5, z6, z7, z8});
    double zmax = std::max({z1, z2, z3, z4, z5, z6, z7, z8});
 
    int NBINS = settings_.NLONGVMBINS() * settings_.NLONGVMBINS();
 
    int zbin1 = NBINS * (zmin + settings_.zlength()) / (2 * settings_.zlength());
    int zbin2 = NBINS * (zmax + settings_.zlength()) / (2 * settings_.zlength());
 
    if (zbin1 >= NBINS)
      return -1;
    if (zbin2 < 0)
      return -1;
 
    if (zbin2 >= NBINS)
      zbin2 = NBINS - 1;
    if (zbin1 < 0)
      zbin1 = 0;
 
    // This is a 10 bit word:
    // xxx|yyy|z|rrr
    // xxx is the delta z window
    // yyy is the z bin
    // z is flag to look in next bin
    // rrr first fine z bin
    // NOTE : this encoding is not efficient z is one if xxx+rrr is greater than 8
    //        and xxx is only 1,2, or 3
    //        should also reject xxx=0 as this means projection is outside range
 
    int value = zbin1 / 8;
    value *= 2;
    if (zbin2 / 8 - zbin1 / 8 > 0)
      value += 1;
    value *= 8;
    value += (zbin1 & 7);
    assert(value / 8 < 15);
    int deltaz = zbin2 - zbin1;
    if (deltaz > 7) {
      deltaz = 7;
    }
    assert(deltaz < 8);
    value += (deltaz << 7);
 
    return value;
 
  } else {
    if (std::abs(z) < 2.0 * z0cut)
      return -1;
 
    double zmean = settings_.zmean(layerdisk - N_LAYER);
    if (z < 0.0)
      zmean = -zmean;
 
    double r1 = (r + 0.5 * dr_) * (zmean + z0cut) / (z + 0.5 * dz_ + z0cut);
    double r2 = (r - 0.5 * dr_) * (zmean - z0cut) / (z + 0.5 * dz_ - z0cut);
    double r3 = (r + 0.5 * dr_) * (zmean + z0cut) / (z - 0.5 * dz_ + z0cut);
    double r4 = (r - 0.5 * dr_) * (zmean - z0cut) / (z - 0.5 * dz_ - z0cut);
    double r5 = (r + 0.5 * dr_) * (zmean - z0cut) / (z + 0.5 * dz_ - z0cut);
    double r6 = (r - 0.5 * dr_) * (zmean + z0cut) / (z + 0.5 * dz_ + z0cut);
    double r7 = (r + 0.5 * dr_) * (zmean - z0cut) / (z - 0.5 * dz_ - z0cut);
    double r8 = (r - 0.5 * dr_) * (zmean + z0cut) / (z - 0.5 * dz_ + z0cut);
 
    double rmin = std::min({r1, r2, r3, r4, r5, r6, r7, r8});
    double rmax = std::max({r1, r2, r3, r4, r5, r6, r7, r8});
 
    int NBINS = settings_.NLONGVMBINS() * settings_.NLONGVMBINS() / 2;
 
    double rmindisk = settings_.rmindiskvm();
    double rmaxdisk = settings_.rmaxdiskvm();
 
    if (iseed == 6)
      rmaxdisk = settings_.rmaxdiskl1overlapvm();
    if (iseed == 7)
      rmindisk = settings_.rmindiskl2overlapvm();
    if (iseed == 10)
      rmaxdisk = settings_.rmaxdisk();
 
    if (rmin > rmaxdisk)
      return -1;
    if (rmax > rmaxdisk)
      rmax = rmaxdisk;
 
    if (rmax < rmindisk)
      return -1;
    if (rmin < rmindisk)
      rmin = rmindisk;
 
    int rbin1 = NBINS * (rmin - settings_.rmindiskvm()) / (settings_.rmaxdiskvm() - settings_.rmindiskvm());
    int rbin2 = NBINS * (rmax - settings_.rmindiskvm()) / (settings_.rmaxdiskvm() - settings_.rmindiskvm());
 
    if (iseed == 10) {
      constexpr double rminspec = 40.0;
      rbin1 = NBINS * (rmin - rminspec) / (settings_.rmaxdisk() - rminspec);
      rbin2 = NBINS * (rmax - rminspec) / (settings_.rmaxdisk() - rminspec);
    }
 
    if (rbin2 >= NBINS)
      rbin2 = NBINS - 1;
    if (rbin1 < 0)
      rbin1 = 0;
 
    // This is a 9 bit word:
    // xxx|yy|z|rrr
    // xxx is the delta r window
    // yy is the r bin yy is three bits for overlaps
    // z is flag to look in next bin
    // rrr fine r bin
    // NOTE : this encoding is not efficient z is one if xxx+rrr is greater than 8
    //        and xxx is only 1,2, or 3
    //        should also reject xxx=0 as this means projection is outside range
 
    bool overlap = iseed == 6 || iseed == 7 || iseed == 10;
 
    int value = rbin1 / 8;
    if (overlap) {
      if (z < 0.0)
        value += 4;
    }
    value *= 2;
    if (rbin2 / 8 - rbin1 / 8 > 0)
      value += 1;
    value *= 8;
    value += (rbin1 & 7);
    assert(value / 8 < 15);
    int deltar = rbin2 - rbin1;
    if (deltar > 7)
      deltar = 7;
    if (overlap) {
      value += (deltar << 7);
    } else {
      value += (deltar << 6);
    }
 
    return value;
  }
}

References funct::abs(), cms::cuda::assert(), dr_, dz_, iseed, SiStripPI::max, min(), trklet::N_LAYER, NBINS, trklet::Settings::NLONGVMBINS(), muon::overlap(), alignCSCRings::r, diffTwoXMLs::r1, diffTwoXMLs::r2, trklet::Settings::rmaxdisk(), trklet::Settings::rmaxdiskl1overlapvm(), trklet::Settings::rmaxdiskvm(), trklet::Settings::rmean(), trklet::Settings::rmindiskl2overlapvm(), trklet::Settings::rmindiskvm(), settings_, relativeConstraints::value, trklet::Settings::z0cut(), testProducerWithPsetDescEmpty_cfi::z2, trklet::Settings::zlength(), SiStripMonitorCluster_cfi::zmax, trklet::Settings::zmean(), and SiStripMonitorCluster_cfi::zmin.

Referenced by init().

init()

void VMRouterTable::init	(	unsigned int	layerdisk,
		std::string const &	name
	)

Definition at line 17 of file VMRouterTable.cc.

                                                                      {
  zbits_ = settings_.vmrlutzbits(layerdisk);
  rbits_ = settings_.vmrlutrbits(layerdisk);
 
  rbins_ = (1 << rbits_);
  zbins_ = (1 << zbits_);
 
  if (layerdisk < N_LAYER) {
    zmin_ = -settings_.zlength();
    zmax_ = settings_.zlength();
    rmin_ = settings_.rmean(layerdisk) - settings_.drmax();
    rmax_ = settings_.rmean(layerdisk) + settings_.drmax();
  } else {
    rmin_ = 0;
    rmax_ = settings_.rmaxdisk();
    zmin_ = settings_.zmean(layerdisk - N_LAYER) - settings_.dzmax();
    zmax_ = settings_.zmean(layerdisk - N_LAYER) + settings_.dzmax();
  }
 
  dr_ = (rmax_ - rmin_) / rbins_;
  dz_ = (zmax_ - zmin_) / zbins_;
 
  int NBINS = settings_.NLONGVMBINS() * settings_.NLONGVMBINS();
 
  for (int izbin = 0; izbin < zbins_; izbin++) {
    for (int irbin = 0; irbin < rbins_; irbin++) {
      double r = rmin_ + (irbin + 0.5) * dr_;
      double z = zmin_ + (izbin + 0.5) * dz_;
 
      if (layerdisk > (N_LAYER - 1) && irbin < 10)  //special case for the tabulated radii in 2S disks
        r = (layerdisk <= 7) ? settings_.rDSSinner(irbin) : settings_.rDSSouter(irbin);
 
      int bin;
      if (layerdisk < N_LAYER) {
        double zproj = z * settings_.rmean(layerdisk) / r;
        bin = NBINS * (zproj + settings_.zlength()) / (2 * settings_.zlength());
      } else {
        double rproj = r * settings_.zmean(layerdisk - N_LAYER) / z;
        bin = NBINS * (rproj - settings_.rmindiskvm()) / (settings_.rmaxdisk() - settings_.rmindiskvm());
      }
      if (bin < 0)
        bin = 0;
      if (bin >= NBINS)
        bin = NBINS - 1;
      vmrtable_.push_back(bin);
 
      if (layerdisk >= N_LAYER) {
        double rproj = r * settings_.zmean(layerdisk - N_LAYER) / z;
        bin = 0.5 * NBINS * (rproj - settings_.rmindiskvm()) / (settings_.rmaxdiskvm() - settings_.rmindiskvm());
        if (bin < 0)
          bin = 0;
        if (bin >= NBINS / 2)
          bin = NBINS / 2 - 1;
        vmrtabletedisk_.push_back(bin);
      }
 
      if (layerdisk == 0 || layerdisk == 2 || layerdisk == 4 || layerdisk == 6 || layerdisk == 8) {
        vmrtableteinner_.push_back(getLookup(layerdisk + 1, z, r));
      }
 
      if (layerdisk == 1) {
        vmrtableteinner_.push_back(getLookup(layerdisk + 1, z, r, 1));
      }
 
      if (layerdisk == 1) {  //projection from L2 to D1 for L2L3D1 seeding
        vmrtableteinnerThird_.push_back(getLookup(6, z, r, 10));
      }
 
      if (layerdisk == 4) {  //projection from L5 to L4 for L5L6L4 seeding
        vmrtableteinnerThird_.push_back(getLookup(3, z, r));
      }
 
      if (layerdisk == 2) {  //projection from L3 to L5 for L3L4L2 seeding
        vmrtableteinnerThird_.push_back(getLookup(1, z, r));
      }
 
      if (layerdisk == 6) {  //projection from D1 to L2 for D1D2L2 seeding
        vmrtableteinnerThird_.push_back(getLookup(1, z, r));
      }
 
      if (layerdisk == 0 || layerdisk == 1) {
        vmrtableteinneroverlap_.push_back(getLookup(6, z, r, layerdisk + 6));
      }
    }
  }
 
  if (settings_.writeTable()) {
    if (not std::filesystem::exists(settings_.tablePath())) {
      int fail = system((string("mkdir -p ") + settings_.tablePath()).c_str());
      if (fail)
        throw cms::Exception("BadDir") << __FILE__ << " " << __LINE__ << " could not create directory "
                                       << settings_.tablePath();
    }
 
    // write finebin tables
    writeVMTable(settings_.tablePath() + name + "_finebin.tab", vmrtable_);
    // write barrel seed teinner tables (L1L2, L2L3, L3L4, L5L6)
    if (layerdisk == 0 || layerdisk == 1 || layerdisk == 2 || layerdisk == 4) {
      std::string fnamesuffix = "L" + to_string(layerdisk + 1) + "L" + std::to_string(layerdisk + 2);
      writeVMTable(settings_.tablePath() + "VMTableInner" + fnamesuffix + ".tab", vmrtableteinner_);
    }
    // write disk seed teinner tables (D1D2, D3D4)
    if (layerdisk == 6 || layerdisk == 8) {
      std::string fnamesuffix = "D" + to_string(layerdisk - N_LAYER + 1) + "D" + to_string(layerdisk - N_LAYER + 2);
      writeVMTable(settings_.tablePath() + "VMTableInner" + fnamesuffix + ".tab", vmrtableteinner_);
    }
    // write overlap seed teinner tables (L1D1, L2D1)
    if (layerdisk == 0 || layerdisk == 1) {
      std::string fnamesuffix = "L" + to_string(layerdisk + 1) + "D1";
      writeVMTable(settings_.tablePath() + "VMTableInner" + fnamesuffix + ".tab", vmrtableteinneroverlap_);
    }
    // write barrel teouter tables (L2, L3, L4, L6, same as finebin tables)
    if (layerdisk == 1 || layerdisk == 2 || layerdisk == 3 || layerdisk == 5) {
      std::string fnamesuffix = "L" + to_string(layerdisk + 1);
      writeVMTable(settings_.tablePath() + "VMTableOuter" + fnamesuffix + ".tab", vmrtable_);
    }
    // write disk teouter tables (D1, D2, D4)
    if (layerdisk == 6 || layerdisk == 7 || layerdisk == 9) {
      std::string fnamesuffix = "D" + to_string(layerdisk - N_LAYER + 1);
      writeVMTable(settings_.tablePath() + "VMTableOuter" + fnamesuffix + ".tab", vmrtabletedisk_);
    }
  }
}

References newFWLiteAna::bin, dr_, trklet::Settings::drmax(), dz_, trklet::Settings::dzmax(), Exception, getLookup(), trklet::N_LAYER, Skims_PA_cff::name, NBINS, trklet::Settings::NLONGVMBINS(), alignCSCRings::r, rbins_, rbits_, trklet::Settings::rDSSinner(), trklet::Settings::rDSSouter(), rmax_, trklet::Settings::rmaxdisk(), trklet::Settings::rmaxdiskvm(), trklet::Settings::rmean(), rmin_, trklet::Settings::rmindiskvm(), settings_, AlCaHLTBitMon_QueryRunRegistry::string, trklet::Settings::tablePath(), trklet::Settings::vmrlutrbits(), trklet::Settings::vmrlutzbits(), vmrtable_, vmrtabletedisk_, vmrtableteinner_, vmrtableteinneroverlap_, vmrtableteinnerThird_, trklet::Settings::writeTable(), writeVMTable(), zbins_, zbits_, trklet::Settings::zlength(), zmax_, trklet::Settings::zmean(), and zmin_.

Referenced by trklet::TrackletProcessor::TrackletProcessor(), trklet::VMRouter::VMRouter(), trklet::VMRouterCM::VMRouterCM(), and VMRouterTable().

lookup()

int VMRouterTable::lookup	(	int	zbin,
		int	rbin
	)

Definition at line 298 of file VMRouterTable.cc.

                                            {
  int index = zbin * rbins_ + rbin;
  assert(index >= 0 && index < (int)vmrtable_.size());
  return vmrtable_[index];
}

References cms::cuda::assert(), rbins_, and vmrtable_.

Referenced by trklet::VMRouterCM::execute(), and trklet::VMRouter::execute().

lookupdisk()

int VMRouterTable::lookupdisk	(	int	zbin,
		int	rbin
	)

Definition at line 304 of file VMRouterTable.cc.

                                                {
  int index = zbin * rbins_ + rbin;
  assert(index >= 0 && index < (int)vmrtabletedisk_.size());
  return vmrtabletedisk_[index];
}

References cms::cuda::assert(), rbins_, and vmrtabletedisk_.

Referenced by trklet::VMRouterCM::execute(), and trklet::VMRouter::execute().

lookupinner()

int VMRouterTable::lookupinner	(	int	zbin,
		int	rbin
	)

Definition at line 310 of file VMRouterTable.cc.

                                                 {
  int index = zbin * rbins_ + rbin;
  assert(index >= 0 && index < (int)vmrtableteinner_.size());
  return vmrtableteinner_[index];
}

References cms::cuda::assert(), rbins_, and vmrtableteinner_.

Referenced by trklet::TrackletProcessor::execute(), and trklet::VMRouter::execute().

lookupinneroverlap()

int VMRouterTable::lookupinneroverlap	(	int	zbin,
		int	rbin
	)

Definition at line 316 of file VMRouterTable.cc.

                                                        {
  int index = zbin * rbins_ + rbin;
  assert(index >= 0 && index < (int)vmrtableteinneroverlap_.size());
  return vmrtableteinneroverlap_[index];
}

References cms::cuda::assert(), rbins_, and vmrtableteinneroverlap_.

Referenced by trklet::TrackletProcessor::execute(), and trklet::VMRouter::execute().

lookupinnerThird()

int VMRouterTable::lookupinnerThird	(	int	zbin,
		int	rbin
	)

Definition at line 322 of file VMRouterTable.cc.

                                                      {
  int index = zbin * rbins_ + rbin;
  assert(index >= 0 && index < (int)vmrtableteinnerThird_.size());
  return vmrtableteinnerThird_[index];
}

References cms::cuda::assert(), rbins_, and vmrtableteinnerThird_.

Referenced by trklet::VMRouter::execute().

writeVMTable()

void VMRouterTable::writeVMTable	(	std::string const &	name,
		std::vector< int > const &	table
	)

Definition at line 328 of file VMRouterTable.cc.

                                                                                   {
  ofstream out(name);
  if (out.fail())
    throw cms::Exception("BadFile") << __FILE__ << " " << __LINE__ << " could not create file " << name;
 
  out << "{" << endl;
  for (unsigned int i = 0; i < table.size(); i++) {
    if (i != 0) {
      out << "," << endl;
    }
    int itable = table[i];
    out << itable;
  }
  out << endl << "};" << endl;
  out.close();
}

References mps_fire::i, Skims_PA_cff::name, MillePedeFileConverter_cfg::out, and TableParser::table.

Referenced by init().

Member Data Documentation

dr_

double trklet::VMRouterTable::dr_

private

Definition at line 44 of file VMRouterTable.h.

Referenced by getLookup(), and init().

dz_

double trklet::VMRouterTable::dz_

private

Definition at line 45 of file VMRouterTable.h.

Referenced by getLookup(), and init().

rbins_

int trklet::VMRouterTable::rbins_

private

Definition at line 51 of file VMRouterTable.h.

Referenced by init(), lookup(), lookupdisk(), lookupinner(), lookupinneroverlap(), and lookupinnerThird().

rbits_

int trklet::VMRouterTable::rbits_

private

Definition at line 48 of file VMRouterTable.h.

Referenced by init().

rmax_

double trklet::VMRouterTable::rmax_

private

Definition at line 39 of file VMRouterTable.h.

Referenced by init().

rmin_

double trklet::VMRouterTable::rmin_

private

Definition at line 38 of file VMRouterTable.h.

Referenced by init().

settings_

Settings const& trklet::VMRouterTable::settings_

private

Definition at line 36 of file VMRouterTable.h.

Referenced by getLookup(), and init().

vmrtable_

std::vector<int> trklet::VMRouterTable::vmrtable_

private

Definition at line 54 of file VMRouterTable.h.

Referenced by init(), and lookup().

vmrtabletedisk_

std::vector<int> trklet::VMRouterTable::vmrtabletedisk_

private

Definition at line 55 of file VMRouterTable.h.

Referenced by init(), and lookupdisk().

vmrtableteinner_

std::vector<int> trklet::VMRouterTable::vmrtableteinner_

private

Definition at line 56 of file VMRouterTable.h.

Referenced by init(), and lookupinner().

vmrtableteinneroverlap_

std::vector<int> trklet::VMRouterTable::vmrtableteinneroverlap_

private

Definition at line 57 of file VMRouterTable.h.

Referenced by init(), and lookupinneroverlap().

vmrtableteinnerThird_

std::vector<int> trklet::VMRouterTable::vmrtableteinnerThird_

private

Definition at line 58 of file VMRouterTable.h.

Referenced by init(), and lookupinnerThird().

zbins_

int trklet::VMRouterTable::zbins_

private

Definition at line 50 of file VMRouterTable.h.

Referenced by init().

zbits_

int trklet::VMRouterTable::zbits_

private

Definition at line 47 of file VMRouterTable.h.

Referenced by init().

zmax_

double trklet::VMRouterTable::zmax_

private

Definition at line 42 of file VMRouterTable.h.

Referenced by init().

zmin_

double trklet::VMRouterTable::zmin_

private

Definition at line 41 of file VMRouterTable.h.

Referenced by init().