CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
L1TkHTMissEmulatorProducer Class Reference
Inheritance diagram for L1TkHTMissEmulatorProducer:
edm::stream::EDProducer<>

Public Member Functions

 L1TkHTMissEmulatorProducer (const edm::ParameterSet &)
 
 ~L1TkHTMissEmulatorProducer () override=default
 
- 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 EDProduceroperator= (const EDProducer &)=delete
 

Private Member Functions

virtual void beginJob ()
 
virtual void endJob ()
 
void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

int aSteps = 8
 
std::vector< l1tmhtemu::MHTphi_tatanLUT_
 
std::vector< l1tmhtemu::phi_tcosLUT_
 
int cosLUTbins
 
bool debug_ = false
 
bool displaced_
 
l1tmhtemu::eta_t jetMaxEta_
 
l1tmhtemu::pt_t jetMinPt_
 
const edm::EDGetTokenT< TkJetWordCollectionjetToken_
 
std::string L1MHTCollectionName_
 
std::vector< l1tmhtemu::Et_tmagNormalisationLUT_
 
l1tmhtemu::ntracks_t minNtracksHighPt_
 
l1tmhtemu::ntracks_t minNtracksLowPt_
 

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 34 of file L1TkHTMissEmulatorProducer.cc.

Constructor & Destructor Documentation

◆ L1TkHTMissEmulatorProducer()

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

Definition at line 66 of file L1TkHTMissEmulatorProducer.cc.

References aSteps, atanLUT_, cosLUT_, cosLUTbins, debug_, displaced_, ALCARECOEcalPhiSym_cff::float, l1tmhtemu::generateaTanLUT(), l1tmhtemu::generateCosLUT(), l1tmhtemu::generatemagNormalisationLUT(), edm::ParameterSet::getParameter(), createfilelist::int, jetMaxEta_, jetMinPt_, l1tmhtemu::kInternalEtaWidth, l1tmhtemu::kInternalPtWidth, l1tmhtemu::kMaxCosLUTPhi, l1tmhtemu::kStepEta, l1tmhtemu::kStepPhi, l1tmhtemu::kStepPt, L1MHTCollectionName_, magNormalisationLUT_, minNtracksHighPt_, minNtracksLowPt_, and AlCaHLTBitMon_QueryRunRegistry::string.

67  : jetToken_(consumes<TkJetWordCollection>(iConfig.getParameter<edm::InputTag>("L1TkJetEmulationInputTag"))) {
68  debug_ = iConfig.getParameter<bool>("debug");
69  displaced_ = iConfig.getParameter<bool>("displaced");
70 
71  jetMinPt_ = l1tmhtemu::digitizeSignedValue<l1tmhtemu::pt_t>(
72  (float)iConfig.getParameter<double>("jet_minPt"), l1tmhtemu::kInternalPtWidth, l1tmhtemu::kStepPt);
73  jetMaxEta_ = l1tmhtemu::digitizeSignedValue<l1tmhtemu::eta_t>(
74  (float)iConfig.getParameter<double>("jet_maxEta"), l1tmhtemu::kInternalEtaWidth, l1tmhtemu::kStepEta);
75  minNtracksHighPt_ = (l1tmhtemu::ntracks_t)iConfig.getParameter<int>("jet_minNtracksHighPt");
76  minNtracksLowPt_ = (l1tmhtemu::ntracks_t)iConfig.getParameter<int>("jet_minNtracksLowPt");
77 
80 
83 
84  // Name of output ED Product
85  L1MHTCollectionName_ = (std::string)iConfig.getParameter<std::string>("L1MHTCollectionName");
86 
87  produces<std::vector<EtSum>>(L1MHTCollectionName_);
88 
89  if (debug_) {
90  edm::LogVerbatim("L1TrackerHTMissEmulatorProducer")
91  << "-------------------------------------------------------------------------\n"
92  << "====BITWIDTHS====\n"
93  << "pt: " << l1t::TkJetWord::TkJetBitWidths::kPtSize << " eta: " << l1t::TkJetWord::TkJetBitWidths::kGlbEtaSize
94  << " phi:" << l1t::TkJetWord::TkJetBitWidths::kGlbPhiSize << "\n"
95  << "====CUT AP_INTS====\n"
96  << "minpt: " << jetMinPt_ << " maxeta: " << jetMaxEta_ << " minNtracksHighPt: " << minNtracksHighPt_
97  << " minNtracksLowPt: " << minNtracksLowPt_ << "\n"
98  << "====CUT AP_INTS TO FLOATS====\n"
99  << "minpt: " << (float)jetMinPt_ * l1tmhtemu::kStepPt << " maxeta: " << (float)jetMaxEta_ * l1tmhtemu::kStepEta
100  << " minNtracksHighPt: " << (int)minNtracksHighPt_ << " minNtracksLowPt: " << (int)minNtracksLowPt_ << "\n"
101  << "-------------------------------------------------------------------------\n";
102  }
103 }
std::vector< l1tmhtemu::MHTphi_t > atanLUT_
Log< level::Info, true > LogVerbatim
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
const edm::EDGetTokenT< TkJetWordCollection > jetToken_
ap_uint< 5 > ntracks_t
std::vector< MHTphi_t > generateaTanLUT(int cordicSteps)
const float kMaxCosLUTPhi
std::vector< l1tmhtemu::Et_t > magNormalisationLUT_
const unsigned int kInternalPtWidth
const unsigned int kInternalEtaWidth
std::vector< l1tmhtemu::phi_t > cosLUT_
std::vector< phi_t > generateCosLUT(unsigned int size)
std::vector< Et_t > generatemagNormalisationLUT(int cordicSteps)

◆ ~L1TkHTMissEmulatorProducer()

L1TkHTMissEmulatorProducer::~L1TkHTMissEmulatorProducer ( )
overridedefault

Member Function Documentation

◆ beginJob()

void L1TkHTMissEmulatorProducer::beginJob ( void  )
privatevirtual

Definition at line 269 of file L1TkHTMissEmulatorProducer.cc.

269 {}

◆ endJob()

void L1TkHTMissEmulatorProducer::endJob ( void  )
privatevirtual

Definition at line 271 of file L1TkHTMissEmulatorProducer.cc.

271 {}

◆ produce()

void L1TkHTMissEmulatorProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivate

Definition at line 105 of file L1TkHTMissEmulatorProducer.cc.

References aSteps, Matriplex::atan2(), atanLUT_, l1tmhtemu::cordicSqrt(), funct::cos(), cosLUT_, cosLUTbins, debug_, displaced_, l1tmhtemu::EtMiss::Et, ALCARECOEcalPhiSym_cff::float, iEvent, createfilelist::int, edm::HandleBase::isValid(), jetMaxEta_, jetMinPt_, jetToken_, l1tmhtemu::kInternalEtaWidth, l1tmhtemu::kInternalPhiWidth, l1tmhtemu::kInternalPtWidth, l1tmhtemu::kMHTPhiBins, L1Analysis::kMissingHt, l1tmhtemu::kStepEta, l1tmhtemu::kStepMHTPhi, l1tmhtemu::kStepPhi, l1tmhtemu::kStepPt, L1MHTCollectionName_, M_PI, magNormalisationLUT_, eostools::move(), or, phi, l1tmhtemu::EtMiss::Phi, funct::sin(), and mathSSE::sqrt().

105  {
106  using namespace edm;
107  std::unique_ptr<std::vector<l1t::EtSum>> MHTCollection(new std::vector<l1t::EtSum>(0));
108 
109  // L1 track-trigger jets
110  edm::Handle<TkJetWordCollection> L1TkJetsHandle;
111  iEvent.getByToken(jetToken_, L1TkJetsHandle);
112  std::vector<TkJetWord>::const_iterator jetIter;
113 
114  if (!L1TkJetsHandle.isValid() && !displaced_) {
115  LogError("TkHTMissEmulatorProducer") << "\nWarning: TkJetCollection not found in the event. Exit\n";
116  return;
117  }
118 
119  if (!L1TkJetsHandle.isValid() && displaced_) {
120  LogError("TkHTMissEmulatorProducer") << "\nWarning: TkJetExtendedCollection not found in the event. Exit\n";
121  return;
122  }
123 
124  // floats used for debugging
125  float sumPx_ = 0;
126  float sumPy_ = 0;
127  float HT_ = 0;
128 
129  l1tmhtemu::Et_t sumPx = 0;
130  l1tmhtemu::Et_t sumPy = 0;
131  l1tmhtemu::MHT_t HT = 0;
132 
133  // loop over jets
134  int jetn = 0;
135 
136  for (jetIter = L1TkJetsHandle->begin(); jetIter != L1TkJetsHandle->end(); ++jetIter) {
137  // floats used for debugging
138  float tmp_jet_px_ = jetIter->pt() * cos(jetIter->glbphi());
139  float tmp_jet_py_ = jetIter->pt() * sin(jetIter->glbphi());
140  //float tmp_jet_et_ = jetIter->pt(); // FIXME Get Et from the emulated jets
141  float tmp_jet_pt_ = jetIter->pt();
142 
143  bool tmp_jet_isDisplaced_ = jetIter->dispflag();
144 
145  l1tmhtemu::pt_t tmp_jet_pt =
146  l1tmhtemu::digitizeSignedValue<l1tmhtemu::pt_t>(jetIter->pt(), l1tmhtemu::kInternalPtWidth, l1tmhtemu::kStepPt);
147  l1tmhtemu::eta_t tmp_jet_eta = l1tmhtemu::digitizeSignedValue<l1tmhtemu::eta_t>(
149  l1tmhtemu::phi_t tmp_jet_phi = l1tmhtemu::digitizeSignedValue<l1tmhtemu::phi_t>(
151  l1tmhtemu::ntracks_t tmp_jet_nt = l1tmhtemu::ntracks_t(jetIter->nt());
152 
153  l1tmhtemu::phi_t tmp_jet_cos_phi = l1tmhtemu::phi_t(-999);
154  l1tmhtemu::phi_t tmp_jet_sin_phi = l1tmhtemu::phi_t(-999);
155 
156  if (tmp_jet_phi >= 0) {
157  tmp_jet_cos_phi = cosLUT_[tmp_jet_phi];
158 
159  if (cosLUTbins / 2 + 1 - tmp_jet_phi >= 0)
160  tmp_jet_sin_phi = cosLUT_[cosLUTbins / 2 + 1 - tmp_jet_phi];
161  else
162  tmp_jet_sin_phi = cosLUT_[-1 * (cosLUTbins / 2 + 1 - tmp_jet_phi)];
163 
164  } else {
165  tmp_jet_cos_phi = cosLUT_[-1 * tmp_jet_phi];
166 
167  if (cosLUTbins / 2 + 1 - (-1 * tmp_jet_phi) >= 0)
168  tmp_jet_sin_phi = -1 * cosLUT_[cosLUTbins / 2 + 1 - (-1 * tmp_jet_phi)];
169  else
170  tmp_jet_sin_phi = -1 * cosLUT_[-1 * (cosLUTbins / 2 + 1 - (-1 * tmp_jet_phi))];
171  }
172 
173  l1tmhtemu::Et_t tmp_jet_px = tmp_jet_pt * tmp_jet_cos_phi;
174  l1tmhtemu::Et_t tmp_jet_py = tmp_jet_pt * tmp_jet_sin_phi;
175 
176  jetn++;
177 
178  if (debug_) {
179  edm::LogVerbatim("L1TrackerHTMissEmulatorProducer")
180  << "****JET EMULATION" << jetn << "****\n"
181  << "FLOATS ORIGINAL\n"
182  << "PT: " << jetIter->pt() << "| ETA: " << jetIter->glbeta() << "| PHI: " << jetIter->glbphi()
183  << "| NTRACKS: " << jetIter->nt() << "| COS(PHI): " << cos(jetIter->glbphi())
184  << "| SIN(PHI): " << sin(jetIter->glbphi()) << "| Px: " << jetIter->pt() * cos(jetIter->glbphi())
185  << "| Py: " << jetIter->pt() * sin(jetIter->glbphi()) << "\n"
186  << "AP_INTS RAW\n"
187  << "PT: " << jetIter->ptWord() << "| ETA: " << jetIter->glbEtaWord() << "| PHI: " << jetIter->glbPhiWord()
188  << "| NTRACKS: " << jetIter->ntWord() << "\n"
189  << "AP_INTS NEW\n"
190  << "PT: " << tmp_jet_pt << "| ETA: " << tmp_jet_eta << "| PHI: " << tmp_jet_phi << "| NTRACKS: " << tmp_jet_nt
191  << "| COS(PHI): " << tmp_jet_cos_phi << "| SIN(PHI): " << tmp_jet_sin_phi << "| Px: " << tmp_jet_px
192  << "| Py: " << tmp_jet_py << "\n"
193  << "AP_INTS NEW TO FLOATS\n"
194  << "PT: " << (float)tmp_jet_pt * l1tmhtemu::kStepPt << "| ETA: " << (float)tmp_jet_eta * l1tmhtemu::kStepEta
195  << "| PHI: " << (float)tmp_jet_phi * l1tmhtemu::kStepPhi << "| NTRACKS: " << (int)tmp_jet_nt
196  << "| COS(PHI): " << (float)tmp_jet_cos_phi * l1tmhtemu::kStepPhi
197  << "| SIN(PHI): " << (float)tmp_jet_sin_phi * l1tmhtemu::kStepPhi
198  << "| Px: " << (float)tmp_jet_px * l1tmhtemu::kStepPt * l1tmhtemu::kStepPhi
199  << "| Py: " << (float)tmp_jet_py * l1tmhtemu::kStepPt * l1tmhtemu::kStepPhi << "\n"
200  << "-------------------------------------------------------------------------\n";
201  }
202 
203  if (tmp_jet_pt < jetMinPt_)
204  continue;
205  if (tmp_jet_eta > jetMaxEta_ or tmp_jet_eta < -1 * jetMaxEta_)
206  continue;
207  if (tmp_jet_nt < minNtracksLowPt_ && tmp_jet_pt > 200)
208  continue;
209  if (tmp_jet_nt < minNtracksHighPt_ && tmp_jet_pt > 400)
210  continue;
211  if (displaced_ && !tmp_jet_isDisplaced_)
212  continue;
213 
214  if (debug_) {
215  sumPx_ += tmp_jet_px_;
216  sumPy_ += tmp_jet_py_;
217  HT_ += tmp_jet_pt_;
218  }
219 
220  sumPx += tmp_jet_pt * tmp_jet_cos_phi;
221  sumPy += tmp_jet_pt * tmp_jet_sin_phi;
222  HT += tmp_jet_pt;
223 
224  } // end jet loop
225 
226  // define missing HT
227 
228  // Perform cordic sqrt, take x,y and converts to polar coordinate r,phi where
229  // r=sqrt(x**2+y**2) and phi = atan(y/x)
231  math::XYZTLorentzVector missingEt(-sumPx, -sumPy, 0, EtMiss.Et);
232 
234 
235  if ((sumPx < 0) && (sumPy < 0))
236  phi = EtMiss.Phi - l1tmhtemu::kMHTPhiBins / 2;
237  else if ((sumPx >= 0) && (sumPy >= 0))
238  phi = (EtMiss.Phi) + l1tmhtemu::kMHTPhiBins / 2;
239  else if ((sumPx >= 0) && (sumPy < 0))
240  phi = EtMiss.Phi - l1tmhtemu::kMHTPhiBins / 2;
241  else if ((sumPx < 0) && (sumPy >= 0))
242  phi = EtMiss.Phi - 3 * l1tmhtemu::kMHTPhiBins / 2;
243 
244  if (debug_) {
245  edm::LogVerbatim("L1TrackerHTMissEmulatorProducer")
246  << "-------------------------------------------------------------------------\n"
247  << "====MHT FLOATS====\n"
248  << "sumPx: " << sumPx_ << "| sumPy: " << sumPy_ << "| ET: " << sqrt(sumPx_ * sumPx_ + sumPy_ * sumPy_)
249  << "| HT: " << HT_ << "| PHI: " << atan2(sumPy_, sumPx_) << "\n"
250  << "====MHT AP_INTS====\n"
251  << "sumPx: " << sumPx << "| sumPy: " << sumPy << "| ET: " << EtMiss.Et << "| HT: " << HT << "| PHI: " << phi
252  << "\n"
253  << "====MHT AP_INTS TO FLOATS====\n"
254  << "sumPx: " << (float)sumPx * l1tmhtemu::kStepPt * l1tmhtemu::kStepPhi
255  << "| sumPy: " << (float)sumPy * l1tmhtemu::kStepPt * l1tmhtemu::kStepPhi << "| ET: " << EtMiss.Et.to_double()
256  << "| HT: " << (float)HT * l1tmhtemu::kStepPt << "| PHI: " << (float)phi * l1tmhtemu::kStepMHTPhi - M_PI << "\n"
257  << "-------------------------------------------------------------------------\n";
258  }
259  //rescale HT to correct output range
260  HT = HT / (int)(1 / l1tmhtemu::kStepPt);
261 
262  EtSum L1HTSum(missingEt, EtSum::EtSumType::kMissingHt, (int)HT.range(), 0, (int)phi, (int)jetn);
263 
264  MHTCollection->push_back(L1HTSum);
266 
267 } //end producer
std::vector< l1tmhtemu::MHTphi_t > atanLUT_
Log< level::Info, true > LogVerbatim
const edm::EDGetTokenT< TkJetWordCollection > jetToken_
ap_uint< 5 > ntracks_t
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
ap_uint< kMHTPhiSize > MHTphi_t
EtMiss cordicSqrt(Et_t x, Et_t y, int cordicSteps, std::vector< l1tmhtemu::MHTphi_t > atanLUT, std::vector< Et_t > magNormalisationLUT)
Log< level::Error, false > LogError
const unsigned int kInternalPhiWidth
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ap_uint< kInternalPtWidth > pt_t
int iEvent
Definition: GenABIO.cc:224
const double kStepMHTPhi
std::vector< l1tmhtemu::Et_t > magNormalisationLUT_
ap_ufixed< kMHTSize, kMHTIntSize > MHT_t
T sqrt(T t)
Definition: SSEVec.h:23
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
std::vector< pat::MHT > MHTCollection
Definition: MHT.h:39
#define M_PI
const unsigned int kInternalPtWidth
ap_int< kInternalEtaWidth > eta_t
const unsigned int kInternalEtaWidth
bool isValid() const
Definition: HandleBase.h:70
std::vector< l1tmhtemu::phi_t > cosLUT_
HLT enums.
const unsigned int kMHTPhiBins
Definition: HT.h:21
ap_int< kInternalPhiWidth > phi_t
ap_int< kHTSize > Et_t
def move(src, dest)
Definition: eostools.py:511
MPlex< T, D1, D2, N > atan2(const MPlex< T, D1, D2, N > &y, const MPlex< T, D1, D2, N > &x)
Definition: Matriplex.h:648

Member Data Documentation

◆ aSteps

int L1TkHTMissEmulatorProducer::aSteps = 8
private

Definition at line 50 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ atanLUT_

std::vector<l1tmhtemu::MHTphi_t> L1TkHTMissEmulatorProducer::atanLUT_
private

Definition at line 58 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ cosLUT_

std::vector<l1tmhtemu::phi_t> L1TkHTMissEmulatorProducer::cosLUT_
private

Definition at line 57 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ cosLUTbins

int L1TkHTMissEmulatorProducer::cosLUTbins
private

Definition at line 49 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ debug_

bool L1TkHTMissEmulatorProducer::debug_ = false
private

Definition at line 46 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ displaced_

bool L1TkHTMissEmulatorProducer::displaced_
private

Definition at line 47 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ jetMaxEta_

l1tmhtemu::eta_t L1TkHTMissEmulatorProducer::jetMaxEta_
private

Definition at line 53 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ jetMinPt_

l1tmhtemu::pt_t L1TkHTMissEmulatorProducer::jetMinPt_
private

Definition at line 52 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ jetToken_

const edm::EDGetTokenT<TkJetWordCollection> L1TkHTMissEmulatorProducer::jetToken_
private

Definition at line 63 of file L1TkHTMissEmulatorProducer.cc.

Referenced by produce().

◆ L1MHTCollectionName_

std::string L1TkHTMissEmulatorProducer::L1MHTCollectionName_
private

Definition at line 61 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ magNormalisationLUT_

std::vector<l1tmhtemu::Et_t> L1TkHTMissEmulatorProducer::magNormalisationLUT_
private

Definition at line 59 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer(), and produce().

◆ minNtracksHighPt_

l1tmhtemu::ntracks_t L1TkHTMissEmulatorProducer::minNtracksHighPt_
private

Definition at line 54 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer().

◆ minNtracksLowPt_

l1tmhtemu::ntracks_t L1TkHTMissEmulatorProducer::minNtracksLowPt_
private

Definition at line 55 of file L1TkHTMissEmulatorProducer.cc.

Referenced by L1TkHTMissEmulatorProducer().