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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Public Attributes
const DTGeometry *	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_

int	df_extended_

double	dT0_correlate_TP_

edm::EDGetTokenT< DTDigiCollection >	dtDigisToken_

bool	dump_

std::shared_ptr< GlobalCoordsObtainer >	globalcoordsobtainer_

std::unique_ptr< MotherGrouping >	grouping_obj_

int	max_index_

std::unique_ptr< MuonPathAnalyzer >	mpathanalyzer_

std::unique_ptr< MuonPathAssociator >	mpathassociator_

std::unique_ptr< MPFilter >	mpathhitsfilter_

std::unique_ptr< MPFilter >	mpathqualityenhancer_

std::unique_ptr< MPFilter >	mpathqualityenhancerbayes_

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<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 72 of file DTTrigPhase2Prod.cc.

Member Typedef Documentation

◆ DTDigiMap

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

private

Definition at line 73 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_const_iterator

typedef DTDigiMap::const_iterator DTTrigPhase2Prod::DTDigiMap_const_iterator

private

Definition at line 75 of file DTTrigPhase2Prod.cc.

◆ DTDigiMap_iterator

typedef DTDigiMap::iterator DTTrigPhase2Prod::DTDigiMap_iterator

private

Definition at line 74 of file DTTrigPhase2Prod.cc.

Constructor & Destructor Documentation

◆ DTTrigPhase2Prod()

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

Constructor.

Definition at line 172 of file DTTrigPhase2Prod.cc.

     : qmap_({{8, 8}, {7, 7}, {6, 6}, {4, 4}, {3, 3}, {2, 2}, {1, 1}}) {
   produces<L1Phase2MuDTPhContainer>();
   produces<L1Phase2MuDTThContainer>();
   produces<L1Phase2MuDTExtPhContainer>();
   produces<L1Phase2MuDTExtThContainer>();
 
   debug_ = pset.getUntrackedParameter<bool>("debug");
   dump_ = pset.getUntrackedParameter<bool>("dump");
 
   scenario_ = pset.getParameter<int>("scenario");
 
   df_extended_ = pset.getParameter<int>("df_extended");
   max_index_ = pset.getParameter<int>("max_primitives") - 1;
 
   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());
   globalcoordsobtainer_ = std::make_shared<GlobalCoordsObtainer>(pset);
   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, globalcoordsobtainer_);
   }
 
   // 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);
   mpathqualityenhancerbayes_ = std::make_unique<MPQualityEnhancerFilterBayes>(pset);
   mpathredundantfilter_ = std::make_unique<MPRedundantFilter>(pset);
   mpathhitsfilter_ = std::make_unique<MPCleanHitsFilter>(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 234 of file DTTrigPhase2Prod.cc.

References debug_, and LogDebug.

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

Member Function Documentation

◆ assignIndex()

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

private

Definition at line 775 of file DTTrigPhase2Prod.cc.

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

Referenced by produce().

                                                                      {
   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)
       if (primitive.index <= max_index_)
         inMPaths.push_back(primitive);
   }
 }

◆ assignIndexPerBX()

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

private

Definition at line 790 of file DTTrigPhase2Prod.cc.

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

Referenced by assignIndex().

                                                                           {
   // 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
 }

◆ assignQualityOrder()

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

private

Definition at line 823 of file DTTrigPhase2Prod.cc.

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

Referenced by assignIndexPerBX().

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

◆ beginRun()

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

override

Create Trigger Units before starting event processing.

Definition at line 239 of file DTTrigPhase2Prod.cc.

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

                                                                                     {
   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
   mpathqualityenhancerbayes_->initialise(iEventSetup);  // Filter object initialisation
   mpathhitsfilter_->initialise(iEventSetup);
   mpathassociator_->initialise(iEventSetup);  // Associator object initialisation
 
   if (auto geom = iEventSetup.getHandle(dtGeomH)) {
     dtGeo_ = &(*geom);
   }
 }

◆ distribDigis()

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

private

Definition at line 830 of file DTTrigPhase2Prod.cc.

References processDigi().

Referenced by produce().

                                                                                                      {
   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;
 }

◆ endRun()

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

override

endRun: finish things

Definition at line 717 of file DTTrigPhase2Prod.cc.

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

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

◆ fillDescriptions()

void DTTrigPhase2Prod::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 879 of file DTTrigPhase2Prod.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), submitPVResolutionJobs::desc, and HLT_2022v15_cff::InputTag.

                                                                                   {
   // dtTriggerPhase2PrimitiveDigis
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("digiTag", edm::InputTag("CalibratedDigis"));
   desc.add<int>("trigger_with_sl", 4);
   desc.add<int>("timeTolerance", 999999);
   desc.add<double>("tanPhiTh", 1.0);
   desc.add<double>("tanPhiThw2max", 1.3);
   desc.add<double>("tanPhiThw2min", 0.5);
   desc.add<double>("tanPhiThw1max", 0.9);
   desc.add<double>("tanPhiThw1min", 0.2);
   desc.add<double>("tanPhiThw0", 0.5);
   desc.add<double>("chi2Th", 0.01);
   desc.add<double>("chi2corTh", 0.1);
   desc.add<bool>("useBX_correlation", false);
   desc.add<double>("dT0_correlate_TP", 25.0);
   desc.add<int>("dBX_correlate_TP", 0);
   desc.add<double>("dTanPsi_correlate_TP", 99999.0);
   desc.add<bool>("clean_chi2_correlation", true);
   desc.add<bool>("allow_confirmation", true);
   desc.add<double>("minx_match_2digis", 1.0);
   desc.add<int>("scenario", 0);
   desc.add<int>("df_extended", 0);
   desc.add<int>("max_primitives", 999);
   desc.add<edm::FileInPath>("ttrig_filename", edm::FileInPath("L1Trigger/DTTriggerPhase2/data/wire_rawId_ttrig.txt"));
   desc.add<edm::FileInPath>("z_filename", edm::FileInPath("L1Trigger/DTTriggerPhase2/data/wire_rawId_z.txt"));
   desc.add<edm::FileInPath>("shift_filename", edm::FileInPath("L1Trigger/DTTriggerPhase2/data/wire_rawId_x.txt"));
   desc.add<edm::FileInPath>("shift_theta_filename", edm::FileInPath("L1Trigger/DTTriggerPhase2/data/theta_shift.txt"));
   desc.add<edm::FileInPath>("global_coords_filename",
                             edm::FileInPath("L1Trigger/DTTriggerPhase2/data/global_coord_perp_x_phi0.txt"));
   desc.add<int>("algo", 0);
   desc.add<int>("minHits4Fit", 3);
   desc.add<bool>("splitPathPerSL", true);
   desc.addUntracked<bool>("debug", false);
   desc.addUntracked<bool>("dump", false);
   desc.add<edm::InputTag>("rpcRecHits", edm::InputTag("rpcRecHits"));
   desc.add<bool>("useRPC", false);
   desc.add<int>("bx_window", 1);
   desc.add<double>("phi_window", 50.0);
   desc.add<int>("max_quality_to_overwrite_t0", 9);
   desc.add<bool>("storeAllRPCHits", false);
   desc.add<bool>("activateBuffer", false);
   desc.add<double>("superCelltimewidth", 400);
   desc.add<int>("superCellspacewidth", 20);
   {
     edm::ParameterSetDescription psd0;
     psd0.addUntracked<bool>("debug", false);
     psd0.add<double>("angletan", 0.3);
     psd0.add<double>("anglebinwidth", 1.0);
     psd0.add<double>("posbinwidth", 2.1);
     psd0.add<double>("maxdeltaAngDeg", 10);
     psd0.add<double>("maxdeltaPos", 10);
     psd0.add<int>("UpperNumber", 6);
     psd0.add<int>("LowerNumber", 4);
     psd0.add<double>("MaxDistanceToWire", 0.03);
     psd0.add<int>("minNLayerHits", 6);
     psd0.add<int>("minSingleSLHitsMax", 3);
     psd0.add<int>("minSingleSLHitsMin", 3);
     psd0.add<bool>("allowUncorrelatedPatterns", true);
     psd0.add<int>("minUncorrelatedHits", 3);
     desc.add<edm::ParameterSetDescription>("HoughGrouping", psd0);
   }
   {
     edm::ParameterSetDescription psd0;
     psd0.add<edm::FileInPath>(
         "pattern_filename", edm::FileInPath("L1Trigger/DTTriggerPhase2/data/PseudoBayesPatterns_uncorrelated_v0.root"));
     psd0.addUntracked<bool>("debug", false);
     psd0.add<int>("minNLayerHits", 3);
     psd0.add<int>("minSingleSLHitsMax", 3);
     psd0.add<int>("minSingleSLHitsMin", 0);
     psd0.add<int>("allowedVariance", 1);
     psd0.add<bool>("allowDuplicates", false);
     psd0.add<bool>("setLateralities", true);
     psd0.add<bool>("allowUncorrelatedPatterns", true);
     psd0.add<int>("minUncorrelatedHits", 3);
     psd0.add<bool>("saveOnPlace", true);
     psd0.add<int>("maxPathsPerMatch", 256);
     desc.add<edm::ParameterSetDescription>("PseudoBayesPattern", psd0);
   }
   descriptions.add("dtTriggerPhase2PrimitiveDigis", desc);
 }

◆ getMinimumQuality()

MP_QUALITY DTTrigPhase2Prod::getMinimumQuality ( void )

◆ hasPosRF()

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

Definition at line 738 of file DTTrigPhase2Prod.cc.

References fileinputsource_cfi::sec.

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

fileinputsource_cfi.sec

sec

Definition: fileinputsource_cfi.py:88

◆ inner()

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

Definition at line 733 of file DTTrigPhase2Prod.cc.

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

Referenced by produce().

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

◆ outer()

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

Definition at line 728 of file DTTrigPhase2Prod.cc.

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

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

◆ printmP()

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

Definition at line 740 of file DTTrigPhase2Prod.cc.

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().

                                                                               {
   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);
 }

◆ printmPC()

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

Definition at line 750 of file DTTrigPhase2Prod.cc.

Referenced by produce().

                                                                                {
   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);
 }

◆ processDigi()

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

private

Definition at line 849 of file DTTrigPhase2Prod.cc.

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

Referenced by distribDigis().

                                                                                           {
   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);
 }

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

References activateBuffer_, algo_, assignIndex(), DTSuperLayerId::chamberId(), DTGeometry::chambers(), cmsdt::CHI2RES_CONV, cmsdt::CHIGHQ, cscdqm::DATA, debug_, df_extended_, 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, SiStripPI::max, cmsdt::MC, eostools::move(), mpathanalyzer_, mpathassociator_, mpathhitsfilter_, 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.

                                                                            {
   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);
   } else {
     mpathhitsfilter_->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->phi_cmssw(),
                                             muonpath->phiB_cmssw(),
                                             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<L1Phase2MuDTExtPhDigi> outExtP2Ph;
   vector<L1Phase2MuDTThDigi> outP2Th;
   vector<L1Phase2MuDTExtThDigi> outExtP2Th;
 
   // 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) {
       if (df_extended_ == 1 || df_extended_ == 2) {
         int pathWireId[8] = {metaPrimitiveIt.wi1,
                              metaPrimitiveIt.wi2,
                              metaPrimitiveIt.wi3,
                              metaPrimitiveIt.wi4,
                              metaPrimitiveIt.wi5,
                              metaPrimitiveIt.wi6,
                              metaPrimitiveIt.wi7,
                              metaPrimitiveIt.wi8};
 
         int pathTDC[8] = {max((int)round(metaPrimitiveIt.tdc1 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc2 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc3 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc4 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc5 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc6 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc7 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc8 - shift_back * LHC_CLK_FREQ), -1)};
 
         int pathLat[8] = {metaPrimitiveIt.lat1,
                           metaPrimitiveIt.lat2,
                           metaPrimitiveIt.lat3,
                           metaPrimitiveIt.lat4,
                           metaPrimitiveIt.lat5,
                           metaPrimitiveIt.lat6,
                           metaPrimitiveIt.lat7,
                           metaPrimitiveIt.lat8};
 
         // phiTP (extended DF)
         outExtP2Ph.emplace_back(
             L1Phase2MuDTExtPhDigi((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  (m_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)
                                   (int)round(metaPrimitiveIt.x * 1000),                        // ux    (m_xLocal)
                                   (int)round(metaPrimitiveIt.tanPhi * 1000),                   // utan  (m_tanPsi)
                                   (int)round(metaPrimitiveIt.phi_cmssw * PHIRES_CONV),    // uphi  (m_phiAngleCMSSW)
                                   (int)round(metaPrimitiveIt.phiB_cmssw * PHIBRES_CONV),  // uphib (m_phiBendingCMSSW)
                                   metaPrimitiveIt.rpcFlag,                                // urpc  (m_rpcFlag)
                                   pathWireId,
                                   pathTDC,
                                   pathLat));
       }
       if (df_extended_ == 0 || df_extended_ == 2) {
         // phiTP (standard DF)
         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 {
       if (df_extended_ == 1 || df_extended_ == 2) {
         int pathWireId[4] = {metaPrimitiveIt.wi1, metaPrimitiveIt.wi2, metaPrimitiveIt.wi3, metaPrimitiveIt.wi4};
 
         int pathTDC[4] = {max((int)round(metaPrimitiveIt.tdc1 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc2 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc3 - shift_back * LHC_CLK_FREQ), -1),
                           max((int)round(metaPrimitiveIt.tdc4 - shift_back * LHC_CLK_FREQ), -1)};
 
         int pathLat[4] = {metaPrimitiveIt.lat1, metaPrimitiveIt.lat2, metaPrimitiveIt.lat3, metaPrimitiveIt.lat4};
 
         // thTP (extended DF)
         outExtP2Th.emplace_back(
             L1Phase2MuDTExtThDigi((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)
                                   (int)round(metaPrimitiveIt.x * 1000),                        // ux    (m_yLocal)
                                   (int)round(metaPrimitiveIt.phi_cmssw * PHIRES_CONV),         // uphi  (m_zCMSSW)
                                   (int)round(metaPrimitiveIt.phiB_cmssw * PHIBRES_CONV),       // uphib (m_kCMSSW)
                                   metaPrimitiveIt.rpcFlag,                                     // urpc  (m_rpcFlag)
                                   pathWireId,
                                   pathTDC,
                                   pathLat));
       }
       if (df_extended_ == 0 || df_extended_ == 2) {
         // thTP (standard DF)
         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);
     }
   }
 
   // Storing Phi results
   if (df_extended_ == 1 || df_extended_ == 2) {
     std::unique_ptr<L1Phase2MuDTExtPhContainer> resultExtP2Ph(new L1Phase2MuDTExtPhContainer);
     resultExtP2Ph->setContainer(outExtP2Ph);
     iEvent.put(std::move(resultExtP2Ph));
   }
   if (df_extended_ == 0 || df_extended_ == 2) {
     std::unique_ptr<L1Phase2MuDTPhContainer> resultP2Ph(new L1Phase2MuDTPhContainer);
     resultP2Ph->setContainer(outP2Ph);
     iEvent.put(std::move(resultP2Ph));
   }
   outExtP2Ph.clear();
   outExtP2Ph.erase(outExtP2Ph.begin(), outExtP2Ph.end());
   outP2Ph.clear();
   outP2Ph.erase(outP2Ph.begin(), outP2Ph.end());
 
   // Storing Theta results
   if (df_extended_ == 1 || df_extended_ == 2) {
     std::unique_ptr<L1Phase2MuDTExtThContainer> resultExtP2Th(new L1Phase2MuDTExtThContainer);
     resultExtP2Th->setContainer(outExtP2Th);
     iEvent.put(std::move(resultExtP2Th));
   }
   if (df_extended_ == 0 || df_extended_ == 2) {
     std::unique_ptr<L1Phase2MuDTThContainer> resultP2Th(new L1Phase2MuDTThContainer);
     resultP2Th->setContainer(outP2Th);
     iEvent.put(std::move(resultP2Th));
   }
   outExtP2Th.clear();
   outExtP2Th.erase(outExtP2Th.begin(), outExtP2Th.end());
   outP2Th.clear();
   outP2Th.erase(outP2Th.begin(), outP2Th.end());
 }

◆ rango()

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

Definition at line 767 of file DTTrigPhase2Prod.cc.

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

Referenced by printmP(), and printmPC().

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

◆ setChiSquareThreshold()

void DTTrigPhase2Prod::setChiSquareThreshold ( float ch2Thr )

◆ setMinimumQuality()

void DTTrigPhase2Prod::setMinimumQuality ( MP_QUALITY q )

Member Data Documentation

◆ activateBuffer_

bool DTTrigPhase2Prod::activateBuffer_

private

Definition at line 146 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ algo_

int DTTrigPhase2Prod::algo_

private

Definition at line 135 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ debug_

bool DTTrigPhase2Prod::debug_

private

Definition at line 123 of file DTTrigPhase2Prod.cc.

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

◆ df_extended_

int DTTrigPhase2Prod::df_extended_

private

Definition at line 127 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ dT0_correlate_TP_

double DTTrigPhase2Prod::dT0_correlate_TP_

private

Definition at line 125 of file DTTrigPhase2Prod.cc.

◆ dtDigisToken_

edm::EDGetTokenT<DTDigiCollection> DTTrigPhase2Prod::dtDigisToken_

private

Definition at line 131 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ dtGeo_

const DTGeometry* DTTrigPhase2Prod::dtGeo_

Definition at line 111 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), and produce().

◆ dtGeomH

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

Definition at line 112 of file DTTrigPhase2Prod.cc.

Referenced by beginRun().

◆ dump_

bool DTTrigPhase2Prod::dump_

private

Definition at line 124 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ globalcoordsobtainer_

std::shared_ptr<GlobalCoordsObtainer> DTTrigPhase2Prod::globalcoordsobtainer_

private

Definition at line 143 of file DTTrigPhase2Prod.cc.

◆ grouping_obj_

std::unique_ptr<MotherGrouping> DTTrigPhase2Prod::grouping_obj_

private

Definition at line 136 of file DTTrigPhase2Prod.cc.

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

◆ max_index_

int DTTrigPhase2Prod::max_index_

private

Definition at line 128 of file DTTrigPhase2Prod.cc.

Referenced by assignIndex().

◆ mpathanalyzer_

std::unique_ptr<MuonPathAnalyzer> DTTrigPhase2Prod::mpathanalyzer_

private

Definition at line 137 of file DTTrigPhase2Prod.cc.

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

◆ mpathassociator_

std::unique_ptr<MuonPathAssociator> DTTrigPhase2Prod::mpathassociator_

private

Definition at line 142 of file DTTrigPhase2Prod.cc.

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

◆ mpathhitsfilter_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathhitsfilter_

private

Definition at line 141 of file DTTrigPhase2Prod.cc.

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

◆ mpathqualityenhancer_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathqualityenhancer_

private

Definition at line 138 of file DTTrigPhase2Prod.cc.

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

◆ mpathqualityenhancerbayes_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathqualityenhancerbayes_

private

Definition at line 139 of file DTTrigPhase2Prod.cc.

Referenced by beginRun(), and endRun().

◆ mpathredundantfilter_

std::unique_ptr<MPFilter> DTTrigPhase2Prod::mpathredundantfilter_

private

Definition at line 140 of file DTTrigPhase2Prod.cc.

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

◆ my_BXoffset_

int DTTrigPhase2Prod::my_BXoffset_

private

Definition at line 120 of file DTTrigPhase2Prod.cc.

◆ my_CCBValid_

bool DTTrigPhase2Prod::my_CCBValid_

private

Definition at line 117 of file DTTrigPhase2Prod.cc.

◆ primitives_

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

Definition at line 113 of file DTTrigPhase2Prod.cc.

◆ qmap_

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

private

Definition at line 161 of file DTTrigPhase2Prod.cc.

Referenced by assignQualityOrder().

◆ rpc_integrator_

std::unique_ptr<RPCIntegrator> DTTrigPhase2Prod::rpc_integrator_

private

Definition at line 154 of file DTTrigPhase2Prod.cc.

Referenced by endRun(), and produce().

◆ rpcRecHitsLabel_

edm::EDGetTokenT<RPCRecHitCollection> DTTrigPhase2Prod::rpcRecHitsLabel_

private

Definition at line 132 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ scenario_

int DTTrigPhase2Prod::scenario_

private

Definition at line 126 of file DTTrigPhase2Prod.cc.

Referenced by produce().

◆ superCellhalfspacewidth_

int DTTrigPhase2Prod::superCellhalfspacewidth_

private

Definition at line 147 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ superCelltimewidth_

float DTTrigPhase2Prod::superCelltimewidth_

private

Definition at line 148 of file DTTrigPhase2Prod.cc.

Referenced by processDigi().

◆ useRPC_

bool DTTrigPhase2Prod::useRPC_

private

Definition at line 155 of file DTTrigPhase2Prod.cc.

Referenced by produce().

Public Member Functions

Static 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()

◆ fillDescriptions()

◆ getMinimumQuality()

◆ hasPosRF()

◆ inner()

◆ outer()

◆ printmP()

◆ printmPC()

◆ processDigi()

◆ produce()

◆ rango()

◆ setChiSquareThreshold()

◆ setMinimumQuality()

Member Data Documentation

◆ activateBuffer_

◆ algo_

◆ debug_

◆ df_extended_

◆ dT0_correlate_TP_

◆ dtDigisToken_

◆ dtGeo_

◆ dtGeomH

◆ dump_

◆ globalcoordsobtainer_

◆ grouping_obj_

◆ max_index_

◆ mpathanalyzer_

◆ mpathassociator_

◆ mpathhitsfilter_

◆ mpathqualityenhancer_

◆ mpathqualityenhancerbayes_

◆ mpathredundantfilter_

◆ my_BXoffset_

◆ my_CCBValid_

◆ primitives_

◆ qmap_

◆ rpc_integrator_

◆ rpcRecHitsLabel_

◆ scenario_

◆ superCellhalfspacewidth_

◆ superCelltimewidth_

◆ useRPC_