trklet::ProjectionRouter Class Reference

#include <ProjectionRouter.h>

Inheritance diagram for trklet::ProjectionRouter:

Public Member Functions
void	addInput (MemoryBase *memory, std::string input) override
void	addOutput (MemoryBase *memory, std::string output) override
void	execute ()
	ProjectionRouter (std::string name, Settings const &settings, Globals *global)
	~ProjectionRouter () override=default
Public Member Functions inherited from trklet::ProcessBase
unsigned int	getISeed (const std::string &name)
std::string const &	getName () const
void	initLayerDisk (unsigned int pos, int &layer, int &disk)
void	initLayerDisk (unsigned int pos, int &layer, int &disk, int &layerdisk)
unsigned int	initLayerDisk (unsigned int pos)
void	initLayerDisksandISeed (unsigned int &layerdisk1, unsigned int &layerdisk2, unsigned int &iSeed)
unsigned int	nbits (unsigned int power)
	ProcessBase (std::string name, Settings const &settings, Globals *global)
virtual	~ProcessBase ()=default

Private Attributes
AllProjectionsMemory *	allproj_
std::vector< TrackletProjectionsMemory * >	inputproj_
unsigned int	layerdisk_
int	nphiderbits_
int	nrbits_
TrackletLUT	rinvbendlut_
std::vector< VMProjectionsMemory * >	vmprojs_

Additional Inherited Members
Protected Attributes inherited from trklet::ProcessBase
Globals *	globals_
std::string	name_
Settings const &	settings_

Detailed Description

Definition at line 16 of file ProjectionRouter.h.

Constructor & Destructor Documentation

ProjectionRouter()

ProjectionRouter::ProjectionRouter	(	std::string	name,
		Settings const &	settings,
		Globals *	global
	)

Definition at line 14 of file ProjectionRouter.cc.

References trklet::IMATH_TrackletCalculator::der_phiD_final, trklet::ProcessBase::initLayerDisk(), trklet::TrackletLUT::initProjectionBend(), trklet::Globals::ITC_L1L2(), trklet::VarBase::K(), layerdisk_, trklet::N_LAYER, nphiderbits_, nrbits_, trklet::Settings::nvmme(), rinvbendlut_, trklet::ProcessBase::settings_, and vmprojs_.

    : ProcessBase(name, settings, global), rinvbendlut_(settings) {
  layerdisk_ = initLayerDisk(3);

  vmprojs_.resize(settings_.nvmme(layerdisk_), nullptr);

  nrbits_ = 5;
  nphiderbits_ = 6;

  if (layerdisk_ >= N_LAYER) {
    rinvbendlut_.initProjectionBend(
        global->ITC_L1L2()->der_phiD_final.K(), layerdisk_ - N_LAYER, nrbits_, nphiderbits_);
  }
}

~ProjectionRouter()

trklet::ProjectionRouter::~ProjectionRouter ( )

overridedefault

Member Function Documentation

addInput()

void ProjectionRouter::addInput ( MemoryBase *  memory, std::string  input  )

overridevirtual

Implements trklet::ProcessBase.

Definition at line 61 of file ProjectionRouter.cc.

References cms::cuda::assert(), Exception, input, inputproj_, mps_setup::memory, trklet::ProcessBase::name_, trklet::ProcessBase::settings_, createJobs::tmp, and trklet::Settings::writetrace().

                                                                {
  if (settings_.writetrace()) {
    edm::LogVerbatim("Tracklet") << "In " << name_ << " adding input from " << memory->getName() << " to input "
                                 << input;
  }
  if (input.substr(0, 4) == "proj" && input.substr(input.size() - 2, 2) == "in") {
    auto* tmp = dynamic_cast<TrackletProjectionsMemory*>(memory);
    assert(tmp != nullptr);
    inputproj_.push_back(tmp);
    return;
  }
  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " could not find input: " << input;
}

addOutput()

void ProjectionRouter::addOutput ( MemoryBase *  memory, std::string  output  )

overridevirtual

Implements trklet::ProcessBase.

Definition at line 29 of file ProjectionRouter.cc.

References allproj_, cms::cuda::assert(), Exception, layerdisk_, mps_setup::memory, trklet::Settings::nallstubs(), Skims_PA_cff::name, trklet::ProcessBase::name_, trklet::Settings::nvmme(), convertSQLitetoXML_cfg::output, trklet::ProcessBase::settings_, AlCaHLTBitMon_QueryRunRegistry::string, createJobs::tmp, cond::impl::to_string(), vmprojs_, and trklet::Settings::writetrace().

                                                                  {
  if (settings_.writetrace()) {
    edm::LogVerbatim("Tracklet") << "In " << name_ << " adding output to " << memory->getName() << " to output "
                                 << output;
  }
  if (output == "allprojout") {
    auto* tmp = dynamic_cast<AllProjectionsMemory*>(memory);
    assert(tmp != nullptr);
    allproj_ = tmp;
    return;
  }

  unsigned int nproj = settings_.nallstubs(layerdisk_);
  unsigned int nprojvm = settings_.nvmme(layerdisk_);

  for (unsigned int iproj = 0; iproj < nproj; iproj++) {
    for (unsigned int iprojvm = 0; iprojvm < nprojvm; iprojvm++) {
      std::string name = "vmprojoutPHI";
      name += char(iproj + 'A');
      name += std::to_string(iproj * nprojvm + iprojvm + 1);
      if (output == name) {
        auto* tmp = dynamic_cast<VMProjectionsMemory*>(memory);
        assert(tmp != nullptr);
        vmprojs_[iprojvm] = tmp;
        return;
      }
    }
  }

  throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__ << " could not find output: " << output;
}

execute()

void ProjectionRouter::execute

(

void

)

Definition at line 75 of file ProjectionRouter.cc.

References trklet::AllProjectionsMemory::addTracklet(), allproj_, trklet::FPGAWord::bits(), trklet::Settings::debugTracklet(), trklet::Projection::fpgaphiproj(), trklet::ProcessBase::getName(), trklet::ProcessBase::globals_, mps_fire::i, inputproj_, layerdisk_, trklet::TrackletLUT::lookup(), trklet::Settings::maxStep(), trklet::N_LAYER, trklet::FPGAWord::nbits(), trklet::Settings::nbitsallstubs(), trklet::Settings::nbitsvmme(), nphiderbits_, nrbits_, trklet::AllProjectionsMemory::nTracklets(), trklet::Globals::ofstream(), MillePedeFileConverter_cfg::out, amptDefault_cfi::proj, trklet::Tracklet::proj(), rinvbendlut_, trklet::ProcessBase::settings_, trklet::Tracklet::TCID(), vmprojs_, and trklet::Settings::writeMonitorData().

                               {
  unsigned int allprojcount = 0;

  //These are just here to test that the order is correct. Does not affect the actual execution

  int lastTCID = -1;

  for (auto& iproj : inputproj_) {
    for (unsigned int i = 0; i < iproj->nTracklets(); i++) {
      if (allprojcount >= settings_.maxStep("PR"))
        continue;

      Tracklet* tracklet = iproj->getTracklet(i);

      FPGAWord fpgaphi;

      if (layerdisk_ < N_LAYER) {
        fpgaphi = tracklet->proj(layerdisk_).fpgaphiproj();
      } else {
        Projection& proj = tracklet->proj(layerdisk_);
        fpgaphi = proj.fpgaphiproj();

        //The next lines looks up the predicted bend based on:
        // 1 - r projections
        // 2 - phi derivative
        // 3 - the sign - i.e. if track is forward or backward

        int rindex = (proj.fpgarzproj().value() >> (proj.fpgarzproj().nbits() - nrbits_)) & ((1 << nrbits_) - 1);

        int phiderindex = (proj.fpgaphiprojder().value() >> (proj.fpgaphiprojder().nbits() - nphiderbits_)) &
                          ((1 << nphiderbits_) - 1);

        int signindex = (proj.fpgarzprojder().value() < 0);

        int bendindex = (signindex << (nphiderbits_ + nrbits_)) + (rindex << (nphiderbits_)) + phiderindex;

        int ibendproj = rinvbendlut_.lookup(bendindex);

        proj.setBendIndex(ibendproj);
      }

      unsigned int iphivm =
          fpgaphi.bits(fpgaphi.nbits() - settings_.nbitsallstubs(layerdisk_) - settings_.nbitsvmme(layerdisk_),
                       settings_.nbitsvmme(layerdisk_));

      //This block of code just checks that the configuration is consistent
      if (lastTCID >= tracklet->TCID()) {
        edm::LogPrint("Tracklet") << "Wrong TCID ordering for projections in " << getName();
      } else {
        lastTCID = tracklet->TCID();
      }

      allproj_->addTracklet(tracklet);

      vmprojs_[iphivm]->addTracklet(tracklet, allprojcount);

      if (settings_.debugTracklet()) {
        edm::LogVerbatim("Tracklet") << getName() << " projection to " << vmprojs_[iphivm]->getName() << " iphivm "
                                     << iphivm;
      }

      allprojcount++;
    }
  }

  if (settings_.writeMonitorData("AP")) {
    globals_->ofstream("allprojections.txt") << getName() << " " << allproj_->nTracklets() << endl;
  }

  if (settings_.writeMonitorData("VMP")) {
    ofstream& out = globals_->ofstream("chisq.txt");
    for (unsigned int i = 0; i < 8; i++) {
      if (vmprojs_[i] != nullptr) {
        out << vmprojs_[i]->getName() << " " << vmprojs_[i]->nTracklets() << endl;
      }
    }
  }
}

Member Data Documentation

allproj_

AllProjectionsMemory* trklet::ProjectionRouter::allproj_

private

Definition at line 38 of file ProjectionRouter.h.

Referenced by addOutput(), and execute().

inputproj_

std::vector<TrackletProjectionsMemory*> trklet::ProjectionRouter::inputproj_

private

Definition at line 36 of file ProjectionRouter.h.

Referenced by addInput(), and execute().

layerdisk_

unsigned int trklet::ProjectionRouter::layerdisk_

private

Definition at line 28 of file ProjectionRouter.h.

Referenced by addOutput(), execute(), and ProjectionRouter().

nphiderbits_

int trklet::ProjectionRouter::nphiderbits_

private

Definition at line 31 of file ProjectionRouter.h.

Referenced by execute(), and ProjectionRouter().

nrbits_

int trklet::ProjectionRouter::nrbits_

private

Definition at line 30 of file ProjectionRouter.h.

Referenced by execute(), and ProjectionRouter().

rinvbendlut_

TrackletLUT trklet::ProjectionRouter::rinvbendlut_

private

Definition at line 34 of file ProjectionRouter.h.

Referenced by execute(), and ProjectionRouter().

vmprojs_

std::vector<VMProjectionsMemory*> trklet::ProjectionRouter::vmprojs_

private

Definition at line 39 of file ProjectionRouter.h.

Referenced by addOutput(), execute(), and ProjectionRouter().