L2TauJetsMerger.cc

Go to the documentation of this file.

#include "RecoTauTag/HLTProducers/interface/L2TauJetsMerger.h"
#include "Math/GenVector/VectorUtil.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "FWCore/Utilities/interface/EDMException.h"
//
// class decleration
//
using namespace reco;
using namespace std;
using namespace edm;

L2TauJetsMerger::L2TauJetsMerger(const edm::ParameterSet& iConfig):
  jetSrc( iConfig.getParameter<vtag>("JetSrc") ),
  mEt_Min( iConfig.getParameter<double>("EtMin") )
{
   for(vtag::const_iterator it = jetSrc.begin(); it != jetSrc.end(); ++it) {
    edm::EDGetTokenT<CaloJetCollection> aToken = consumes<CaloJetCollection>(*it);
    jetSrc_token.push_back(aToken);
  }
  
  produces<CaloJetCollection>();
}

L2TauJetsMerger::~L2TauJetsMerger(){ }

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

 using namespace edm;
 using namespace std;
 using namespace reco;

 //Getting the Collections of L2ReconstructedJets from L1Seeds
 //and removing the collinear jets
 CaloJetCollection myTmpJets;

 int iL1Jet = 0;
 for( vtoken_cjets::const_iterator s = jetSrc_token.begin(); s != jetSrc_token.end(); ++ s ) {
   edm::Handle<CaloJetCollection> tauJets;
   iEvent.getByToken( * s, tauJets );
   for(CaloJetCollection::const_iterator iTau = tauJets->begin(); iTau !=tauJets->end(); ++iTau)
     { 
       if(iTau->et() > mEt_Min) {

     //Add the Pdg Id here 
     CaloJet myJet = *iTau;
     myJet.setPdgId(15);
     myTmpJets.push_back(myJet);
       }
     }
   iL1Jet++;
 }

 std::auto_ptr<CaloJetCollection> tauL2jets(new CaloJetCollection); 

 //Removing the collinear jets correctly!
 
 //First sort the jets you have merged
 SorterByPt sorter;
 std::sort(myTmpJets.begin(),myTmpJets.end(),sorter);
 
 //Remove Collinear Jets by prefering the highest ones!
 while(myTmpJets.size()>0) {
   tauL2jets->push_back(myTmpJets[0]);
   CaloJetCollection tmp;
   for(unsigned int i=1 ;i<myTmpJets.size();++i) {
     double DR = ROOT::Math::VectorUtil::DeltaR(myTmpJets[0].p4(),myTmpJets[i].p4());
     if(DR>0.1) 
       tmp.push_back(myTmpJets[i]);
   }
   myTmpJets.swap(tmp);
 }
 

  iEvent.put(tauL2jets);

}

void L2TauJetsMerger::fillDescriptions(edm::ConfigurationDescriptions& descriptions) 
{

  edm::ParameterSetDescription desc;
  std::vector<edm::InputTag> inputTags;  
  inputTags.push_back( edm::InputTag("hltAkIsoTau1Regional") );
  inputTags.push_back( edm::InputTag("hltAkIsoTau2Regional") );
  inputTags.push_back( edm::InputTag("hltAkIsoTau3Regional") );
  inputTags.push_back( edm::InputTag("hltAkIsoTau4Regional") );
  desc.add< std::vector<edm::InputTag> >("JetSrc",inputTags)->setComment("CaloJet collections to merge");
  desc.add<double>("EtMin",20.0)->setComment("Minimal ET of jet to merge");
  descriptions.setComment("Merges CaloJet collections removing duplicates");
  descriptions.add("L2TauJetsMerger",desc);

}