L1HLTTauMatching Class Reference

#include <L1HLTTauMatching.h>

Inheritance diagram for L1HLTTauMatching:

Public Member Functions
	L1HLTTauMatching (const edm::ParameterSet &)
virtual void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
	~L1HLTTauMatching ()
Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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

Private Attributes
const edm::EDGetTokenT< reco::PFTauCollection >	jetSrc
const double	mEt_Min
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs >	tauTrigger

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

Detailed Description

Definition at line 22 of file L1HLTTauMatching.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 of file L1HLTTauMatching.cc.

                                                                :
  jetSrc( consumes<PFTauCollection>(iConfig.getParameter<InputTag>("JetSrc") ) ),
  tauTrigger( consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<InputTag>("L1TauTrigger") ) ),
  mEt_Min( iConfig.getParameter<double>("EtMin") )
{  
  produces<PFTauCollection>();
}

L1HLTTauMatching::~L1HLTTauMatching

(

)

Definition at line 24 of file L1HLTTauMatching.cc.

{ }

Member Function Documentation

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

static

Definition at line 107 of file L1HLTTauMatching.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::setComment(), and edm::ParameterDescriptionNode::setComment().

{
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("L1TauTrigger",edm::InputTag("hltL1sDoubleIsoTau40er"))->setComment("Name of trigger filter");
  desc.add<edm::InputTag>("JetSrc",edm::InputTag("hltSelectedPFTausTrackPt1MediumIsolationReg"))->setComment("Input collection of PFTaus");
  desc.add<double>("EtMin",0.0)->setComment("Minimal pT of PFTau to match");
  descriptions.setComment("This module produces collection of PFTaus matched to L1ExtraTaus/Jets passing a HLT filter (Only p4 and vertex of returned PFTaus are set).");
  descriptions.add("L1HLTJetsMatching",desc);
}

void L1HLTTauMatching::produce ( edm::StreamID  iSId, edm::Event &  iEvent, const edm::EventSetup &  iES  ) const

overridevirtual

Definition at line 26 of file L1HLTTauMatching.cc.

References a, deltaR(), HiRegitMuonDetachedTripletStep_cff::DeltaR, edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), edm::Ptr< T >::isNonnull(), jetSrc, reco::PFTau::leadPFChargedHadrCand(), mEt_Min, eostools::move(), reco::LeafCandidate::p4(), reco::LeafCandidate::pt(), edm::Event::put(), tauTrigger, trigger::TriggerL1CenJet, trigger::TriggerL1TauJet, and reco::PFCandidate::vertex().

{
    
  using namespace edm;
  using namespace std;
  using namespace reco;
  using namespace trigger;
  using namespace l1extra;
    
  unique_ptr<PFTauCollection> tauL2jets(new PFTauCollection);
    
  double deltaR = 1.0;
  double matchingR = 0.5;
  //Getting HLT jets to be matched
  edm::Handle<PFTauCollection > tauJets;
  iEvent.getByToken( jetSrc, tauJets );

  //        std::cout <<"Size of input jet collection "<<tauJets->size()<<std::endl;
        
  edm::Handle<trigger::TriggerFilterObjectWithRefs> l1TriggeredTaus;
  iEvent.getByToken(tauTrigger,l1TriggeredTaus);
                
  vector<l1extra::L1JetParticleRef> tauCandRefVec;                                                
  vector<l1extra::L1JetParticleRef> jetCandRefVec;
  l1TriggeredTaus->getObjects( trigger::TriggerL1TauJet,tauCandRefVec);
  l1TriggeredTaus->getObjects( trigger::TriggerL1CenJet,jetCandRefVec);

  math::XYZPoint a(0.,0.,0.);
        
  for( unsigned int iL1Tau=0; iL1Tau <tauCandRefVec.size();iL1Tau++)
    {  
      for(unsigned int iJet=0;iJet<tauJets->size();iJet++)
    {
      //Find the relative L2TauJets, to see if it has been reconstructed
      const PFTau &  myJet = (*tauJets)[iJet];
      deltaR = ROOT::Math::VectorUtil::DeltaR(myJet.p4().Vect(), (tauCandRefVec[iL1Tau]->p4()).Vect());
      if(deltaR < matchingR ) {
        //       LeafCandidate myLC(myJet);
        if(myJet.leadPFChargedHadrCand().isNonnull()){
          a =  myJet.leadPFChargedHadrCand()->vertex();  
        }
        PFTau myPFTau(std::numeric_limits<int>::quiet_NaN(), myJet.p4(), a);
        if(myJet.pt() > mEt_Min) {
          //          tauL2LC->push_back(myLC);
          tauL2jets->push_back(myPFTau);
        }
        break;
      }
    }
    }  
  
  for(unsigned int iL1Tau=0; iL1Tau <jetCandRefVec.size();iL1Tau++)
    {  
      for(unsigned int iJet=0;iJet<tauJets->size();iJet++)
    {
      const PFTau &  myJet = (*tauJets)[iJet];
      //Find the relative L2TauJets, to see if it has been reconstructed
      deltaR = ROOT::Math::VectorUtil::DeltaR(myJet.p4().Vect(), (jetCandRefVec[iL1Tau]->p4()).Vect());
      if(deltaR < matchingR ) {
        //       LeafCandidate myLC(myJet);
        if(myJet.leadPFChargedHadrCand().isNonnull()){
          a =  myJet.leadPFChargedHadrCand()->vertex();  
        }
            
        PFTau myPFTau(std::numeric_limits<int>::quiet_NaN(), myJet.p4(),a);
        if(myJet.pt() > mEt_Min) {
          //          tauL2LC->push_back(myLC);
          tauL2jets->push_back(myPFTau);
        }
        break;
      }
    }
    }
  

  //std::cout <<"Size of L1HLT matched jets "<<tauL2jets->size()<<std::endl;

  iEvent.put(std::move(tauL2jets));
  // iEvent.put(std::move(tauL2LC));
}

Member Data Documentation

const edm::EDGetTokenT<reco::PFTauCollection> L1HLTTauMatching::jetSrc

private

Definition at line 31 of file L1HLTTauMatching.h.

Referenced by produce().

const double L1HLTTauMatching::mEt_Min

private

Definition at line 33 of file L1HLTTauMatching.h.

Referenced by produce().

const edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> L1HLTTauMatching::tauTrigger

private

Definition at line 32 of file L1HLTTauMatching.h.

Referenced by produce().