HLTNVFilter.cc

Go to the documentation of this file.

#include "HLTrigger/JetMET/interface/HLTNVFilter.h"

#include "DataFormats/Common/interface/Handle.h"


#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/METReco/interface/CaloMET.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"


//
// constructors and destructor
//
HLTNVFilter::HLTNVFilter(const edm::ParameterSet& iConfig) : HLTFilter(iConfig)
{
   inputJetTag_ = iConfig.getParameter< edm::InputTag > ("inputJetTag");
   inputMETTag_ = iConfig.getParameter< edm::InputTag > ("inputMETTag");
   minNV_   = iConfig.getParameter<double> ("minNV");
   minEtjet1_= iConfig.getParameter<double> ("minEtJet1");
   minEtjet2_ = iConfig.getParameter<double> ("minEtJet2");
   m_theJetToken = consumes<reco::CaloJetCollection>(inputJetTag_);
   m_theMETToken = consumes<trigger::TriggerFilterObjectWithRefs>(inputMETTag_);
}

HLTNVFilter::~HLTNVFilter(){}


void HLTNVFilter::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  makeHLTFilterDescription(desc);
  desc.add<edm::InputTag>("inputJetTag",edm::InputTag("iterativeCone5CaloJets"));
  desc.add<edm::InputTag>("inputMETTag",edm::InputTag("hlt1MET60"));
  desc.add<double>("minEtJet2",20.);
  desc.add<double>("minEtJet1",80.);
  desc.add<double>("minNV",0.1);
  descriptions.add("hltNVFilter",desc);
}

// ------------ method called to produce the data  ------------
bool
HLTNVFilter::hltFilter(edm::Event& iEvent, const edm::EventSetup& iSetup, trigger::TriggerFilterObjectWithRefs & filterproduct) const
{
  using namespace std;
  using namespace edm;
  using namespace reco;
  using namespace trigger;

  // The filter object
  if (saveTags()) {
    filterproduct.addCollectionTag(inputJetTag_);
    filterproduct.addCollectionTag(inputMETTag_);
  }

  Handle<CaloJetCollection> recocalojets;
  iEvent.getByToken(m_theJetToken,recocalojets);
  Handle<trigger::TriggerFilterObjectWithRefs> metcal;
  iEvent.getByToken(m_theMETToken,metcal);


  // look at all candidates,  check cuts and add to filter object
  int n(0);

  if(recocalojets->size() > 1){
    // events with two or more jets

    double etjet1=0.;
    double etjet2=0.;
    double etmiss=0.;
    int countjets =0;

    VRcalomet vrefMET;
    metcal->getObjects(TriggerMET,vrefMET);
    CaloMETRef metRef=vrefMET.at(0);
    etmiss=vrefMET.at(0)->et();

    CaloJetRef ref1,ref2;
    for (CaloJetCollection::const_iterator recocalojet = recocalojets->begin();
     recocalojet<=(recocalojets->begin()+1); recocalojet++) {

      if(countjets==0) {
    etjet1 = recocalojet->et();
                ref1  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));
      }
      if(countjets==1) {
    etjet2 = recocalojet->et();
                ref2  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));
      }
      countjets++;
    }

    double NV = (etmiss*etmiss-(etjet1-etjet2)*(etjet1-etjet2))/(etjet2*etjet2);
    if(etjet1>minEtjet1_  && etjet2>minEtjet2_ && NV>minNV_){
      filterproduct.addObject(TriggerMET,metRef);
      filterproduct.addObject(TriggerJet,ref1);
      filterproduct.addObject(TriggerJet,ref2);
      n++;
    }

  } // events with two or more jets



  // filter decision
  bool accept(n>=1);

  return accept;
}