#include <HGCalTowerMap2DImpl.h>

Public Member Functions
void	buildTowerMap2D (const std::vector< edm::Ptr< l1t::HGCalTriggerCell >> &triggerCellsPtrs, l1t::HGCalTowerMapBxCollection &towermaps)

void	eventSetup (const edm::EventSetup &es)

	HGCalTowerMap2DImpl (const edm::ParameterSet &conf)

void	resetTowerMaps ()

Private Member Functions
std::vector< l1t::HGCalTowerMap >	newTowerMaps ()

Private Attributes
std::vector< double >	etaBins_

std::vector< double >	layerWeights_

int	nEtaBins_

int	nPhiBins_

std::vector< double >	phiBins_

HGCalTriggerTools	triggerTools_

bool	useLayerWeights_

Detailed Description

Author: : Thomas Strebler

Description: first iteration of HGCal Tower Maps

Definition at line 13 of file HGCalTowerMap2DImpl.h.

Constructor & Destructor Documentation

HGCalTowerMap2DImpl::HGCalTowerMap2DImpl ( const edm::ParameterSet & conf )

Definition at line 14 of file HGCalTowerMap2DImpl.cc.

References edm::errors::Configuration, PVValHelper::eta, etaBins_, Exception, createfilelist::int, nEtaBins_, newTowerMaps(), nPhiBins_, phi, and phiBins_.

                                                                       :
   nEtaBins_(conf.getParameter<int>("nEtaBins")),
   nPhiBins_(conf.getParameter<int>("nPhiBins")),
   etaBins_(conf.getParameter<std::vector<double> >("etaBins")),
   phiBins_(conf.getParameter<std::vector<double> >("phiBins")),
   useLayerWeights_(conf.getParameter<bool>("useLayerWeights")),
   layerWeights_(conf.getParameter< std::vector<double> >("layerWeights"))
 {
   edm::LogInfo("HGCalTowerMap2DImpl") << "Number of eta bins for the tower maps: " << nEtaBins_<<endl;
   edm::LogInfo("HGCalTowerMap2DImpl") << "Number of phi bins for the tower maps: " << nPhiBins_<<endl;
 
   if(!etaBins_.empty() && int(etaBins_.size())!=nEtaBins_+1){
     throw edm::Exception(edm::errors::Configuration, "Configuration")
       << "HGCalTowerMap2DImpl nEtaBins for the tower map not consistent with etaBins size"<<endl;
   }
   if(!phiBins_.empty() && int(phiBins_.size())!=nPhiBins_+1){
     throw edm::Exception(edm::errors::Configuration, "Configuration")
       << "HGCalTowerMap2DImpl nPhiBins for the tower map not consistent with phiBins size"<<endl;
   }
 
 
   std::vector<l1t::HGCalTowerMap> towerMaps = newTowerMaps();
   edm::LogInfo("HGCalTowerMap2DImpl") << "Eta bins for the tower maps: {";
   for(auto eta : towerMaps[0].etaBins()) edm::LogInfo("HGCalTowerMap2DImpl") << eta << ",";
   edm::LogInfo("HGCalTowerMap2DImpl") << "}" <<endl;
   edm::LogInfo("HGCalTowerMap2DImpl") << "Phi bins for the tower maps: {";
   for(auto phi : towerMaps[0].phiBins()) edm::LogInfo("HGCalTowerMap2DImpl") << phi << ",";
   edm::LogInfo("HGCalTowerMap2DImpl") << "}" <<endl;  
 
 }

Member Function Documentation

void HGCalTowerMap2DImpl::buildTowerMap2D	(	const std::vector< edm::Ptr< l1t::HGCalTriggerCell >> &	triggerCellsPtrs,
		l1t::HGCalTowerMapBxCollection &	towermaps
	)

Definition at line 68 of file HGCalTowerMap2DImpl.cc.

References mps_fire::i, HGCalTriggerTools::lastLayerBH(), HGCalTriggerTools::lastLayerEE(), layerWeights_, HGCalTriggerTools::layerWithOffset(), newTowerMaps(), p4, BXVector< T >::resize(), BXVector< T >::set(), l1t::HGCalTower::setEtEm(), l1t::HGCalTower::setEtHad(), l1t::L1Candidate::setHwEta(), l1t::L1Candidate::setHwPhi(), reco::LeafCandidate::setP4(), triggerTools_, and useLayerWeights_.

                         {
 
   std::vector<l1t::HGCalTowerMap> towerMapsTmp = newTowerMaps();
 
   for(const auto& tc : triggerCellsPtrs){
 
     unsigned layer = triggerTools_.layerWithOffset(tc->detId());
     int iEta = towerMapsTmp[layer-1].iEta(tc->eta());
     int iPhi = towerMapsTmp[layer-1].iPhi(tc->phi());
 
     double calibPt = tc->pt();
     if(useLayerWeights_) calibPt = layerWeights_[layer]*(tc->mipPt());
     math::PtEtaPhiMLorentzVector p4(calibPt,
                     tc->eta(),
                     tc->phi(),
                     0. );
 
     double etEm = layer<=triggerTools_.lastLayerEE() ? calibPt : 0;
     double etHad = layer>triggerTools_.lastLayerEE() ? calibPt : 0;
 
     l1t::HGCalTower tower;
     tower.setP4(p4);
     tower.setEtEm(etEm);
     tower.setEtHad(etHad);
     tower.setHwEta(iEta);
     tower.setHwPhi(iPhi);
 
     towerMapsTmp[layer-1].addTower(iEta,iPhi,tower);
 
   }
 
   /* store towerMaps in the persistent collection */
   towerMaps.resize(0, triggerTools_.lastLayerBH());
   int i=0;
   for(auto towerMap : towerMapsTmp){
     towerMaps.set( 0, i, towerMap);
     i++;
   }
 
 
 }

void HGCalTowerMap2DImpl::eventSetup ( const edm::EventSetup & es )

inline

Definition at line 26 of file HGCalTowerMap2DImpl.h.

References HGCalTriggerTools::eventSetup(), and triggerTools_.

     {
         triggerTools_.eventSetup(es);
     }

std::vector< l1t::HGCalTowerMap > HGCalTowerMap2DImpl::newTowerMaps ( )

private

Definition at line 46 of file HGCalTowerMap2DImpl.cc.

References etaBins_, HGCalTriggerTools::lastLayerBH(), genParticles_cff::map, SiStripPI::max, nEtaBins_, nPhiBins_, phiBins_, and triggerTools_.

Referenced by buildTowerMap2D(), and HGCalTowerMap2DImpl().

                                                              {
 
   //If no custom binning specified, assume uniform one
   l1t::HGCalTowerMap map;
   if(etaBins_.empty() || phiBins_.empty()){
     l1t::HGCalTowerMap mapTmp(nEtaBins_,nPhiBins_);
     map = mapTmp;
   }
   else{
     l1t::HGCalTowerMap mapTmp(etaBins_,phiBins_);
     map = mapTmp;
   }
 
   std::vector<l1t::HGCalTowerMap> towerMaps(max(triggerTools_.lastLayerBH(),unsigned(1)),map); //Always create at least 1 towerMap
   for(unsigned layer=0; layer<triggerTools_.lastLayerBH(); layer++) towerMaps[layer].setLayer(layer+1);
   return towerMaps;
 
 }