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	cmssw_for_global_

bool	debug_

bool	do_correlation_

double	dT0_correlate_TP_

edm::EDGetTokenT< DTDigiCollection >	dtDigisToken_

bool	dump_

std::string	geometry_tag_

std::shared_ptr< GlobalCoordsObtainer >	globalcoordsobtainer_

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_

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 66 of file DTTrigPhase2Prod.cc.

Member Typedef Documentation

◆ DTDigiMap

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

private

Definition at line 67 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_const_iterator

typedef DTDigiMap::const_iterator DTTrigPhase2Prod::DTDigiMap_const_iterator

private

Definition at line 69 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_iterator

typedef DTDigiMap::iterator DTTrigPhase2Prod::DTDigiMap_iterator

private

Definition at line 68 of file DTTrigPhase2Prod.cc.

Constructor & Destructor Documentation

◆ DTTrigPhase2Prod()

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

Constructor.

Definition at line 163 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>();
   produces<L1Phase2MuDTThContainer>();
  
   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");
  
   // Local to global coordinates approach
   cmssw_for_global_ = pset.getUntrackedParameter<bool>("cmssw_for_global", true);
   geometry_tag_ = pset.getUntrackedParameter<std::string>("geometry_tag", "");
  
   edm::ConsumesCollector consumesColl(consumesCollector());
   globalcoordsobtainer_ = std::make_shared<GlobalCoordsObtainer>(pset);
   if (!cmssw_for_global_)
     globalcoordsobtainer_->generate_luts();
  
   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<MuonPathAnalyticAnalyzer>(pset, consumesColl, globalcoordsobtainer_);
   } 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, globalcoordsobtainer_);
   rpc_integrator_ = std::make_unique<RPCIntegrator>(pset, consumesColl);
  
   dtGeomH = esConsumes<DTGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
 }

◆ ~DTTrigPhase2Prod()

DTTrigPhase2Prod::~DTTrigPhase2Prod ( )

override

Destructor.

Definition at line 224 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 646 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 660 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 693 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 229 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
  
   edm::ESHandle<DTGeometry> geom;
   iEventSetup.get<MuonGeometryRecord>().get(geometry_tag_, geom);
   dtGeo_ = &(*geom);
 }

References debug_, dtGeo_, relativeConstraints::geom, geometry_tag_, edm::EventSetup::get(), 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 700 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 590 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 609 of file DTTrigPhase2Prod.cc.

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

References fileinputsource_cfi::sec.

◆ inner()

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

Definition at line 604 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 599 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 611 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 621 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 719 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 246 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 = 400;
  
   // 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;
   vector<L1Phase2MuDTThDigi> outP2Th;
  
   // Assigning index value
   assignIndex(correlatedMetaPrimitives);
   for (const auto& metaPrimitiveIt : correlatedMetaPrimitives) {
     DTChamberId chId(metaPrimitiveIt.rawId);
     DTSuperLayerId slId(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";
  
     if (slId.superLayer() != 2) {
       //phiTP
       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)
           ));
     } else {
       //thTP
       outP2Th.push_back(L1Phase2MuDTThDigi(
           (int)round(metaPrimitiveIt.t0 / (float)LHC_CLK_FREQ) - shift_back,
           chId.wheel(),                                                // uwh (m_wheel)
           sectorTP,                                                    // usc (m_sector)
           chId.station(),                                              // ust (m_station)
           (int)round(metaPrimitiveIt.phi * ZRES_CONV),                 // uz (m_zGlobal)
           (int)round(metaPrimitiveIt.phiB * KRES_CONV),                // uk (m_kSlope)
           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());
  
   std::unique_ptr<L1Phase2MuDTThContainer> resultP2Th(new L1Phase2MuDTThContainer);
   resultP2Th->setContainer(outP2Th);
   iEvent.put(std::move(resultP2Th));
   outP2Th.clear();
   outP2Th.erase(outP2Th.begin(), outP2Th.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::KRES_CONV, 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, jetUpdater_cfi::sort, contentValuesCheck::ss, cmsdt::Standard, DTChamberId::station(), DTSuperLayerId::superLayer(), DTLayerId::superlayerId(), useRPC_, DTChamberId::wheel(), and cmsdt::ZRES_CONV.

◆ rango()

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

Definition at line 638 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 137 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ algo_

int DTTrigPhase2Prod::algo_

private

Definition at line 128 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ cmssw_for_global_

bool DTTrigPhase2Prod::cmssw_for_global_

private

Definition at line 120 of file DTTrigPhase2Prod.cc.

◆ debug_

bool DTTrigPhase2Prod::debug_

private

Definition at line 115 of file DTTrigPhase2Prod.cc.

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

◆ do_correlation_

bool DTTrigPhase2Prod::do_correlation_

private

Definition at line 118 of file DTTrigPhase2Prod.cc.

◆ dT0_correlate_TP_

double DTTrigPhase2Prod::dT0_correlate_TP_

private

Definition at line 117 of file DTTrigPhase2Prod.cc.

◆ dtDigisToken_

edm::EDGetTokenT<DTDigiCollection> DTTrigPhase2Prod::dtDigisToken_

private

Definition at line 124 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ dtGeo_

DTGeometry const* DTTrigPhase2Prod::dtGeo_

Definition at line 103 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), and produce().

◆ dtGeomH

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

Definition at line 104 of file DTTrigPhase2Prod.cc.

◆ dump_

bool DTTrigPhase2Prod::dump_

private

Definition at line 116 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ geometry_tag_

std::string DTTrigPhase2Prod::geometry_tag_

private

Definition at line 121 of file DTTrigPhase2Prod.cc.

Referenced by beginRun().

◆ globalcoordsobtainer_

std::shared_ptr<GlobalCoordsObtainer> DTTrigPhase2Prod::globalcoordsobtainer_

private

Definition at line 134 of file DTTrigPhase2Prod.cc.

◆ grouping_obj_

std::unique_ptr<MotherGrouping> DTTrigPhase2Prod::grouping_obj_

private

Definition at line 129 of file DTTrigPhase2Prod.cc.

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

◆ mpathanalyzer_

std::unique_ptr<MuonPathAnalyzer> DTTrigPhase2Prod::mpathanalyzer_

private

Definition at line 130 of file DTTrigPhase2Prod.cc.

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

◆ mpathassociator_

std::unique_ptr<MuonPathAssociator> DTTrigPhase2Prod::mpathassociator_

private

Definition at line 133 of file DTTrigPhase2Prod.cc.

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

◆ mpathqualityenhancer_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathqualityenhancer_

private

Definition at line 131 of file DTTrigPhase2Prod.cc.

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

◆ mpathredundantfilter_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathredundantfilter_

private

Definition at line 132 of file DTTrigPhase2Prod.cc.

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

◆ my_BXoffset_

int DTTrigPhase2Prod::my_BXoffset_

private

Definition at line 112 of file DTTrigPhase2Prod.cc.

◆ my_CCBValid_

bool DTTrigPhase2Prod::my_CCBValid_

private

Definition at line 109 of file DTTrigPhase2Prod.cc.

◆ primitives_

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

Definition at line 105 of file DTTrigPhase2Prod.cc.

◆ qmap_

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

private

Definition at line 152 of file DTTrigPhase2Prod.cc.

Referenced by assignQualityOrder().

◆ rpc_integrator_

std::unique_ptr<RPCIntegrator> DTTrigPhase2Prod::rpc_integrator_

private

Definition at line 145 of file DTTrigPhase2Prod.cc.

Referenced by endRun(), and produce().

◆ rpcRecHitsLabel_

edm::EDGetTokenT<RPCRecHitCollection> DTTrigPhase2Prod::rpcRecHitsLabel_

private

Definition at line 125 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ scenario_

int DTTrigPhase2Prod::scenario_

private

Definition at line 119 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ superCellhalfspacewidth_

int DTTrigPhase2Prod::superCellhalfspacewidth_

private

Definition at line 138 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ superCelltimewidth_

float DTTrigPhase2Prod::superCelltimewidth_

private

Definition at line 139 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ useRPC_

bool DTTrigPhase2Prod::useRPC_

private

Definition at line 146 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_

◆ cmssw_for_global_

◆ debug_

◆ do_correlation_

◆ dT0_correlate_TP_

◆ dtDigisToken_

◆ dtGeo_

◆ dtGeomH

◆ dump_

◆ geometry_tag_

◆ globalcoordsobtainer_

◆ grouping_obj_

◆ mpathanalyzer_

◆ mpathassociator_

◆ mpathqualityenhancer_

◆ mpathredundantfilter_

◆ my_BXoffset_

◆ my_CCBValid_

◆ primitives_

◆ qmap_

◆ rpc_integrator_

◆ rpcRecHitsLabel_

◆ scenario_

◆ superCellhalfspacewidth_

◆ superCelltimewidth_

◆ useRPC_