l1t::PhysicalEtAdder Class Reference

#include <PhysicalEtAdder.h>

Inheritance diagram for l1t::PhysicalEtAdder:

Public Member Functions
	PhysicalEtAdder (const edm::ParameterSet &ps)
	~PhysicalEtAdder ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	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
	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 (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) 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	endJob () override
virtual void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::EDGetToken	EGammaToken_
edm::EDGetToken	EtSumToken_
edm::EDGetToken	JetToken_
edm::EDGetToken	TauToken_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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

Definition at line 39 of file PhysicalEtAdder.h.

Constructor & Destructor Documentation

l1t::PhysicalEtAdder::PhysicalEtAdder ( const edm::ParameterSet &  ps )

explicit

Definition at line 24 of file PhysicalEtAdder.cc.

References EGammaToken_, EtSumToken_, edm::ParameterSet::getParameter(), JetToken_, and TauToken_.

                                                             {

  produces<l1t::EGammaBxCollection>();
  produces<l1t::TauBxCollection>();
  produces<l1t::JetBxCollection>();
  produces<l1t::EtSumBxCollection>();

  EGammaToken_ = consumes<l1t::EGammaBxCollection>(ps.getParameter<edm::InputTag>("InputCollection"));
  TauToken_ = consumes<l1t::TauBxCollection>(ps.getParameter<edm::InputTag>("InputCollection"));
  JetToken_ = consumes<l1t::JetBxCollection>(ps.getParameter<edm::InputTag>("InputCollection"));
  EtSumToken_ = consumes<l1t::EtSumBxCollection>(ps.getParameter<edm::InputTag>("InputCollection"));
}

l1t::PhysicalEtAdder::~PhysicalEtAdder

(

)

Definition at line 37 of file PhysicalEtAdder.cc.

                                     {

}

Member Function Documentation

void l1t::PhysicalEtAdder::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 171 of file PhysicalEtAdder.cc.

{
}

void l1t::PhysicalEtAdder::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 177 of file PhysicalEtAdder.cc.

                           {
}

void l1t::PhysicalEtAdder::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 216 of file PhysicalEtAdder.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 l1t::PhysicalEtAdder::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 43 of file PhysicalEtAdder.cc.

References eta(), edm::EventSetup::get(), edm::Event::getByToken(), getPhysicalEta(), getPhysicalPhi(), l1t::L1Candidate::hwPt(), metsig::jet, p4, phi, RecoTauCleanerPlugins::pt, edm::Event::put(), and metsig::tau.

{
  // store new collections which include physical quantities
  std::auto_ptr<l1t::EGammaBxCollection> new_egammas (new l1t::EGammaBxCollection);
  std::auto_ptr<l1t::TauBxCollection> new_taus (new l1t::TauBxCollection);
  std::auto_ptr<l1t::JetBxCollection> new_jets (new l1t::JetBxCollection);
  std::auto_ptr<l1t::EtSumBxCollection> new_etsums (new l1t::EtSumBxCollection);

  edm::Handle<l1t::EGammaBxCollection> old_egammas;
  edm::Handle<l1t::TauBxCollection> old_taus;
  edm::Handle<l1t::JetBxCollection> old_jets;
  edm::Handle<l1t::EtSumBxCollection> old_etsums;

  iEvent.getByToken(EGammaToken_, old_egammas);
  iEvent.getByToken(TauToken_, old_taus);
  iEvent.getByToken(JetToken_, old_jets);
  iEvent.getByToken(EtSumToken_, old_etsums);


  //get the proper scales for conversion to physical et
  edm::ESHandle< L1CaloEtScale > emScale ;
  iSetup.get< L1EmEtScaleRcd >().get( emScale ) ;

  edm::ESHandle< L1CaloEtScale > jetScale ;
  iSetup.get< L1JetEtScaleRcd >().get( jetScale ) ;

  edm::ESHandle< L1CaloEtScale > hwForJetScale ;
  iSetup.get< L1JetEtScaleRcd >().get( hwForJetScale ) ;

  edm::ESHandle< L1CaloEtScale > htMissScale ;
  std::vector< bool > htMissMatched ;
  iSetup.get< L1HtMissScaleRcd >().get( htMissScale ) ;

  int firstBX = old_egammas->getFirstBX();
  int lastBX = old_egammas->getLastBX();

  new_egammas->setBXRange(firstBX, lastBX);
  new_taus->setBXRange(firstBX, lastBX);
  new_jets->setBXRange(firstBX, lastBX);
  new_taus->setBXRange(firstBX, lastBX);

  for(int bx = firstBX; bx <= lastBX; ++bx)
  {
    for(l1t::EGammaBxCollection::const_iterator itEGamma = old_egammas->begin(bx);
    itEGamma != old_egammas->end(bx); ++itEGamma)
    {
      const double pt = itEGamma->hwPt() * emScale->linearLsb();
      const double eta = getPhysicalEta(itEGamma->hwEta());
      const double phi = getPhysicalPhi(itEGamma->hwPhi());
      math::PtEtaPhiMLorentzVector *p4 = new math::PtEtaPhiMLorentzVector(pt, eta, phi, 0);

      l1t::EGamma *eg = new l1t::EGamma(*p4, itEGamma->hwPt(),
                       itEGamma->hwEta(), itEGamma->hwPhi(),
                       itEGamma->hwQual(), itEGamma->hwIso());
      new_egammas->push_back(bx, *eg);


    }

    for(l1t::TauBxCollection::const_iterator itTau = old_taus->begin(bx);
    itTau != old_taus->end(bx); ++itTau)
    {
      // use the full-circle conversion to match l1extra, accounts for linearLsb and max value automatically
      //const uint16_t rankPt = jetScale->rank((uint16_t)itTau->hwPt());
      //const double et = jetScale->et( rankPt ) ;

      // or use the emScale to get finer-grained et
      const double et = itTau->hwPt() * emScale->linearLsb();
      
      const double eta = getPhysicalEta(itTau->hwEta());
      const double phi = getPhysicalPhi(itTau->hwPhi());
      math::PtEtaPhiMLorentzVector *p4 = new math::PtEtaPhiMLorentzVector(et, eta, phi, 0);

      l1t::Tau *tau = new l1t::Tau(*p4, itTau->hwPt(),
                   itTau->hwEta(), itTau->hwPhi(),
                   itTau->hwQual());
      new_taus->push_back(bx, *tau);

    }

    for(l1t::JetBxCollection::const_iterator itJet = old_jets->begin(bx);
    itJet != old_jets->end(bx); ++itJet)
    {
      // use the full-circle conversion to match l1extra, accounts for linearLsb and max value automatically
      //const uint16_t rankPt = jetScale->rank((uint16_t)itJet->hwPt());
      //const double et = jetScale->et( rankPt ) ;
      
      // or use the emScale to get finer-grained et
      const double et = itJet->hwPt() * emScale->linearLsb();
      
      const double eta = getPhysicalEta(itJet->hwEta());
      const double phi = getPhysicalPhi(itJet->hwPhi());
      math::PtEtaPhiMLorentzVector *p4 = new math::PtEtaPhiMLorentzVector(et, eta, phi, 0);

      l1t::Jet *jet = new l1t::Jet(*p4, itJet->hwPt(),
                   itJet->hwEta(), itJet->hwPhi(),
                   itJet->hwQual());
      new_jets->push_back(bx, *jet);

    }

    for(l1t::EtSumBxCollection::const_iterator itEtSum = old_etsums->begin(bx);
    itEtSum != old_etsums->end(bx); ++itEtSum)
    {
      const double pt = itEtSum->hwPt() * emScale->linearLsb();
      const double eta = getPhysicalEta(itEtSum->hwEta());
      const double phi = getPhysicalPhi(itEtSum->hwPhi());
      const l1t::EtSum::EtSumType sumType = itEtSum->getType();
      //std::cout << "SumType: " << sumType << std::endl;
      math::PtEtaPhiMLorentzVector *p4 = new math::PtEtaPhiMLorentzVector(pt, eta, phi, 0);

      l1t::EtSum *eg = new l1t::EtSum(*p4, sumType, itEtSum->hwPt(),
                      itEtSum->hwEta(), itEtSum->hwPhi(),
                      itEtSum->hwQual());
      new_etsums->push_back(bx, *eg);

    }
  }

  iEvent.put(new_egammas);
  iEvent.put(new_taus);
  iEvent.put(new_jets);
  iEvent.put(new_etsums);

}

Member Data Documentation

edm::EDGetToken l1t::PhysicalEtAdder::EGammaToken_

private

Definition at line 58 of file PhysicalEtAdder.h.

Referenced by PhysicalEtAdder().

edm::EDGetToken l1t::PhysicalEtAdder::EtSumToken_

private

Definition at line 61 of file PhysicalEtAdder.h.

Referenced by PhysicalEtAdder().

edm::EDGetToken l1t::PhysicalEtAdder::JetToken_

private

Definition at line 60 of file PhysicalEtAdder.h.

Referenced by PhysicalEtAdder().

edm::EDGetToken l1t::PhysicalEtAdder::TauToken_

private

Definition at line 59 of file PhysicalEtAdder.h.

Referenced by PhysicalEtAdder().