L1TCaloLayer1 Class Reference

#include <L1Trigger/L1TCaloLayer1/plugins/L1TCaloLayer1.cc>

Inheritance diagram for L1TCaloLayer1:

Public Member Functions
	L1TCaloLayer1 (const edm::ParameterSet &)
	~L1TCaloLayer1 ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
virtual void	beginJob () override
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override
virtual void	endJob () override
virtual void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::vector< std::vector < std::vector< uint32_t > > >	ecalLUT
edm::EDGetTokenT < EcalTrigPrimDigiCollection >	ecalTPSource
std::string	ecalTPSourceLabel
std::vector< std::vector < std::vector< uint32_t > > >	hcalLUT
edm::EDGetTokenT < HcalTrigPrimDigiCollection >	hcalTPSource
std::string	hcalTPSourceLabel
std::vector< std::vector < uint32_t > >	hfLUT
UCTLayer1 *	layer1
std::vector< UCTTower * >	twrList
bool	unpackEcalMask
bool	unpackHcalMask
bool	useCalib
bool	useECALLUT
bool	useHCALLUT
bool	useHFLUT
bool	useLSB
bool	verbose

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 59 of file L1TCaloLayer1.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 115 of file L1TCaloLayer1.cc.

References a, b, layer1, HLT_25ns10e33_v2_cff::towers, and twrList.

                                                           :
  ecalTPSource(consumes<EcalTrigPrimDigiCollection>(iConfig.getParameter<edm::InputTag>("ecalToken"))),
  ecalTPSourceLabel(iConfig.getParameter<edm::InputTag>("ecalToken").label()),
  hcalTPSource(consumes<HcalTrigPrimDigiCollection>(iConfig.getParameter<edm::InputTag>("hcalToken"))),
  hcalTPSourceLabel(iConfig.getParameter<edm::InputTag>("hcalToken").label()),
  ecalLUT(28, std::vector< std::vector<uint32_t> >(2, std::vector<uint32_t>(256))),
  hcalLUT(28, std::vector< std::vector<uint32_t> >(2, std::vector<uint32_t>(256))),
  hfLUT(12, std::vector < uint32_t >(256)),
  useLSB(iConfig.getParameter<bool>("useLSB")),
  useCalib(iConfig.getParameter<bool>("useCalib")),
  useECALLUT(iConfig.getParameter<bool>("useECALLUT")),
  useHCALLUT(iConfig.getParameter<bool>("useHCALLUT")),
  useHFLUT(iConfig.getParameter<bool>("useHFLUT")),
  verbose(iConfig.getParameter<bool>("verbose")), 
  unpackHcalMask(iConfig.getParameter<bool>("unpackHcalMask")),
  unpackEcalMask(iConfig.getParameter<bool>("unpackEcalMask"))
{
  produces<CaloTowerBxCollection>();
  produces<L1CaloRegionCollection>();
  layer1 = new UCTLayer1;
  vector<UCTCrate*> crates = layer1->getCrates();
  for(uint32_t crt = 0; crt < crates.size(); crt++) {
    vector<UCTCard*> cards = crates[crt]->getCards();
    for(uint32_t crd = 0; crd < cards.size(); crd++) {
      vector<UCTRegion*> regions = cards[crd]->getRegions();
      for(uint32_t rgn = 0; rgn < regions.size(); rgn++) {
    vector<UCTTower*> towers = regions[rgn]->getTowers();
    for(uint32_t twr = 0; twr < towers.size(); twr++) {
      twrList.push_back(towers[twr]);
    }
      }
    }
  }

  // This sort corresponds to the sort condition on
  // the output CaloTowerBxCollection
  std::sort(twrList.begin(), twrList.end(), [](UCTTower* a, UCTTower* b) {
      return CaloTools::caloTowerHash(a->caloEta(), a->caloPhi()) < CaloTools::caloTowerHash(b->caloEta(), b->caloPhi());
      });
}

L1TCaloLayer1::~L1TCaloLayer1

(

)

Definition at line 156 of file L1TCaloLayer1.cc.

References layer1.

                              {
  if(layer1 != 0) delete layer1;
}

Member Function Documentation

void L1TCaloLayer1::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 288 of file L1TCaloLayer1.cc.

{
}

void L1TCaloLayer1::beginRun ( edm::Run const &  iRun, edm::EventSetup const &  iSetup  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 299 of file L1TCaloLayer1.cc.

References ecalLUT, hcalLUT, hfLUT, L1TCaloLayer1FetchLUTs(), LOG_ERROR, twrList, useCalib, useECALLUT, useHCALLUT, useHFLUT, and useLSB.

{
  if(!L1TCaloLayer1FetchLUTs(iSetup, ecalLUT, hcalLUT, hfLUT, useLSB, useCalib, useECALLUT, useHCALLUT, useHFLUT)) {
    LOG_ERROR << "L1TCaloLayer1::beginRun: failed to fetch LUTS - using unity" << std::endl;
  }
  for(uint32_t twr = 0; twr < twrList.size(); twr++) {
    twrList[twr]->setECALLUT(&ecalLUT);
    twrList[twr]->setHCALLUT(&hcalLUT);
    twrList[twr]->setHFLUT(&hfLUT);
  }
}

void L1TCaloLayer1::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 294 of file L1TCaloLayer1.cc.

                      {
}

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

static

Definition at line 338 of file L1TCaloLayer1.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                          {
  //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);
}

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 166 of file L1TCaloLayer1.cc.

References funct::abs(), ecalTPSource, g, edm::Event::getByToken(), hcalTPSource, layer1, LOG_ERROR, edm::Event::put(), HLT_25ns10e33_v2_cff::region, l1t::L1Candidate::setHwEta(), l1t::CaloTower::setHwEtEm(), l1t::CaloTower::setHwEtHad(), l1t::CaloTower::setHwEtRatio(), l1t::L1Candidate::setHwPhi(), l1t::L1Candidate::setHwPt(), l1t::L1Candidate::setHwQual(), lumiQTWidget::t, twrList, unpackEcalMask, and unpackHcalMask.

{
  using namespace edm;

  edm::Handle<EcalTrigPrimDigiCollection> ecalTPs;
  iEvent.getByToken(ecalTPSource, ecalTPs);
  edm::Handle<HcalTrigPrimDigiCollection> hcalTPs;
  iEvent.getByToken(hcalTPSource, hcalTPs);

  std::auto_ptr<CaloTowerBxCollection> towersColl (new CaloTowerBxCollection);
  std::auto_ptr<L1CaloRegionCollection> rgnCollection (new L1CaloRegionCollection);

  uint32_t expectedTotalET = 0;
  if(!layer1->clearEvent()) {
    LOG_ERROR << "UCT: Failed to clear event" << std::endl;
    return;
  }

  for ( const auto& ecalTp : *ecalTPs ) {
    if ( unpackEcalMask && ((ecalTp.sample(0).raw()>>13) & 0x1) ) continue;
    int caloEta = ecalTp.id().ieta();
    int caloPhi = ecalTp.id().iphi();
    int et = ecalTp.compressedEt();
    bool fgVeto = ecalTp.fineGrain();
    UCTTowerIndex t = UCTTowerIndex(caloEta, caloPhi);
    if(!layer1->setECALData(t,fgVeto,et)) {
      LOG_ERROR << "UCT: Failed loading an ECAL tower" << std::endl;
      return;
    }
    expectedTotalET += et;
  }

  for ( const auto& hcalTp : *hcalTPs ) {
    if ( unpackHcalMask && ((hcalTp.sample(0).raw()>>13) & 0x1) ) continue;
    int caloEta = hcalTp.id().ieta();
    uint32_t absCaloEta = abs(caloEta);
    // Tower 29 is not used by Layer-1
    if(absCaloEta == 29) {
      continue;
    }
    // Prevent usage of HF TPs with Layer-1 emulator if HCAL TPs are old style
    else if(hcalTp.id().version() == 0 && absCaloEta > 29) {
      continue;
    }
    else if(absCaloEta <= 41) {
      int caloPhi = hcalTp.id().iphi();
      int et = hcalTp.SOI_compressedEt();
      bool fg  = hcalTp.t0().fineGrain(0);
      bool fg2 = hcalTp.t0().fineGrain(1);
      if(caloPhi <= 72) {
        UCTTowerIndex t = UCTTowerIndex(caloEta, caloPhi);
        uint32_t featureBits = 0;
        if(fg)  featureBits |= 0b01;
        // fg2 should only be set for HF
        if(absCaloEta > 29 && fg2) featureBits |= 0b10;
        if(!layer1->setHCALData(t, featureBits, et)) {
          LOG_ERROR << "caloEta = " << caloEta << "; caloPhi =" << caloPhi << std::endl;
          LOG_ERROR << "UCT: Failed loading an HCAL tower" << std::endl;
          return;
        }
        expectedTotalET += et;
      }
      else {
    LOG_ERROR << "Illegal Tower: caloEta = " << caloEta << "; caloPhi =" << caloPhi << "; et = " << et << std::endl;    
      }
    }
    else {
      LOG_ERROR << "Illegal Tower: caloEta = " << caloEta << std::endl;
    }
  }
  
   //Process
  if(!layer1->process()) {
    LOG_ERROR << "UCT: Failed to process layer 1" << std::endl;
  }

  int theBX = 0; // Currently we only read and process the "hit" BX only

  for(uint32_t twr = 0; twr < twrList.size(); twr++) {
    CaloTower caloTower;
    caloTower.setHwPt(twrList[twr]->et());               // Bits 0-8 of the 16-bit word per the interface protocol document
    caloTower.setHwEtRatio(twrList[twr]->er());          // Bits 9-11 of the 16-bit word per the interface protocol document
    caloTower.setHwQual(twrList[twr]->miscBits());       // Bits 12-15 of the 16-bit word per the interface protocol document
    caloTower.setHwEta(twrList[twr]->caloEta());         // caloEta = 1-28 and 30-41
    caloTower.setHwPhi(twrList[twr]->caloPhi());         // caloPhi = 1-72
    caloTower.setHwEtEm(twrList[twr]->getEcalET());      // This is provided as a courtesy - not available to hardware
    caloTower.setHwEtHad(twrList[twr]->getHcalET());     // This is provided as a courtesy - not available to hardware
    towersColl->push_back(theBX, caloTower);
  }

  iEvent.put(towersColl);

  UCTGeometry g;
  vector<UCTCrate*> crates = layer1->getCrates();
  for(uint32_t crt = 0; crt < crates.size(); crt++) {
    vector<UCTCard*> cards = crates[crt]->getCards();
    for(uint32_t crd = 0; crd < cards.size(); crd++) {
      vector<UCTRegion*> regions = cards[crd]->getRegions();
      for(uint32_t rgn = 0; rgn < regions.size(); rgn++) {
    uint32_t rawData = regions[rgn]->rawData();
    uint32_t regionData = rawData & 0x0000FFFF;
    uint32_t crate = regions[rgn]->getCrate();
    uint32_t card = regions[rgn]->getCard();
    uint32_t region = regions[rgn]->getRegion();
    bool negativeEta = regions[rgn]->isNegativeEta();
    uint32_t rPhi = g.getUCTRegionPhiIndex(crate, card);
    if(region < NRegionsInCard) { // We only store the Barrel and Endcap - HF has changed in the upgrade
      uint32_t rEta = 10 - region; // UCT region is 0-6 for B/E but GCT eta goes 0-21, 0-3 -HF, 4-10 -B/E, 11-17 +B/E, 18-21 +HF
      if(!negativeEta) rEta = 11 + region; // Positive eta portion is offset by 11
      rgnCollection->push_back(L1CaloRegion((uint16_t) regionData, (unsigned) rEta, (unsigned) rPhi, (int16_t) 0));
    }
      }
    }
  }  
  iEvent.put(rgnCollection);

}

Member Data Documentation

std::vector< std::vector< std::vector < uint32_t > > > L1TCaloLayer1::ecalLUT

private

Definition at line 84 of file L1TCaloLayer1.cc.

Referenced by beginRun().

edm::EDGetTokenT<EcalTrigPrimDigiCollection> L1TCaloLayer1::ecalTPSource

private

Definition at line 79 of file L1TCaloLayer1.cc.

Referenced by produce().

std::string L1TCaloLayer1::ecalTPSourceLabel

private

Definition at line 80 of file L1TCaloLayer1.cc.

std::vector< std::vector< std::vector < uint32_t > > > L1TCaloLayer1::hcalLUT

private

Definition at line 85 of file L1TCaloLayer1.cc.

Referenced by beginRun().

edm::EDGetTokenT<HcalTrigPrimDigiCollection> L1TCaloLayer1::hcalTPSource

private

Definition at line 81 of file L1TCaloLayer1.cc.

Referenced by produce().

std::string L1TCaloLayer1::hcalTPSourceLabel

private

Definition at line 82 of file L1TCaloLayer1.cc.

std::vector< std::vector< uint32_t > > L1TCaloLayer1::hfLUT

private

Definition at line 86 of file L1TCaloLayer1.cc.

Referenced by beginRun().

UCTLayer1* L1TCaloLayer1::layer1

private

Definition at line 99 of file L1TCaloLayer1.cc.

Referenced by L1TCaloLayer1(), produce(), and ~L1TCaloLayer1().

std::vector< UCTTower* > L1TCaloLayer1::twrList

private

Definition at line 88 of file L1TCaloLayer1.cc.

Referenced by beginRun(), L1TCaloLayer1(), and produce().

bool L1TCaloLayer1::unpackEcalMask

private

Definition at line 97 of file L1TCaloLayer1.cc.

Referenced by produce().

bool L1TCaloLayer1::unpackHcalMask

private

Definition at line 96 of file L1TCaloLayer1.cc.

Referenced by produce().

bool L1TCaloLayer1::useCalib

private

Definition at line 91 of file L1TCaloLayer1.cc.

Referenced by beginRun().

bool L1TCaloLayer1::useECALLUT

private

Definition at line 92 of file L1TCaloLayer1.cc.

Referenced by beginRun().

bool L1TCaloLayer1::useHCALLUT

private

Definition at line 93 of file L1TCaloLayer1.cc.

Referenced by beginRun().

bool L1TCaloLayer1::useHFLUT

private

Definition at line 94 of file L1TCaloLayer1.cc.

Referenced by beginRun().

bool L1TCaloLayer1::useLSB

private

Definition at line 90 of file L1TCaloLayer1.cc.

Referenced by beginRun().

bool L1TCaloLayer1::verbose

private

Definition at line 95 of file L1TCaloLayer1.cc.

Referenced by python.diff_provenance.difference::list_diff(), python.diffProv.difference::list_diff(), python.diff_provenance.difference::module_diff(), python.diffProv.difference::module_diff(), python.diff_provenance.difference::onefilemodules(), python.diffProv.difference::onefilemodules(), core.TriggerMatchAnalyzer.TriggerMatchAnalyzer::process(), core.SkimAnalyzerCount.SkimAnalyzerCount::process(), and objects.VertexAnalyzer.VertexAnalyzer::process().