Phase2L1CaloPFClusterEmulator Class Reference

#include <L1Trigger/L1CaloTrigger/plugins/Phase2L1CaloPFClusterEmulator.cc>

Inheritance diagram for Phase2L1CaloPFClusterEmulator:

Public Member Functions
	Phase2L1CaloPFClusterEmulator (const edm::ParameterSet &)
	~Phase2L1CaloPFClusterEmulator () 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 EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const edm::EDGetTokenT< l1tp2::CaloTowerCollection >	caloTowerToken_

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

Description: Creates 3x3 PF clusters from GCTintTowers to be sent to correlator. Follows firmware logic, creates 8 clusters per (2+17+2)x(2+4+2).

Implementation: To be run together with Phase2L1CaloEGammaEmulator.

Definition at line 48 of file Phase2L1CaloPFClusterEmulator.cc.

Constructor & Destructor Documentation

Phase2L1CaloPFClusterEmulator()

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

explicit

Definition at line 65 of file Phase2L1CaloPFClusterEmulator.cc.

    : caloTowerToken_(consumes<l1tp2::CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("gctFullTowers"))) {
  produces<l1tp2::CaloPFClusterCollection>("GCTPFCluster");
}

~Phase2L1CaloPFClusterEmulator()

Phase2L1CaloPFClusterEmulator::~Phase2L1CaloPFClusterEmulator ( )

overridedefault

Member Function Documentation

fillDescriptions()

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

static

Definition at line 170 of file Phase2L1CaloPFClusterEmulator.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                               {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("gctFullTowers", edm::InputTag("l1tPhase2L1CaloEGammaEmulator", "GCTFullTowers"));
  descriptions.addWithDefaultLabel(desc);
}

produce()

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

overrideprivate

Definition at line 75 of file Phase2L1CaloPFClusterEmulator.cc.

References caloTowerToken_, gctpf::GCTpfcluster_t::et, gctpf::GCTpfcluster_t::eta, gctpf::PFcluster_t::GCTpfclusters, mps_fire::i, hcalRecHitTable_cff::ieta, iEvent, hcalRecHitTable_cff::iphi, edm::HandleBase::isValid(), dqmdumpme::k, eostools::move(), nPFClusterSLR, nSLR, nTowerEta, nTowerEtaSLR, nTowerPhi, nTowerPhiSLR, gctpf::pfcluster(), gctpf::GCTpfcluster_t::phi, l1tp2::CaloPFCluster::setClusterEt(), l1tp2::CaloPFCluster::setClusterEta(), l1tp2::CaloPFCluster::setClusterIEta(), l1tp2::CaloPFCluster::setClusterIPhi(), and l1tp2::CaloPFCluster::setClusterPhi().

                                                                                         {
  using namespace edm;
  std::unique_ptr<l1tp2::CaloPFClusterCollection> pfclusterCands(make_unique<l1tp2::CaloPFClusterCollection>());

  edm::Handle<std::vector<l1tp2::CaloTower>> caloTowerCollection;
  iEvent.getByToken(caloTowerToken_, caloTowerCollection);
  if (!caloTowerCollection.isValid())
    cms::Exception("Phase2L1CaloPFClusterEmulator") << "Failed to get towers from caloTowerCollection!";

  float GCTintTowers[nTowerEta][nTowerPhi];
  float realEta[nTowerEta][nTowerPhi];
  float realPhi[nTowerEta][nTowerPhi];
  for (const l1tp2::CaloTower& i : *caloTowerCollection) {
    int ieta = i.towerIEta();
    int iphi = i.towerIPhi();
    GCTintTowers[ieta][iphi] = i.ecalTowerEt();
    realEta[ieta][iphi] = i.towerEta();
    realPhi[ieta][iphi] = i.towerPhi();
  }

  float regions[nSLR][nTowerEtaSLR][nTowerPhiSLR];

  for (int ieta = 0; ieta < nTowerEtaSLR; ieta++) {
    for (int iphi = 0; iphi < nTowerPhiSLR; iphi++) {
      for (int k = 0; k < nSLR; k++) {
        regions[k][ieta][iphi] = 0.;
      }
    }
  }

  //Assign ETs to 36 21x8 regions

  for (int ieta = 0; ieta < nTowerEtaSLR; ieta++) {
    for (int iphi = 0; iphi < nTowerPhiSLR; iphi++) {
      if (ieta > 1) {
        if (iphi > 1)
          regions[0][ieta][iphi] = GCTintTowers[ieta - 2][iphi - 2];
        for (int k = 1; k < 17; k++) {
          regions[k * 2][ieta][iphi] = GCTintTowers[ieta - 2][iphi + k * 4 - 2];
        }
        if (iphi < 6)
          regions[34][ieta][iphi] = GCTintTowers[ieta - 2][iphi + 66];
      }
      if (ieta < 19) {
        if (iphi > 1)
          regions[1][ieta][iphi] = GCTintTowers[ieta + 15][iphi - 2];
        for (int k = 1; k < 17; k++) {
          regions[k * 2 + 1][ieta][iphi] = GCTintTowers[ieta + 15][iphi + k * 4 - 2];
        }
        if (iphi < 6)
          regions[35][ieta][iphi] = GCTintTowers[ieta + 15][iphi + 66];
      }
    }
  }

  float temporary[nTowerEtaSLR][nTowerPhiSLR];
  int etaoffset = 0;
  int phioffset = 0;

  //Use same code from firmware for finding clusters
  for (int k = 0; k < nSLR; k++) {
    for (int ieta = 0; ieta < nTowerEtaSLR; ieta++) {
      for (int iphi = 0; iphi < nTowerPhiSLR; iphi++) {
        temporary[ieta][iphi] = regions[k][ieta][iphi];
      }
    }
    if (k % 2 == 0)
      etaoffset = 0;
    else
      etaoffset = 17;
    if (k > 1 && k % 2 == 0)
      phioffset = phioffset + 4;

    gctpf::PFcluster_t tempPfclusters;
    tempPfclusters = gctpf::pfcluster(temporary, etaoffset, phioffset);

    for (int i = 0; i < nPFClusterSLR; i++) {
      int gcteta = tempPfclusters.GCTpfclusters[i].eta;
      int gctphi = tempPfclusters.GCTpfclusters[i].phi;
      float towereta = realEta[gcteta][gctphi];
      float towerphi = realPhi[gcteta][gctphi];
      l1tp2::CaloPFCluster l1CaloPFCluster;
      l1CaloPFCluster.setClusterEt(tempPfclusters.GCTpfclusters[i].et);
      l1CaloPFCluster.setClusterIEta(gcteta);
      l1CaloPFCluster.setClusterIPhi(gctphi);
      l1CaloPFCluster.setClusterEta(towereta);
      l1CaloPFCluster.setClusterPhi(towerphi);
      pfclusterCands->push_back(l1CaloPFCluster);
    }
  }

  iEvent.put(std::move(pfclusterCands), "GCTPFCluster");
}

Member Data Documentation

caloTowerToken_

const edm::EDGetTokenT<l1tp2::CaloTowerCollection> Phase2L1CaloPFClusterEmulator::caloTowerToken_

private

Definition at line 59 of file Phase2L1CaloPFClusterEmulator.cc.

Referenced by produce().