CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Static Private Member Functions | Private Attributes | Static Private Attributes
L1TMicroGMTInputProducerFromGen Class Reference

#include <L1Trigger/L1TGlobalMuon/plugins/L1TMicroGMTInputProducerFromGen.cc>

Inheritance diagram for L1TMicroGMTInputProducerFromGen:
edm::stream::EDProducer<>

Public Member Functions

 L1TMicroGMTInputProducerFromGen (const edm::ParameterSet &)
 
 ~L1TMicroGMTInputProducerFromGen () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Member Functions

void beginLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
 
void beginRun (const edm::Run &, edm::EventSetup const &) override
 
void endLuminosityBlock (const edm::LuminosityBlock &, edm::EventSetup const &) override
 
void endRun (const edm::Run &, edm::EventSetup const &) override
 
void produce (edm::Event &, const edm::EventSetup &) override
 

Static Private Member Functions

static bool compareMuons (const RegionalMuonCand &, const RegionalMuonCand &)
 

Private Attributes

edm::EDGetTokenT< reco::GenParticleCollection > genParticlesToken
 
int m_currEvt
 
TRandom3 m_rnd
 

Static Private Attributes

static const int m_maxMuons = 108
 

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

Description: takes generated muons and fills them in the expected collections for the MicroGMT

Implementation: [Notes on implementation]

Definition at line 51 of file L1TMicroGMTInputProducerFromGen.cc.

Constructor & Destructor Documentation

L1TMicroGMTInputProducerFromGen::L1TMicroGMTInputProducerFromGen ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 87 of file L1TMicroGMTInputProducerFromGen.cc.

References genParticlesToken, and AlCaHLTBitMon_QueryRunRegistry::string.

87  :
88  m_currEvt(0), m_rnd(0)
89 {
90  //register your inputs:
91  genParticlesToken = consumes <reco::GenParticleCollection> (std::string("genParticles"));
92  //register your products
93  produces<RegionalMuonCandBxCollection>("BarrelTFMuons");
94  produces<RegionalMuonCandBxCollection>("OverlapTFMuons");
95  produces<RegionalMuonCandBxCollection>("ForwardTFMuons");
96  produces<MuonCaloSumBxCollection>("TriggerTowerSums");
97 }
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
L1TMicroGMTInputProducerFromGen::genParticlesToken
edm::EDGetTokenT< reco::GenParticleCollection > genParticlesToken
Definition: L1TMicroGMTInputProducerFromGen.cc:69
L1TMicroGMTInputProducerFromGen::~L1TMicroGMTInputProducerFromGen ( )
override

Definition at line 100 of file L1TMicroGMTInputProducerFromGen.cc.

101 {
102  // do anything here that needs to be done at desctruction time
103  // (e.g. close files, deallocate resources etc.)
104 }

Member Function Documentation

void L1TMicroGMTInputProducerFromGen::beginLuminosityBlock ( const edm::LuminosityBlock &  ,
edm::EventSetup const &   
)
overrideprivate

Definition at line 269 of file L1TMicroGMTInputProducerFromGen.cc.

270 {
271 }
void L1TMicroGMTInputProducerFromGen::beginRun ( const edm::Run &  ,
edm::EventSetup const &   
)
overrideprivate

Definition at line 257 of file L1TMicroGMTInputProducerFromGen.cc.

258 {
259 }
bool L1TMicroGMTInputProducerFromGen::compareMuons ( const RegionalMuonCand &  mu1,
const RegionalMuonCand &  mu2 
)
staticprivate

Definition at line 112 of file L1TMicroGMTInputProducerFromGen.cc.

References l1t::RegionalMuonCand::processor().

Referenced by produce().

113 {
114  return mu1.processor() < mu2.processor();
115 }
l1t::RegionalMuonCand::processor
const int processor() const
Get processor ID on which the candidate was found (0..5 for OMTF/EMTF; 0..11 for BMTF) ...
Definition: RegionalMuonCand.h:123
void L1TMicroGMTInputProducerFromGen::endLuminosityBlock ( const edm::LuminosityBlock &  ,
edm::EventSetup const &   
)
overrideprivate

Definition at line 275 of file L1TMicroGMTInputProducerFromGen.cc.

276 {
277 }
void L1TMicroGMTInputProducerFromGen::endRun ( const edm::Run &  ,
edm::EventSetup const &   
)
overrideprivate

Definition at line 263 of file L1TMicroGMTInputProducerFromGen.cc.

264 {
265 }
void L1TMicroGMTInputProducerFromGen::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 281 of file L1TMicroGMTInputProducerFromGen.cc.

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

281  {
282  //The following says we do not know what parameters are allowed so do no validation
283  // Please change this to state exactly what you do use, even if it is no parameters
284  edm::ParameterSetDescription desc;
285  desc.setUnknown();
286  descriptions.addDefault(desc);
287 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
void L1TMicroGMTInputProducerFromGen::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivate

Definition at line 119 of file L1TMicroGMTInputProducerFromGen.cc.

References funct::abs(), l1t::bmtf, reco::Candidate::charge(), compareMuons(), l1t::emtf_neg, l1t::emtf_pos, PVValHelper::eta, reco::Candidate::eta(), GenHFHadronMatcher_cfi::genParticles, genParticlesToken, edm::Event::getByToken(), mps_fire::i, createfilelist::int, LogTrace, m_currEvt, m_maxMuons, m_rnd, MuonErrorMatrixAnalyzer_cfi::maxPt, eostools::move(), RPCpg::mu, l1t::omtf_neg, l1t::omtf_pos, reco::Candidate::pdgId(), phi, reco::Candidate::phi(), Pi, reco::Candidate::pt(), edm::Event::put(), l1t::RegionalMuonCand::setHwEta(), l1t::RegionalMuonCand::setHwPhi(), l1t::RegionalMuonCand::setHwPt(), l1t::RegionalMuonCand::setHwQual(), l1t::RegionalMuonCand::setHwSign(), l1t::RegionalMuonCand::setHwSignValid(), l1t::RegionalMuonCand::setTFIdentifiers(), jetUpdater_cfi::sort, reco::Candidate::status(), and Geom::twoPi().

120 {
121  using namespace edm;
122 
123  std::unique_ptr<RegionalMuonCandBxCollection> barrelMuons (new RegionalMuonCandBxCollection());
124  std::unique_ptr<RegionalMuonCandBxCollection> overlapMuons (new RegionalMuonCandBxCollection());
125  std::unique_ptr<RegionalMuonCandBxCollection> endcapMuons (new RegionalMuonCandBxCollection());
126  std::unique_ptr<MuonCaloSumBxCollection> towerSums (new MuonCaloSumBxCollection());
127 
128  std::vector<RegionalMuonCand> bmMuons;
129  std::vector<RegionalMuonCand> omMuons;
130  std::vector<RegionalMuonCand> emMuons;
131 
132 
133  std::vector<int> muIndices;
134  edm::Handle<reco::GenParticleCollection> genParticles;
135  // Make sure that you can get genParticles
136  if (iEvent.getByToken(genParticlesToken, genParticles)) {
137  int cntr = 0;
138  for (auto it = genParticles->cbegin(); it != genParticles->cend(); ++it) {
139  const reco::Candidate& mcParticle = *it;
140  if( abs(mcParticle.pdgId()) == 13 && mcParticle.status() == 1 ) muIndices.push_back(cntr);
141  cntr++;
142  }
143  }
144  else {
145  LogTrace("GlobalMuon") << " GenParticleCollection not found." << std::endl;
146  }
147 
148  RegionalMuonCand mu;
149  MuonCaloSum tSum;
150  // alternative scale (using full phi bit-width): 163.4521265553765f;
151  const float phiToInt = 91.67324722093171f;
152  // alternative scale: 100.0f;
153  const float etaToInt = 90.9090909090f;
154  const int maxPt = (1 << 9)-1;
155  int muCntr = 0;
156 
157  double twoPi = TMath::Pi() * 2.;
158 
159  for (auto it = muIndices.begin(); it != muIndices.end(); ++it) {
160  // don't really care which muons are taken...
161  // guess there ain't 108 generated anyways
162  if (muCntr == m_maxMuons) break;
163  int gen_idx = *it;
164  const reco::Candidate& mcMuon = genParticles->at(gen_idx);
165  double eta = mcMuon.eta();
166  if (fabs(eta) > 2.45) continue; // out of acceptance
167  int hwPt = int(mcMuon.pt() * 2);
168  hwPt = (hwPt < maxPt ? hwPt : maxPt);
169  int hwEta = int(eta * etaToInt);
170  double phi = mcMuon.phi();
171  if (phi < 0) phi += twoPi; // add 2*pi
172  int hwPhi = (int(phi * phiToInt))%576;
173  int hwQual = 8;
174  int hwCharge = (mcMuon.charge() > 0) ? 0 : 1;
175  int hwChargeValid = 1;
176 
177  mu.setHwPt(hwPt);
178 
179 
180 
181  tftype tf(tftype::bmtf);
182  int globalWedgePhi = (hwPhi+24)%576; // this sets CMS phi = 0 to -15 deg
183  int localPhi = globalWedgePhi%48;
184  int processor = globalWedgePhi / 48 + 1;
185  int globalSectorPhi = (hwPhi-24); // this sets CMS phi = 0 to +15 deg
186  if (globalSectorPhi < 0) {
187  globalSectorPhi += 576;
188  }
189 
190 
191  if (fabs(eta) > 0.8) {
192  if (fabs(eta) < 1.2) {
193  tf = (eta > 0 ? tftype::omtf_pos : tftype::omtf_neg);
194  processor = globalSectorPhi / 96 + 1;
195  localPhi = globalSectorPhi%96;
196  } else {
197  tf = (eta > 0 ? tftype::emtf_pos : tftype::emtf_neg);
198  processor = globalSectorPhi / 96 + 1;
199  localPhi = globalSectorPhi%96;
200  }
201  }
202  mu.setHwPhi(localPhi);
203  mu.setTFIdentifiers(processor, tf);
204 
205  mu.setHwEta(hwEta);
206  mu.setHwSign(hwCharge);
207  mu.setHwSignValid(hwChargeValid);
208  mu.setHwQual(hwQual);
209 
210  if (fabs(eta) < 0.8 && bmMuons.size() < 36) {
211  bmMuons.push_back(mu);
212  muCntr++;
213  } else if (fabs(eta) < 1.2 && omMuons.size() < 36) {
214  omMuons.push_back(mu);
215  muCntr++;
216  } else if (emMuons.size() < 36) {
217  emMuons.push_back(mu);
218  muCntr++;
219  }
220  }
221 
222  std::sort(bmMuons.begin(), bmMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);
223  std::sort(omMuons.begin(), omMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);
224  std::sort(emMuons.begin(), emMuons.end(), L1TMicroGMTInputProducerFromGen::compareMuons);
225 
226  for (const auto& mu:bmMuons) {
227  barrelMuons->push_back(0, mu);
228  }
229 
230  for (const auto& mu:omMuons) {
231  overlapMuons->push_back(0, mu);
232  }
233 
234  for (const auto& mu:emMuons) {
235  endcapMuons->push_back(0, mu);
236  }
237 
238  for (int i = 0; i < 1008; ++i) {
239  // from where could I take the tower energies?
240  int energy = int(m_rnd.Gaus(12, 6));
241  if (energy < 0) energy = 0;
242  if (energy > 31) energy = 31;
243  MuonCaloSum sum(energy, i/28, i%28, i);
244  towerSums->push_back(0, sum);
245  }
246 
247  iEvent.put(std::move(barrelMuons), "BarrelTFMuons");
248  iEvent.put(std::move(overlapMuons), "OverlapTFMuons");
249  iEvent.put(std::move(endcapMuons), "ForwardTFMuons");
250  iEvent.put(std::move(towerSums), "TriggerTowerSums");
251  m_currEvt++;
252 
253 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:29
l1t::RegionalMuonCand::setHwPhi
void setHwPhi(int bits)
Set compressed relative phi as transmitted by hardware LSB = 2*pi/576 (8 bits)
Definition: RegionalMuonCand.h:67
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:125
mps_fire.i
i
Definition: mps_fire.py:338
l1t::RegionalMuonCandBxCollection
BXVector< RegionalMuonCand > RegionalMuonCandBxCollection
Definition: RegionalMuonCandFwd.h:11
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
l1t::MuonCaloSumBxCollection
BXVector< MuonCaloSum > MuonCaloSumBxCollection
Definition: MuonCaloSumFwd.h:7
reco::Candidate::status
virtual int status() const =0
status word
l1t::RegionalMuonCand::setTFIdentifiers
void setTFIdentifiers(int processor, tftype trackFinder)
Set the processor ID, track-finder type. From these two, the link is set.
Definition: RegionalMuonCand.cc:6
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
L1TMicroGMTInputProducerFromGen::genParticlesToken
edm::EDGetTokenT< reco::GenParticleCollection > genParticlesToken
Definition: L1TMicroGMTInputProducerFromGen.cc:69
RPCpg::mu
const int mu
Definition: Constants.h:22
GenHFHadronMatcher_cfi.genParticles
genParticles
Definition: GenHFHadronMatcher_cfi.py:4
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
L1TMicroGMTInputProducerFromGen::compareMuons
static bool compareMuons(const RegionalMuonCand &, const RegionalMuonCand &)
Definition: L1TMicroGMTInputProducerFromGen.cc:112
l1t::RegionalMuonCand::setHwQual
void setHwQual(int bits)
Set compressed quality code as transmitted by hardware (4 bits)
Definition: RegionalMuonCand.h:75
l1t::RegionalMuonCand::setHwPt
void setHwPt(int bits)
Set compressed pT as transmitted by hardware LSB = 0.5 (9 bits)
Definition: RegionalMuonCand.h:61
l1t::RegionalMuonCand::setHwEta
void setHwEta(int bits)
Set compressed eta as transmitted by hardware LSB = 0.010875 (9 bits)
Definition: RegionalMuonCand.h:69
reco::Candidate::eta
virtual double eta() const =0
momentum pseudorapidity
reco::Candidate::pt
virtual double pt() const =0
transverse momentum
reco::Candidate::charge
virtual int charge() const =0
electric charge
edm
HLT enums.
Definition: AlignableModifier.h:17
MuonErrorMatrixAnalyzer_cfi.maxPt
maxPt
Definition: MuonErrorMatrixAnalyzer_cfi.py:19
l1t::RegionalMuonCand::setHwSignValid
void setHwSignValid(int bits)
Set whether charge measurement is valid (0 for high pT muons)
Definition: RegionalMuonCand.h:73
Geom::twoPi
constexpr double twoPi()
Definition: Pi.h:32
reco::Candidate::phi
virtual double phi() const =0
momentum azimuthal angle
eostools.move
def move(src, dest)
Definition: eostools.py:511
l1t::RegionalMuonCand::setHwSign
void setHwSign(int bits)
Set charge sign bit (charge = (-1)^(sign))
Definition: RegionalMuonCand.h:71

Member Data Documentation

edm::EDGetTokenT<reco::GenParticleCollection> L1TMicroGMTInputProducerFromGen::genParticlesToken
private

Definition at line 69 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by L1TMicroGMTInputProducerFromGen(), and produce().

int L1TMicroGMTInputProducerFromGen::m_currEvt
private

Definition at line 70 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

const int L1TMicroGMTInputProducerFromGen::m_maxMuons = 108
staticprivate

Definition at line 71 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

TRandom3 L1TMicroGMTInputProducerFromGen::m_rnd
private

Definition at line 72 of file L1TMicroGMTInputProducerFromGen.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11