L1TPhase2GMTEndcapStubProcessor Class Reference

#include <L1TPhase2GMTEndcapStubProcessor.h>

Public Member Functions
	L1TPhase2GMTEndcapStubProcessor ()
	L1TPhase2GMTEndcapStubProcessor (const edm::ParameterSet &)
l1t::MuonStubCollection	makeStubs (const MuonDigiCollection< CSCDetId, CSCCorrelatedLCTDigi > &csc, const MuonDigiCollection< RPCDetId, RPCDigi > &rpc, const L1TMuon::GeometryTranslator *t, const edm::EventSetup &iSetup)
	~L1TPhase2GMTEndcapStubProcessor ()

Private Member Functions
l1t::MuonStub	buildCSCOnlyStub (const CSCDetId &, const CSCCorrelatedLCTDigi &, const L1TMuon::GeometryTranslator *)
l1t::MuonStub	buildRPCOnlyStub (const RPCDetId &, const RPCDigi &, const L1TMuon::GeometryTranslator *)
l1t::MuonStubCollection	combineStubs (const l1t::MuonStubCollection &, const l1t::MuonStubCollection &)

Private Attributes
double	coord1LSB_
double	coord2LSB_
double	eta1LSB_
double	eta2LSB_
double	etaMatch_
int	maxBX_
int	minBX_
double	phiMatch_
bool	verbose_

Detailed Description

Definition at line 16 of file L1TPhase2GMTEndcapStubProcessor.h.

Constructor & Destructor Documentation

L1TPhase2GMTEndcapStubProcessor() [1/2]

L1TPhase2GMTEndcapStubProcessor::L1TPhase2GMTEndcapStubProcessor

(

)

Definition at line 8 of file L1TPhase2GMTEndcapStubProcessor.cc.

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

L1TPhase2GMTEndcapStubProcessor() [2/2]

L1TPhase2GMTEndcapStubProcessor::L1TPhase2GMTEndcapStubProcessor

(

const edm::ParameterSet &

iConfig

)

Definition at line 10 of file L1TPhase2GMTEndcapStubProcessor.cc.

    : minBX_(iConfig.getParameter<int>("minBX")),
      maxBX_(iConfig.getParameter<int>("maxBX")),
      coord1LSB_(iConfig.getParameter<double>("coord1LSB")),
      coord2LSB_(iConfig.getParameter<double>("coord2LSB")),
      eta1LSB_(iConfig.getParameter<double>("eta1LSB")),
      eta2LSB_(iConfig.getParameter<double>("eta2LSB")),
      etaMatch_(iConfig.getParameter<double>("etaMatch")),
      phiMatch_(iConfig.getParameter<double>("phiMatch")),
      verbose_(iConfig.getParameter<unsigned int>("verbose")) {}

~L1TPhase2GMTEndcapStubProcessor()

L1TPhase2GMTEndcapStubProcessor::~L1TPhase2GMTEndcapStubProcessor

(

)

Definition at line 21 of file L1TPhase2GMTEndcapStubProcessor.cc.

{}

Member Function Documentation

buildCSCOnlyStub()

l1t::MuonStub L1TPhase2GMTEndcapStubProcessor::buildCSCOnlyStub ( const CSCDetId &  detid, const CSCCorrelatedLCTDigi &  digi, const L1TMuon::GeometryTranslator *  translator  )

private

Definition at line 23 of file L1TPhase2GMTEndcapStubProcessor.cc.

References simKBmtfDigis_cfi::bx, relativeConstraints::chamber, CSCDetId::chamber(), coord1LSB_, CSCDetId::endcap(), makeMuonMisalignmentScenario::endcap, HLT_2022v15_cff::eta1, eta1LSB_, CSCCorrelatedLCTDigi::getBX(), L1TMuon::GeometryTranslator::getGlobalPoint(), runTauDisplay::gp, createfilelist::int, phi, funct::primitive(), quality, CSCDetId::ring(), relativeConstraints::ring, hgcalTBTopologyTester_cfi::sector, l1t::MuonStub::setOfflineQuantities(), Validation_hcalonly_cfi::sign, relativeConstraints::station, CSCDetId::station(), parallelization::uint, and makeMuonMisalignmentScenario::wheel.

Referenced by makeStubs().

                                                                                                             {
  int endcap = detid.endcap();
  int station = detid.station();
  int chamber = detid.chamber();
  int ring = detid.ring();

  L1TMuon::TriggerPrimitive primitive(detid, digi);

  const GlobalPoint& gp = translator->getGlobalPoint(primitive);

  int phi = int(gp.phi().value() / coord1LSB_);
  int eta1 = int(gp.eta() / eta1LSB_);

  int wheel = 0;
  int sign = endcap == 1 ? -1 : 1;

  if (ring == 3)
    wheel = sign * 3;
  else if (ring == 2)
    wheel = sign * 4;
  else if (ring == 1)
    wheel = sign * 5;

  int sector = fabs(chamber);

  int bx = digi.getBX() - 8;
  int quality = 1;

  uint tfLayer = 0;
  if ((ring == 3 || ring == 2) && station == 1)  //ME1/3
    tfLayer = 4;
  else if (ring == 1 && station == 1)  //ME1/3
    tfLayer = 0;
  else if (station == 2)  //ME2/2
    tfLayer = 2;
  else if (station == 3)  //ME3/2
    tfLayer = 1;
  else if (station == 4)  //ME4/2
    tfLayer = 3;

  l1t::MuonStub stub(wheel, sector, station, tfLayer, phi, 0, 0, bx, quality, eta1, 0, 1, 0);

  stub.setOfflineQuantities(gp.phi().value(), 0.0, gp.eta(), 0.0);
  return stub;
}

buildRPCOnlyStub()

l1t::MuonStub L1TPhase2GMTEndcapStubProcessor::buildRPCOnlyStub ( const RPCDetId &  detid, const RPCDigi &  digi, const L1TMuon::GeometryTranslator *  translator  )

private

Definition at line 71 of file L1TPhase2GMTEndcapStubProcessor.cc.

References RPCDigi::bx(), simKBmtfDigis_cfi::bx, coord2LSB_, HLT_2022v15_cff::eta2, eta2LSB_, L1TMuon::GeometryTranslator::getGlobalPoint(), runTauDisplay::gp, createfilelist::int, funct::primitive(), quality, RPCDetId::region(), RPCDetId::ring(), relativeConstraints::ring, hgcalTBTopologyTester_cfi::sector, RPCDetId::sector(), l1t::MuonStub::setOfflineQuantities(), relativeConstraints::station, RPCDetId::station(), RPCDetId::subsector(), makeGlobalPositionRcd_cfg::tag, RPCDetId::trIndex(), parallelization::uint, and makeMuonMisalignmentScenario::wheel.

Referenced by makeStubs().

                                                                                                             {
  L1TMuon::TriggerPrimitive primitive(detid, digi);
  const GlobalPoint& gp = translator->getGlobalPoint(primitive);

  int phi2 = int(gp.phi().value() / coord2LSB_);
  int eta2 = int(gp.eta() / eta2LSB_);

  int wheel = -(6 - detid.ring()) * detid.region();
  int sector = (detid.sector() - 1) * 6 + detid.subsector();
  int station = detid.station();
  bool tag = detid.trIndex();
  int bx = digi.bx();
  int quality = 2;

  int ring = detid.ring();

  uint tfLayer = 0;

  if ((ring == 3 || ring == 2) && station == 1)  //ME1/3
    tfLayer = 4;
  //  else if (ring==1 && station==1) //ME1/3
  //    tfLayer=0;
  else if (station == 2)  //ME2/2
    tfLayer = 2;
  else if (station == 3)  //ME3/2
    tfLayer = 1;
  else if (station == 4)  //ME4/2
    tfLayer = 3;

  l1t::MuonStub stub(wheel, sector, station, tfLayer, 0, phi2, tag, bx, quality, 0, eta2, 2, 0);
  stub.setOfflineQuantities(gp.phi().value(), gp.phi().value(), gp.eta(), gp.eta());
  return stub;
}

combineStubs()

l1t::MuonStubCollection L1TPhase2GMTEndcapStubProcessor::combineStubs ( const l1t::MuonStubCollection &  cscStubs, const l1t::MuonStubCollection &  rpcStubs  )

private

Definition at line 107 of file L1TPhase2GMTEndcapStubProcessor.cc.

References SiPixelRawToDigiRegional_cfi::deltaPhi, l1t::MuonStub::depthRegion(), PVValHelper::eta, l1t::MuonStub::eta1(), etaMatch_, l1t::MuonStub::etaRegion(), mps_fire::i, runTheMatrix::keep, l1t::MuonStub::offline_coord1(), l1t::MuonStub::offline_eta1(), MillePedeFileConverter_cfg::out, phi, phiMatch_, omtf::DataWord64::rpc, l1t::MuonStub::setEta(), l1t::MuonStub::setOfflineQuantities(), and parallelization::uint.

Referenced by makeStubs().

                                                                                                             {
  l1t::MuonStubCollection out;
  l1t::MuonStubCollection usedRPC;
  l1t::MuonStubCollection cleanedRPC;

  l1t::MuonStubCollection cleanedCSC;

  //clean ME11 ambiguities
  l1t::MuonStubCollection allCSC = cscStubs;

  while (!allCSC.empty()) {
    l1t::MuonStub stub = allCSC[0];
    l1t::MuonStubCollection freeCSC;
    for (uint i = 1; i < allCSC.size(); ++i) {
      if ((stub.etaRegion() == allCSC[i].etaRegion()) && (stub.depthRegion() == allCSC[i].depthRegion()) &&
          (fabs(deltaPhi(stub.offline_coord1(), allCSC[i].offline_coord1())) < 0.001)) {
        if (fabs(stub.offline_eta1() - allCSC[i].offline_eta1()) > 0.001) {
          stub.setEta(stub.eta1(), allCSC[i].eta1(), 3);
          stub.setOfflineQuantities(stub.offline_coord1(), 0.0, stub.offline_eta1(), allCSC[i].offline_eta1());
        }
      } else {
        freeCSC.push_back(allCSC[i]);
      }
    }
    cleanedCSC.push_back(stub);
    allCSC = freeCSC;
  }

  for (const auto& csc : cleanedCSC) {
    int nRPC = 0;
    float phiF = 0.0;
    float etaF = 0.0;
    int phi = 0;
    int eta = 0;
    for (const auto& rpc : rpcStubs) {
      if (csc.depthRegion() != rpc.depthRegion())
        continue;
      if (fabs(deltaPhi(csc.offline_coord1(), rpc.offline_coord2())) < phiMatch_ &&
          fabs(csc.offline_eta1() - rpc.offline_eta2()) < etaMatch_ && csc.bxNum() == rpc.bxNum()) {
        phiF += rpc.offline_coord2();
        etaF += rpc.offline_eta2();
        phi += rpc.coord2();
        eta += rpc.eta2();
        nRPC++;
        usedRPC.push_back(rpc);
      }
    }

    int finalRPCPhi = 0;
    int finalRPCEta = 0;
    double offline_finalRPCPhi = 0;
    double offline_finalRPCEta = 0;
    if (nRPC != 0) {
      finalRPCPhi = phi / nRPC;
      finalRPCEta = eta / nRPC;
      offline_finalRPCPhi = phiF / nRPC;
      offline_finalRPCEta = etaF / nRPC;
      l1t::MuonStub stub(csc.etaRegion(),
                         csc.phiRegion(),
                         csc.depthRegion(),
                         csc.tfLayer(),
                         csc.coord1(),
                         finalRPCPhi,
                         0,
                         csc.bxNum(),
                         3,
                         csc.eta1(),
                         finalRPCEta,
                         3,
                         0);
      stub.setOfflineQuantities(csc.offline_coord1(), offline_finalRPCPhi, csc.offline_eta1(), offline_finalRPCEta);
      out.push_back(stub);
    } else {
      out.push_back(csc);
    }
  }

  //clean the RPC from the used ones

  for (const auto& rpc : rpcStubs) {
    bool keep = true;
    for (const auto& rpc2 : usedRPC) {
      if (rpc == rpc2) {
        keep = false;
        break;
      }
    }
    if (keep)
      cleanedRPC.push_back(rpc);
  }

  while (!cleanedRPC.empty()) {
    l1t::MuonStubCollection freeRPC;

    int nRPC = 1;
    float phiF = cleanedRPC[0].offline_coord2();
    float etaF = cleanedRPC[0].offline_eta2();
    int phi = cleanedRPC[0].coord2();
    int eta = cleanedRPC[0].eta2();

    for (unsigned i = 1; i < cleanedRPC.size(); ++i) {
      if (fabs(deltaPhi(cleanedRPC[0].offline_coord2(), cleanedRPC[i].offline_coord2())) < phiMatch_ &&
          cleanedRPC[0].depthRegion() == cleanedRPC[i].depthRegion() &&
          fabs(cleanedRPC[0].offline_eta2() - cleanedRPC[i].offline_eta2()) < etaMatch_ &&
          cleanedRPC[0].bxNum() == cleanedRPC[i].bxNum()) {
        phiF += cleanedRPC[i].offline_coord2();
        etaF += cleanedRPC[i].offline_eta2();
        phi += cleanedRPC[i].coord2();
        eta += cleanedRPC[i].eta2();
        nRPC++;
      } else {
        freeRPC.push_back(cleanedRPC[i]);
      }
    }
    l1t::MuonStub stub(cleanedRPC[0].etaRegion(),
                       cleanedRPC[0].phiRegion(),
                       cleanedRPC[0].depthRegion(),
                       cleanedRPC[0].tfLayer(),
                       0,
                       phi / nRPC,
                       0,
                       cleanedRPC[0].bxNum(),
                       2,
                       0,
                       eta / nRPC,
                       2,
                       0);
    stub.setOfflineQuantities(phiF / nRPC, phiF / nRPC, etaF / nRPC, etaF / nRPC);
    out.push_back(stub);
    cleanedRPC = freeRPC;
  };
  return out;
}

makeStubs()

l1t::MuonStubCollection L1TPhase2GMTEndcapStubProcessor::makeStubs	(	const MuonDigiCollection< CSCDetId, CSCCorrelatedLCTDigi > &	csc,
		const MuonDigiCollection< RPCDetId, RPCDigi > &	rpc,
		const L1TMuon::GeometryTranslator *	t,
		const edm::EventSetup &	iSetup
	)

Definition at line 242 of file L1TPhase2GMTEndcapStubProcessor.cc.

References MuonDigiCollection< IndexType, DigiType >::begin(), buildCSCOnlyStub(), buildRPCOnlyStub(), l1t::MuonStub::bxNum(), relativeConstraints::chamber, combineStubs(), MuonDigiCollection< IndexType, DigiType >::end(), maxBX_, minBX_, submitPVValidationJobs::t, and verbose_.

Referenced by Phase2L1TGMTStubProducer::produce().

                                 {
  l1t::MuonStubCollection cscStubs;
  auto chamber = csc.begin();
  auto chend = csc.end();
  for (; chamber != chend; ++chamber) {
    auto digi = (*chamber).second.first;
    auto dend = (*chamber).second.second;
    for (; digi != dend; ++digi) {
      l1t::MuonStub stub = buildCSCOnlyStub((*chamber).first, *digi, t);
      if (stub.bxNum() >= minBX_ && stub.bxNum() <= maxBX_)
        cscStubs.push_back(stub);
    }
  }

  l1t::MuonStubCollection rpcStubs;

  auto rpcchamber = cleaned.begin();
  auto rpcchend = cleaned.end();
  for (; rpcchamber != rpcchend; ++rpcchamber) {
    if ((*rpcchamber).first.region() == 0)
      continue;
    auto digi = (*rpcchamber).second.first;
    auto dend = (*rpcchamber).second.second;
    for (; digi != dend; ++digi) {
      l1t::MuonStub stub = buildRPCOnlyStub((*rpcchamber).first, *digi, t);
      if (stub.bxNum() >= minBX_ && stub.bxNum() <= maxBX_)
        rpcStubs.push_back(stub);
    }
  }

  l1t::MuonStubCollection combinedStubs = combineStubs(cscStubs, rpcStubs);

  if (verbose_) {
    printf("CSC Stubs\n");
    for (const auto& stub : cscStubs)
      printf(
          "CSC Stub bx=%d TF=%d etaRegion=%d phiRegion=%d depthRegion=%d  coord1=%f,%d coord2=%f,%d eta1=%f,%d "
          "eta2=%f,%d quality=%d etaQuality=%d\n",
          stub.bxNum(),
          stub.tfLayer(),
          stub.etaRegion(),
          stub.phiRegion(),
          stub.depthRegion(),
          stub.offline_coord1(),
          stub.coord1(),
          stub.offline_coord2(),
          stub.coord2(),
          stub.offline_eta1(),
          stub.eta1(),
          stub.offline_eta2(),
          stub.eta2(),
          stub.quality(),
          stub.etaQuality());
    printf("RPC Stubs\n");
    for (const auto& stub : rpcStubs)
      printf(
          "RPC Stub bx=%d TF=%d etaRegion=%d phiRegion=%d depthRegion=%d  coord1=%f,%d coord2=%f,%d eta1=%f,%d "
          "eta2=%f,%d quality=%d etaQuality=%d\n",
          stub.bxNum(),
          stub.tfLayer(),
          stub.etaRegion(),
          stub.phiRegion(),
          stub.depthRegion(),
          stub.offline_coord1(),
          stub.coord1(),
          stub.offline_coord2(),
          stub.coord2(),
          stub.offline_eta1(),
          stub.eta1(),
          stub.offline_eta2(),
          stub.eta2(),
          stub.quality(),
          stub.etaQuality());
    for (const auto& stub : combinedStubs)
      printf(
          "Combined Stub bx=%d TF=%d etaRegion=%d phiRegion=%d depthRegion=%d  coord1=%f,%d coord2=%f,%d eta1=%f,%d "
          "eta2=%f,%d quality=%d etaQuality=%d\n",
          stub.bxNum(),
          stub.tfLayer(),
          stub.etaRegion(),
          stub.phiRegion(),
          stub.depthRegion(),
          stub.offline_coord1(),
          stub.coord1(),
          stub.offline_coord2(),
          stub.coord2(),
          stub.offline_eta1(),
          stub.eta1(),
          stub.offline_eta2(),
          stub.eta2(),
          stub.quality(),
          stub.etaQuality());
  }

  return combinedStubs;
}

Member Data Documentation

coord1LSB_

double L1TPhase2GMTEndcapStubProcessor::coord1LSB_

private

Definition at line 34 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by buildCSCOnlyStub().

coord2LSB_

double L1TPhase2GMTEndcapStubProcessor::coord2LSB_

private

Definition at line 35 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by buildRPCOnlyStub().

eta1LSB_

double L1TPhase2GMTEndcapStubProcessor::eta1LSB_

private

Definition at line 36 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by buildCSCOnlyStub().

eta2LSB_

double L1TPhase2GMTEndcapStubProcessor::eta2LSB_

private

Definition at line 37 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by buildRPCOnlyStub().

etaMatch_

double L1TPhase2GMTEndcapStubProcessor::etaMatch_

private

Definition at line 38 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by combineStubs().

maxBX_

int L1TPhase2GMTEndcapStubProcessor::maxBX_

private

Definition at line 33 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by makeStubs().

minBX_

int L1TPhase2GMTEndcapStubProcessor::minBX_

private

Definition at line 32 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by makeStubs().

phiMatch_

double L1TPhase2GMTEndcapStubProcessor::phiMatch_

private

Definition at line 39 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by combineStubs().

verbose_

bool L1TPhase2GMTEndcapStubProcessor::verbose_

private

Definition at line 40 of file L1TPhase2GMTEndcapStubProcessor.h.

Referenced by makeStubs().