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

#include <L1TrackerEtMissEmulatorProducer.cc>

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

Public Types

typedef std::vector< L1TTTrackType > L1TTTrackCollectionType
 
typedef edm::RefVector< L1TTTrackCollectionType > L1TTTrackRefCollectionType
 
typedef TTTrack< Ref_Phase2TrackerDigi_ > L1TTTrackType
 
typedef l1t::VertexWordCollection L1VertexCollectionType
 
typedef l1t::VertexWord L1VertexType
 
- 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
 

Public Member Functions

 L1TrackerEtMissEmulatorProducer (const edm::ParameterSet &)
 
 ~L1TrackerEtMissEmulatorProducer () override
 
- 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)
 

Private Member Functions

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

Private Attributes

bool cordicDebug_ = false
 
int cordicSteps_
 
std::vector< l1tmetemu::global_phi_t > cosLUT_
 
int debug_
 
l1tmetemu::z_t deltaZ0_ = 0
 
bool GTTinput_ = false
 
std::string L1MetCollectionName_
 
std::vector< l1tmetemu::global_phi_t > phiQuadrants_
 
std::vector< l1tmetemu::global_phi_t > phiShifts_
 
const edm::EDGetTokenT< L1VertexCollectionType > pvToken_
 
const edm::EDGetTokenT< L1TTTrackRefCollectionType > trackToken_
 
L1TkEtMissEmuTrackTransform TrackTransform
 
const edm::EDGetTokenT< L1TTTrackRefCollectionType > vtxAssocTrackToken_
 

Detailed Description

L1Trigger/L1TTrackMatch/plugins/L1TrackerEtMissEmulatorProducer.cc

Description: Takes L1TTTracks and performs a integer emulation of Track-based missing Et, outputting a collection of EtSum

Definition at line 41 of file L1TrackerEtMissEmulatorProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

typedef std::vector<L1TTTrackType> L1TrackerEtMissEmulatorProducer::L1TTTrackCollectionType

Definition at line 44 of file L1TrackerEtMissEmulatorProducer.cc.

◆ L1TTTrackRefCollectionType

typedef edm::RefVector<L1TTTrackCollectionType> L1TrackerEtMissEmulatorProducer::L1TTTrackRefCollectionType

Definition at line 45 of file L1TrackerEtMissEmulatorProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackerEtMissEmulatorProducer::L1TTTrackType

Definition at line 43 of file L1TrackerEtMissEmulatorProducer.cc.

◆ L1VertexCollectionType

typedef l1t::VertexWordCollection L1TrackerEtMissEmulatorProducer::L1VertexCollectionType

Definition at line 46 of file L1TrackerEtMissEmulatorProducer.cc.

◆ L1VertexType

typedef l1t::VertexWord L1TrackerEtMissEmulatorProducer::L1VertexType

Definition at line 47 of file L1TrackerEtMissEmulatorProducer.cc.

Constructor & Destructor Documentation

◆ L1TrackerEtMissEmulatorProducer()

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

Definition at line 83 of file L1TrackerEtMissEmulatorProducer.cc.

References electrons_cff::bool, cordicDebug_, cordicSteps_, cosLUT_, debug_, l1tmetemu::generateCosLUT(), L1TkEtMissEmuTrackTransform::generateLUTs(), edm::ParameterSet::getParameter(), L1TkEtMissEmuTrackTransform::getPhiQuad(), L1TkEtMissEmuTrackTransform::getPhiShift(), GTTinput_, createfilelist::int, l1tmetemu::kMaxCosLUTPhi, l1tmetemu::kStepPhi, L1MetCollectionName_, phiQuadrants_, phiShifts_, L1TkEtMissEmuTrackTransform::setGTTinput(), AlCaHLTBitMon_QueryRunRegistry::string, and TrackTransform.

84  : pvToken_(consumes<L1VertexCollectionType>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))),
85  trackToken_(consumes<L1TTTrackRefCollectionType>(iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))),
86  vtxAssocTrackToken_(
87  consumes<L1TTTrackRefCollectionType>(iConfig.getParameter<edm::InputTag>("L1TrackAssociatedInputTag"))) {
88  // Setup LUTs
89  TrackTransform.generateLUTs();
90  phiQuadrants_ = TrackTransform.getPhiQuad();
91  phiShifts_ = TrackTransform.getPhiShift();
92 
93  // Get Emulator config parameters
94  cordicSteps_ = (int)iConfig.getParameter<int>("nCordicSteps");
95  debug_ = (int)iConfig.getParameter<int>("debug");
96 
97  GTTinput_ = (bool)iConfig.getParameter<bool>("useGTTinput");
98 
99  TrackTransform.setGTTinput(GTTinput_);
100 
101  // Name of output ED Product
102  L1MetCollectionName_ = (std::string)iConfig.getParameter<std::string>("L1MetCollectionName");
103 
104  if (debug_ == 5) {
105  cordicDebug_ = true;
106  }
107 
108  // To have same bin spacing between 0 and pi/2 as between original phi
109  // granularity
110  int cosLUTbins = floor(l1tmetemu::kMaxCosLUTPhi / l1tmetemu::kStepPhi);
111 
112  // Compute LUTs
113  cosLUT_ = l1tmetemu::generateCosLUT(cosLUTbins);
114 
115  produces<std::vector<EtSum>>(L1MetCollectionName_);
116 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
L1TrackerEtMissEmulatorProducer::phiQuadrants_
std::vector< l1tmetemu::global_phi_t > phiQuadrants_
Definition: L1TrackerEtMissEmulatorProducer.cc:62
L1TkEtMissEmuTrackTransform::getPhiShift
std::vector< l1tmetemu::global_phi_t > getPhiShift() const
Definition: L1TkEtMissEmuTrackTransform.h:61
L1TrackerEtMissEmulatorProducer::pvToken_
const edm::EDGetTokenT< L1VertexCollectionType > pvToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:77
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
L1TkEtMissEmuTrackTransform::getPhiQuad
std::vector< l1tmetemu::global_phi_t > getPhiQuad() const
Definition: L1TkEtMissEmuTrackTransform.h:60
L1TrackerEtMissEmulatorProducer::vtxAssocTrackToken_
const edm::EDGetTokenT< L1TTTrackRefCollectionType > vtxAssocTrackToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:79
L1TrackerEtMissEmulatorProducer::trackToken_
const edm::EDGetTokenT< L1TTTrackRefCollectionType > trackToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:78
electrons_cff.bool
bool
Definition: electrons_cff.py:381
L1TrackerEtMissEmulatorProducer::phiShifts_
std::vector< l1tmetemu::global_phi_t > phiShifts_
Definition: L1TrackerEtMissEmulatorProducer.cc:63
l1tmetemu::generateCosLUT
std::vector< global_phi_t > generateCosLUT(unsigned int size)
Definition: L1TkEtMissEmuAlgo.cc:6
l1tmetemu::kMaxCosLUTPhi
const float kMaxCosLUTPhi
Definition: L1TkEtMissEmuAlgo.h:78
L1TrackerEtMissEmulatorProducer::cosLUT_
std::vector< l1tmetemu::global_phi_t > cosLUT_
Definition: L1TrackerEtMissEmulatorProducer.cc:61
l1tmetemu::kStepPhi
const double kStepPhi
Definition: L1TkEtMissEmuAlgo.h:71

◆ ~L1TrackerEtMissEmulatorProducer()

L1TrackerEtMissEmulatorProducer::~L1TrackerEtMissEmulatorProducer ( )
override

Definition at line 118 of file L1TrackerEtMissEmulatorProducer.cc.

118 {}

Member Function Documentation

◆ beginJob()

void L1TrackerEtMissEmulatorProducer::beginJob ( void  )
privatevirtual

Definition at line 356 of file L1TrackerEtMissEmulatorProducer.cc.

356 {}

◆ endJob()

void L1TrackerEtMissEmulatorProducer::endJob ( void  )
privatevirtual

Definition at line 358 of file L1TrackerEtMissEmulatorProducer.cc.

358 {}

◆ fillDescriptions()

void L1TrackerEtMissEmulatorProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 361 of file L1TrackerEtMissEmulatorProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

361  {
362  //The following says we do not know what parameters are allowed so do no validation
363  // Please change this to state exactly what you do use, even if it is no parameters
364  edm::ParameterSetDescription desc;
365  desc.setUnknown();
366  descriptions.addDefault(desc);
367 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99

◆ produce()

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

Definition at line 120 of file L1TrackerEtMissEmulatorProducer.cc.

References InternalEtWord::bendChi2, InternalEtWord::chi2rphidof, InternalEtWord::chi2rzdof, cordicDebug_, l1tmhtemu::cordicSqrt(), cordicSteps_, funct::cos(), cosLUT_, debug_, l1tmetemu::EtMiss::Et, InternalEtWord::eta, InternalEtWord::EtaSector, spr::find(), dqmMemoryStats::float, InternalEtWord::globalPhi, mps_fire::i, iEvent, edm::HandleBase::isValid(), l1tmetemu::kGlobalPhiBins, l1tmetemu::kMETPhiBins, l1tmetemu::kMETSize, L1Analysis::kMissingEt, l1tmetemu::kNSector, l1tmetemu::kStepMET, l1tmetemu::kStepMETPhi, l1tmetemu::kStepPt, L1MetCollectionName_, M_PI, eostools::move(), InternalEtWord::nstubs, InternalEtWord::phi, l1tmetemu::EtMiss::Phi, phiQuadrants_, InternalEtWord::pt, InternalEtWord::pV, pvToken_, InternalEtWord::Sector, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, cond::impl::to_string(), HLT_2022v14_cff::track, trackToken_, TrackTransform, L1TkEtMissEmuTrackTransform::transformTrack(), extraflags_cff::vtx, vtxAssocTrackToken_, and InternalEtWord::z0.

120  {
121  using namespace edm;
122 
123  std::unique_ptr<std::vector<l1t::EtSum>> METCollection(new std::vector<l1t::EtSum>(0));
124 
125  edm::Handle<L1VertexCollectionType> L1VertexHandle;
126  iEvent.getByToken(pvToken_, L1VertexHandle);
127 
128  edm::Handle<L1TTTrackRefCollectionType> L1TTTrackHandle;
129  iEvent.getByToken(trackToken_, L1TTTrackHandle);
130 
131  edm::Handle<L1TTTrackRefCollectionType> L1TTTrackAssociatedHandle;
132  iEvent.getByToken(vtxAssocTrackToken_, L1TTTrackAssociatedHandle);
133 
134  // Initialize cordic class
135  Cordic cordicSqrt(l1tmetemu::kMETPhiBins, l1tmetemu::kMETSize, cordicSteps_, cordicDebug_);
136 
137  if (!L1VertexHandle.isValid()) {
138  LogError("L1TrackerEtMissEmulatorProducer") << "\nWarning: VertexCollection not found in the event. Exit\n";
139  return;
140  }
141 
142  if (!L1TTTrackHandle.isValid()) {
143  LogError("L1TrackerEtMissEmulatorProducer") << "\nWarning: L1TTTrackCollection not found in the event. Exit\n";
144  return;
145  }
146 
147  if (!L1TTTrackAssociatedHandle.isValid()) {
148  LogError("L1TrackerEtMissEmulatorProducer")
149  << "\nWarning: L1TTTrackAssociatedCollection not found in the event. Exit\n";
150  return;
151  }
152 
153  // Initialize sector sums, need 0 initialization in case a sector has no
154  // tracks
155  l1tmetemu::Et_t sumPx[l1tmetemu::kNSector * 2] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
156  l1tmetemu::Et_t sumPy[l1tmetemu::kNSector * 2] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
157  int sector_totals[l1tmetemu::kNSector] = {0, 0, 0, 0, 0, 0, 0, 0, 0};
158 
159  // Track counters
160  int num_quality_tracks{0};
161  int num_assoc_tracks{0};
162 
163  // Get reference to first vertex in event vertex collection
164  L1VertexType& vtx = const_cast<L1VertexType&>(L1VertexHandle->at(0));
165 
166  for (const auto& track : *L1TTTrackHandle) {
167  num_quality_tracks++;
168  L1TTTrackType& track_ref = const_cast<L1TTTrackType&>(*track); // Get Reference to track to pass to TrackTransform
169 
170  // Convert to internal track representation
171  InternalEtWord EtTrack = TrackTransform.transformTrack<L1TTTrackType, L1VertexType>(track_ref, vtx);
172 
173  if (std::find(L1TTTrackAssociatedHandle->begin(), L1TTTrackAssociatedHandle->end(), track) !=
174  L1TTTrackAssociatedHandle->end()) {
175  num_assoc_tracks++;
176  if (debug_ == 7) {
177  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
178  << "Track to Vertex ID: " << num_quality_tracks << "\n"
179  << "Phi Sector: " << EtTrack.Sector << " pT: " << EtTrack.pt << " Phi: " << EtTrack.globalPhi
180  << " TanL: " << EtTrack.eta << " Z0: " << EtTrack.z0 << " Nstub: " << EtTrack.nstubs
181  << " Chi2rphi: " << EtTrack.chi2rphidof << " Chi2rz: " << EtTrack.chi2rzdof
182  << " bendChi2: " << EtTrack.bendChi2 << " PV: " << EtTrack.pV << "\n"
183  << "--------------------------------------------------------------\n";
184  }
185 
186  if (debug_ == 2) {
187  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
188  << "========================Phi debug=================================\n"
189  << "Int pT: " << EtTrack.pt << "\n"
190  << "Int Phi: " << EtTrack.globalPhi << " Float Phi: " << EtTrack.phi
191  << " Actual Float Cos(Phi): " << cos(EtTrack.phi) << " Actual Float Sin(Phi): " << sin(EtTrack.phi) << "\n";
192  }
193  l1tmetemu::Et_t temppx = 0;
194  l1tmetemu::Et_t temppy = 0;
195 
196  // Split tracks in phi quadrants and access cosLUT_, backwards iteration
197  // through cosLUT_ gives sin Sum sector Et -ve when cos or sin phi are -ve
198  sector_totals[EtTrack.Sector] += 1;
199  if (EtTrack.globalPhi >= phiQuadrants_[0] && EtTrack.globalPhi < phiQuadrants_[1]) {
200  temppx = (EtTrack.pt * cosLUT_[EtTrack.globalPhi]);
201  temppy = (EtTrack.pt * cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi]);
202 
203  if (debug_ == 2) {
204  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
205  << "Sector: " << EtTrack.Sector << " Quadrant: " << 1 << "\n"
206  << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): " << cosLUT_[EtTrack.globalPhi]
207  << " Int Sin(Phi): " << cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi] << "\n"
208 
209  << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins
210  << " Float Cos(Phi): " << (float)cosLUT_[EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
211  << " Float Sin(Phi): "
212  << (float)cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins << "\n";
213  }
214  } else if (EtTrack.globalPhi >= phiQuadrants_[1] && EtTrack.globalPhi < phiQuadrants_[2]) {
215  temppx = -(EtTrack.pt * cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi]);
216  temppy = (EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]]);
217 
218  if (debug_ == 2) {
219  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
220  << "Sector: " << EtTrack.Sector << " Quadrant: " << 2 << "\n"
221  << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): -"
222  << cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi]
223  << " Int Sin(Phi): " << cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]] << "\n"
224 
225  << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): -"
226  << (float)cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
227  << " Float Sin(Phi): " << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]] / l1tmetemu::kGlobalPhiBins
228  << "\n";
229  }
230  } else if (EtTrack.globalPhi >= phiQuadrants_[2] && EtTrack.globalPhi < phiQuadrants_[3]) {
231  temppx = -(EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]]);
232  temppy = -(EtTrack.pt * cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi]);
233 
234  if (debug_ == 2) {
235  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
236  << "Sector: " << EtTrack.Sector << " Quadrant: " << 3 << "\n"
237  << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): -" << cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]]
238  << " Int Sin(Phi): -" << cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi] << "\n"
239 
240  << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): -"
241  << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]] / l1tmetemu::kGlobalPhiBins
242  << " Float Sin(Phi): -"
243  << (float)cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins << "\n";
244  }
245 
246  } else if (EtTrack.globalPhi >= phiQuadrants_[3] && EtTrack.globalPhi < phiQuadrants_[4]) {
247  temppx = (EtTrack.pt * cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi]);
248  temppy = -(EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]]);
249 
250  if (debug_ == 2) {
251  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
252  << "Sector: " << EtTrack.Sector << " Quadrant: " << 4 << "\n"
253  << "Int Phi: " << EtTrack.globalPhi
254  << " Int Cos(Phi): " << cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi] << " Int Sin(Phi): -"
255  << cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]] << "\n"
256 
257  << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): "
258  << (float)cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
259  << " Float Sin(Phi): -"
260  << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]] / l1tmetemu::kGlobalPhiBins << "\n";
261  }
262  } else {
263  temppx = 0;
264  temppy = 0;
265  }
266  if (EtTrack.EtaSector) {
267  sumPx[EtTrack.Sector] = sumPx[EtTrack.Sector] + temppx;
268  sumPy[EtTrack.Sector] = sumPy[EtTrack.Sector] + temppy;
269  } else {
270  sumPx[l1tmetemu::kNSector + EtTrack.Sector] = sumPx[l1tmetemu::kNSector + EtTrack.Sector] + temppx;
271  sumPy[l1tmetemu::kNSector + EtTrack.Sector] = sumPy[l1tmetemu::kNSector + EtTrack.Sector] + temppy;
272  }
273  }
274 
275  } // end loop over tracks
276 
277  l1tmetemu::Et_t GlobalPx = 0;
278  l1tmetemu::Et_t GlobalPy = 0;
279 
280  float tempsumPx = 0;
281  float tempsumPy = 0;
282 
283  // Global Et sum as floats to emulate rounding in HW
284  for (unsigned int i = 0; i < l1tmetemu::kNSector * 2; i++) {
285  tempsumPx += floor((float)sumPx[i] / (float)l1tmetemu::kGlobalPhiBins);
286  tempsumPy += floor((float)sumPy[i] / (float)l1tmetemu::kGlobalPhiBins);
287  }
288 
289  // Recast rounded temporary sums into Et_t datatype
290  GlobalPx = tempsumPx;
291  GlobalPy = tempsumPy;
292 
293  // Perform cordic sqrt, take x,y and converts to polar coordinate r,phi where
294  // r=sqrt(x**2+y**2) and phi = atan(y/x)
295  l1tmetemu::EtMiss EtMiss = cordicSqrt.toPolar(GlobalPx, GlobalPy);
296 
297  // Recentre phi
298  l1tmetemu::METphi_t tempPhi = 0;
299 
300  if ((GlobalPx < 0) && (GlobalPy < 0))
301  tempPhi = EtMiss.Phi - l1tmetemu::kMETPhiBins / 2;
302  else if ((GlobalPx >= 0) && (GlobalPy >= 0))
303  tempPhi = (EtMiss.Phi) + l1tmetemu::kMETPhiBins / 2;
304  else if ((GlobalPx >= 0) && (GlobalPy < 0))
305  tempPhi = EtMiss.Phi - l1tmetemu::kMETPhiBins / 2;
306  else if ((GlobalPx < 0) && (GlobalPy >= 0))
307  tempPhi = EtMiss.Phi - 3 * l1tmetemu::kMETPhiBins / 2;
308 
309  if (debug_ == 6) {
310  std::string flpxarray[l1tmetemu::kNSector * 2];
311  std::string flpyarray[l1tmetemu::kNSector * 2];
312 
313  std::string intpxarray[l1tmetemu::kNSector * 2];
314  std::string intpyarray[l1tmetemu::kNSector * 2];
315 
316  std::string totalsarray[l1tmetemu::kNSector * 2];
317 
318  for (unsigned int i = 0; i < l1tmetemu::kNSector * 2; i++) {
319  flpxarray[i] = to_string(sumPx[i] * l1tmetemu::kStepPt) + "|";
320  flpyarray[i] = to_string(sumPy[i] * l1tmetemu::kStepPt) + "|";
321  intpxarray[i] = to_string(floor((float)sumPx[i] / (float)l1tmetemu::kGlobalPhiBins)) + "|";
322  intpyarray[i] = to_string(floor((float)sumPy[i] / (float)l1tmetemu::kGlobalPhiBins)) + "|";
323  totalsarray[i] = to_string(sector_totals[i]) + "|";
324  }
325 
326  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
327  << "====Sector Pt====\n"
328  << "Float Px: " << flpxarray << "Float Py: " << flpyarray << "Integer Px: " << intpyarray
329  << "Integer Px: " << intpyarray << "Sector Totals: " << totalsarray
330 
331  << "====Global Pt====\n"
332  << "Integer Global Px: " << GlobalPx << "| Integer Global Py: " << GlobalPy << "\n"
333  << "Float Global Px: " << GlobalPx * l1tmetemu::kStepPt
334  << "| Float Global Py: " << GlobalPy * l1tmetemu::kStepPt << "\n";
335  }
336 
337  if (debug_ == 6) {
338  edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
339  << "====MET===\n"
340  << "Integer MET: " << EtMiss.Et << "| Integer MET phi: " << EtMiss.Phi << "\n"
341  << "Float MET: " << (EtMiss.Et) * l1tmetemu::kStepMET
342  << "| Float MET phi: " << (float)tempPhi * l1tmetemu::kStepMETPhi - M_PI << "\n"
343  << "# Tracks after Quality Cuts: " << num_quality_tracks << "\n"
344  << "# Tracks Associated to Vertex: " << num_assoc_tracks << "\n"
345  << "========================================================\n";
346  }
347 
348  math::XYZTLorentzVector missingEt(-GlobalPx, -GlobalPy, 0, EtMiss.Et);
349  EtSum L1EtSum(missingEt, EtSum::EtSumType::kMissingEt, (int)EtMiss.Et, 0, (int)tempPhi, (int)num_assoc_tracks);
350 
351  METCollection->push_back(L1EtSum);
352 
353  iEvent.put(std::move(METCollection), L1MetCollectionName_);
354 } // end producer
l1tmetemu::kGlobalPhiBins
const unsigned int kGlobalPhiBins
Definition: L1TkEtMissEmuAlgo.h:53
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
l1tmetemu::METphi_t
ap_uint< kMETPhiSize > METphi_t
Definition: L1TkEtMissEmuAlgo.h:51
mps_fire.i
i
Definition: mps_fire.py:429
L1TrackerEtMissEmulatorProducer::L1TTTrackType
TTTrack< Ref_Phase2TrackerDigi_ > L1TTTrackType
Definition: L1TrackerEtMissEmulatorProducer.cc:43
InternalEtWord::chi2rphidof
TTTrack_TrackWord::chi2rphi_t chi2rphidof
Definition: L1TkEtMissEmuTrackTransform.h:32
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:19
L1TrackerEtMissEmulatorProducer::phiQuadrants_
std::vector< l1tmetemu::global_phi_t > phiQuadrants_
Definition: L1TrackerEtMissEmulatorProducer.cc:62
HLT_2022v14_cff.track
track
Definition: HLT_2022v14_cff.py:9541
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
cond::impl::to_string
std::string to_string(const V &value)
Definition: OMSAccess.h:71
l1tmhtemu::cordicSqrt
EtMiss cordicSqrt(Et_t x, Et_t y, int cordicSteps, std::vector< l1tmhtemu::MHTphi_t > atanLUT, std::vector< Et_t > magNormalisationLUT)
Definition: L1TkHTMissEmulatorProducer.h:104
l1tmetemu::kStepPt
const double kStepPt
Definition: L1TkEtMissEmuAlgo.h:67
l1tmetemu::kStepMET
const double kStepMET
Definition: L1TkEtMissEmuAlgo.h:73
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
reco::METCollection
std::vector< reco::MET > METCollection
collection of MET objects
Definition: METCollection.h:22
InternalEtWord::chi2rzdof
TTTrack_TrackWord::chi2rz_t chi2rzdof
Definition: L1TkEtMissEmuTrackTransform.h:33
L1TrackerEtMissEmulatorProducer::pvToken_
const edm::EDGetTokenT< L1VertexCollectionType > pvToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:77
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
l1tmetemu::kNSector
const unsigned int kNSector
Definition: L1TkEtMissEmuAlgo.h:58
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
iEvent
int iEvent
Definition: GenABIO.cc:224
L1TrackerEtMissEmulatorProducer::vtxAssocTrackToken_
const edm::EDGetTokenT< L1TTTrackRefCollectionType > vtxAssocTrackToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:79
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
L1TrackerEtMissEmulatorProducer::trackToken_
const edm::EDGetTokenT< L1TTTrackRefCollectionType > trackToken_
Definition: L1TrackerEtMissEmulatorProducer.cc:78
l1tmetemu::kMETPhiBins
const unsigned int kMETPhiBins
Definition: L1TkEtMissEmuAlgo.h:55
Cordic
Definition: Cordic.h:7
l1tmetemu::kMETSize
const unsigned int kMETSize
Definition: L1TkEtMissEmuAlgo.h:45
METCollection
Collection of MET.
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
InternalEtWord::nstubs
l1tmetemu::nstub_t nstubs
Definition: L1TkEtMissEmuTrackTransform.h:28
TTTrack
Class to store the L1 Track Trigger tracks.
Definition: TTTrack.h:29
InternalEtWord::globalPhi
l1tmetemu::global_phi_t globalPhi
Definition: L1TkEtMissEmuTrackTransform.h:27
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
L1TrackerEtMissEmulatorProducer::cosLUT_
std::vector< l1tmetemu::global_phi_t > cosLUT_
Definition: L1TrackerEtMissEmulatorProducer.cc:61
edm
HLT enums.
Definition: AlignableModifier.h:19
InternalEtWord::bendChi2
TTTrack_TrackWord::bendChi2_t bendChi2
Definition: L1TkEtMissEmuTrackTransform.h:31
L1TkEtMissEmuTrackTransform::transformTrack
InternalEtWord transformTrack(track &track_ref, vertex &PV)
Definition: L1TkEtMissEmuTrackTransform.h:74
eostools.move
def move(src, dest)
Definition: eostools.py:511
l1tmetemu::Et_t
ap_int< kInternalPtWidth+kEtExtra > Et_t
Definition: L1TkEtMissEmuAlgo.h:40
l1tmetemu::kStepMETPhi
const double kStepMETPhi
Definition: L1TkEtMissEmuAlgo.h:74

Member Data Documentation

◆ cordicDebug_

bool L1TrackerEtMissEmulatorProducer::cordicDebug_ = false
private

Definition at line 69 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ cordicSteps_

int L1TrackerEtMissEmulatorProducer::cordicSteps_
private

Definition at line 67 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ cosLUT_

std::vector<l1tmetemu::global_phi_t> L1TrackerEtMissEmulatorProducer::cosLUT_
private

Definition at line 61 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ debug_

int L1TrackerEtMissEmulatorProducer::debug_
private

Definition at line 68 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ deltaZ0_

l1tmetemu::z_t L1TrackerEtMissEmulatorProducer::deltaZ0_ = 0
private

Definition at line 65 of file L1TrackerEtMissEmulatorProducer.cc.

◆ GTTinput_

bool L1TrackerEtMissEmulatorProducer::GTTinput_ = false
private

Definition at line 71 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().

◆ L1MetCollectionName_

std::string L1TrackerEtMissEmulatorProducer::L1MetCollectionName_
private

Definition at line 75 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ phiQuadrants_

std::vector<l1tmetemu::global_phi_t> L1TrackerEtMissEmulatorProducer::phiQuadrants_
private

Definition at line 62 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ phiShifts_

std::vector<l1tmetemu::global_phi_t> L1TrackerEtMissEmulatorProducer::phiShifts_
private

Definition at line 63 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().

◆ pvToken_

const edm::EDGetTokenT<L1VertexCollectionType> L1TrackerEtMissEmulatorProducer::pvToken_
private

Definition at line 77 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

◆ trackToken_

const edm::EDGetTokenT<L1TTTrackRefCollectionType> L1TrackerEtMissEmulatorProducer::trackToken_
private

Definition at line 78 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

◆ TrackTransform

L1TkEtMissEmuTrackTransform L1TrackerEtMissEmulatorProducer::TrackTransform
private

Definition at line 73 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

◆ vtxAssocTrackToken_

const edm::EDGetTokenT<L1TTTrackRefCollectionType> L1TrackerEtMissEmulatorProducer::vtxAssocTrackToken_
private

Definition at line 79 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.14