Inheritance diagram for DTTrigPhase2Prod:

Public Member Functions
void	beginRun (edm::Run const &iRun, const edm::EventSetup &iEventSetup) override
	Create Trigger Units before starting event processing. More...

	DTTrigPhase2Prod (const edm::ParameterSet &pset)
	Constructor. More...

void	endRun (edm::Run const &iRun, const edm::EventSetup &iEventSetup) override
	endRun: finish things More...

MP_QUALITY	getMinimumQuality (void)

bool	hasPosRF (int wh, int sec) const

bool	inner (const metaPrimitive &mp) const

bool	outer (const metaPrimitive &mp) const

void	printmP (const std::string &ss, const metaPrimitive &mP) const

void	printmPC (const std::string &ss, const metaPrimitive &mP) const

void	produce (edm::Event &iEvent, const edm::EventSetup &iEventSetup) override
	Producer: process every event and generates trigger data. More...

int	rango (const metaPrimitive &mp) const

void	setChiSquareThreshold (float ch2Thr)

void	setMinimumQuality (MP_QUALITY q)

	~DTTrigPhase2Prod () override
	Destructor. More...

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Public Attributes
DTGeometry const *	dtGeo_

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	dtGeomH

std::vector< std::pair< int, MuonPath > >	primitives_

Private Types
typedef std::map< DTChamberId, DTDigiCollection, std::less< DTChamberId > >	DTDigiMap

typedef DTDigiMap::const_iterator	DTDigiMap_const_iterator

typedef DTDigiMap::iterator	DTDigiMap_iterator

Private Member Functions
void	assignIndex (std::vector< metaPrimitive > &inMPaths)

void	assignIndexPerBX (std::vector< metaPrimitive > &inMPaths)

int	assignQualityOrder (const metaPrimitive &mP) const

std::vector< DTDigiCollection * >	distribDigis (std::queue< std::pair< DTLayerId, DTDigi >> &inQ)

void	processDigi (std::queue< std::pair< DTLayerId, DTDigi >> &inQ, std::vector< std::queue< std::pair< DTLayerId, DTDigi >> * > &vec)

Private Attributes
bool	activateBuffer_

int	algo_

bool	debug_

bool	do_correlation_

double	dT0_correlate_TP_

edm::EDGetTokenT< DTDigiCollection >	dtDigisToken_

bool	dump_

std::unique_ptr< MotherGrouping >	grouping_obj_

std::unique_ptr< MuonPathAnalyzer >	mpathanalyzer_

std::unique_ptr< MuonPathAssociator >	mpathassociator_

std::unique_ptr< MPFilter >	mpathqualityenhancer_

std::unique_ptr< MPFilter >	mpathredundantfilter_

int	my_BXoffset_

bool	my_CCBValid_

const std::unordered_map< int, int >	qmap_

std::unique_ptr< RPCIntegrator >	rpc_integrator_

edm::EDGetTokenT< RPCRecHitCollection >	rpcRecHitsLabel_

int	scenario_

edm::FileInPath	shift_filename_

std::map< int, float >	shiftinfo_

int	superCellhalfspacewidth_

float	superCelltimewidth_

bool	useRPC_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 63 of file DTTrigPhase2Prod.cc.

Member Typedef Documentation

◆ DTDigiMap

typedef std::map<DTChamberId, DTDigiCollection, std::less<DTChamberId> > DTTrigPhase2Prod::DTDigiMap

private

Definition at line 64 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_const_iterator

typedef DTDigiMap::const_iterator DTTrigPhase2Prod::DTDigiMap_const_iterator

private

Definition at line 66 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_iterator

typedef DTDigiMap::iterator DTTrigPhase2Prod::DTDigiMap_iterator

private

Definition at line 65 of file DTTrigPhase2Prod.cc.

Constructor & Destructor Documentation

◆ DTTrigPhase2Prod()

DTTrigPhase2Prod::DTTrigPhase2Prod ( const edm::ParameterSet & pset )

Constructor.

Definition at line 161 of file DTTrigPhase2Prod.cc.

     : qmap_({{9, 9}, {8, 8}, {7, 6}, {6, 7}, {5, 3}, {4, 5}, {3, 4}, {2, 2}, {1, 1}}) {
   produces<L1Phase2MuDTPhContainer>();
  
   debug_ = pset.getUntrackedParameter<bool>("debug");
   dump_ = pset.getUntrackedParameter<bool>("dump");
  
   do_correlation_ = pset.getParameter<bool>("do_correlation");
   scenario_ = pset.getParameter<int>("scenario");
  
   dtDigisToken_ = consumes<DTDigiCollection>(pset.getParameter<edm::InputTag>("digiTag"));
  
   rpcRecHitsLabel_ = consumes<RPCRecHitCollection>(pset.getParameter<edm::InputTag>("rpcRecHits"));
   useRPC_ = pset.getParameter<bool>("useRPC");
  
   // Choosing grouping scheme:
   algo_ = pset.getParameter<int>("algo");
  
   edm::ConsumesCollector consumesColl(consumesCollector());
  
   if (algo_ == PseudoBayes) {
     grouping_obj_ =
         std::make_unique<PseudoBayesGrouping>(pset.getParameter<edm::ParameterSet>("PseudoBayesPattern"), consumesColl);
   } else if (algo_ == HoughTrans) {
     grouping_obj_ =
         std::make_unique<HoughGrouping>(pset.getParameter<edm::ParameterSet>("HoughGrouping"), consumesColl);
   } else {
     grouping_obj_ = std::make_unique<InitialGrouping>(pset, consumesColl);
   }
  
   if (algo_ == Standard) {
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "DTp2:constructor: JM analyzer";
     mpathanalyzer_ = std::make_unique<MuonPathAnalyzerPerSL>(pset, consumesColl);
   } else {
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "DTp2:constructor: Full chamber analyzer";
     mpathanalyzer_ = std::make_unique<MuonPathAnalyzerInChamber>(pset, consumesColl);
   }
  
   // Getting buffer option
   activateBuffer_ = pset.getParameter<bool>("activateBuffer");
   superCellhalfspacewidth_ = pset.getParameter<int>("superCellspacewidth") / 2;
   superCelltimewidth_ = pset.getParameter<double>("superCelltimewidth");
  
   mpathqualityenhancer_ = std::make_unique<MPQualityEnhancerFilter>(pset);
   mpathredundantfilter_ = std::make_unique<MPRedundantFilter>(pset);
   mpathassociator_ = std::make_unique<MuonPathAssociator>(pset, consumesColl);
   rpc_integrator_ = std::make_unique<RPCIntegrator>(pset, consumesColl);
  
   dtGeomH = esConsumes<DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
 }

◆ ~DTTrigPhase2Prod()

DTTrigPhase2Prod::~DTTrigPhase2Prod ( )

override

Destructor.

Definition at line 214 of file DTTrigPhase2Prod.cc.

                                     {
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "DTp2: calling destructor" << std::endl;
 }

References debug_, and LogDebug.

Member Function Documentation

◆ assignIndex()

void DTTrigPhase2Prod::assignIndex ( std::vector< metaPrimitive > & inMPaths )

private

Definition at line 608 of file DTTrigPhase2Prod.cc.

                                                                      {
   std::map<int, std::vector<metaPrimitive>> primsPerBX;
   for (const auto& metaPrimitive : inMPaths) {
     int BX = round(metaPrimitive.t0 / 25.);
     primsPerBX[BX].push_back(metaPrimitive);
   }
   inMPaths.clear();
   for (auto& prims : primsPerBX) {
     assignIndexPerBX(prims.second);
     for (const auto& primitive : prims.second)
       inMPaths.push_back(primitive);
   }
 }

References assignIndexPerBX(), L1TStage2uGTEmulatorClient_cff::BX, funct::primitive(), and cmsdt::metaPrimitive::t0.

Referenced by produce().

◆ assignIndexPerBX()

void DTTrigPhase2Prod::assignIndexPerBX ( std::vector< metaPrimitive > & inMPaths )

private

Definition at line 622 of file DTTrigPhase2Prod.cc.

                                                                           {
   // First we asociate a new index to the metaprimitive depending on quality or phiB;
   uint32_t rawId = -1;
   int numP = -1;
   for (auto& metaPrimitiveIt : inMPaths) {
     numP++;
     rawId = metaPrimitiveIt.rawId;
     int iOrder = assignQualityOrder(metaPrimitiveIt);
     int inf = 0;
     int numP2 = -1;
     for (auto& metaPrimitiveItN : inMPaths) {
       int nOrder = assignQualityOrder(metaPrimitiveItN);
       numP2++;
       if (rawId != metaPrimitiveItN.rawId)
         continue;
       if (numP2 == numP) {
         metaPrimitiveIt.index = inf;
         break;
       } else if (iOrder < nOrder) {
         inf++;
       } else if (iOrder > nOrder) {
         metaPrimitiveItN.index++;
       } else if (iOrder == nOrder) {
         if (std::abs(metaPrimitiveIt.phiB) >= std::abs(metaPrimitiveItN.phiB)) {
           inf++;
         } else if (std::abs(metaPrimitiveIt.phiB) < std::abs(metaPrimitiveItN.phiB)) {
           metaPrimitiveItN.index++;
         }
       }
     }  // ending second for
   }    // ending first for
 }

References funct::abs(), assignQualityOrder(), and dqmiodatasetharvest::inf.

Referenced by assignIndex().

◆ assignQualityOrder()

int DTTrigPhase2Prod::assignQualityOrder ( const metaPrimitive & mP ) const

private

Definition at line 655 of file DTTrigPhase2Prod.cc.

                                                                       {
   if (mP.quality > 9 || mP.quality < 1)
     return -1;
  
   return qmap_.find(mP.quality)->second;
 }

References qmap_, and cmsdt::metaPrimitive::quality.

Referenced by assignIndexPerBX().

◆ beginRun()

void DTTrigPhase2Prod::beginRun	(	edm::Run const &	iRun,
		const edm::EventSetup &	iEventSetup
	)

override

Create Trigger Units before starting event processing.

Definition at line 219 of file DTTrigPhase2Prod.cc.

                                                                                     {
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "beginRun " << iRun.id().run();
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "beginRun: getting DT geometry";
  
   grouping_obj_->initialise(iEventSetup);          // Grouping object initialisation
   mpathanalyzer_->initialise(iEventSetup);         // Analyzer object initialisation
   mpathqualityenhancer_->initialise(iEventSetup);  // Filter object initialisation
   mpathredundantfilter_->initialise(iEventSetup);  // Filter object initialisation
   mpathassociator_->initialise(iEventSetup);       // Associator object initialisation
  
   const MuonGeometryRecord& geom = iEventSetup.get<MuonGeometryRecord>();
   dtGeo_ = &geom.get(dtGeomH);
 }

References debug_, dtGeo_, dtGeomH, relativeConstraints::geom, edm::EventSetup::get(), grouping_obj_, edm::RunBase::id(), LogDebug, mpathanalyzer_, mpathassociator_, mpathqualityenhancer_, mpathredundantfilter_, and edm::RunID::run().

◆ distribDigis()

std::vector< DTDigiCollection * > DTTrigPhase2Prod::distribDigis ( std::queue< std::pair< DTLayerId, DTDigi >> & inQ )

private

Definition at line 662 of file DTTrigPhase2Prod.cc.

                                                                                                      {
   std::vector<std::queue<std::pair<DTLayerId, DTDigi>>*> tmpVector;
   tmpVector.clear();
   std::vector<DTDigiCollection*> collVector;
   collVector.clear();
   while (!inQ.empty()) {
     processDigi(inQ, tmpVector);
   }
   for (auto& sQ : tmpVector) {
     DTDigiCollection tmpColl;
     while (!sQ->empty()) {
       tmpColl.insertDigi((sQ->front().first), (sQ->front().second));
       sQ->pop();
     }
     collVector.push_back(&tmpColl);
   }
   return collVector;
 }

References processDigi().

Referenced by produce().

◆ endRun()

void DTTrigPhase2Prod::endRun	(	edm::Run const &	iRun,
		const edm::EventSetup &	iEventSetup
	)

override

endRun: finish things

Definition at line 552 of file DTTrigPhase2Prod.cc.

                                                                                   {
   grouping_obj_->finish();
   mpathanalyzer_->finish();
   mpathqualityenhancer_->finish();
   mpathredundantfilter_->finish();
   mpathassociator_->finish();
   rpc_integrator_->finish();
 };

References grouping_obj_, mpathanalyzer_, mpathassociator_, mpathqualityenhancer_, mpathredundantfilter_, and rpc_integrator_.

◆ getMinimumQuality()

MP_QUALITY DTTrigPhase2Prod::getMinimumQuality ( void )

◆ hasPosRF()

bool DTTrigPhase2Prod::hasPosRF	(	int	wh,
		int	sec
	)		const

Definition at line 571 of file DTTrigPhase2Prod.cc.

571 { return wh > 0 || (wh == 0 && sec % 4 > 1); }

References fileinputsource_cfi::sec.

◆ inner()

bool DTTrigPhase2Prod::inner ( const metaPrimitive & mp ) const

Definition at line 566 of file DTTrigPhase2Prod.cc.

                                                           {
   int counter = (mp.wi1 != -1) + (mp.wi2 != -1) + (mp.wi3 != -1) + (mp.wi4 != -1);
   return (counter > 2);
 }

References cmsdt::metaPrimitive::wi1, cmsdt::metaPrimitive::wi2, cmsdt::metaPrimitive::wi3, and cmsdt::metaPrimitive::wi4.

Referenced by produce().

◆ outer()

bool DTTrigPhase2Prod::outer ( const metaPrimitive & mp ) const

Definition at line 561 of file DTTrigPhase2Prod.cc.

                                                           {
   int counter = (mp.wi5 != -1) + (mp.wi6 != -1) + (mp.wi7 != -1) + (mp.wi8 != -1);
   return (counter > 2);
 }

References cmsdt::metaPrimitive::wi5, cmsdt::metaPrimitive::wi6, cmsdt::metaPrimitive::wi7, and cmsdt::metaPrimitive::wi8.

◆ printmP()

void DTTrigPhase2Prod::printmP	(	const std::string &	ss,
		const metaPrimitive &	mP
	)		const

Definition at line 573 of file DTTrigPhase2Prod.cc.

                                                                               {
   DTSuperLayerId slId(mP.rawId);
   LogInfo("DTTrigPhase2Prod") << ss << (int)slId << "\t " << setw(2) << left << mP.wi1 << " " << setw(2) << left
                               << mP.wi2 << " " << setw(2) << left << mP.wi3 << " " << setw(2) << left << mP.wi4 << " "
                               << setw(5) << left << mP.tdc1 << " " << setw(5) << left << mP.tdc2 << " " << setw(5)
                               << left << mP.tdc3 << " " << setw(5) << left << mP.tdc4 << " " << setw(10) << right
                               << mP.x << " " << setw(9) << left << mP.tanPhi << " " << setw(5) << left << mP.t0 << " "
                               << setw(13) << left << mP.chi2 << " r:" << rango(mP);
 }

References cmsdt::metaPrimitive::chi2, createfilelist::int, rango(), cmsdt::metaPrimitive::rawId, contentValuesCheck::ss, cmsdt::metaPrimitive::t0, cmsdt::metaPrimitive::tanPhi, cmsdt::metaPrimitive::tdc1, cmsdt::metaPrimitive::tdc2, cmsdt::metaPrimitive::tdc3, cmsdt::metaPrimitive::tdc4, cmsdt::metaPrimitive::wi1, cmsdt::metaPrimitive::wi2, cmsdt::metaPrimitive::wi3, cmsdt::metaPrimitive::wi4, and cmsdt::metaPrimitive::x.

Referenced by produce().

◆ printmPC()

void DTTrigPhase2Prod::printmPC	(	const std::string &	ss,
		const metaPrimitive &	mP
	)		const

Definition at line 583 of file DTTrigPhase2Prod.cc.

                                                                                {
   DTChamberId ChId(mP.rawId);
   LogInfo("DTTrigPhase2Prod") << ss << (int)ChId << "\t  " << setw(2) << left << mP.wi1 << " " << setw(2) << left
                               << mP.wi2 << " " << setw(2) << left << mP.wi3 << " " << setw(2) << left << mP.wi4 << " "
                               << setw(2) << left << mP.wi5 << " " << setw(2) << left << mP.wi6 << " " << setw(2) << left
                               << mP.wi7 << " " << setw(2) << left << mP.wi8 << " " << setw(5) << left << mP.tdc1 << " "
                               << setw(5) << left << mP.tdc2 << " " << setw(5) << left << mP.tdc3 << " " << setw(5)
                               << left << mP.tdc4 << " " << setw(5) << left << mP.tdc5 << " " << setw(5) << left
                               << mP.tdc6 << " " << setw(5) << left << mP.tdc7 << " " << setw(5) << left << mP.tdc8
                               << " " << setw(2) << left << mP.lat1 << " " << setw(2) << left << mP.lat2 << " "
                               << setw(2) << left << mP.lat3 << " " << setw(2) << left << mP.lat4 << " " << setw(2)
                               << left << mP.lat5 << " " << setw(2) << left << mP.lat6 << " " << setw(2) << left
                               << mP.lat7 << " " << setw(2) << left << mP.lat8 << " " << setw(10) << right << mP.x << " "
                               << setw(9) << left << mP.tanPhi << " " << setw(5) << left << mP.t0 << " " << setw(13)
                               << left << mP.chi2 << " r:" << rango(mP);
 }

Referenced by produce().

◆ processDigi()

void DTTrigPhase2Prod::processDigi	(	std::queue< std::pair< DTLayerId, DTDigi >> &	inQ,
		std::vector< std::queue< std::pair< DTLayerId, DTDigi >> * > &	vec
	)

private

Definition at line 681 of file DTTrigPhase2Prod.cc.

                                                                                           {
   bool classified = false;
   if (!vec.empty()) {
     for (auto& sC : vec) {  // Conditions for entering a super cell.
       if ((sC->front().second.time() + superCelltimewidth_) > inQ.front().second.time()) {
         // Time requirement
         if (TMath::Abs(sC->front().second.wire() - inQ.front().second.wire()) <= superCellhalfspacewidth_) {
           // Spatial requirement
           sC->push(std::move(inQ.front()));
           classified = true;
         }
       }
     }
   }
   if (classified) {
     inQ.pop();
     return;
   }
  
   std::queue<std::pair<DTLayerId, DTDigi>> newQueue;
  
   std::pair<DTLayerId, DTDigi> tmpPair;
   tmpPair = std::move(inQ.front());
   newQueue.push(tmpPair);
   inQ.pop();
  
   vec.push_back(&newQueue);
 }

References Abs(), eostools::move(), superCellhalfspacewidth_, and superCelltimewidth_.

Referenced by distribDigis().

◆ produce()

void DTTrigPhase2Prod::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iEventSetup
	)

override

Producer: process every event and generates trigger data.

FITTING SECTION;

STORING RESULTs

Definition at line 235 of file DTTrigPhase2Prod.cc.

                                                                            {
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "produce";
   edm::Handle<DTDigiCollection> dtdigis;
   iEvent.getByToken(dtDigisToken_, dtdigis);
  
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "\t Getting the RPC RecHits" << std::endl;
   edm::Handle<RPCRecHitCollection> rpcRecHits;
   iEvent.getByToken(rpcRecHitsLabel_, rpcRecHits);
  
   // GROUPING CODE:
   DTDigiMap digiMap;
   DTDigiCollection::DigiRangeIterator detUnitIt;
   for (const auto& detUnitIt : *dtdigis) {
     const DTLayerId& layId = detUnitIt.first;
     const DTChamberId chambId = layId.superlayerId().chamberId();
     const DTDigiCollection::Range& range = detUnitIt.second;
     digiMap[chambId].put(range, layId);
   }
  
   // generate a list muon paths for each event!!!
   if (debug_ && activateBuffer_)
     LogDebug("DTTrigPhase2Prod") << "produce - Getting and grouping digis per chamber using a buffer and super cells.";
   else if (debug_)
     LogDebug("DTTrigPhase2Prod") << "produce - Getting and grouping digis per chamber.";
  
   MuonPathPtrs muonpaths;
   for (const auto& ich : dtGeo_->chambers()) {
     // The code inside this for loop would ideally later fit inside a trigger unit (in principle, a DT station) of the future Phase 2 DT Trigger.
     const DTChamber* chamb = ich;
     DTChamberId chid = chamb->id();
     DTDigiMap_iterator dmit = digiMap.find(chid);
  
     if (dmit == digiMap.end())
       continue;
  
     if (activateBuffer_) {  // Use buffering (per chamber) or not
       // Import digis from the station
       std::vector<std::pair<DTLayerId, DTDigi>> tmpvec;
       tmpvec.clear();
  
       for (const auto& dtLayerIdIt : (*dmit).second) {
         for (DTDigiCollection::const_iterator digiIt = (dtLayerIdIt.second).first;
              digiIt != (dtLayerIdIt.second).second;
              digiIt++) {
           tmpvec.emplace_back(dtLayerIdIt.first, *digiIt);
         }
       }
  
       // Check to enhance CPU time usage
       if (tmpvec.empty())
         continue;
  
       // Order digis depending on TDC time and insert them into a queue (FIFO buffer). TODO: adapt for MC simulations.
       std::sort(tmpvec.begin(), tmpvec.end(), DigiTimeOrdering);
       std::queue<std::pair<DTLayerId, DTDigi>> timequeue;
  
       for (const auto& elem : tmpvec)
         timequeue.emplace(elem);
       tmpvec.clear();
  
       // Distribute the digis from the queue into supercells
       std::vector<DTDigiCollection*> superCells;
       superCells = distribDigis(timequeue);
  
       // Process each supercell & collect the resulting muonpaths (as the muonpaths std::vector is only enlarged each time
       // the groupings access it, it's not needed to "collect" the final products).
       while (!superCells.empty()) {
         grouping_obj_->run(iEvent, iEventSetup, *(superCells.back()), muonpaths);
         superCells.pop_back();
       }
     } else {
       grouping_obj_->run(iEvent, iEventSetup, (*dmit).second, muonpaths);
     }
   }
   digiMap.clear();
  
   if (dump_) {
     for (unsigned int i = 0; i < muonpaths.size(); i++) {
       stringstream ss;
       ss << iEvent.id().event() << "      mpath " << i << ": ";
       for (int lay = 0; lay < muonpaths.at(i)->nprimitives(); lay++)
         ss << muonpaths.at(i)->primitive(lay)->channelId() << " ";
       for (int lay = 0; lay < muonpaths.at(i)->nprimitives(); lay++)
         ss << muonpaths.at(i)->primitive(lay)->tdcTimeStamp() << " ";
       for (int lay = 0; lay < muonpaths.at(i)->nprimitives(); lay++)
         ss << muonpaths.at(i)->primitive(lay)->laterality() << " ";
       LogInfo("DTTrigPhase2Prod") << ss.str();
     }
   }
  
   // FILTER GROUPING
   MuonPathPtrs filteredmuonpaths;
   if (algo_ == Standard) {
     mpathredundantfilter_->run(iEvent, iEventSetup, muonpaths, filteredmuonpaths);
   }
  
   if (dump_) {
     for (unsigned int i = 0; i < filteredmuonpaths.size(); i++) {
       stringstream ss;
       ss << iEvent.id().event() << " filt. mpath " << i << ": ";
       for (int lay = 0; lay < filteredmuonpaths.at(i)->nprimitives(); lay++)
         ss << filteredmuonpaths.at(i)->primitive(lay)->channelId() << " ";
       for (int lay = 0; lay < filteredmuonpaths.at(i)->nprimitives(); lay++)
         ss << filteredmuonpaths.at(i)->primitive(lay)->tdcTimeStamp() << " ";
       LogInfo("DTTrigPhase2Prod") << ss.str();
     }
   }
  
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "MUON PATHS found: " << muonpaths.size() << " (" << filteredmuonpaths.size()
                                  << ") in event " << iEvent.id().event();
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "filling NmetaPrimtives" << std::endl;
   std::vector<metaPrimitive> metaPrimitives;
   MuonPathPtrs outmpaths;
   if (algo_ == Standard) {
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "Fitting 1SL ";
     mpathanalyzer_->run(iEvent, iEventSetup, filteredmuonpaths, metaPrimitives);
   } else {
     // implementation for advanced (2SL) grouping, no filter required..
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "Fitting 2SL at once ";
     mpathanalyzer_->run(iEvent, iEventSetup, muonpaths, outmpaths);
   }
  
   if (dump_) {
     for (unsigned int i = 0; i < outmpaths.size(); i++) {
       LogInfo("DTTrigPhase2Prod") << iEvent.id().event() << " mp " << i << ": " << outmpaths.at(i)->bxTimeValue() << " "
                                   << outmpaths.at(i)->horizPos() << " " << outmpaths.at(i)->tanPhi() << " "
                                   << outmpaths.at(i)->phi() << " " << outmpaths.at(i)->phiB() << " "
                                   << outmpaths.at(i)->quality() << " " << outmpaths.at(i)->chiSquare();
     }
     for (unsigned int i = 0; i < metaPrimitives.size(); i++) {
       stringstream ss;
       ss << iEvent.id().event() << " mp " << i << ": ";
       printmP(ss.str(), metaPrimitives.at(i));
     }
   }
  
   muonpaths.clear();
   filteredmuonpaths.clear();
  
   //  FILTER SECTIONS:
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "declaring new vector for filtered" << std::endl;
  
   std::vector<metaPrimitive> filteredMetaPrimitives;
   if (algo_ == Standard)
     mpathqualityenhancer_->run(iEvent, iEventSetup, metaPrimitives, filteredMetaPrimitives);
  
   if (dump_) {
     for (unsigned int i = 0; i < filteredMetaPrimitives.size(); i++) {
       stringstream ss;
       ss << iEvent.id().event() << " filtered mp " << i << ": ";
       printmP(ss.str(), filteredMetaPrimitives.at(i));
     }
   }
  
   metaPrimitives.clear();
   metaPrimitives.erase(metaPrimitives.begin(), metaPrimitives.end());
  
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "DTp2 in event:" << iEvent.id().event() << " we found "
                                  << filteredMetaPrimitives.size() << " filteredMetaPrimitives (superlayer)"
                                  << std::endl;
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "filteredMetaPrimitives: starting correlations" << std::endl;
  
   std::vector<metaPrimitive> correlatedMetaPrimitives;
   if (algo_ == Standard)
     mpathassociator_->run(iEvent, iEventSetup, dtdigis, filteredMetaPrimitives, correlatedMetaPrimitives);
   else {
     for (const auto& muonpath : outmpaths) {
       correlatedMetaPrimitives.emplace_back(muonpath->rawId(),
                                             (double)muonpath->bxTimeValue(),
                                             muonpath->horizPos(),
                                             muonpath->tanPhi(),
                                             muonpath->phi(),
                                             muonpath->phiB(),
                                             muonpath->chiSquare(),
                                             (int)muonpath->quality(),
                                             muonpath->primitive(0)->channelId(),
                                             muonpath->primitive(0)->tdcTimeStamp(),
                                             muonpath->primitive(0)->laterality(),
                                             muonpath->primitive(1)->channelId(),
                                             muonpath->primitive(1)->tdcTimeStamp(),
                                             muonpath->primitive(1)->laterality(),
                                             muonpath->primitive(2)->channelId(),
                                             muonpath->primitive(2)->tdcTimeStamp(),
                                             muonpath->primitive(2)->laterality(),
                                             muonpath->primitive(3)->channelId(),
                                             muonpath->primitive(3)->tdcTimeStamp(),
                                             muonpath->primitive(3)->laterality(),
                                             muonpath->primitive(4)->channelId(),
                                             muonpath->primitive(4)->tdcTimeStamp(),
                                             muonpath->primitive(4)->laterality(),
                                             muonpath->primitive(5)->channelId(),
                                             muonpath->primitive(5)->tdcTimeStamp(),
                                             muonpath->primitive(5)->laterality(),
                                             muonpath->primitive(6)->channelId(),
                                             muonpath->primitive(6)->tdcTimeStamp(),
                                             muonpath->primitive(6)->laterality(),
                                             muonpath->primitive(7)->channelId(),
                                             muonpath->primitive(7)->tdcTimeStamp(),
                                             muonpath->primitive(7)->laterality());
     }
   }
   filteredMetaPrimitives.clear();
  
   if (debug_)
     LogDebug("DTTrigPhase2Prod") << "DTp2 in event:" << iEvent.id().event() << " we found "
                                  << correlatedMetaPrimitives.size() << " correlatedMetPrimitives (chamber)";
  
   if (dump_) {
     LogInfo("DTTrigPhase2Prod") << "DTp2 in event:" << iEvent.id().event() << " we found "
                                 << correlatedMetaPrimitives.size() << " correlatedMetPrimitives (chamber)";
  
     for (unsigned int i = 0; i < correlatedMetaPrimitives.size(); i++) {
       stringstream ss;
       ss << iEvent.id().event() << " correlated mp " << i << ": ";
       printmPC(ss.str(), correlatedMetaPrimitives.at(i));
     }
   }
  
   double shift_back = 0;
   if (scenario_ == MC)  //scope for MC
     shift_back = 400;
   else if (scenario_ == DATA)  //scope for data
     shift_back = 0;
   else if (scenario_ == SLICE_TEST)  //scope for slice test
     shift_back = 0;
  
   // RPC integration
   if (useRPC_) {
     rpc_integrator_->initialise(iEventSetup, shift_back);
     rpc_integrator_->prepareMetaPrimitives(rpcRecHits);
     rpc_integrator_->matchWithDTAndUseRPCTime(correlatedMetaPrimitives);
     rpc_integrator_->makeRPCOnlySegments();
     rpc_integrator_->storeRPCSingleHits();
     rpc_integrator_->removeRPCHitsUsed();
   }
  
   vector<L1Phase2MuDTPhDigi> outP2Ph;
  
   // Assigning index value
   assignIndex(correlatedMetaPrimitives);
   for (const auto& metaPrimitiveIt : correlatedMetaPrimitives) {
     DTChamberId chId(metaPrimitiveIt.rawId);
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "looping in final vector: SuperLayerId" << chId << " x=" << metaPrimitiveIt.x
                                    << " quality=" << metaPrimitiveIt.quality
                                    << " BX=" << round(metaPrimitiveIt.t0 / 25.) << " index=" << metaPrimitiveIt.index;
  
     int sectorTP = chId.sector();
     //sectors 13 and 14 exist only for the outermost stations for sectors 4 and 10 respectively
     //due to the larger MB4 that are divided into two.
     if (sectorTP == 13)
       sectorTP = 4;
     if (sectorTP == 14)
       sectorTP = 10;
     sectorTP = sectorTP - 1;
     int sl = 0;
     if (metaPrimitiveIt.quality < LOWLOWQ || metaPrimitiveIt.quality == CHIGHQ) {
       if (inner(metaPrimitiveIt))
         sl = 1;
       else
         sl = 3;
     }
  
     if (debug_)
       LogDebug("DTTrigPhase2Prod") << "pushing back phase-2 dataformat carlo-federica dataformat";
     outP2Ph.push_back(L1Phase2MuDTPhDigi(
         (int)round(metaPrimitiveIt.t0 / (float)LHC_CLK_FREQ) - shift_back,
         chId.wheel(),                                                // uwh (m_wheel)
         sectorTP,                                                    // usc (m_sector)
         chId.station(),                                              // ust (m_station)
         sl,                                                          // ust (m_station)
         (int)round(metaPrimitiveIt.phi * PHIRES_CONV),               // uphi (_phiAngle)
         (int)round(metaPrimitiveIt.phiB * PHIBRES_CONV),             // uphib (m_phiBending)
         metaPrimitiveIt.quality,                                     // uqua (m_qualityCode)
         metaPrimitiveIt.index,                                       // uind (m_segmentIndex)
         (int)round(metaPrimitiveIt.t0) - shift_back * LHC_CLK_FREQ,  // ut0 (m_t0Segment)
         (int)round(metaPrimitiveIt.chi2 * CHI2RES_CONV),             // uchi2 (m_chi2Segment)
         metaPrimitiveIt.rpcFlag                                      // urpc (m_rpcFlag)
         ));
   }
  
   // Storing RPC hits that were not used elsewhere
   if (useRPC_) {
     for (auto rpc_dt_digi = rpc_integrator_->rpcRecHits_translated_.begin();
          rpc_dt_digi != rpc_integrator_->rpcRecHits_translated_.end();
          rpc_dt_digi++) {
       outP2Ph.push_back(*rpc_dt_digi);
     }
   }
  
   auto resultP2Ph = std::make_unique<L1Phase2MuDTPhContainer>();
   resultP2Ph->setContainer(outP2Ph);
   iEvent.put(std::move(resultP2Ph));
   outP2Ph.clear();
   outP2Ph.erase(outP2Ph.begin(), outP2Ph.end());
 }

References activateBuffer_, algo_, assignIndex(), DTSuperLayerId::chamberId(), DTGeometry::chambers(), cmsdt::CHI2RES_CONV, cmsdt::CHIGHQ, cscdqm::DATA, debug_, distribDigis(), dtDigisToken_, dtGeo_, dump_, grouping_obj_, mps_fire::i, DTChamber::id(), iEvent, inner(), createfilelist::int, cmsdt::LHC_CLK_FREQ, LogDebug, cmsdt::LOWLOWQ, cmsdt::MC, eostools::move(), mpathanalyzer_, mpathassociator_, mpathqualityenhancer_, mpathredundantfilter_, cmsdt::PHIBRES_CONV, cmsdt::PHIRES_CONV, printmP(), printmPC(), FastTimerService_cff::range, rpc_integrator_, dtTriggerPhase2PrimitiveDigis_cfi::rpcRecHits, rpcRecHitsLabel_, scenario_, DTChamberId::sector(), cmsdt::SLICE_TEST, contentValuesCheck::ss, cmsdt::Standard, DTChamberId::station(), DTLayerId::superlayerId(), useRPC_, and DTChamberId::wheel().

◆ rango()

int DTTrigPhase2Prod::rango ( const metaPrimitive & mp ) const

Definition at line 600 of file DTTrigPhase2Prod.cc.

                                                          {
   if (mp.quality == 1 or mp.quality == 2)
     return 3;
   if (mp.quality == 3 or mp.quality == 4)
     return 4;
   return mp.quality;
 }

References or, and cmsdt::metaPrimitive::quality.

Referenced by printmP(), and printmPC().

◆ setChiSquareThreshold()

void DTTrigPhase2Prod::setChiSquareThreshold ( float ch2Thr )

◆ setMinimumQuality()

void DTTrigPhase2Prod::setMinimumQuality ( MP_QUALITY q )

Member Data Documentation

◆ activateBuffer_

bool DTTrigPhase2Prod::activateBuffer_

private

Definition at line 135 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ algo_

int DTTrigPhase2Prod::algo_

private

Definition at line 127 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ debug_

bool DTTrigPhase2Prod::debug_

private

Definition at line 112 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), produce(), and ~DTTrigPhase2Prod().

◆ do_correlation_

bool DTTrigPhase2Prod::do_correlation_

private

Definition at line 115 of file DTTrigPhase2Prod.cc.

◆ dT0_correlate_TP_

double DTTrigPhase2Prod::dT0_correlate_TP_

private

Definition at line 114 of file DTTrigPhase2Prod.cc.

◆ dtDigisToken_

edm::EDGetTokenT<DTDigiCollection> DTTrigPhase2Prod::dtDigisToken_

private

Definition at line 123 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ dtGeo_

DTGeometry const* DTTrigPhase2Prod::dtGeo_

Definition at line 100 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), and produce().

◆ dtGeomH

edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTTrigPhase2Prod::dtGeomH

Definition at line 101 of file DTTrigPhase2Prod.cc.

Referenced by beginRun().

◆ dump_

bool DTTrigPhase2Prod::dump_

private

Definition at line 113 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ grouping_obj_

std::unique_ptr<MotherGrouping> DTTrigPhase2Prod::grouping_obj_

private

Definition at line 128 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), endRun(), and produce().

◆ mpathanalyzer_

std::unique_ptr<MuonPathAnalyzer> DTTrigPhase2Prod::mpathanalyzer_

private

Definition at line 129 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), endRun(), and produce().

◆ mpathassociator_

std::unique_ptr<MuonPathAssociator> DTTrigPhase2Prod::mpathassociator_

private

Definition at line 132 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), endRun(), and produce().

◆ mpathqualityenhancer_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathqualityenhancer_

private

Definition at line 130 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), endRun(), and produce().

◆ mpathredundantfilter_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathredundantfilter_

private

Definition at line 131 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), endRun(), and produce().

◆ my_BXoffset_

int DTTrigPhase2Prod::my_BXoffset_

private

Definition at line 109 of file DTTrigPhase2Prod.cc.

◆ my_CCBValid_

bool DTTrigPhase2Prod::my_CCBValid_

private

Definition at line 106 of file DTTrigPhase2Prod.cc.

◆ primitives_

std::vector<std::pair<int, MuonPath> > DTTrigPhase2Prod::primitives_

Definition at line 102 of file DTTrigPhase2Prod.cc.

◆ qmap_

const std::unordered_map<int, int> DTTrigPhase2Prod::qmap_

private

Definition at line 150 of file DTTrigPhase2Prod.cc.

Referenced by assignQualityOrder().

◆ rpc_integrator_

std::unique_ptr<RPCIntegrator> DTTrigPhase2Prod::rpc_integrator_

private

Definition at line 143 of file DTTrigPhase2Prod.cc.

Referenced by endRun(), and produce().

◆ rpcRecHitsLabel_

edm::EDGetTokenT<RPCRecHitCollection> DTTrigPhase2Prod::rpcRecHitsLabel_

private

Definition at line 124 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ scenario_

int DTTrigPhase2Prod::scenario_

private

Definition at line 116 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ shift_filename_

edm::FileInPath DTTrigPhase2Prod::shift_filename_

private

Definition at line 119 of file DTTrigPhase2Prod.cc.

◆ shiftinfo_

std::map<int, float> DTTrigPhase2Prod::shiftinfo_

private

Definition at line 120 of file DTTrigPhase2Prod.cc.

◆ superCellhalfspacewidth_

int DTTrigPhase2Prod::superCellhalfspacewidth_

private

Definition at line 136 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ superCelltimewidth_

float DTTrigPhase2Prod::superCelltimewidth_

private

Definition at line 137 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ useRPC_

bool DTTrigPhase2Prod::useRPC_

private

Definition at line 144 of file DTTrigPhase2Prod.cc.

Referenced by produce().

Public Member Functions

Public Attributes

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DTDigiMap

◆ DTDigiMap_const_iterator

◆ DTDigiMap_iterator

Constructor & Destructor Documentation

◆ DTTrigPhase2Prod()

◆ ~DTTrigPhase2Prod()

Member Function Documentation

◆ assignIndex()

◆ assignIndexPerBX()

◆ assignQualityOrder()

◆ beginRun()

◆ distribDigis()

◆ endRun()

◆ getMinimumQuality()

◆ hasPosRF()

◆ inner()

◆ outer()

◆ printmP()

◆ printmPC()

◆ processDigi()

◆ produce()

◆ rango()

◆ setChiSquareThreshold()

◆ setMinimumQuality()

Member Data Documentation

◆ activateBuffer_

◆ algo_

◆ debug_

◆ do_correlation_

◆ dT0_correlate_TP_

◆ dtDigisToken_

◆ dtGeo_

◆ dtGeomH

◆ dump_

◆ grouping_obj_

◆ mpathanalyzer_

◆ mpathassociator_

◆ mpathqualityenhancer_

◆ mpathredundantfilter_

◆ my_BXoffset_

◆ my_CCBValid_

◆ primitives_

◆ qmap_

◆ rpc_integrator_

◆ rpcRecHitsLabel_

◆ scenario_

◆ shift_filename_

◆ shiftinfo_

◆ superCellhalfspacewidth_

◆ superCelltimewidth_

◆ useRPC_