PFTauL1TJetsMatching.cc

Go to the documentation of this file.

#include "RecoTauTag/HLTProducers/interface/PFTauL1TJetsMatching.h"
#include "Math/GenVector/VectorUtil.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "DataFormats/TauReco/interface/PFTau.h"
#include "DataFormats/Math/interface/deltaR.h"

//
// class declaration
//
PFTauL1TJetsMatching::PFTauL1TJetsMatching(const edm::ParameterSet& iConfig):
    tauSrc_    ( consumes<reco::PFTauCollection>(iConfig.getParameter<edm::InputTag>("TauSrc"      ) ) ),
    L1JetSrc_  ( consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<edm::InputTag>("L1JetSrc") ) ),
    matchingR2_ ( iConfig.getParameter<double>("MatchingdR")*iConfig.getParameter<double>("MatchingdR") ),
    minTauPt_ (iConfig.getParameter<double>("MinTauPt") ),
    minL1TPt_ (iConfig.getParameter<double>("MinL1TPt") )
{  
    produces<reco::PFTauCollection>();
}
PFTauL1TJetsMatching::~PFTauL1TJetsMatching(){ }

void PFTauL1TJetsMatching::produce(edm::StreamID iSId, edm::Event& iEvent, const edm::EventSetup& iES) const
{
    std::unique_ptr<reco::PFTauCollection> L1TmatchedPFTau(new reco::PFTauCollection);
    
    edm::Handle<reco::PFTauCollection> taus;
    iEvent.getByToken(tauSrc_, taus);
  
    edm::Handle<trigger::TriggerFilterObjectWithRefs> L1Jets;
    iEvent.getByToken(L1JetSrc_,L1Jets);
                  
    l1t::JetVectorRef jetCandRefVec;
    L1Jets->getObjects(trigger::TriggerL1Jet,jetCandRefVec);

    /* Loop over taus that must pass a certain minTauPt_ cut */ 
    /* then loop over L1T jets that must pass minL1TPt_ */
    /* and check whether they match, if yes -> include the taus in */ 
    /* the new L1T matched PFTau collection */

    for(unsigned int iTau = 0; iTau < taus->size(); iTau++){  
        bool isMatched = false;
        if ((*taus)[iTau].pt() > minTauPt_){
            for (unsigned int iJet = 0; iJet < jetCandRefVec.size(); iJet++) {
                if (jetCandRefVec[iJet]->pt() > minL1TPt_){

                    if(reco::deltaR2((*taus)[iTau].p4(), jetCandRefVec[iJet]->p4()) < matchingR2_){
                        isMatched = true;
                        break;
                    }
                }
            }
        }
        if(isMatched) L1TmatchedPFTau->push_back((*taus)[iTau]);
    }
    iEvent.put(std::move(L1TmatchedPFTau));
}

void PFTauL1TJetsMatching::fillDescriptions(edm::ConfigurationDescriptions& descriptions) 
{
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag>("L1JetSrc", edm::InputTag("hltL1VBFDiJetOR"))->setComment("Input filter objects passing L1 seed"    );
    desc.add<edm::InputTag>("TauSrc", edm::InputTag("hltSelectedPFTausTrackFindingLooseChargedIsolationAgainstMuon"))->setComment("Input collection of PFTaus");
    desc.add<double>       ("MatchingdR",0.5)->setComment("Maximum dR for matching between PFTaus and L1 filter jets");
    desc.add<double>       ("MinTauPt",20.0)->setComment("PFTaus above this pt will be considered");
    desc.add<double>       ("MinL1TPt",115.0)->setComment("L1T Objects above this pt will be considered");
    descriptions.setComment("This module produces a collection of PFTaus matched to the leading jet passing the L1 seed filter.");
    descriptions.add       ("PFTauL1TJetsMatching",desc);
}