trklet::ProjectionRouterBendTable Class Reference

#include <ProjectionRouterBendTable.h>

Public Member Functions
int	bendLoookup (int diskindex, int bendindex)
void	init (Settings const &settings, Globals *globals, unsigned int nrbits, unsigned int nphiderbits)
	ProjectionRouterBendTable ()
	~ProjectionRouterBendTable ()=default

Private Attributes
std::vector< int >	bendtable_ [N_DISK]

Detailed Description

Definition at line 16 of file ProjectionRouterBendTable.h.

Constructor & Destructor Documentation

ProjectionRouterBendTable()

trklet::ProjectionRouterBendTable::ProjectionRouterBendTable ( )

inline

Definition at line 18 of file ProjectionRouterBendTable.h.

{}

~ProjectionRouterBendTable()

trklet::ProjectionRouterBendTable::~ProjectionRouterBendTable ( )

default

Member Function Documentation

bendLoookup()

int ProjectionRouterBendTable::bendLoookup	(	int	diskindex,
		int	bendindex
	)

Definition at line 54 of file ProjectionRouterBendTable.cc.

{ return bendtable_[diskindex][bendindex]; }

References bendtable_.

Referenced by trklet::ProjectionRouter::execute(), and trklet::MatchProcessor::execute().

init()

void ProjectionRouterBendTable::init	(	Settings const &	settings,
		Globals *	globals,
		unsigned int	nrbits,
		unsigned int	nphiderbits
	)

Definition at line 8 of file ProjectionRouterBendTable.cc.

                                                               {
  for (unsigned int idisk = 0; idisk < N_DISK; idisk++) {
    unsigned int nsignbins = 2;
    unsigned int nrbins = 1 << (nrbits);
    unsigned int nphiderbins = 1 << (nphiderbits);
 
    for (unsigned int isignbin = 0; isignbin < nsignbins; isignbin++) {
      for (unsigned int irbin = 0; irbin < nrbins; irbin++) {
        int ir = irbin;
        if (ir > (1 << (nrbits - 1)))
          ir -= (1 << nrbits);
        ir = ir << (settings.nrbitsstub(N_LAYER) - nrbits);
        for (unsigned int iphiderbin = 0; iphiderbin < nphiderbins; iphiderbin++) {
          int iphider = iphiderbin;
          if (iphider > (1 << (nphiderbits - 1)))
            iphider -= (1 << nphiderbits);
          iphider = iphider << (settings.nbitsphiprojderL123() - nphiderbits);
 
          double rproj = ir * settings.krprojshiftdisk();
          double phider = iphider * globals->ITC_L1L2()->der_phiD_final.K();
          double t = settings.zmean(idisk) / rproj;
 
          if (isignbin)
            t = -t;
 
          double rinv = -phider * (2.0 * t);
 
          double stripPitch = (rproj < settings.rcrit()) ? settings.stripPitch(true) : settings.stripPitch(false);
          double bendproj = 0.5 * bend(rproj, rinv, stripPitch);
 
          int ibendproj = 2.0 * bendproj + 15.5;
          if (ibendproj < 0)
            ibendproj = 0;
          if (ibendproj > 31)
            ibendproj = 31;
 
          bendtable_[idisk].push_back(ibendproj);
        }
      }
    }
  }
}

References trklet::bend(), bendtable_, trklet::IMATH_TrackletCalculator::der_phiD_final, trklet::Globals::ITC_L1L2(), trklet::VarBase::K(), trklet::Settings::krprojshiftdisk(), trklet::N_DISK, trklet::N_LAYER, trklet::Settings::nbitsphiprojderL123(), trklet::Settings::nrbitsstub(), trklet::Settings::rcrit(), trklet::rinv(), trklet::Settings::stripPitch(), OrderedSet::t, and trklet::Settings::zmean().

Member Data Documentation

bendtable_

std::vector<int> trklet::ProjectionRouterBendTable::bendtable_[N_DISK]

private

Definition at line 27 of file ProjectionRouterBendTable.h.

Referenced by bendLoookup(), and init().