L1TPhase2GMTBarrelStubProcessor Class Reference

#include <L1TPhase2GMTBarrelStubProcessor.h>

Public Member Functions
	L1TPhase2GMTBarrelStubProcessor ()
	L1TPhase2GMTBarrelStubProcessor (const edm::ParameterSet &)
l1t::MuonStubCollection	makeStubs (const L1Phase2MuDTPhContainer *, const L1MuDTChambThContainer *)
	~L1TPhase2GMTBarrelStubProcessor ()

Private Member Functions
l1t::MuonStub	buildStub (const L1Phase2MuDTPhDigi &, const L1MuDTChambThDigi *)
l1t::MuonStub	buildStubNoEta (const L1Phase2MuDTPhDigi &)
int	calculateEta (uint, int, uint, uint)

Private Attributes
std::vector< int >	coarseEta1_
std::vector< int >	coarseEta2_
std::vector< int >	coarseEta3_
std::vector< int >	coarseEta4_
std::vector< int >	eta1_
std::vector< int >	eta2_
std::vector< int >	eta3_
double	etaLSB_
int	maxBX_
int	minBX_
int	minPhiQuality_
int	phiBFactor_
double	phiLSB_
std::vector< int >	phiOffset_
int	verbose_

Detailed Description

Definition at line 13 of file L1TPhase2GMTBarrelStubProcessor.h.

Constructor & Destructor Documentation

L1TPhase2GMTBarrelStubProcessor() [1/2]

L1TPhase2GMTBarrelStubProcessor::L1TPhase2GMTBarrelStubProcessor

(

)

Definition at line 9 of file L1TPhase2GMTBarrelStubProcessor.cc.

: minPhiQuality_(0), minBX_(-3), maxBX_(3) {}

L1TPhase2GMTBarrelStubProcessor() [2/2]

L1TPhase2GMTBarrelStubProcessor::L1TPhase2GMTBarrelStubProcessor

(

const edm::ParameterSet &

iConfig

)

Definition at line 11 of file L1TPhase2GMTBarrelStubProcessor.cc.

    : minPhiQuality_(iConfig.getParameter<int>("minPhiQuality")),
      minBX_(iConfig.getParameter<int>("minBX")),
      maxBX_(iConfig.getParameter<int>("maxBX")),
      eta1_(iConfig.getParameter<std::vector<int> >("eta_1")),
      eta2_(iConfig.getParameter<std::vector<int> >("eta_2")),
      eta3_(iConfig.getParameter<std::vector<int> >("eta_3")),
      coarseEta1_(iConfig.getParameter<std::vector<int> >("coarseEta_1")),
      coarseEta2_(iConfig.getParameter<std::vector<int> >("coarseEta_2")),
      coarseEta3_(iConfig.getParameter<std::vector<int> >("coarseEta_3")),
      coarseEta4_(iConfig.getParameter<std::vector<int> >("coarseEta_4")),
      phiOffset_(iConfig.getParameter<std::vector<int> >("phiOffset")),
      phiBFactor_(iConfig.getParameter<int>("phiBDivider")),
      verbose_(iConfig.getParameter<int>("verbose")),
      phiLSB_(iConfig.getParameter<double>("phiLSB")),
      etaLSB_(iConfig.getParameter<double>("etaLSB")) {}

~L1TPhase2GMTBarrelStubProcessor()

L1TPhase2GMTBarrelStubProcessor::~L1TPhase2GMTBarrelStubProcessor

(

)

Definition at line 28 of file L1TPhase2GMTBarrelStubProcessor.cc.

{}

Member Function Documentation

buildStub()

l1t::MuonStub L1TPhase2GMTBarrelStubProcessor::buildStub ( const L1Phase2MuDTPhDigi &  phiS, const L1MuDTChambThDigi *  etaS  )

private

Definition at line 30 of file L1TPhase2GMTBarrelStubProcessor.cc.

References buildStubNoEta(), calculateEta(), HLT_2024v13_cff::eta1, HLT_2024v13_cff::eta2, etaLSB_, mps_fire::i, l1t::MuonStub::offline_coord1(), l1t::MuonStub::offline_coord2(), L1MuDTChambThDigi::position(), L1MuDTChambThDigi::quality(), L1MuDTChambThDigi::scNum(), l1t::MuonStub::setEta(), l1t::MuonStub::setOfflineQuantities(), L1MuDTChambThDigi::stNum(), parallelization::uint, and L1MuDTChambThDigi::whNum().

Referenced by makeStubs().

                                                                                        {
  l1t::MuonStub stub = buildStubNoEta(phiS);

  //Now full eta
  int qeta1 = -16384;
  int qeta2 = -16384;
  int eta1 = -16384;
  int eta2 = -16384;

  bool hasEta = false;
  for (uint i = 0; i < 7; ++i) {
    if (etaS->position(i) == 0)
      continue;
    if (!hasEta) {
      hasEta = true;
      eta1 = calculateEta(i, etaS->whNum(), etaS->scNum(), etaS->stNum());

      if (etaS->quality(i) == 1)
        qeta1 = 2;
      else
        qeta1 = 1;
    } else {
      eta2 = calculateEta(i, etaS->whNum(), etaS->scNum(), etaS->stNum());
      if (etaS->quality(i) == 1)
        qeta2 = 2;
      else
        qeta2 = 1;
    }
  }

  if (qeta2 > 0) {  //both stubs->average
    stub.setEta(eta1, eta2, 3);
    stub.setOfflineQuantities(stub.offline_coord1(), stub.offline_coord2(), eta1 * etaLSB_, eta2 * etaLSB_);

  } else if (qeta1 > 0) {  //Good single stub
    stub.setEta(eta1, 0, 1);
    stub.setOfflineQuantities(stub.offline_coord1(), stub.offline_coord2(), eta1 * etaLSB_, 0.0);
  }

  return stub;
}

buildStubNoEta()

l1t::MuonStub L1TPhase2GMTBarrelStubProcessor::buildStubNoEta ( const L1Phase2MuDTPhDigi &  phiS )

private

Definition at line 73 of file L1TPhase2GMTBarrelStubProcessor.cc.

References nano_mu_digi_cff::bx, L1Phase2MuDTPhDigi::bxNum(), coarseEta1_, coarseEta2_, coarseEta3_, coarseEta4_, PVValHelper::eta, etaLSB_, L1Phase2MuDTPhDigi::index(), createfilelist::int, M_PI, phi, L1Phase2MuDTPhDigi::phi(), L1Phase2MuDTPhDigi::phiBend(), L1Phase2MuDTPhDigi::quality(), quality, L1Phase2MuDTPhDigi::scNum(), nano_mu_digi_cff::sector, l1t::MuonStub::setOfflineQuantities(), Validation_hcalonly_cfi::sign, relativeConstraints::station, L1Phase2MuDTPhDigi::stNum(), makeGlobalPositionRcd_cfg::tag, parallelization::uint, makeMuonMisalignmentScenario::wheel, and L1Phase2MuDTPhDigi::whNum().

Referenced by buildStub(), and makeStubs().

                                                                                          {
  int wheel = phiS.whNum();
  int abswheel = fabs(phiS.whNum());
  int sign = wheel > 0 ? 1 : -1;
  int sector = phiS.scNum();
  int station = phiS.stNum();

  ap_uint<18> normalization0 = sector * ap_uint<15>(21845);
  ap_int<18> normalization1 = ap_int<18>(ap_int<17>(phiS.phi()) * ap_ufixed<8, 0>(0.3183));
  ap_int<18> kmtf_phi = ap_int<18>(normalization0 + normalization1);
  int phi = int(kmtf_phi);
  float globalPhi = phi * M_PI / (1 << 17);

  //  double globalPhi = (sector * 30) + phiS.phi() * 30. / 65535.;
  int tag = phiS.index();

  int bx = phiS.bxNum() - 20;
  int quality = phiS.quality();
  uint tfLayer = phiS.stNum() - 1;
  int eta = -16384;
  if (station == 1) {
    eta = coarseEta1_[abswheel];
  } else if (station == 2) {
    eta = coarseEta2_[abswheel];
  } else if (station == 3) {
    eta = coarseEta3_[abswheel];
  } else if (station == 4) {
    eta = coarseEta4_[abswheel];
  }

  //override!!!
  //  eta=abswheel;

  //Now full eta

  eta = eta * sign;
  l1t::MuonStub stub(wheel, sector, station, tfLayer, phi, phiS.phiBend(), tag, bx, quality, eta, 0, 0, 1);

  stub.setOfflineQuantities(globalPhi, float(phiS.phiBend() * 0.49e-3), eta * etaLSB_, 0.0);
  return stub;
}

calculateEta()

int L1TPhase2GMTBarrelStubProcessor::calculateEta ( uint  i, int  wheel, uint  sector, uint  station  )

private

Definition at line 182 of file L1TPhase2GMTBarrelStubProcessor.cc.

References PVValHelper::eta, eta1_, eta2_, eta3_, mps_fire::i, nano_mu_digi_cff::sector, relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

Referenced by buildStub().

                                                                                              {
  int eta = 0;
  if (wheel > 0) {
    eta = 7 * wheel + 3 - i;
  } else if (wheel < 0) {
    eta = 7 * wheel + i - 3;
  } else {
    if (sector == 0 || sector == 3 || sector == 4 || sector == 7 || sector == 8 || sector == 11)
      eta = i - 3;
    else
      eta = 3 - i;
  }

  if (station == 1)
    eta = eta1_[eta + 17];
  else if (station == 2)
    eta = eta2_[eta + 17];
  else
    eta = eta3_[eta + 17];

  return eta;
}

makeStubs()

l1t::MuonStubCollection L1TPhase2GMTBarrelStubProcessor::makeStubs	(	const L1Phase2MuDTPhContainer *	phiContainer,
		const L1MuDTChambThContainer *	etaContainer
	)

Definition at line 115 of file L1TPhase2GMTBarrelStubProcessor.cc.

References buildStub(), buildStubNoEta(), nano_mu_digi_cff::bx, L1MuDTChambThContainer::chThetaSegm(), TauDecayModes::dec, L1Phase2MuDTPhContainer::getContainer(), maxBX_, minBX_, minPhiQuality_, l1tGTMenu_BTagSeeds_cff::os, MillePedeFileConverter_cfg::out, nano_mu_digi_cff::sector, relativeConstraints::station, verbose_, and makeMuonMisalignmentScenario::wheel.

Referenced by Phase2L1TGMTStubProducer::produce().

                                                                                                               {
  l1t::MuonStubCollection out;
  for (int bx = minBX_; bx <= maxBX_; bx++) {
    ostringstream os;
    if (verbose_ == 2)
      os << "PATTERN ";
    for (int wheel = -2; wheel <= 2; wheel++) {
      for (int sector = 0; sector < 12; sector++) {
        for (int station = 1; station < 5; station++) {
          bool hasEta = false;
          const L1MuDTChambThDigi* tseta = etaContainer->chThetaSegm(wheel, station, sector, bx);
          if (tseta != nullptr) {
            hasEta = true;
          }

          for (const auto& phiDigi : *phiContainer->getContainer()) {
            if ((phiDigi.bxNum() - 20) != bx || phiDigi.whNum() != wheel || phiDigi.scNum() != sector ||
                phiDigi.stNum() != station)
              continue;
            if (phiDigi.quality() < minPhiQuality_)
              continue;

            if (verbose_ == 2) {
              ap_uint<64> wphi = ap_uint<17>(phiDigi.phi());
              ap_uint<64> wphib = ap_uint<13>(phiDigi.phiBend());
              ap_uint<64> wr1 = ap_uint<21>(0);
              ap_uint<64> wq = ap_uint<4>(phiDigi.quality());
              ap_uint<64> wr2 = ap_uint<9>(0);
              ap_uint<64> sN = 0;
              sN = sN | wphi;
              sN = sN | (wphib << 17);
              sN = sN | (wr1 << 30);
              sN = sN | (wq << 51);
              sN = sN | (wr2 << 55);
              os << std::setw(0) << std::dec << sector << " " << wheel << " " << station << " ";
              os << std::uppercase << std::setfill('0') << std::setw(16) << std::hex << uint64_t(sN) << " ";
            }

            if (hasEta) {
              out.push_back(buildStub(phiDigi, tseta));
            } else {
              out.push_back(buildStubNoEta(phiDigi));
            }
          }
        }
      }
    }
    if (verbose_ == 2)
      edm::LogInfo("BarrelStub") << os.str() << std::endl;
  }

  if (verbose_) {
    edm::LogInfo("BarrelStub") << "Barrel Stubs";
    for (const auto& stub : out)
      edm::LogInfo("BarrelStub") << "Barrel Stub bx=" << stub.bxNum() << " TF=" << stub.tfLayer()
                                 << " etaRegion=" << stub.etaRegion() << " phiRegion=" << stub.phiRegion()
                                 << " depthRegion=" << stub.depthRegion() << "  coord1=" << stub.offline_coord1() << ","
                                 << stub.coord1() << " coord2=" << stub.offline_coord2() << "," << stub.coord2()
                                 << " eta1=" << stub.offline_eta1() << "," << stub.eta1()
                                 << " eta2=" << stub.offline_eta2() << "," << stub.eta2()
                                 << " quality=" << stub.quality() << " etaQuality=" << stub.etaQuality();
  }

  return out;
}

Member Data Documentation

coarseEta1_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::coarseEta1_

private

Definition at line 35 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by buildStubNoEta().

coarseEta2_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::coarseEta2_

private

Definition at line 36 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by buildStubNoEta().

coarseEta3_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::coarseEta3_

private

Definition at line 37 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by buildStubNoEta().

coarseEta4_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::coarseEta4_

private

Definition at line 38 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by buildStubNoEta().

eta1_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::eta1_

private

Definition at line 32 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by calculateEta().

eta2_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::eta2_

private

Definition at line 33 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by calculateEta().

eta3_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::eta3_

private

Definition at line 34 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by calculateEta().

etaLSB_

double L1TPhase2GMTBarrelStubProcessor::etaLSB_

private

Definition at line 43 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by buildStub(), and buildStubNoEta().

maxBX_

int L1TPhase2GMTBarrelStubProcessor::maxBX_

private

Definition at line 30 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by makeStubs().

minBX_

int L1TPhase2GMTBarrelStubProcessor::minBX_

private

Definition at line 29 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by makeStubs().

minPhiQuality_

int L1TPhase2GMTBarrelStubProcessor::minPhiQuality_

private

Definition at line 27 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by makeStubs().

phiBFactor_

int L1TPhase2GMTBarrelStubProcessor::phiBFactor_

private

Definition at line 40 of file L1TPhase2GMTBarrelStubProcessor.h.

phiLSB_

double L1TPhase2GMTBarrelStubProcessor::phiLSB_

private

Definition at line 42 of file L1TPhase2GMTBarrelStubProcessor.h.

phiOffset_

std::vector<int> L1TPhase2GMTBarrelStubProcessor::phiOffset_

private

Definition at line 39 of file L1TPhase2GMTBarrelStubProcessor.h.

verbose_

int L1TPhase2GMTBarrelStubProcessor::verbose_

private

Definition at line 41 of file L1TPhase2GMTBarrelStubProcessor.h.

Referenced by makeStubs().