L1TrackerEtMissEmulatorProducer Class Reference

#include <L1TrackerEtMissEmulatorProducer.cc>

Inheritance diagram for L1TrackerEtMissEmulatorProducer:

Public Types
typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType
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
TTTrack_TrackWord::bendChi2_t	bendChi2Max_
TTTrack_TrackWord::chi2rphi_t	chi2rphiMax_
TTTrack_TrackWord::chi2rz_t	chi2rzMax_
bool	cordicDebug_ = false
int	cordicSteps_
std::vector< l1tmetemu::global_phi_t >	cosLUT_
int	debug_
l1tmetemu::z_t	deltaZ0_ = 0
std::vector< l1tmetemu::z_t >	DeltaZLUT_
vector< double >	etaRegions_
std::vector< l1tmetemu::eta_t >	EtaRegionsLUT_
bool	GTTinput_ = false
std::string	L1MetCollectionName_
l1tmetemu::eta_t	maxEta_
l1tmetemu::z_t	maxZ0_
l1tmetemu::pt_t	minPt_
l1tmetemu::z_t	minZ0_
l1tmetemu::nstub_t	nStubsmin_
std::vector< l1tmetemu::global_phi_t >	phiQuadrants_
std::vector< l1tmetemu::global_phi_t >	phiShifts_
const edm::EDGetTokenT< L1VertexCollectionType >	pvToken_
const edm::EDGetTokenT< L1TTTrackCollectionType >	trackToken_
L1TkEtMissEmuTrackTransform	TrackTransform
vector< double >	z0Thresholds_

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.

L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TrackerEtMissEmulatorProducer::L1TTTrackType

Definition at line 43 of file L1TrackerEtMissEmulatorProducer.cc.

L1VertexCollectionType

typedef l1t::VertexWordCollection L1TrackerEtMissEmulatorProducer::L1VertexCollectionType

Definition at line 45 of file L1TrackerEtMissEmulatorProducer.cc.

L1VertexType

typedef l1t::VertexWord L1TrackerEtMissEmulatorProducer::L1VertexType

Definition at line 46 of file L1TrackerEtMissEmulatorProducer.cc.

Constructor & Destructor Documentation

L1TrackerEtMissEmulatorProducer()

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

explicit

Definition at line 95 of file L1TrackerEtMissEmulatorProducer.cc.

References TTTrack_TrackWord::bendChi2Bins, bendChi2Max_, electrons_cff::bool, TTTrack_TrackWord::chi2RPhiBins, chi2rphiMax_, TTTrack_TrackWord::chi2RZBins, chi2rzMax_, cordicDebug_, cordicSteps_, cosLUT_, debug_, DeltaZLUT_, etaRegions_, EtaRegionsLUT_, l1tmetemu::generateCosLUT(), l1tmetemu::generateDeltaZLUT(), l1tmetemu::generateEtaRegionLUT(), L1TkEtMissEmuTrackTransform::generateLUTs(), l1tmetemu::getBin(), edm::ParameterSet::getParameter(), L1TkEtMissEmuTrackTransform::getPhiQuad(), L1TkEtMissEmuTrackTransform::getPhiShift(), GTTinput_, createfilelist::int, l1tmetemu::kInternalEtaWidth, l1tmetemu::kInternalPtWidth, l1tmetemu::kInternalVTXWidth, l1tmetemu::kMaxCosLUTPhi, l1tmetemu::kStepEta, l1tmetemu::kStepPhi, l1tmetemu::kStepPt, l1tmetemu::kStepZ0, L1MetCollectionName_, maxEta_, maxZ0_, minPt_, minZ0_, nStubsmin_, phiQuadrants_, phiShifts_, L1TkEtMissEmuTrackTransform::setGTTinput(), AlCaHLTBitMon_QueryRunRegistry::string, TrackTransform, and z0Thresholds_.

    : pvToken_(consumes<L1VertexCollectionType>(iConfig.getParameter<edm::InputTag>("L1VertexInputTag"))),
      trackToken_(consumes<L1TTTrackCollectionType>(iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))) {
  // Setup LUTs
  TrackTransform.generateLUTs();
  phiQuadrants_ = TrackTransform.getPhiQuad();
  phiShifts_ = TrackTransform.getPhiShift();

  // Get Emulator config parameters
  cordicSteps_ = (int)iConfig.getParameter<int>("nCordicSteps");
  debug_ = (int)iConfig.getParameter<int>("debug");

  GTTinput_ = (bool)iConfig.getParameter<bool>("useGTTinput");

  TrackTransform.setGTTinput(GTTinput_);

  // Name of output ED Product
  L1MetCollectionName_ = (std::string)iConfig.getParameter<std::string>("L1MetCollectionName");

  // Input parameter cuts and convert to correct integer representations
  maxZ0_ = l1tmetemu::digitizeSignedValue<l1tmetemu::z_t>(
      (double)iConfig.getParameter<double>("maxZ0"), l1tmetemu::kInternalVTXWidth, l1tmetemu::kStepZ0);
  minZ0_ = (1 << TTTrack_TrackWord::TrackBitWidths::kZ0Size) - maxZ0_;

  maxEta_ = l1tmetemu::digitizeSignedValue<l1tmetemu::eta_t>(
      (double)iConfig.getParameter<double>("maxEta"), l1tmetemu::kInternalEtaWidth, l1tmetemu::kStepEta);

  chi2rphiMax_ =
      l1tmetemu::getBin((double)iConfig.getParameter<double>("chi2rphidofMax"), TTTrack_TrackWord::chi2RPhiBins);
  chi2rzMax_ = l1tmetemu::getBin((double)iConfig.getParameter<double>("chi2rzdofMax"), TTTrack_TrackWord::chi2RZBins);
  bendChi2Max_ =
      l1tmetemu::getBin((double)iConfig.getParameter<double>("bendChi2Max"), TTTrack_TrackWord::bendChi2Bins);

  minPt_ = l1tmetemu::digitizeSignedValue<l1tmetemu::pt_t>(
      (double)iConfig.getParameter<double>("minPt"), l1tmetemu::kInternalPtWidth, l1tmetemu::kStepPt);

  nStubsmin_ = (l1tmetemu::nstub_t)iConfig.getParameter<int>("nStubsmin");

  z0Thresholds_ = iConfig.getParameter<std::vector<double>>("z0Thresholds");
  etaRegions_ = iConfig.getParameter<std::vector<double>>("etaRegions");

  if (debug_ == 5) {
    cordicDebug_ = true;
  }

  // To have same bin spacing between 0 and pi/2 as between original phi
  // granularity
  int cosLUTbins = floor(l1tmetemu::kMaxCosLUTPhi / l1tmetemu::kStepPhi);

  // Compute LUTs
  cosLUT_ = l1tmetemu::generateCosLUT(cosLUTbins);
  EtaRegionsLUT_ = l1tmetemu::generateEtaRegionLUT(etaRegions_);
  DeltaZLUT_ = l1tmetemu::generateDeltaZLUT(z0Thresholds_);

  produces<std::vector<EtSum>>(L1MetCollectionName_);
}

~L1TrackerEtMissEmulatorProducer()

L1TrackerEtMissEmulatorProducer::~L1TrackerEtMissEmulatorProducer ( )

override

Definition at line 152 of file L1TrackerEtMissEmulatorProducer.cc.

{}

Member Function Documentation

beginJob()

void L1TrackerEtMissEmulatorProducer::beginJob ( void  )

privatevirtual

Definition at line 427 of file L1TrackerEtMissEmulatorProducer.cc.

{}

endJob()

void L1TrackerEtMissEmulatorProducer::endJob ( void  )

privatevirtual

Definition at line 429 of file L1TrackerEtMissEmulatorProducer.cc.

{}

fillDescriptions()

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

static

Definition at line 432 of file L1TrackerEtMissEmulatorProducer.cc.

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

                                                                                                 {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

produce()

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

overrideprivate

Definition at line 154 of file L1TrackerEtMissEmulatorProducer.cc.

References funct::abs(), InternalEtWord::bendChi2, bendChi2Max_, InternalEtWord::chi2rphidof, chi2rphiMax_, InternalEtWord::chi2rzdof, chi2rzMax_, cordicDebug_, l1tmhtemu::cordicSqrt(), cordicSteps_, funct::cos(), cosLUT_, debug_, deltaZ0_, DeltaZLUT_, l1tmetemu::EtMiss::Et, InternalEtWord::eta, EtaRegionsLUT_, InternalEtWord::EtaSector, dqmMemoryStats::float, InternalEtWord::globalPhi, mps_fire::i, iEvent, edm::HandleBase::isValid(), l1tmetemu::kGlobalPhiBins, l1tmetemu::kInternalVTXWidth, l1tmetemu::kMETPhiBins, l1tmetemu::kMETSize, L1Analysis::kMissingEt, l1tmetemu::kNEtaRegion, l1tmetemu::kNSector, l1tmetemu::kStepMET, l1tmetemu::kStepMETPhi, l1tmetemu::kStepPt, L1MetCollectionName_, M_PI, maxEta_, maxZ0_, minPt_, minZ0_, eostools::move(), InternalEtWord::nstubs, nStubsmin_, InternalEtWord::phi, l1tmetemu::EtMiss::Phi, phiQuadrants_, InternalEtWord::pt, InternalEtWord::pV, pvToken_, InternalEtWord::Sector, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, cond::impl::to_string(), HLT_2022v11_cff::track, trackToken_, TrackTransform, L1TkEtMissEmuTrackTransform::transformTrack(), l1tmetemu::unpackSignedValue(), extraflags_cff::vtx, and InternalEtWord::z0.

                                                                                           {
  using namespace edm;

  std::unique_ptr<std::vector<l1t::EtSum>> METCollection(new std::vector<l1t::EtSum>(0));

  edm::Handle<L1VertexCollectionType> L1VertexHandle;
  iEvent.getByToken(pvToken_, L1VertexHandle);

  edm::Handle<L1TTTrackCollectionType> L1TTTrackHandle;
  iEvent.getByToken(trackToken_, L1TTTrackHandle);

  // Initialize cordic class
  Cordic cordicSqrt(l1tmetemu::kMETPhiBins, l1tmetemu::kMETSize, cordicSteps_, cordicDebug_);

  if (!L1VertexHandle.isValid()) {
    LogError("L1TrackerEtMissEmulatorProducer") << "\nWarning: VertexCollection not found in the event. Exit\n";
    return;
  }

  if (!L1TTTrackHandle.isValid()) {
    LogError("L1TrackerEtMissEmulatorProducer") << "\nWarning: L1TTTrackCollection not found in the event. Exit\n";
    return;
  }

  // Initialize sector sums, need 0 initialization in case a sector has no
  // tracks
  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};
  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};
  int sector_totals[l1tmetemu::kNSector] = {0, 0, 0, 0, 0, 0, 0, 0, 0};

  // Track counters
  int num_tracks{0};
  int num_assoc_tracks{0};
  int num_quality_tracks{0};

  // Get reference to first vertex in event vertex collection
  L1VertexType& vtx = const_cast<L1VertexType&>(L1VertexHandle->at(0));

  for (const auto& track : *L1TTTrackHandle) {
    num_tracks++;
    L1TTTrackType& track_ref = const_cast<L1TTTrackType&>(track);  // Get Reference to track to pass to TrackTransform

    // Convert to internal track representation
    InternalEtWord EtTrack = TrackTransform.transformTrack<L1TTTrackType, L1VertexType>(track_ref, vtx);

    // Parameter cuts
    if (EtTrack.pt < minPt_)
      continue;

    // Z signed so double bound
    if (EtTrack.z0 & (1 << (l1tmetemu::kInternalVTXWidth - 1))) {
      // if negative
      if (EtTrack.z0 <= maxZ0_)
        continue;
    } else {
      if (EtTrack.z0 > minZ0_)
        continue;
    }
    if (EtTrack.eta > maxEta_)
      continue;

    // Quality Cuts
    if (EtTrack.chi2rphidof >= chi2rphiMax_)
      continue;

    if (EtTrack.chi2rzdof >= chi2rzMax_)
      continue;

    if (EtTrack.bendChi2 >= bendChi2Max_)
      continue;

    if (EtTrack.nstubs < nStubsmin_)
      continue;

    num_quality_tracks++;

    // Temporary int representation to get the difference
    int tempz = l1tmetemu::unpackSignedValue(EtTrack.z0, l1tmetemu::kInternalVTXWidth);
    int temppv = l1tmetemu::unpackSignedValue(EtTrack.pV, l1tmetemu::kInternalVTXWidth);

    l1tmetemu::z_t z_diff = abs(tempz - temppv);

    // Track to vertex association, adaptive z window based on eta region
    for (unsigned int reg = 0; reg < l1tmetemu::kNEtaRegion; reg++) {
      if (EtTrack.eta >= EtaRegionsLUT_[reg] && EtTrack.eta < EtaRegionsLUT_[reg + 1]) {
        deltaZ0_ = DeltaZLUT_[reg];
        break;
      }
    }

    if (z_diff <= deltaZ0_) {
      num_assoc_tracks++;
      if (debug_ == 7) {
        edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
            << "Track to Vertex ID: " << num_tracks << "\n"
            << "Phi Sector: " << EtTrack.Sector << " pT: " << EtTrack.pt << " Phi: " << EtTrack.globalPhi
            << " TanL: " << EtTrack.eta << " Z0: " << EtTrack.z0 << " Nstub: " << EtTrack.nstubs
            << " Chi2rphi: " << EtTrack.chi2rphidof << " Chi2rz: " << EtTrack.chi2rzdof
            << " bendChi2: " << EtTrack.bendChi2 << " PV: " << EtTrack.pV << "\n"
            << "--------------------------------------------------------------\n";
      }

      if (debug_ == 2) {
        edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
            << "========================Phi debug=================================\n"
            << "Int pT: " << EtTrack.pt << "\n"
            << "Int Phi: " << EtTrack.globalPhi << " Float Phi: " << EtTrack.phi
            << " Actual Float Cos(Phi): " << cos(EtTrack.phi) << " Actual Float Sin(Phi): " << sin(EtTrack.phi) << "\n";
      }
      l1tmetemu::Et_t temppx = 0;
      l1tmetemu::Et_t temppy = 0;

      // Split tracks in phi quadrants and access cosLUT_, backwards iteration
      // through cosLUT_ gives sin Sum sector Et -ve when cos or sin phi are -ve
      sector_totals[EtTrack.Sector] += 1;
      if (EtTrack.globalPhi >= phiQuadrants_[0] && EtTrack.globalPhi < phiQuadrants_[1]) {
        temppx = (EtTrack.pt * cosLUT_[EtTrack.globalPhi]);
        temppy = (EtTrack.pt * cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi]);

        if (debug_ == 2) {
          edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
              << "Sector: " << EtTrack.Sector << " Quadrant: " << 1 << "\n"
              << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): " << cosLUT_[EtTrack.globalPhi]
              << " Int Sin(Phi): " << cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi] << "\n"

              << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins
              << " Float Cos(Phi): " << (float)cosLUT_[EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
              << " Float Sin(Phi): "
              << (float)cosLUT_[phiQuadrants_[1] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins << "\n";
        }
      } else if (EtTrack.globalPhi >= phiQuadrants_[1] && EtTrack.globalPhi < phiQuadrants_[2]) {
        temppx = -(EtTrack.pt * cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi]);
        temppy = (EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]]);

        if (debug_ == 2) {
          edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
              << "Sector: " << EtTrack.Sector << " Quadrant: " << 2 << "\n"
              << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): -"
              << cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi]
              << " Int Sin(Phi): " << cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]] << "\n"

              << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): -"
              << (float)cosLUT_[phiQuadrants_[2] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
              << " Float Sin(Phi): " << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[1]] / l1tmetemu::kGlobalPhiBins
              << "\n";
        }
      } else if (EtTrack.globalPhi >= phiQuadrants_[2] && EtTrack.globalPhi < phiQuadrants_[3]) {
        temppx = -(EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]]);
        temppy = -(EtTrack.pt * cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi]);

        if (debug_ == 2) {
          edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
              << "Sector: " << EtTrack.Sector << " Quadrant: " << 3 << "\n"
              << "Int Phi: " << EtTrack.globalPhi << " Int Cos(Phi): -" << cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]]
              << " Int Sin(Phi): -" << cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi] << "\n"

              << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): -"
              << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[2]] / l1tmetemu::kGlobalPhiBins
              << " Float Sin(Phi): -"
              << (float)cosLUT_[phiQuadrants_[3] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins << "\n";
        }

      } else if (EtTrack.globalPhi >= phiQuadrants_[3] && EtTrack.globalPhi < phiQuadrants_[4]) {
        temppx = (EtTrack.pt * cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi]);
        temppy = -(EtTrack.pt * cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]]);

        if (debug_ == 2) {
          edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
              << "Sector: " << EtTrack.Sector << " Quadrant: " << 4 << "\n"
              << "Int Phi: " << EtTrack.globalPhi
              << " Int Cos(Phi): " << cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi] << " Int Sin(Phi): -"
              << cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]] << "\n"

              << "Float Phi: " << (float)EtTrack.globalPhi / l1tmetemu::kGlobalPhiBins << " Float Cos(Phi): "
              << (float)cosLUT_[phiQuadrants_[4] - 1 - EtTrack.globalPhi] / l1tmetemu::kGlobalPhiBins
              << " Float Sin(Phi): -"
              << (float)cosLUT_[EtTrack.globalPhi - phiQuadrants_[3]] / l1tmetemu::kGlobalPhiBins << "\n";
        }
      } else {
        temppx = 0;
        temppy = 0;
      }
      if (EtTrack.EtaSector) {
        sumPx[EtTrack.Sector] = sumPx[EtTrack.Sector] + temppx;
        sumPy[EtTrack.Sector] = sumPy[EtTrack.Sector] + temppy;
      } else {
        sumPx[l1tmetemu::kNSector + EtTrack.Sector] = sumPx[l1tmetemu::kNSector + EtTrack.Sector] + temppx;
        sumPy[l1tmetemu::kNSector + EtTrack.Sector] = sumPy[l1tmetemu::kNSector + EtTrack.Sector] + temppy;
      }
    }

  }  // end loop over tracks

  l1tmetemu::Et_t GlobalPx = 0;
  l1tmetemu::Et_t GlobalPy = 0;

  float tempsumPx = 0;
  float tempsumPy = 0;

  // Global Et sum as floats to emulate rounding in HW
  for (unsigned int i = 0; i < l1tmetemu::kNSector * 2; i++) {
    tempsumPx += floor((float)sumPx[i] / (float)l1tmetemu::kGlobalPhiBins);
    tempsumPy += floor((float)sumPy[i] / (float)l1tmetemu::kGlobalPhiBins);
  }

  // Recast rounded temporary sums into Et_t datatype
  GlobalPx = tempsumPx;
  GlobalPy = tempsumPy;

  // Perform cordic sqrt, take x,y and converts to polar coordinate r,phi where
  // r=sqrt(x**2+y**2) and phi = atan(y/x)
  l1tmetemu::EtMiss EtMiss = cordicSqrt.toPolar(GlobalPx, GlobalPy);

  // Recentre phi
  l1tmetemu::METphi_t tempPhi = 0;

  if ((GlobalPx < 0) && (GlobalPy < 0))
    tempPhi = EtMiss.Phi - l1tmetemu::kMETPhiBins / 2;
  else if ((GlobalPx >= 0) && (GlobalPy >= 0))
    tempPhi = (EtMiss.Phi) + l1tmetemu::kMETPhiBins / 2;
  else if ((GlobalPx >= 0) && (GlobalPy < 0))
    tempPhi = EtMiss.Phi - l1tmetemu::kMETPhiBins / 2;
  else if ((GlobalPx < 0) && (GlobalPy >= 0))
    tempPhi = EtMiss.Phi - 3 * l1tmetemu::kMETPhiBins / 2;

  if (debug_ == 6) {
    std::string flpxarray[l1tmetemu::kNSector * 2];
    std::string flpyarray[l1tmetemu::kNSector * 2];

    std::string intpxarray[l1tmetemu::kNSector * 2];
    std::string intpyarray[l1tmetemu::kNSector * 2];

    std::string totalsarray[l1tmetemu::kNSector * 2];

    for (unsigned int i = 0; i < l1tmetemu::kNSector * 2; i++) {
      flpxarray[i] = to_string(sumPx[i] * l1tmetemu::kStepPt) + "|";
      flpyarray[i] = to_string(sumPy[i] * l1tmetemu::kStepPt) + "|";
      intpxarray[i] = to_string(floor((float)sumPx[i] / (float)l1tmetemu::kGlobalPhiBins)) + "|";
      intpyarray[i] = to_string(floor((float)sumPy[i] / (float)l1tmetemu::kGlobalPhiBins)) + "|";
      totalsarray[i] = to_string(sector_totals[i]) + "|";
    }

    edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
        << "====Sector Pt====\n"
        << "Float Px: " << flpxarray << "Float Py: " << flpyarray << "Integer Px: " << intpyarray
        << "Integer Px: " << intpyarray << "Sector Totals: " << totalsarray

        << "====Global Pt====\n"
        << "Integer Global Px: " << GlobalPx << "| Integer Global Py: " << GlobalPy << "\n"
        << "Float Global Px: " << GlobalPx * l1tmetemu::kStepPt
        << "| Float Global Py: " << GlobalPy * l1tmetemu::kStepPt << "\n";
  }

  if (debug_ == 6) {
    edm::LogVerbatim("L1TrackerEtMissEmulatorProducer")
        << "====MET===\n"
        << "Integer MET: " << EtMiss.Et << "| Integer MET phi: " << EtMiss.Phi << "\n"
        << "Float MET: " << (EtMiss.Et) * l1tmetemu::kStepMET
        << "| Float MET phi: " << (float)tempPhi * l1tmetemu::kStepMETPhi - M_PI << "\n"
        << "# Intial Tracks: " << num_tracks << "\n"
        << "# Tracks after Quality Cuts: " << num_quality_tracks << "\n"
        << "# Tracks Associated to Vertex: " << num_assoc_tracks << "\n"
        << "========================================================\n";
  }

  math::XYZTLorentzVector missingEt(-GlobalPx, -GlobalPy, 0, EtMiss.Et);
  EtSum L1EtSum(missingEt, EtSum::EtSumType::kMissingEt, (int)EtMiss.Et, 0, (int)tempPhi, (int)num_assoc_tracks);

  METCollection->push_back(L1EtSum);

  iEvent.put(std::move(METCollection), L1MetCollectionName_);
}  // end producer

Member Data Documentation

bendChi2Max_

TTTrack_TrackWord::bendChi2_t L1TrackerEtMissEmulatorProducer::bendChi2Max_

private

Definition at line 71 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

chi2rphiMax_

TTTrack_TrackWord::chi2rphi_t L1TrackerEtMissEmulatorProducer::chi2rphiMax_

private

Definition at line 69 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

chi2rzMax_

TTTrack_TrackWord::chi2rz_t L1TrackerEtMissEmulatorProducer::chi2rzMax_

private

Definition at line 70 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

cordicDebug_

bool L1TrackerEtMissEmulatorProducer::cordicDebug_ = false

private

Definition at line 82 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

cordicSteps_

int L1TrackerEtMissEmulatorProducer::cordicSteps_

private

Definition at line 80 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

cosLUT_

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

private

Definition at line 60 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

debug_

int L1TrackerEtMissEmulatorProducer::debug_

private

Definition at line 81 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

deltaZ0_

l1tmetemu::z_t L1TrackerEtMissEmulatorProducer::deltaZ0_ = 0

private

Definition at line 78 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

DeltaZLUT_

std::vector<l1tmetemu::z_t> L1TrackerEtMissEmulatorProducer::DeltaZLUT_

private

Definition at line 62 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

etaRegions_

vector<double> L1TrackerEtMissEmulatorProducer::etaRegions_

private

Definition at line 76 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().

EtaRegionsLUT_

std::vector<l1tmetemu::eta_t> L1TrackerEtMissEmulatorProducer::EtaRegionsLUT_

private

Definition at line 61 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

GTTinput_

bool L1TrackerEtMissEmulatorProducer::GTTinput_ = false

private

Definition at line 84 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().

L1MetCollectionName_

std::string L1TrackerEtMissEmulatorProducer::L1MetCollectionName_

private

Definition at line 88 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

maxEta_

l1tmetemu::eta_t L1TrackerEtMissEmulatorProducer::maxEta_

private

Definition at line 68 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

maxZ0_

l1tmetemu::z_t L1TrackerEtMissEmulatorProducer::maxZ0_

private

Definition at line 67 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

minPt_

l1tmetemu::pt_t L1TrackerEtMissEmulatorProducer::minPt_

private

Definition at line 72 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

minZ0_

l1tmetemu::z_t L1TrackerEtMissEmulatorProducer::minZ0_

private

Definition at line 66 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

nStubsmin_

l1tmetemu::nstub_t L1TrackerEtMissEmulatorProducer::nStubsmin_

private

Definition at line 73 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

phiQuadrants_

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

private

Definition at line 63 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

phiShifts_

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

private

Definition at line 64 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().

pvToken_

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

private

Definition at line 90 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

trackToken_

const edm::EDGetTokenT<L1TTTrackCollectionType> L1TrackerEtMissEmulatorProducer::trackToken_

private

Definition at line 91 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by produce().

TrackTransform

L1TkEtMissEmuTrackTransform L1TrackerEtMissEmulatorProducer::TrackTransform

private

Definition at line 86 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer(), and produce().

z0Thresholds_

vector<double> L1TrackerEtMissEmulatorProducer::z0Thresholds_

private

Definition at line 75 of file L1TrackerEtMissEmulatorProducer.cc.

Referenced by L1TrackerEtMissEmulatorProducer().