CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Static Private Member Functions | Private Attributes
L1TMuonProducer Class Reference

#include <L1Trigger/L1TMuon/src/L1TMuonProducer.cc>

Inheritance diagram for L1TMuonProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 L1TMuonProducer (const edm::ParameterSet &)
 
 ~L1TMuonProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 
registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void addMuonsToCollections (MicroGMTConfiguration::InterMuonList &coll, MicroGMTConfiguration::InterMuonList &interout, std::auto_ptr< MuonBxCollection > &out, int bx) const
 
virtual void beginJob ()
 
virtual void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void beginRun (edm::Run const &, edm::EventSetup const &)
 
void calculateRank (MicroGMTConfiguration::InterMuonList &muons) const
 
void convertMuons (edm::Handle< MicroGMTConfiguration::InputCollection > const &in, MicroGMTConfiguration::InterMuonList &out, GMTInternalWedges &wedges, int bx) const
 
virtual void endJob ()
 
virtual void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void endRun (edm::Run const &, edm::EventSetup const &)
 
virtual void produce (edm::Event &, const edm::EventSetup &)
 
void sortMuons (MicroGMTConfiguration::InterMuonList &, unsigned) const
 
void splitAndConvertMuons (edm::Handle< MicroGMTConfiguration::InputCollection > const &in, MicroGMTConfiguration::InterMuonList &out_pos, MicroGMTConfiguration::InterMuonList &out_neg, GMTInternalWedges &wedges_pos, GMTInternalWedges &wedges_neg, int bx) const
 

Static Private Member Functions

static bool compareMuons (const std::shared_ptr< MicroGMTConfiguration::InterMuon > &mu1, const std::shared_ptr< MicroGMTConfiguration::InterMuon > &mu2)
 

Private Attributes

edm::InputTag m_barrelTfInputTag
 
edm::EDGetTokenT
< MicroGMTConfiguration::InputCollection
m_barrelTfInputToken
 
int m_bxMax
 
int m_bxMin
 
edm::EDGetTokenT
< MicroGMTConfiguration::CaloInputCollection
m_caloTowerInputToken
 
MicroGMTCancelOutUnit m_cancelOutUnit
 
std::ofstream m_debugOut
 
edm::InputTag m_endcapTfInputTag
 
edm::EDGetTokenT
< MicroGMTConfiguration::InputCollection
m_endcapTfInputToken
 
MicroGMTIsolationUnit m_isolationUnit
 
edm::InputTag m_overlapTfInputTag
 
edm::EDGetTokenT
< MicroGMTConfiguration::InputCollection
m_overlapTfInputToken
 
std::shared_ptr
< MicroGMTRankPtQualLUT
m_rankPtQualityLUT
 
edm::InputTag m_trigTowerTag
 
std::unique_ptr
< L1TMuonGlobalParams
microGMTParams
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 
TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Description: Takes txt-file input and produces barrel- / overlap- / forward TF muons

Implementation: [Notes on implementation]

Definition at line 57 of file L1TMuonProducer.cc.

Constructor & Destructor Documentation

L1TMuonProducer::L1TMuonProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 130 of file L1TMuonProducer.cc.

References edm::ParameterSet::getParameter(), m_barrelTfInputTag, m_barrelTfInputToken, m_caloTowerInputToken, m_endcapTfInputTag, m_endcapTfInputToken, m_overlapTfInputTag, m_overlapTfInputToken, and m_trigTowerTag.

130  : m_debugOut("test/debug/iso_debug.dat")
131 {
132  // edm::InputTag barrelTfInputTag = iConfig.getParameter<edm::InputTag>("barrelTFInput");
133  // edm::InputTag overlapTfInputTag = iConfig.getParameter<edm::InputTag>("overlapTFInput");
134  // edm::InputTag forwardTfInputTag = iConfig.getParameter<edm::InputTag>("forwardTFInput");
135 
136  m_barrelTfInputTag = iConfig.getParameter<edm::InputTag>("barrelTFInput");
137  m_overlapTfInputTag = iConfig.getParameter<edm::InputTag>("overlapTFInput");
138  m_endcapTfInputTag = iConfig.getParameter<edm::InputTag>("forwardTFInput");
139  m_trigTowerTag = iConfig.getParameter<edm::InputTag>("triggerTowerInput");
140 
141  m_barrelTfInputToken = consumes<MicroGMTConfiguration::InputCollection>(m_barrelTfInputTag);
142  m_overlapTfInputToken = consumes<MicroGMTConfiguration::InputCollection>(m_overlapTfInputTag);
143  m_endcapTfInputToken = consumes<MicroGMTConfiguration::InputCollection>(m_endcapTfInputTag);
144  m_caloTowerInputToken = consumes<MicroGMTConfiguration::CaloInputCollection>(m_trigTowerTag);
145 
146  //register your products
147  produces<MuonBxCollection>();
148  produces<MuonBxCollection>("imdMuonsBMTF");
149  produces<MuonBxCollection>("imdMuonsEMTFPos");
150  produces<MuonBxCollection>("imdMuonsEMTFNeg");
151  produces<MuonBxCollection>("imdMuonsOMTFPos");
152  produces<MuonBxCollection>("imdMuonsOMTFNeg");
153 }
T getParameter(std::string const &) const
std::ofstream m_debugOut
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_endcapTfInputToken
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_overlapTfInputToken
edm::EDGetTokenT< MicroGMTConfiguration::CaloInputCollection > m_caloTowerInputToken
edm::InputTag m_overlapTfInputTag
edm::InputTag m_barrelTfInputTag
edm::InputTag m_endcapTfInputTag
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_barrelTfInputToken
edm::InputTag m_trigTowerTag
L1TMuonProducer::~L1TMuonProducer ( )

Definition at line 155 of file L1TMuonProducer.cc.

References m_debugOut.

156 {
157  m_debugOut.close();
158 }
std::ofstream m_debugOut

Member Function Documentation

void L1TMuonProducer::addMuonsToCollections ( MicroGMTConfiguration::InterMuonList coll,
MicroGMTConfiguration::InterMuonList interout,
std::auto_ptr< MuonBxCollection > &  out,
int  bx 
) const
private

Definition at line 355 of file L1TMuonProducer.cc.

References RPCpg::mu.

Referenced by produce().

358 {
359  for (auto& mu : coll) {
360  interout.push_back(mu);
361  math::PtEtaPhiMLorentzVector vec{(mu->hwPt()-1)*0.5, mu->hwEta()*0.010875, mu->hwGlobalPhi()*0.010908, 0.0};
362  // FIXME: once we debugged the change global -> local: Change hwLocalPhi -> hwGlobalPhi to test offsets
363  // FIXME: set tfMuonIndex to actual index
364  Muon outMu{vec, mu->hwPt(), mu->hwEta(), mu->hwGlobalPhi(), mu->hwQual(), mu->hwSign(), mu->hwSignValid(), -1, -1, 0, true, -1, mu->hwDPhi(), mu->hwDEta(), mu->hwRank()};
365  if (mu->hwSignValid()) {
366  outMu.setCharge(1 - 2 * mu->hwSign());
367  } else {
368  outMu.setCharge(0);
369  }
370 
371  out->push_back(bx, outMu);
372  }
373 }
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
const int mu
Definition: Constants.h:22
Definition: Muon.h:16
JetCorrectorParametersCollection coll
Definition: classes.h:10
void L1TMuonProducer::beginJob ( void  )
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 432 of file L1TMuonProducer.cc.

433 {
434 }
void L1TMuonProducer::beginLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 470 of file L1TMuonProducer.cc.

471 {
472 }
void L1TMuonProducer::beginRun ( edm::Run const &  run,
edm::EventSetup const &  iSetup 
)
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 443 of file L1TMuonProducer.cc.

References l1t::MicroGMTRankPtQualLUTFactory::create(), edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), l1t::MicroGMTIsolationUnit::initialise(), l1t::MicroGMTCancelOutUnit::initialise(), m_bxMax, m_bxMin, m_cancelOutUnit, m_isolationUnit, m_rankPtQualityLUT, and microGMTParams.

444 {
445  const L1TMuonGlobalParamsRcd& microGMTParamsRcd = iSetup.get<L1TMuonGlobalParamsRcd>();
446  edm::ESHandle<L1TMuonGlobalParams> microGMTParamsHandle;
447  microGMTParamsRcd.get(microGMTParamsHandle);
448 
449  microGMTParams = std::unique_ptr<L1TMuonGlobalParams>(new L1TMuonGlobalParams(*microGMTParamsHandle.product()));
450  if (!microGMTParams) {
451  edm::LogError("L1TMuonProducer") << "Could not retrieve parameters from Event Setup" << std::endl;
452  }
453 
454  //microGMTParams->print(std::cout);
455  m_bxMin = microGMTParams->bxMin();
456  m_bxMax = microGMTParams->bxMax();
457  m_rankPtQualityLUT = l1t::MicroGMTRankPtQualLUTFactory::create(microGMTParams->sortRankLUTPath(), microGMTParams->fwVersion(), microGMTParams->sortRankLUTPtFactor(), microGMTParams->sortRankLUTQualFactor());
460 }
void initialise(L1TMuonGlobalParams *)
Initialisation from ES record.
MicroGMTIsolationUnit m_isolationUnit
std::shared_ptr< MicroGMTRankPtQualLUT > m_rankPtQualityLUT
static ReturnType create(const std::string &filename, const int fwVersion, const unsigned ptFactor, const unsigned qualFactor)
void get(HolderT &iHolder) const
void initialise(L1TMuonGlobalParams *)
Initialisation from ES record.
std::unique_ptr< L1TMuonGlobalParams > microGMTParams
MicroGMTCancelOutUnit m_cancelOutUnit
void L1TMuonProducer::calculateRank ( MicroGMTConfiguration::InterMuonList muons) const
private

Definition at line 345 of file L1TMuonProducer.cc.

References m_rankPtQualityLUT.

Referenced by produce().

346 {
347  for (auto& mu1 : muons) {
348  int rank = m_rankPtQualityLUT->lookup(mu1->hwPt(), mu1->hwQual());
349  mu1->setHwRank(rank);
350  }
351 }
std::shared_ptr< MicroGMTRankPtQualLUT > m_rankPtQualityLUT
tuple muons
Definition: patZpeak.py:38
bool L1TMuonProducer::compareMuons ( const std::shared_ptr< MicroGMTConfiguration::InterMuon > &  mu1,
const std::shared_ptr< MicroGMTConfiguration::InterMuon > &  mu2 
)
staticprivate

Definition at line 293 of file L1TMuonProducer.cc.

Referenced by sortMuons().

294  {
295  return (mu1->hwWins() >= mu2->hwWins());
296 }
void L1TMuonProducer::convertMuons ( edm::Handle< MicroGMTConfiguration::InputCollection > const &  in,
MicroGMTConfiguration::InterMuonList out,
GMTInternalWedges wedges,
int  bx 
) const
private

Definition at line 409 of file L1TMuonProducer.cc.

References i, and findQualityFiles::size.

Referenced by produce().

412 {
413  // initialize the wedge collection:
414  for (int i = 0; i < 12; ++i) {
415  wedges[i] = std::vector<std::shared_ptr<GMTInternalMuon>>();
416  wedges[i].reserve(3);
417  }
418  if (bx < in->getFirstBX() || bx > in->getLastBX()) return;
419  for (size_t i = 0; i < in->size(bx); ++i) {
420  int gPhi = MicroGMTConfiguration::calcGlobalPhi(in->at(bx, i).hwPhi(), in->at(bx, i).trackFinderType(), in->at(bx, i).processor());
421  std::shared_ptr<GMTInternalMuon> outMu = std::make_shared<GMTInternalMuon>(in->at(bx, i), gPhi);
422  out.emplace_back(outMu);
423  wedges[in->at(bx, i).processor()].push_back(outMu);
424  }
425  for (int i = 0; i < 12; ++i) {
426  if(wedges[i].size() > 3) edm::LogWarning("Input Mismatch") << " too many inputs per processor for barrel. Wedge " << i << ": Size " << wedges[i].size() << std::endl;
427  }
428 }
int i
Definition: DBlmapReader.cc:9
tuple size
Write out results.
void L1TMuonProducer::endJob ( void  )
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 438 of file L1TMuonProducer.cc.

438  {
439 }
void L1TMuonProducer::endLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 476 of file L1TMuonProducer.cc.

477 {
478 }
void L1TMuonProducer::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 464 of file L1TMuonProducer.cc.

465 {
466 }
void L1TMuonProducer::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 482 of file L1TMuonProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

482  {
483  //The following says we do not know what parameters are allowed so do no validation
484  // Please change this to state exactly what you do use, even if it is no parameters
486  desc.setUnknown();
487  descriptions.addDefault(desc);
488 }
void addDefault(ParameterSetDescription const &psetDescription)
void L1TMuonProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
privatevirtual

Implements edm::EDProducer.

Definition at line 169 of file L1TMuonProducer.cc.

References addMuonsToCollections(), l1t::bmtf, calculateRank(), convertMuons(), l1t::coordinate, l1t::emtf_neg, l1t::emtf_pos, l1t::MicroGMTIsolationUnit::extrapolateMuons(), edm::Event::getByToken(), l1t::MicroGMTIsolationUnit::isolatePreSummed(), m_barrelTfInputToken, m_bxMax, m_bxMin, m_caloTowerInputToken, m_cancelOutUnit, m_debugOut, m_endcapTfInputToken, m_isolationUnit, m_overlapTfInputToken, RPCpg::mu, l1t::omtf_neg, l1t::omtf_pos, edm::Event::put(), l1t::MicroGMTCancelOutUnit::setCancelOutBits(), l1t::MicroGMTCancelOutUnit::setCancelOutBitsOverlapBarrel(), l1t::MicroGMTCancelOutUnit::setCancelOutBitsOverlapEndcap(), l1t::MicroGMTIsolationUnit::setTowerSums(), sortMuons(), splitAndConvertMuons(), and testEve_cfg::tracks.

170 {
171  using namespace edm;
172  std::auto_ptr<MuonBxCollection> outMuons (new MuonBxCollection());
173  std::auto_ptr<MuonBxCollection> imdMuonsBMTF (new MuonBxCollection());
174  std::auto_ptr<MuonBxCollection> imdMuonsEMTFPos (new MuonBxCollection());
175  std::auto_ptr<MuonBxCollection> imdMuonsEMTFNeg (new MuonBxCollection());
176  std::auto_ptr<MuonBxCollection> imdMuonsOMTFPos (new MuonBxCollection());
177  std::auto_ptr<MuonBxCollection> imdMuonsOMTFNeg (new MuonBxCollection());
178 
179 
184 
185  iEvent.getByToken(m_barrelTfInputToken, bmtfMuons);
186  iEvent.getByToken(m_endcapTfInputToken, emtfMuons);
187  iEvent.getByToken(m_overlapTfInputToken, omtfMuons);
188  iEvent.getByToken(m_caloTowerInputToken, trigTowers);
189 
190  // set BX range for outputs
191  outMuons->setBXRange(m_bxMin, m_bxMax);
192  imdMuonsBMTF->setBXRange(m_bxMin, m_bxMax);
193  imdMuonsEMTFPos->setBXRange(m_bxMin, m_bxMax);
194  imdMuonsEMTFNeg->setBXRange(m_bxMin, m_bxMax);
195  imdMuonsOMTFPos->setBXRange(m_bxMin, m_bxMax);
196  imdMuonsOMTFNeg->setBXRange(m_bxMin, m_bxMax);
197 
198  for (int bx = m_bxMin; bx <= m_bxMax; ++bx) {
199  m_isolationUnit.setTowerSums(*trigTowers, bx);
200  MicroGMTConfiguration::InterMuonList internMuonsBmtf;
201  MicroGMTConfiguration::InterMuonList internMuonsEmtfPos;
202  MicroGMTConfiguration::InterMuonList internMuonsEmtfNeg;
203  MicroGMTConfiguration::InterMuonList internMuonsOmtfPos;
204  MicroGMTConfiguration::InterMuonList internMuonsOmtfNeg;
205 
206  // These wedges contain shared pointers to the ones in the InterMuonList
207  GMTInternalWedges omtfNegWedges;
208  GMTInternalWedges bmtfWedges;
209  GMTInternalWedges emtfPosWedges;
210  GMTInternalWedges emtfNegWedges;
211  GMTInternalWedges omtfPosWedges;
212 
213  // this converts the InputMuon type to the InternalMuon type and splits them into
214  // positive / negative eta collections necessary as LUTs may differ for pos / neg.
215  convertMuons(bmtfMuons, internMuonsBmtf, bmtfWedges, bx);
216  splitAndConvertMuons(emtfMuons, internMuonsEmtfPos, internMuonsEmtfNeg, emtfPosWedges, emtfNegWedges, bx);
217  splitAndConvertMuons(omtfMuons, internMuonsOmtfPos, internMuonsOmtfNeg, omtfPosWedges, omtfNegWedges, bx);
218 
219  // cancel out within the track finders:
223  // cancel-out for endcap will be done in the sorter
226 
227  // cancel out between track finder acceptance overlaps:
232 
233  m_isolationUnit.extrapolateMuons(internMuonsBmtf);
234  m_isolationUnit.extrapolateMuons(internMuonsEmtfNeg);
235  m_isolationUnit.extrapolateMuons(internMuonsEmtfPos);
236  m_isolationUnit.extrapolateMuons(internMuonsOmtfNeg);
237  m_isolationUnit.extrapolateMuons(internMuonsOmtfPos);
238 
239  // the rank calculated here is used in the sort below
240  calculateRank(internMuonsBmtf);
241  calculateRank(internMuonsEmtfNeg);
242  calculateRank(internMuonsEmtfPos);
243  calculateRank(internMuonsOmtfNeg);
244  calculateRank(internMuonsOmtfPos);
245 
246  // The sort function both sorts and removes all but best "nSurvivors"
247  sortMuons(internMuonsBmtf, 8);
248  sortMuons(internMuonsOmtfPos, 4);
249  sortMuons(internMuonsOmtfNeg, 4);
250  sortMuons(internMuonsEmtfPos, 4);
251  sortMuons(internMuonsEmtfNeg, 4);
252 
253  // This combines the 5 streams into one InternalMuon collection for
254  // the final global sort.
256  addMuonsToCollections(internMuonsEmtfPos, internalMuons, imdMuonsEMTFPos, bx);
257  addMuonsToCollections(internMuonsOmtfPos, internalMuons, imdMuonsOMTFPos, bx);
258  addMuonsToCollections(internMuonsBmtf, internalMuons, imdMuonsBMTF, bx);
259  addMuonsToCollections(internMuonsOmtfNeg, internalMuons, imdMuonsOMTFNeg, bx);
260  addMuonsToCollections(internMuonsEmtfNeg, internalMuons, imdMuonsEMTFNeg, bx);
261 
262  // sort internal muons and delete all but best 8
263  sortMuons(internalMuons, 8);
264 
265  m_isolationUnit.isolatePreSummed(internalMuons);
266  // copy muons to output collection...
267  for (const auto& mu : internalMuons) {
268  if (mu->hwPt() > 0) {
269  math::PtEtaPhiMLorentzVector vec{(mu->hwPt()-1)*0.5, mu->hwEta()*0.010875, mu->hwGlobalPhi()*0.010908, 0.0};
270  int iso = mu->hwAbsIso() + (mu->hwRelIso() << 1);
271  Muon outMu{vec, mu->hwPt(), mu->hwEta(), mu->hwGlobalPhi(), mu->hwQual(), mu->hwSign(), mu->hwSignValid(), iso, -1, 0, true, mu->hwIsoSum(), mu->hwDPhi(), mu->hwDEta(), mu->hwRank()};
272  if (mu->hwSignValid()) {
273  outMu.setCharge(1 - 2 * mu->hwSign());
274  } else {
275  outMu.setCharge(0);
276  }
277  m_debugOut << mu->hwCaloPhi() << " " << mu->hwCaloEta() << std::endl;
278  outMuons->push_back(bx, outMu);
279  }
280  }
281  }
282 
283  iEvent.put(outMuons);
284  iEvent.put(imdMuonsBMTF, "imdMuonsBMTF");
285  iEvent.put(imdMuonsEMTFPos, "imdMuonsEMTFPos");
286  iEvent.put(imdMuonsEMTFNeg, "imdMuonsEMTFNeg");
287  iEvent.put(imdMuonsOMTFPos, "imdMuonsOMTFPos");
288  iEvent.put(imdMuonsOMTFNeg, "imdMuonsOMTFNeg");
289 }
void extrapolateMuons(MicroGMTConfiguration::InterMuonList &) const
void setCancelOutBitsOverlapEndcap(GMTInternalWedges &, GMTInternalWedges &, cancelmode)
Cancel-out between overlap and endcap track finders.
std::ofstream m_debugOut
void setCancelOutBitsOverlapBarrel(GMTInternalWedges &, GMTInternalWedges &, cancelmode)
Cancel-out between overlap and barrel track finders.
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_endcapTfInputToken
void convertMuons(edm::Handle< MicroGMTConfiguration::InputCollection > const &in, MicroGMTConfiguration::InterMuonList &out, GMTInternalWedges &wedges, int bx) const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
void setTowerSums(const MicroGMTConfiguration::CaloInputCollection &inputs, int bx)
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_overlapTfInputToken
edm::EDGetTokenT< MicroGMTConfiguration::CaloInputCollection > m_caloTowerInputToken
MicroGMTIsolationUnit m_isolationUnit
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
Definition: LorentzVector.h:25
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:121
void calculateRank(MicroGMTConfiguration::InterMuonList &muons) const
void sortMuons(MicroGMTConfiguration::InterMuonList &, unsigned) const
const int mu
Definition: Constants.h:22
GMTInternalMuonList InterMuonList
Definition: Muon.h:16
void addMuonsToCollections(MicroGMTConfiguration::InterMuonList &coll, MicroGMTConfiguration::InterMuonList &interout, std::auto_ptr< MuonBxCollection > &out, int bx) const
BXVector< Muon > MuonBxCollection
Definition: Muon.h:10
void setCancelOutBits(GMTInternalWedges &, tftype, cancelmode)
Cancel out between sectors/wedges in one track finder.
tuple tracks
Definition: testEve_cfg.py:39
edm::EDGetTokenT< MicroGMTConfiguration::InputCollection > m_barrelTfInputToken
void splitAndConvertMuons(edm::Handle< MicroGMTConfiguration::InputCollection > const &in, MicroGMTConfiguration::InterMuonList &out_pos, MicroGMTConfiguration::InterMuonList &out_neg, GMTInternalWedges &wedges_pos, GMTInternalWedges &wedges_neg, int bx) const
std::map< int, std::vector< std::shared_ptr< GMTInternalMuon > > > GMTInternalWedges
MicroGMTCancelOutUnit m_cancelOutUnit
void isolatePreSummed(MicroGMTConfiguration::InterMuonList &muons) const
void L1TMuonProducer::sortMuons ( MicroGMTConfiguration::InterMuonList muons,
unsigned  nSurvivors 
) const
private

Definition at line 299 of file L1TMuonProducer.cc.

References compareMuons().

Referenced by produce().

299  {
300  MicroGMTConfiguration::InterMuonList::iterator mu1;
301  // reset from previous sort stage
302  for (mu1 = muons.begin(); mu1 != muons.end(); ++mu1) {
303  (*mu1)->setHwWins(0);
304  }
305 
306  int nCancelled = 0;
307  for (mu1 = muons.begin(); mu1 != muons.end(); ++mu1) {
308  int mu1CancelBit = (*mu1)->hwCancelBit();
309  nCancelled += mu1CancelBit;
310  auto mu2 = mu1;
311  mu2++;
312  for ( ; mu2 != muons.end(); ++mu2) {
313  if (mu1CancelBit != 1 && (*mu2)->hwCancelBit() != 1) {
314  if ((*mu1)->hwRank() >= (*mu2)->hwRank()) {
315  (*mu1)->increaseWins();
316  } else {
317  (*mu2)->increaseWins();
318  }
319  } else if (mu1CancelBit != 1) {
320  (*mu1)->increaseWins();
321  } else if ((*mu2)->hwCancelBit() != 1) {
322  (*mu2)->increaseWins();
323  }
324  }
325  }
326 
327  size_t nMuonsBefore = muons.size();
328  int minWins = nMuonsBefore - nCancelled - nSurvivors;
329 
330  // remove all muons that were cancelled or that do not have sufficient rank
331  // (reduces the container size to nSurvivors)
332  mu1 = muons.begin();
333  while (mu1 != muons.end()) {
334  if ((*mu1)->hwWins() < minWins || (*mu1)->hwCancelBit() == 1) {
335  muons.erase(mu1);
336  }
337  ++mu1;
338  }
340 }
tuple muons
Definition: patZpeak.py:38
static bool compareMuons(const std::shared_ptr< MicroGMTConfiguration::InterMuon > &mu1, const std::shared_ptr< MicroGMTConfiguration::InterMuon > &mu2)
void L1TMuonProducer::splitAndConvertMuons ( edm::Handle< MicroGMTConfiguration::InputCollection > const &  in,
MicroGMTConfiguration::InterMuonList out_pos,
MicroGMTConfiguration::InterMuonList out_neg,
GMTInternalWedges wedges_pos,
GMTInternalWedges wedges_neg,
int  bx 
) const
private

Definition at line 376 of file L1TMuonProducer.cc.

References i, GenerateHcalLaserBadRunList::out, and findQualityFiles::size.

Referenced by produce().

382 {
383  // initialize the wedge collections:
384  for (int i = 0; i < 6; ++i) {
385  wedges_pos[i] = std::vector<std::shared_ptr<GMTInternalMuon>>();
386  wedges_pos[i].reserve(3);
387  wedges_neg[i] = std::vector<std::shared_ptr<GMTInternalMuon>>();
388  wedges_neg[i].reserve(3);
389  }
390  if (bx < in->getFirstBX() || bx > in->getLastBX()) return;
391  for (size_t i = 0; i < in->size(bx); ++i) {
392  int gPhi = MicroGMTConfiguration::calcGlobalPhi(in->at(bx, i).hwPhi(), in->at(bx, i).trackFinderType(), in->at(bx, i).processor());
393  std::shared_ptr<GMTInternalMuon> out = std::make_shared<GMTInternalMuon>(in->at(bx, i), gPhi);
394  if(in->at(bx, i).hwEta() > 0) {
395  out_pos.push_back(out);
396  wedges_pos[in->at(bx, i).processor()].push_back(out);
397  } else {
398  out_neg.emplace_back(out);
399  wedges_neg[in->at(bx, i).processor()].push_back(out);
400  }
401  }
402  for (int i = 0; i < 6; ++i) {
403  if(wedges_pos[i].size() > 3) edm::LogWarning("Input Mismatch") << " too many inputs per processor for emtf+ / omtf+. Wedge " << i << ": Size " << wedges_pos[i].size() << std::endl;
404  if(wedges_neg[i].size() > 3) edm::LogWarning("Input Mismatch") << " too many inputs per processor for emtf- / omtf-. Wedge " << i << ": Size " << wedges_neg[i].size() << std::endl;
405  }
406 }
int i
Definition: DBlmapReader.cc:9
tuple size
Write out results.

Member Data Documentation

edm::InputTag L1TMuonProducer::m_barrelTfInputTag
private

Definition at line 102 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer().

edm::EDGetTokenT<MicroGMTConfiguration::InputCollection> L1TMuonProducer::m_barrelTfInputToken
private

Definition at line 111 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer(), and produce().

int L1TMuonProducer::m_bxMax
private

Definition at line 100 of file L1TMuonProducer.cc.

Referenced by beginRun(), and produce().

int L1TMuonProducer::m_bxMin
private

Definition at line 99 of file L1TMuonProducer.cc.

Referenced by beginRun(), and produce().

edm::EDGetTokenT<MicroGMTConfiguration::CaloInputCollection> L1TMuonProducer::m_caloTowerInputToken
private

Definition at line 114 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer(), and produce().

MicroGMTCancelOutUnit L1TMuonProducer::m_cancelOutUnit
private

Definition at line 108 of file L1TMuonProducer.cc.

Referenced by beginRun(), and produce().

std::ofstream L1TMuonProducer::m_debugOut
private

Definition at line 109 of file L1TMuonProducer.cc.

Referenced by produce(), and ~L1TMuonProducer().

edm::InputTag L1TMuonProducer::m_endcapTfInputTag
private

Definition at line 104 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer().

edm::EDGetTokenT<MicroGMTConfiguration::InputCollection> L1TMuonProducer::m_endcapTfInputToken
private

Definition at line 113 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer(), and produce().

MicroGMTIsolationUnit L1TMuonProducer::m_isolationUnit
private

Definition at line 107 of file L1TMuonProducer.cc.

Referenced by beginRun(), and produce().

edm::InputTag L1TMuonProducer::m_overlapTfInputTag
private

Definition at line 103 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer().

edm::EDGetTokenT<MicroGMTConfiguration::InputCollection> L1TMuonProducer::m_overlapTfInputToken
private

Definition at line 112 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer(), and produce().

std::shared_ptr<MicroGMTRankPtQualLUT> L1TMuonProducer::m_rankPtQualityLUT
private

Definition at line 106 of file L1TMuonProducer.cc.

Referenced by beginRun(), and calculateRank().

edm::InputTag L1TMuonProducer::m_trigTowerTag
private

Definition at line 105 of file L1TMuonProducer.cc.

Referenced by L1TMuonProducer().

std::unique_ptr<L1TMuonGlobalParams> L1TMuonProducer::microGMTParams
private

Definition at line 101 of file L1TMuonProducer.cc.

Referenced by beginRun().