#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	(	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_);
   }
 }

trklet::ProjectionRouter::~ProjectionRouter ( )

overridedefault

Member Function Documentation

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, trklet::MemoryBase::getName(), input, inputproj_, 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;
 }

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, trklet::MemoryBase::getName(), layerdisk_, trklet::Settings::nallstubs(), mergeVDriftHistosByStation::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;
 }

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::Projection::fpgaphiprojder(), trklet::Projection::fpgarzproj(), trklet::Projection::fpgarzprojder(), 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(), submitPVResolutionJobs::out, trklet::Tracklet::proj(), rinvbendlut_, trklet::Projection::setBendIndex(), trklet::ProcessBase::settings_, trklet::Tracklet::TCID(), trklet::FPGAWord::value(), 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;
       }
     }
   }
 }