#include <L1ExtraParticlesProd.h>

Inheritance diagram for L1ExtraParticlesProd:

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

	~L1ExtraParticlesProd ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

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

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

Static Private Attributes
static double	muonMassGeV_ = 0.105658369

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

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 42 of file L1ExtraParticlesProd.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 68 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" ) ),
      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< L1MuonParticleCollection >() ;
    produces< L1EtMissParticleCollection >( "MET" ) ;
    produces< L1EtMissParticleCollection >( "MHT" ) ;
    produces< L1HFRingsCollection >() ;
 
    //now do what ever other initialization is needed
 }

L1ExtraParticlesProd::~L1ExtraParticlesProd ( )

Definition at line 117 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

void L1ExtraParticlesProd::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 1092 of file L1ExtraParticlesProd.cc.

1093 {

1094 }

void L1ExtraParticlesProd::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 1098 of file L1ExtraParticlesProd.cc.

1098 {

1099 }

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

private

Definition at line 1064 of file L1ExtraParticlesProd.cc.

References L1CaloGeometry::emJetPhiBinCenter(), eta(), L1CaloGeometry::etaBinCenter(), L1GctCand::etaIndex(), phi, and L1GctCand::phiIndex().

Referenced by produce().

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

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

privatevirtual

Implements edm::EDProducer.

Definition at line 132 of file L1ExtraParticlesProd.cc.

References cenJetSource_, centralBxOnly_, eta(), etHadSource_, etMissSource_, L1CaloGeometry::etSumPhiBinCenter(), etTotSource_, forJetSource_, gctLorentzVector(), edm::EventSetup::get(), edm::Event::getByLabel(), hfRingBitCountsSource_, hfRingEtSumsSource_, htMissSource_, L1CaloGeometry::htSumPhiBinCenter(), i, ignoreHtMiss_, isoEmSource_, edm::HandleBase::isValid(), L1GctEtHad::kEtHadMaxValue, L1GctEtMiss::kEtMissMaxValue, L1GctEtTotal::kEtTotalMaxValue, LogDebug, muonMassGeV_, muonSource_, nonIsoEmSource_, p4, phi, produceCaloParticles_, produceMuonParticles_, edm::Event::put(), tauJetSource_, and funct::true.

 {
    using namespace edm ;
    using namespace l1extra ;
    using namespace std ;
 
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    // ~~~~~~~~~~~~~~~~~~~~ Muons ~~~~~~~~~~~~~~~~~~~~
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
    auto_ptr< L1MuonParticleCollection > muColl(
      new L1MuonParticleCollection );
 
    if( produceMuonParticles_ )
    {
       ESHandle< L1MuTriggerScales > muScales ;
       iSetup.get< L1MuTriggerScalesRcd >().get( muScales ) ;
 
       ESHandle< L1MuTriggerPtScale > muPtScale ;
       iSetup.get< L1MuTriggerPtScaleRcd >().get( muPtScale ) ;
 
       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( muColl );
 
 
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    // ~~~~~~~~~~~~~~~~~~~~ Calorimeter ~~~~~~~~~~~~~~~~~~~~
    // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 
    auto_ptr< L1EmParticleCollection > isoEmColl(
       new L1EmParticleCollection );
 
    auto_ptr< L1EmParticleCollection > nonIsoEmColl(
       new L1EmParticleCollection );
 
    auto_ptr< L1JetParticleCollection > cenJetColl(
       new L1JetParticleCollection );
 
    auto_ptr< L1JetParticleCollection > forJetColl(
       new L1JetParticleCollection );
 
    auto_ptr< L1JetParticleCollection > tauJetColl(
       new L1JetParticleCollection );
 
    auto_ptr< L1EtMissParticleCollection > etMissColl(
       new L1EtMissParticleCollection );
 
    auto_ptr< L1EtMissParticleCollection > htMissColl(
       new L1EtMissParticleCollection );
 
    auto_ptr< L1HFRingsCollection > hfRingsColl( new L1HFRingsCollection );
 
    if( produceCaloParticles_ )
    {
       // ~~~~~~~~~~~~~~~~~~~~ Geometry ~~~~~~~~~~~~~~~~~~~~
 
       ESHandle< L1CaloGeometry > caloGeomESH ;
       iSetup.get< L1CaloGeometryRecord >().get( caloGeomESH ) ;
       const L1CaloGeometry* caloGeom = &( *caloGeomESH ) ;
 
       // ~~~~~~~~~~~~~~~~~~~~ EM ~~~~~~~~~~~~~~~~~~~~
 
       ESHandle< L1CaloEtScale > emScale ;
       iSetup.get< L1EmEtScaleRcd >().get( emScale ) ;
 
       // 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.get< L1JetEtScaleRcd >().get( jetScale ) ;
 
       // 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() ) ) ;
         }
         }
     }
 
       // ~~~~~~~~~~~~~~~~~~~~ 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.get< L1GctJetFinderParamsRcd >().get( jetFinderParams ) ;
       double htSumLSB = jetFinderParams->getHtLsbGeV();
 
       ESHandle< L1CaloEtScale > htMissScale ;
       std::vector< bool > htMissMatched ;
       if( !ignoreHtMiss_ )
     {
       iSetup.get< L1HtMissScaleRcd >().get( htMissScale ) ;
 
       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.get< L1HfRingEtScaleRcd >().get( hfRingEtScale ) ;
 
       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( isoEmColl, "Isolated" ) ;
 
    OrphanHandle< L1EmParticleCollection > nonIsoEmHandle =
      iEvent.put( nonIsoEmColl, "NonIsolated" ) ;
 
    OrphanHandle< L1JetParticleCollection > cenJetHandle =
      iEvent.put( cenJetColl, "Central" ) ;
 
    OrphanHandle< L1JetParticleCollection > forJetHandle =
      iEvent.put( forJetColl, "Forward" ) ;
 
    OrphanHandle< L1JetParticleCollection > tauJetHandle =
      iEvent.put( tauJetColl, "Tau" ) ;
 
    OrphanHandle< L1EtMissParticleCollection > etMissCollHandle =
      iEvent.put( etMissColl, "MET" ) ;
 
    OrphanHandle< L1EtMissParticleCollection > htMissCollHandle =
      iEvent.put( htMissColl, "MHT" ) ;
 
    OrphanHandle< L1HFRingsCollection > hfRingsCollHandle =
      iEvent.put( hfRingsColl ) ;
 }

Member Data Documentation

edm::InputTag L1ExtraParticlesProd::cenJetSource_

private

Definition at line 65 of file L1ExtraParticlesProd.h.

Referenced by produce().

bool L1ExtraParticlesProd::centralBxOnly_

private

Definition at line 77 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::etHadSource_

private

Definition at line 69 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::etMissSource_

private

Definition at line 70 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::etTotSource_

private

Definition at line 68 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::forJetSource_

private

Definition at line 66 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::hfRingBitCountsSource_

private

Definition at line 73 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::hfRingEtSumsSource_

private

Definition at line 72 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::htMissSource_

private

Definition at line 71 of file L1ExtraParticlesProd.h.

Referenced by produce().

bool L1ExtraParticlesProd::ignoreHtMiss_

private

Definition at line 81 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::isoEmSource_

private

Definition at line 63 of file L1ExtraParticlesProd.h.

Referenced by produce().

double L1ExtraParticlesProd::muonMassGeV_ = 0.105658369

staticprivate

Definition at line 75 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::muonSource_

private

Definition at line 60 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::nonIsoEmSource_

private

Definition at line 64 of file L1ExtraParticlesProd.h.

Referenced by produce().

bool L1ExtraParticlesProd::produceCaloParticles_

private

Definition at line 62 of file L1ExtraParticlesProd.h.

Referenced by produce().

bool L1ExtraParticlesProd::produceMuonParticles_

private

Definition at line 59 of file L1ExtraParticlesProd.h.

Referenced by produce().

edm::InputTag L1ExtraParticlesProd::tauJetSource_

private

Definition at line 67 of file L1ExtraParticlesProd.h.

Referenced by produce().

Private Attributes
edm::InputTag	cenJetSource_

bool	centralBxOnly_

edm::InputTag	etHadSource_

edm::InputTag	etMissSource_

edm::InputTag	etTotSource_

edm::InputTag	forJetSource_

edm::InputTag	hfRingBitCountsSource_

edm::InputTag	hfRingEtSumsSource_

edm::InputTag	htMissSource_

bool	ignoreHtMiss_

edm::InputTag	isoEmSource_

edm::InputTag	muonSource_

edm::InputTag	nonIsoEmSource_

bool	produceCaloParticles_

bool	produceMuonParticles_

edm::InputTag	tauJetSource_

Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation