#include <L1ExtraParticlesProd.h>

Inheritance diagram for L1ExtraParticlesProd:

Public Member Functions
	L1ExtraParticlesProd (const edm::ParameterSet &)

	~L1ExtraParticlesProd () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
math::PtEtaPhiMLorentzVector	gctLorentzVector (const double &et, const L1GctCand &cand, const L1CaloGeometry *geom, bool central)

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

Private Attributes
edm::ESGetToken< L1CaloGeometry, L1CaloGeometryRecord >	caloGeomToken_

edm::InputTag	cenJetSource_

bool	centralBxOnly_

edm::ESGetToken< L1CaloEtScale, L1EmEtScaleRcd >	emScaleToken_

edm::InputTag	etHadSource_

edm::InputTag	etMissSource_

edm::InputTag	etTotSource_

edm::InputTag	forJetSource_

edm::InputTag	hfRingBitCountsSource_

edm::ESGetToken< L1CaloEtScale, L1HfRingEtScaleRcd >	hfRingEtScaleToken_

edm::InputTag	hfRingEtSumsSource_

edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcd >	htMissScaleToken_

edm::InputTag	htMissSource_

bool	ignoreHtMiss_

edm::InputTag	isoEmSource_

edm::InputTag	isoTauJetSource_

edm::ESGetToken< L1GctJetFinderParams, L1GctJetFinderParamsRcd >	jetFinderParamsToken_

edm::ESGetToken< L1CaloEtScale, L1JetEtScaleRcd >	jetScaleToken_

edm::InputTag	muonSource_

edm::ESGetToken< L1MuTriggerPtScale, L1MuTriggerPtScaleRcd >	muPtScaleToken_

edm::ESGetToken< L1MuTriggerScales, L1MuTriggerScalesRcd >	muScalesToken_

edm::InputTag	nonIsoEmSource_

bool	produceCaloParticles_

bool	produceMuonParticles_

edm::InputTag	tauJetSource_

Static Private Attributes
static const double	muonMassGeV_ = 0.105658369

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

L1Trigger/L1ExtraFromDigis/interface/L1ExtraParticlesProd.h

Author: Werner Sun

Description: producer of L1Extra particle objects from Level-1 hardware objects.

src/L1ExtraParticlesProd/src/L1ExtraParticlesProd.cc

Definition at line 56 of file L1ExtraParticlesProd.h.

Constructor & Destructor Documentation

◆ L1ExtraParticlesProd()

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

explicit

Definition at line 51 of file L1ExtraParticlesProd.cc.

     : produceMuonParticles_(iConfig.getParameter<bool>("produceMuonParticles")),
       muonSource_(iConfig.getParameter<edm::InputTag>("muonSource")),
       produceCaloParticles_(iConfig.getParameter<bool>("produceCaloParticles")),
       isoEmSource_(iConfig.getParameter<edm::InputTag>("isolatedEmSource")),
       nonIsoEmSource_(iConfig.getParameter<edm::InputTag>("nonIsolatedEmSource")),
       cenJetSource_(iConfig.getParameter<edm::InputTag>("centralJetSource")),
       forJetSource_(iConfig.getParameter<edm::InputTag>("forwardJetSource")),
       tauJetSource_(iConfig.getParameter<edm::InputTag>("tauJetSource")),
       isoTauJetSource_(iConfig.getParameter<edm::InputTag>("isoTauJetSource")),
       etTotSource_(iConfig.getParameter<edm::InputTag>("etTotalSource")),
       etHadSource_(iConfig.getParameter<edm::InputTag>("etHadSource")),
       etMissSource_(iConfig.getParameter<edm::InputTag>("etMissSource")),
       htMissSource_(iConfig.getParameter<edm::InputTag>("htMissSource")),
       hfRingEtSumsSource_(iConfig.getParameter<edm::InputTag>("hfRingEtSumsSource")),
       hfRingBitCountsSource_(iConfig.getParameter<edm::InputTag>("hfRingBitCountsSource")),
       centralBxOnly_(iConfig.getParameter<bool>("centralBxOnly")),
       ignoreHtMiss_(iConfig.getParameter<bool>("ignoreHtMiss")) {
   using namespace l1extra;
  
   // register your products
   produces<L1EmParticleCollection>("Isolated");
   produces<L1EmParticleCollection>("NonIsolated");
   produces<L1JetParticleCollection>("Central");
   produces<L1JetParticleCollection>("Forward");
   produces<L1JetParticleCollection>("Tau");
   produces<L1JetParticleCollection>("IsoTau");
   produces<L1MuonParticleCollection>();
   produces<L1EtMissParticleCollection>("MET");
   produces<L1EtMissParticleCollection>("MHT");
   produces<L1HFRingsCollection>();
  
   // now do what ever other initialization is needed
   consumes<L1MuGMTReadoutCollection>(muonSource_);
   consumes<L1GctEmCandCollection>(isoEmSource_);
   consumes<L1GctEmCandCollection>(nonIsoEmSource_);
   consumes<L1GctJetCandCollection>(cenJetSource_);
   consumes<L1GctJetCandCollection>(forJetSource_);
   consumes<L1GctJetCandCollection>(tauJetSource_);
   consumes<L1GctJetCandCollection>(isoTauJetSource_);
   consumes<L1GctEtTotalCollection>(etTotSource_);
   consumes<L1GctEtMissCollection>(etMissSource_);
   consumes<L1GctEtHadCollection>(etHadSource_);
   consumes<L1GctHtMissCollection>(htMissSource_);
   consumes<L1GctHFRingEtSumsCollection>(hfRingEtSumsSource_);
   consumes<L1GctHFBitCountsCollection>(hfRingBitCountsSource_);
  
   if (produceMuonParticles_) {
     muScalesToken_ = esConsumes<L1MuTriggerScales, L1MuTriggerScalesRcd>();
     muPtScaleToken_ = esConsumes<L1MuTriggerPtScale, L1MuTriggerPtScaleRcd>();
   }
   if (produceCaloParticles_) {
     caloGeomToken_ = esConsumes<L1CaloGeometry, L1CaloGeometryRecord>();
     emScaleToken_ = esConsumes<L1CaloEtScale, L1EmEtScaleRcd>();
     jetScaleToken_ = esConsumes<L1CaloEtScale, L1JetEtScaleRcd>();
     jetFinderParamsToken_ = esConsumes<L1GctJetFinderParams, L1GctJetFinderParamsRcd>();
     hfRingEtScaleToken_ = esConsumes<L1CaloEtScale, L1HfRingEtScaleRcd>();
     if (!ignoreHtMiss_) {
       htMissScaleToken_ = esConsumes<L1CaloEtScale, L1HtMissScaleRcd>();
     }
   }
 }

References caloGeomToken_, cenJetSource_, emScaleToken_, etHadSource_, etMissSource_, etTotSource_, forJetSource_, hfRingBitCountsSource_, hfRingEtScaleToken_, hfRingEtSumsSource_, htMissScaleToken_, htMissSource_, ignoreHtMiss_, isoEmSource_, isoTauJetSource_, jetFinderParamsToken_, jetScaleToken_, muonSource_, muPtScaleToken_, muScalesToken_, nonIsoEmSource_, produceCaloParticles_, produceMuonParticles_, and tauJetSource_.

◆ ~L1ExtraParticlesProd()

L1ExtraParticlesProd::~L1ExtraParticlesProd ( )

override

Definition at line 114 of file L1ExtraParticlesProd.cc.

                                             {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

◆ gctLorentzVector()

math::PtEtaPhiMLorentzVector L1ExtraParticlesProd::gctLorentzVector	(	const double &	et,
		const L1GctCand &	cand,
		const L1CaloGeometry *	geom,
		bool	central
	)

private

Definition at line 870 of file L1ExtraParticlesProd.cc.

                                                                                   {
   // To keep x and y components non-zero.
   double etCorr = et + 1.e-6;  // protect against roundoff, not only for et=0
  
   double eta = geom->etaBinCenter(cand.etaIndex(), central);
  
   //    double tanThOver2 = exp( -eta ) ;
   //    double ez = etCorr * ( 1. - tanThOver2 * tanThOver2 ) / ( 2. *
   //    tanThOver2 ); double e  = etCorr * ( 1. + tanThOver2 * tanThOver2 ) /
   //    ( 2. * tanThOver2 );
  
   double phi = geom->emJetPhiBinCenter(cand.phiIndex());
  
   //    return math::XYZTLorentzVector( etCorr * cos( phi ),
   //                   etCorr * sin( phi ),
   //                   ez,
   //                   e ) ;
   return math::PtEtaPhiMLorentzVector(etCorr, eta, phi, 0.);
 }

References central, EgHLTOffHistBins_cfi::et, PVValHelper::eta, relativeConstraints::geom, and phi.

Referenced by produce().

◆ produce()

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

overrideprivate

Definition at line 124 of file L1ExtraParticlesProd.cc.

                                                                                 {
   using namespace edm;
   using namespace l1extra;
   using namespace std;
  
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~ Muons ~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  
   unique_ptr<L1MuonParticleCollection> muColl(new L1MuonParticleCollection);
  
   if (produceMuonParticles_) {
     ESHandle<L1MuTriggerScales> muScales = iSetup.getHandle(muScalesToken_);
  
     ESHandle<L1MuTriggerPtScale> muPtScale = iSetup.getHandle(muPtScaleToken_);
  
     Handle<L1MuGMTReadoutCollection> hwMuCollection;
     iEvent.getByLabel(muonSource_, hwMuCollection);
  
     vector<L1MuGMTExtendedCand> hwMuCands;
  
     if (!hwMuCollection.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1MuGMTReadoutCollection with " << muonSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       if (centralBxOnly_) {
         // Get GMT candidates from central bunch crossing only
         hwMuCands = hwMuCollection->getRecord().getGMTCands();
       } else {
         // Get GMT candidates from all bunch crossings
         vector<L1MuGMTReadoutRecord> records = hwMuCollection->getRecords();
         vector<L1MuGMTReadoutRecord>::const_iterator rItr = records.begin();
         vector<L1MuGMTReadoutRecord>::const_iterator rEnd = records.end();
  
         for (; rItr != rEnd; ++rItr) {
           vector<L1MuGMTExtendedCand> tmpCands = rItr->getGMTCands();
  
           hwMuCands.insert(hwMuCands.end(), tmpCands.begin(), tmpCands.end());
         }
       }
  
       //       cout << "HW muons" << endl ;
       vector<L1MuGMTExtendedCand>::const_iterator muItr = hwMuCands.begin();
       vector<L1MuGMTExtendedCand>::const_iterator muEnd = hwMuCands.end();
       for (int i = 0; muItr != muEnd; ++muItr, ++i) {
         //   cout << "#" << i
         //        << " name " << muItr->name()
         //        << " empty " << muItr->empty()
         //        << " pt " << muItr->ptIndex()
         //        << " eta " << muItr->etaIndex()
         //        << " phi " << muItr->phiIndex()
         //        << " iso " << muItr->isol()
         //        << " mip " << muItr->mip()
         //        << " bx " << muItr->bx()
         //        << endl ;
  
         if (!muItr->empty()) {
           // keep x and y components non-zero and protect against roundoff.
           double pt = muPtScale->getPtScale()->getLowEdge(muItr->ptIndex()) + 1.e-6;
  
           //        cout << "L1Extra pt " << pt << endl ;
  
           double eta = muScales->getGMTEtaScale()->getCenter(muItr->etaIndex());
  
           double phi = muScales->getPhiScale()->getLowEdge(muItr->phiIndex());
  
           math::PtEtaPhiMLorentzVector p4(pt, eta, phi, muonMassGeV_);
  
           muColl->push_back(L1MuonParticle(muItr->charge(), p4, *muItr, muItr->bx()));
         }
       }
     }
   }
  
   OrphanHandle<L1MuonParticleCollection> muHandle = iEvent.put(std::move(muColl));
  
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~ Calorimeter ~~~~~~~~~~~~~~~~~~~~
   // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  
   unique_ptr<L1EmParticleCollection> isoEmColl(new L1EmParticleCollection);
  
   unique_ptr<L1EmParticleCollection> nonIsoEmColl(new L1EmParticleCollection);
  
   unique_ptr<L1JetParticleCollection> cenJetColl(new L1JetParticleCollection);
  
   unique_ptr<L1JetParticleCollection> forJetColl(new L1JetParticleCollection);
  
   unique_ptr<L1JetParticleCollection> tauJetColl(new L1JetParticleCollection);
  
   unique_ptr<L1JetParticleCollection> isoTauJetColl(new L1JetParticleCollection);
  
   unique_ptr<L1EtMissParticleCollection> etMissColl(new L1EtMissParticleCollection);
  
   unique_ptr<L1EtMissParticleCollection> htMissColl(new L1EtMissParticleCollection);
  
   unique_ptr<L1HFRingsCollection> hfRingsColl(new L1HFRingsCollection);
  
   if (produceCaloParticles_) {
     // ~~~~~~~~~~~~~~~~~~~~ Geometry ~~~~~~~~~~~~~~~~~~~~
  
     ESHandle<L1CaloGeometry> caloGeomESH = iSetup.getHandle(caloGeomToken_);
     const L1CaloGeometry *caloGeom = &(*caloGeomESH);
  
     // ~~~~~~~~~~~~~~~~~~~~ EM ~~~~~~~~~~~~~~~~~~~~
  
     ESHandle<L1CaloEtScale> emScale = iSetup.getHandle(emScaleToken_);
  
     // Isolated EM
     Handle<L1GctEmCandCollection> hwIsoEmCands;
     iEvent.getByLabel(isoEmSource_, hwIsoEmCands);
  
     if (!hwIsoEmCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEmCandCollection with " << isoEmSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       //       cout << "HW iso EM" << endl ;
  
       L1GctEmCandCollection::const_iterator emItr = hwIsoEmCands->begin();
       L1GctEmCandCollection::const_iterator emEnd = hwIsoEmCands->end();
       for (int i = 0; emItr != emEnd; ++emItr, ++i) {
         //   cout << "#" << i
         //        << " name " << emItr->name()
         //        << " empty " << emItr->empty()
         //        << " rank " << emItr->rank()
         //        << " eta " << emItr->etaIndex()
         //        << " sign " << emItr->etaSign()
         //        << " phi " << emItr->phiIndex()
         //        << " iso " << emItr->isolated()
         //        << " bx " << emItr->bx()
         //        << endl ;
  
         if (!emItr->empty() && (!centralBxOnly_ || emItr->bx() == 0)) {
           double et = emScale->et(emItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           isoEmColl->push_back(L1EmParticle(
               gctLorentzVector(et, *emItr, caloGeom, true), Ref<L1GctEmCandCollection>(hwIsoEmCands, i), emItr->bx()));
         }
       }
     }
  
     // Non-isolated EM
     Handle<L1GctEmCandCollection> hwNonIsoEmCands;
     iEvent.getByLabel(nonIsoEmSource_, hwNonIsoEmCands);
  
     if (!hwNonIsoEmCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEmCandCollection with " << nonIsoEmSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       //       cout << "HW non-iso EM" << endl ;
       L1GctEmCandCollection::const_iterator emItr = hwNonIsoEmCands->begin();
       L1GctEmCandCollection::const_iterator emEnd = hwNonIsoEmCands->end();
       for (int i = 0; emItr != emEnd; ++emItr, ++i) {
         //   cout << "#" << i
         //        << " name " << emItr->name()
         //        << " empty " << emItr->empty()
         //        << " rank " << emItr->rank()
         //        << " eta " << emItr->etaIndex()
         //        << " sign " << emItr->etaSign()
         //        << " phi " << emItr->phiIndex()
         //        << " iso " << emItr->isolated()
         //        << " bx " << emItr->bx()
         //        << endl ;
  
         if (!emItr->empty() && (!centralBxOnly_ || emItr->bx() == 0)) {
           double et = emScale->et(emItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           nonIsoEmColl->push_back(L1EmParticle(gctLorentzVector(et, *emItr, caloGeom, true),
                                                Ref<L1GctEmCandCollection>(hwNonIsoEmCands, i),
                                                emItr->bx()));
         }
       }
     }
  
     // ~~~~~~~~~~~~~~~~~~~~ Jets ~~~~~~~~~~~~~~~~~~~~
  
     ESHandle<L1CaloEtScale> jetScale = iSetup.getHandle(jetScaleToken_);
  
     // Central jets.
     Handle<L1GctJetCandCollection> hwCenJetCands;
     iEvent.getByLabel(cenJetSource_, hwCenJetCands);
  
     if (!hwCenJetCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctJetCandCollection with " << cenJetSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       //       cout << "HW central jets" << endl ;
       L1GctJetCandCollection::const_iterator jetItr = hwCenJetCands->begin();
       L1GctJetCandCollection::const_iterator jetEnd = hwCenJetCands->end();
       for (int i = 0; jetItr != jetEnd; ++jetItr, ++i) {
         //   cout << "#" << i
         //        << " name " << jetItr->name()
         //        << " empty " << jetItr->empty()
         //        << " rank " << jetItr->rank()
         //        << " eta " << jetItr->etaIndex()
         //        << " sign " << jetItr->etaSign()
         //        << " phi " << jetItr->phiIndex()
         //        << " cen " << jetItr->isCentral()
         //        << " for " << jetItr->isForward()
         //        << " tau " << jetItr->isTau()
         //        << " bx " << jetItr->bx()
         //        << endl ;
  
         if (!jetItr->empty() && (!centralBxOnly_ || jetItr->bx() == 0)) {
           double et = jetScale->et(jetItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           cenJetColl->push_back(L1JetParticle(gctLorentzVector(et, *jetItr, caloGeom, true),
                                               Ref<L1GctJetCandCollection>(hwCenJetCands, i),
                                               jetItr->bx()));
         }
       }
     }
  
     // Forward jets.
     Handle<L1GctJetCandCollection> hwForJetCands;
     iEvent.getByLabel(forJetSource_, hwForJetCands);
  
     if (!hwForJetCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEmCandCollection with " << forJetSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       //       cout << "HW forward jets" << endl ;
       L1GctJetCandCollection::const_iterator jetItr = hwForJetCands->begin();
       L1GctJetCandCollection::const_iterator jetEnd = hwForJetCands->end();
       for (int i = 0; jetItr != jetEnd; ++jetItr, ++i) {
         //   cout << "#" << i
         //        << " name " << jetItr->name()
         //        << " empty " << jetItr->empty()
         //        << " rank " << jetItr->rank()
         //        << " eta " << jetItr->etaIndex()
         //        << " sign " << jetItr->etaSign()
         //        << " phi " << jetItr->phiIndex()
         //        << " cen " << jetItr->isCentral()
         //        << " for " << jetItr->isForward()
         //        << " tau " << jetItr->isTau()
         //        << " bx " << jetItr->bx()
         //        << endl ;
  
         if (!jetItr->empty() && (!centralBxOnly_ || jetItr->bx() == 0)) {
           double et = jetScale->et(jetItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           forJetColl->push_back(L1JetParticle(gctLorentzVector(et, *jetItr, caloGeom, false),
                                               Ref<L1GctJetCandCollection>(hwForJetCands, i),
                                               jetItr->bx()));
         }
       }
     }
  
     // Tau jets.
     //       cout << "HW tau jets" << endl ;
     Handle<L1GctJetCandCollection> hwTauJetCands;
     iEvent.getByLabel(tauJetSource_, hwTauJetCands);
  
     if (!hwTauJetCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctJetCandCollection with " << tauJetSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       L1GctJetCandCollection::const_iterator jetItr = hwTauJetCands->begin();
       L1GctJetCandCollection::const_iterator jetEnd = hwTauJetCands->end();
       for (int i = 0; jetItr != jetEnd; ++jetItr, ++i) {
         //   cout << "#" << i
         //        << " name " << jetItr->name()
         //        << " empty " << jetItr->empty()
         //        << " rank " << jetItr->rank()
         //        << " eta " << jetItr->etaIndex()
         //        << " sign " << jetItr->etaSign()
         //        << " phi " << jetItr->phiIndex()
         //        << " cen " << jetItr->isCentral()
         //        << " for " << jetItr->isForward()
         //        << " tau " << jetItr->isTau()
         //        << " bx " << jetItr->bx()
         //        << endl ;
  
         if (!jetItr->empty() && (!centralBxOnly_ || jetItr->bx() == 0)) {
           double et = jetScale->et(jetItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           tauJetColl->push_back(L1JetParticle(gctLorentzVector(et, *jetItr, caloGeom, true),
                                               Ref<L1GctJetCandCollection>(hwTauJetCands, i),
                                               jetItr->bx()));
         }
       }
     }
  
     // Isolated Tau jets.
     //       cout << "HW tau jets" << endl ;
     Handle<L1GctJetCandCollection> hwIsoTauJetCands;
     iEvent.getByLabel(isoTauJetSource_, hwIsoTauJetCands);
  
     if (!hwIsoTauJetCands.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctJetCandCollection with " << isoTauJetSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       L1GctJetCandCollection::const_iterator jetItr = hwIsoTauJetCands->begin();
       L1GctJetCandCollection::const_iterator jetEnd = hwIsoTauJetCands->end();
       for (int i = 0; jetItr != jetEnd; ++jetItr, ++i) {
         //   cout << "#" << i
         //        << " name " << jetItr->name()
         //        << " empty " << jetItr->empty()
         //        << " rank " << jetItr->rank()
         //        << " eta " << jetItr->etaIndex()
         //        << " sign " << jetItr->etaSign()
         //        << " phi " << jetItr->phiIndex()
         //        << " cen " << jetItr->isCentral()
         //        << " for " << jetItr->isForward()
         //        << " tau " << jetItr->isTau()
         //        << " bx " << jetItr->bx()
         //        << endl ;
  
         if (!jetItr->empty() && (!centralBxOnly_ || jetItr->bx() == 0)) {
           double et = jetScale->et(jetItr->rank());
  
           //        cout << "L1Extra et " << et << endl ;
  
           isoTauJetColl->push_back(L1JetParticle(gctLorentzVector(et, *jetItr, caloGeom, true),
                                                  Ref<L1GctJetCandCollection>(hwIsoTauJetCands, i),
                                                  jetItr->bx()));
         }
       }
     }
  
     // ~~~~~~~~~~~~~~~~~~~~ ET Sums ~~~~~~~~~~~~~~~~~~~~
  
     double etSumLSB = jetScale->linearLsb();
  
     Handle<L1GctEtTotalCollection> hwEtTotColl;
     iEvent.getByLabel(etTotSource_, hwEtTotColl);
  
     Handle<L1GctEtMissCollection> hwEtMissColl;
     iEvent.getByLabel(etMissSource_, hwEtMissColl);
  
     if (!hwEtTotColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEtTotalCollection with " << etTotSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else if (!hwEtMissColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEtMissCollection with " << etMissSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       // Make a L1EtMissParticle even if either L1GctEtTotal or L1GctEtMiss
       // is missing for a given bx.  Keep track of which L1GctEtMiss objects
       // have a corresponding L1GctEtTotal object.
       std::vector<bool> etMissMatched;
  
       L1GctEtMissCollection::const_iterator hwEtMissItr = hwEtMissColl->begin();
       L1GctEtMissCollection::const_iterator hwEtMissEnd = hwEtMissColl->end();
       for (; hwEtMissItr != hwEtMissEnd; ++hwEtMissItr) {
         etMissMatched.push_back(false);
       }
  
       // Collate energy sums by bx
       L1GctEtTotalCollection::const_iterator hwEtTotItr = hwEtTotColl->begin();
       L1GctEtTotalCollection::const_iterator hwEtTotEnd = hwEtTotColl->end();
  
       int iTot = 0;
       for (; hwEtTotItr != hwEtTotEnd; ++hwEtTotItr, ++iTot) {
         int bx = hwEtTotItr->bx();
  
         if (!centralBxOnly_ || bx == 0) {
           // ET bin low edge
           double etTot =
               (hwEtTotItr->overFlow() ? (double)L1GctEtTotal::kEtTotalMaxValue : (double)hwEtTotItr->et()) * etSumLSB +
               1.e-6;
  
           int iMiss = 0;
           hwEtMissItr = hwEtMissColl->begin();
           hwEtMissEnd = hwEtMissColl->end();
           for (; hwEtMissItr != hwEtMissEnd; ++hwEtMissItr, ++iMiss) {
             if (hwEtMissItr->bx() == bx) {
               etMissMatched[iMiss] = true;
               break;
             }
           }
  
           double etMiss = 0.;
           double phi = 0.;
           math::PtEtaPhiMLorentzVector p4;
           Ref<L1GctEtMissCollection> metRef;
  
           // If a L1GctEtMiss with the right bx is not found, itr == end.
           if (hwEtMissItr != hwEtMissEnd) {
             // ET bin low edge
             etMiss = (hwEtMissItr->overFlow() ? (double)L1GctEtMiss::kEtMissMaxValue : (double)hwEtMissItr->et()) *
                          etSumLSB +
                      1.e-6;
             // keep x and y components non-zero and
             // protect against roundoff.
  
             phi = caloGeom->etSumPhiBinCenter(hwEtMissItr->phi());
  
             p4 = math::PtEtaPhiMLorentzVector(etMiss, 0., phi, 0.);
  
             metRef = Ref<L1GctEtMissCollection>(hwEtMissColl, iMiss);
  
             //         cout << "HW ET Sums " << endl
             //          << "MET: phi " << hwEtMissItr->phi() << " = " << phi
             //          << " et " << hwEtMissItr->et() << " = " << etMiss
             //          << " EtTot " << hwEtTotItr->et() << " = " << etTot
             //          << " bx " << bx
             //          << endl ;
           }
           //       else
           //         {
           //           cout << "HW ET Sums " << endl
           //            << "MET: phi " << phi
           //            << " et "<< etMiss
           //            << " EtTot " << hwEtTotItr->et() << " = " << etTot
           //            << " bx " << bx
           //            << endl ;
           //         }
  
           etMissColl->push_back(L1EtMissParticle(p4,
                                                  L1EtMissParticle::kMET,
                                                  etTot,
                                                  metRef,
                                                  Ref<L1GctEtTotalCollection>(hwEtTotColl, iTot),
                                                  Ref<L1GctHtMissCollection>(),
                                                  Ref<L1GctEtHadCollection>(),
                                                  bx));
         }
       }
  
       if (!centralBxOnly_) {
         // Make L1EtMissParticles for those L1GctEtMiss objects without
         // a matched L1GctEtTotal object.
  
         double etTot = 0.;
  
         hwEtMissItr = hwEtMissColl->begin();
         hwEtMissEnd = hwEtMissColl->end();
         int iMiss = 0;
         for (; hwEtMissItr != hwEtMissEnd; ++hwEtMissItr, ++iMiss) {
           if (!etMissMatched[iMiss]) {
             int bx = hwEtMissItr->bx();
  
             // ET bin low edge
             double etMiss =
                 (hwEtMissItr->overFlow() ? (double)L1GctEtMiss::kEtMissMaxValue : (double)hwEtMissItr->et()) *
                     etSumLSB +
                 1.e-6;
             // keep x and y components non-zero and
             // protect against roundoff.
  
             double phi = caloGeom->etSumPhiBinCenter(hwEtMissItr->phi());
  
             math::PtEtaPhiMLorentzVector p4(etMiss, 0., phi, 0.);
  
             //        cout << "HW ET Sums " << endl
             //             << "MET: phi " << hwEtMissItr->phi() << " = " <<
             // phi
             //             << " et " << hwEtMissItr->et() << " = " << etMiss
             //             << " EtTot " << etTot
             //             << " bx " << bx
             //             << endl ;
  
             etMissColl->push_back(L1EtMissParticle(p4,
                                                    L1EtMissParticle::kMET,
                                                    etTot,
                                                    Ref<L1GctEtMissCollection>(hwEtMissColl, iMiss),
                                                    Ref<L1GctEtTotalCollection>(),
                                                    Ref<L1GctHtMissCollection>(),
                                                    Ref<L1GctEtHadCollection>(),
                                                    bx));
           }
         }
       }
     }
  
     // ~~~~~~~~~~~~~~~~~~~~ HT Sums ~~~~~~~~~~~~~~~~~~~~
  
     Handle<L1GctEtHadCollection> hwEtHadColl;
     iEvent.getByLabel(etHadSource_, hwEtHadColl);
  
     Handle<L1GctHtMissCollection> hwHtMissColl;
     if (!ignoreHtMiss_) {
       iEvent.getByLabel(htMissSource_, hwHtMissColl);
     }
  
     ESHandle<L1GctJetFinderParams> jetFinderParams = iSetup.getHandle(jetFinderParamsToken_);
     double htSumLSB = jetFinderParams->getHtLsbGeV();
  
     ESHandle<L1CaloEtScale> htMissScale;
     std::vector<bool> htMissMatched;
     if (!ignoreHtMiss_) {
       htMissScale = iSetup.getHandle(htMissScaleToken_);
  
       if (!hwEtHadColl.isValid()) {
         LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEtHadCollection with " << etHadSource_
                                          << "\nrequested in configuration, but not found in the event." << std::endl;
       } else if (!hwHtMissColl.isValid()) {
         LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctHtMissCollection with " << htMissSource_
                                          << "\nrequested in configuration, but not found in the event." << std::endl;
       } else {
         // Make a L1EtMissParticle even if either L1GctEtHad or L1GctHtMiss
         // is missing for a given bx. Keep track of which L1GctHtMiss objects
         // have a corresponding L1GctHtTotal object.
         L1GctHtMissCollection::const_iterator hwHtMissItr = hwHtMissColl->begin();
         L1GctHtMissCollection::const_iterator hwHtMissEnd = hwHtMissColl->end();
         for (; hwHtMissItr != hwHtMissEnd; ++hwHtMissItr) {
           htMissMatched.push_back(false);
         }
       }
     }
  
     if (!hwEtHadColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctEtHadCollection with " << etHadSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else if (!hwHtMissColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctHtMissCollection with " << htMissSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       L1GctEtHadCollection::const_iterator hwEtHadItr = hwEtHadColl->begin();
       L1GctEtHadCollection::const_iterator hwEtHadEnd = hwEtHadColl->end();
  
       int iHad = 0;
       for (; hwEtHadItr != hwEtHadEnd; ++hwEtHadItr, ++iHad) {
         int bx = hwEtHadItr->bx();
  
         if (!centralBxOnly_ || bx == 0) {
           // HT bin low edge
           double htTot =
               (hwEtHadItr->overFlow() ? (double)L1GctEtHad::kEtHadMaxValue : (double)hwEtHadItr->et()) * htSumLSB +
               1.e-6;
  
           double htMiss = 0.;
           double phi = 0.;
           math::PtEtaPhiMLorentzVector p4;
           Ref<L1GctHtMissCollection> mhtRef;
  
           if (!ignoreHtMiss_) {
             L1GctHtMissCollection::const_iterator hwHtMissItr = hwHtMissColl->begin();
             L1GctHtMissCollection::const_iterator hwHtMissEnd = hwHtMissColl->end();
  
             int iMiss = 0;
             for (; hwHtMissItr != hwHtMissEnd; ++hwHtMissItr, ++iMiss) {
               if (hwHtMissItr->bx() == bx) {
                 htMissMatched[iMiss] = true;
                 break;
               }
             }
  
             // If a L1GctHtMiss with the right bx is not found, itr == end.
             if (hwHtMissItr != hwHtMissEnd) {
               // HT bin low edge
               htMiss =
                   htMissScale->et(hwHtMissItr->overFlow() ? htMissScale->rankScaleMax() : hwHtMissItr->et()) + 1.e-6;
               // keep x and y components non-zero and
               // protect against roundoff.
  
               phi = caloGeom->htSumPhiBinCenter(hwHtMissItr->phi());
  
               p4 = math::PtEtaPhiMLorentzVector(htMiss, 0., phi, 0.);
  
               mhtRef = Ref<L1GctHtMissCollection>(hwHtMissColl, iMiss);
  
               //           cout << "HW HT Sums " << endl
               //            << "MHT: phi " << hwHtMissItr->phi() <<
               // " = " << phi
               //            << " ht " << hwHtMissItr->et() << " = "
               // << htMiss
               //            << " HtTot " << hwEtHadItr->et() << " =
               // "
               // << htTot
               //            << " bx " << bx
               //            << endl ;
             }
             //         else
             //       {
             //         cout << "HW HT Sums " << endl
             //          << "MHT: phi " << phi
             //          << " ht " << htMiss
             //          << " HtTot " << hwEtHadItr->et() << " = " <<
             // htTot
             //          << " bx " << bx
             //          << endl ;
             //       }
           }
  
           htMissColl->push_back(L1EtMissParticle(p4,
                                                  L1EtMissParticle::kMHT,
                                                  htTot,
                                                  Ref<L1GctEtMissCollection>(),
                                                  Ref<L1GctEtTotalCollection>(),
                                                  mhtRef,
                                                  Ref<L1GctEtHadCollection>(hwEtHadColl, iHad),
                                                  bx));
         }
       }
  
       if (!centralBxOnly_ && !ignoreHtMiss_) {
         // Make L1EtMissParticles for those L1GctHtMiss objects without
         // a matched L1GctHtTotal object.
         double htTot = 0.;
  
         L1GctHtMissCollection::const_iterator hwHtMissItr = hwHtMissColl->begin();
         L1GctHtMissCollection::const_iterator hwHtMissEnd = hwHtMissColl->end();
  
         int iMiss = 0;
         for (; hwHtMissItr != hwHtMissEnd; ++hwHtMissItr, ++iMiss) {
           if (!htMissMatched[iMiss]) {
             int bx = hwHtMissItr->bx();
  
             // HT bin low edge
             double htMiss =
                 htMissScale->et(hwHtMissItr->overFlow() ? htMissScale->rankScaleMax() : hwHtMissItr->et()) + 1.e-6;
             // keep x and y components non-zero and
             // protect against roundoff.
  
             double phi = caloGeom->htSumPhiBinCenter(hwHtMissItr->phi());
  
             math::PtEtaPhiMLorentzVector p4(htMiss, 0., phi, 0.);
  
             //        cout << "HW HT Sums " << endl
             //             << "MHT: phi " << hwHtMissItr->phi() << " = " <<
             // phi
             //             << " ht " << hwHtMissItr->et() << " = " << htMiss
             //             << " HtTot " << htTot
             //             << " bx " << bx
             //             << endl ;
  
             htMissColl->push_back(L1EtMissParticle(p4,
                                                    L1EtMissParticle::kMHT,
                                                    htTot,
                                                    Ref<L1GctEtMissCollection>(),
                                                    Ref<L1GctEtTotalCollection>(),
                                                    Ref<L1GctHtMissCollection>(hwHtMissColl, iMiss),
                                                    Ref<L1GctEtHadCollection>(),
                                                    bx));
           }
         }
       }
     }
  
     // ~~~~~~~~~~~~~~~~~~~~ HF Rings ~~~~~~~~~~~~~~~~~~~~
  
     Handle<L1GctHFRingEtSumsCollection> hwHFEtSumsColl;
     iEvent.getByLabel(hfRingEtSumsSource_, hwHFEtSumsColl);
  
     Handle<L1GctHFBitCountsCollection> hwHFBitCountsColl;
     iEvent.getByLabel(hfRingBitCountsSource_, hwHFBitCountsColl);
  
     ESHandle<L1CaloEtScale> hfRingEtScale = iSetup.getHandle(hfRingEtScaleToken_);
  
     if (!hwHFEtSumsColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctHFRingEtSumsCollection with " << hfRingEtSumsSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else if (!hwHFBitCountsColl.isValid()) {
       LogDebug("L1ExtraParticlesProd") << "\nWarning: L1GctHFBitCountsCollection with " << hfRingBitCountsSource_
                                        << "\nrequested in configuration, but not found in the event." << std::endl;
     } else {
       L1GctHFRingEtSumsCollection::const_iterator hwHFEtSumsItr = hwHFEtSumsColl->begin();
       L1GctHFRingEtSumsCollection::const_iterator hwHFEtSumsEnd = hwHFEtSumsColl->end();
  
       int iEtSums = 0;
       for (; hwHFEtSumsItr != hwHFEtSumsEnd; ++hwHFEtSumsItr, ++iEtSums) {
         int bx = hwHFEtSumsItr->bx();
  
         if (!centralBxOnly_ || bx == 0) {
           L1GctHFBitCountsCollection::const_iterator hwHFBitCountsItr = hwHFBitCountsColl->begin();
           L1GctHFBitCountsCollection::const_iterator hwHFBitCountsEnd = hwHFBitCountsColl->end();
  
           int iBitCounts = 0;
           for (; hwHFBitCountsItr != hwHFBitCountsEnd; ++hwHFBitCountsItr, ++iBitCounts) {
             if (hwHFBitCountsItr->bx() == bx) {
               break;
             }
           }
  
           // If a L1GctHFBitCounts with the right bx is not found, itr == end.
           if (hwHFBitCountsItr != hwHFBitCountsEnd) {
             // Construct L1HFRings only if both HW objects are present.
  
             // cout << "HF Rings " << endl
  
             // HF Et sums
             double etSums[L1HFRings::kNumRings];
             etSums[L1HFRings::kRing1PosEta] = hfRingEtScale->et(hwHFEtSumsItr->etSum(0)) + 1.e-6;
             etSums[L1HFRings::kRing1NegEta] = hfRingEtScale->et(hwHFEtSumsItr->etSum(1)) + 1.e-6;
             etSums[L1HFRings::kRing2PosEta] = hfRingEtScale->et(hwHFEtSumsItr->etSum(2)) + 1.e-6;
             etSums[L1HFRings::kRing2NegEta] = hfRingEtScale->et(hwHFEtSumsItr->etSum(3)) + 1.e-6;
             // protect against roundoff.
  
             //         cout << "HF Et Sums "
             //          << hwHFEtSumsItr->etSum( 0 ) << " = "
             //          << etSums[ L1HFRings::kRing1PosEta ] << " "
             //          << hwHFEtSumsItr->etSum( 1 ) << " = "
             //          << etSums[ L1HFRings::kRing1NegEta ] << " "
             //          << hwHFEtSumsItr->etSum( 2 ) << " = "
             //          << etSums[ L1HFRings::kRing2PosEta ] << " "
             //          << hwHFEtSumsItr->etSum( 3 ) << " = "
             //          << etSums[ L1HFRings::kRing2NegEta ]
             //          << endl ;
  
             // HF bit counts
             int bitCounts[L1HFRings::kNumRings];
             bitCounts[L1HFRings::kRing1PosEta] = hwHFBitCountsItr->bitCount(0);
             bitCounts[L1HFRings::kRing1NegEta] = hwHFBitCountsItr->bitCount(1);
             bitCounts[L1HFRings::kRing2PosEta] = hwHFBitCountsItr->bitCount(2);
             bitCounts[L1HFRings::kRing2NegEta] = hwHFBitCountsItr->bitCount(3);
  
             //         cout << "HF bit counts "
             //          << hwHFBitCountsItr->bitCount( 0 ) << " "
             //          << hwHFBitCountsItr->bitCount( 1 ) << " "
             //          << hwHFBitCountsItr->bitCount( 2 ) << " "
             //          << hwHFBitCountsItr->bitCount( 3 ) << " "
             //          << endl ;
  
             hfRingsColl->push_back(L1HFRings(etSums,
                                              bitCounts,
                                              Ref<L1GctHFRingEtSumsCollection>(hwHFEtSumsColl, iEtSums),
                                              Ref<L1GctHFBitCountsCollection>(hwHFBitCountsColl, iBitCounts),
                                              bx));
           }
         }
       }
     }
   }
  
   OrphanHandle<L1EmParticleCollection> isoEmHandle = iEvent.put(std::move(isoEmColl), "Isolated");
  
   OrphanHandle<L1EmParticleCollection> nonIsoEmHandle = iEvent.put(std::move(nonIsoEmColl), "NonIsolated");
  
   OrphanHandle<L1JetParticleCollection> cenJetHandle = iEvent.put(std::move(cenJetColl), "Central");
  
   OrphanHandle<L1JetParticleCollection> forJetHandle = iEvent.put(std::move(forJetColl), "Forward");
  
   OrphanHandle<L1JetParticleCollection> tauJetHandle = iEvent.put(std::move(tauJetColl), "Tau");
  
   OrphanHandle<L1JetParticleCollection> IsoTauJetHandle = iEvent.put(std::move(isoTauJetColl), "IsoTau");
  
   OrphanHandle<L1EtMissParticleCollection> etMissCollHandle = iEvent.put(std::move(etMissColl), "MET");
  
   OrphanHandle<L1EtMissParticleCollection> htMissCollHandle = iEvent.put(std::move(htMissColl), "MHT");
  
   OrphanHandle<L1HFRingsCollection> hfRingsCollHandle = iEvent.put(std::move(hfRingsColl));
 }