HLTAcoFilter.cc

Go to the documentation of this file.

#include <string>
#include "HLTrigger/JetMET/interface/HLTAcoFilter.h"

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

#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"

#include "FWCore/MessageLogger/interface/MessageLogger.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"
#include "FWCore/Utilities/interface/InputTag.h"


//
// constructors and destructor
//
HLTAcoFilter::HLTAcoFilter(const edm::ParameterSet& iConfig) : HLTFilter(iConfig)
{
   inputJetTag_ = iConfig.getParameter< edm::InputTag > ("inputJetTag");
   inputMETTag_ = iConfig.getParameter< edm::InputTag > ("inputMETTag");
   minDPhi_     = iConfig.getParameter<double> ("minDeltaPhi");
   maxDPhi_     = iConfig.getParameter<double> ("maxDeltaPhi");
   minEtjet1_   = iConfig.getParameter<double> ("minEtJet1");
   minEtjet2_   = iConfig.getParameter<double> ("minEtJet2");
   AcoString_   = iConfig.getParameter<std::string> ("Acoplanar");

   m_theJetToken = consumes<reco::CaloJetCollection>(inputJetTag_);
   m_theMETToken = consumes<trigger::TriggerFilterObjectWithRefs>(inputMETTag_);
}

HLTAcoFilter::~HLTAcoFilter()= default;

void
HLTAcoFilter::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  makeHLTFilterDescription(desc);
  desc.add<edm::InputTag>("inputJetTag",edm::InputTag("IterativeCone5CaloJets"));
  desc.add<edm::InputTag>("inputMETTag",edm::InputTag("MET"));
  desc.add<double>("minDeltaPhi",0.0);
  desc.add<double>("maxDeltaPhi",2.0);
  desc.add<double>("minEtJet1",20.0);
  desc.add<double>("minEtJet2",20.0);
  desc.add<std::string>("Acoplanar","Jet1Jet2");
  descriptions.add("hltAcoFilter",desc);
}

// ------------ method called to produce the data  ------------
bool
HLTAcoFilter::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);
  int JetNum = recocalojets->size();

  // events with two or more jets

  double etjet1=0.;
  double etjet2=0.;
  double phijet1=0.;
  double phijet2=0.;
  //double etmiss=0.;
  double phimiss=0.;

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

  CaloJetRef ref1,ref2;

  if (JetNum>0) {
    auto recocalojet = recocalojets->begin();
    
    etjet1 = recocalojet->et();
    phijet1 = recocalojet->phi();
    ref1  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));

    if(JetNum>1) {
      recocalojet++;
      etjet2 = recocalojet->et();
      phijet2 = recocalojet->phi();
      ref2  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));
    }
    double Dphi= -1.;
    int JetSel = 0;

    if (AcoString_ == "Jet2Met") {
      Dphi = std::abs(phimiss-phijet2);
      if (JetNum>=2 && etjet1>minEtjet1_  && etjet2>minEtjet2_) {JetSel=1;}
    }
    if (AcoString_ == "Jet1Jet2") {
      Dphi = std::abs(phijet1-phijet2);
      if (JetNum>=2 && etjet1>minEtjet1_  && etjet2>minEtjet2_) {JetSel=1;}
    }
    if (AcoString_ == "Jet1Met") {
      Dphi = std::abs(phimiss-phijet1);
      if (JetNum>=1 && etjet1>minEtjet1_ ) {JetSel=1;}
    }


    if (Dphi>M_PI) {Dphi=2.0*M_PI-Dphi;}
    if(JetSel>0 && Dphi>=minDPhi_ && Dphi<=maxDPhi_){

      if (AcoString_=="Jet2Met" || AcoString_=="Jet1Met")  {filterproduct.addObject(TriggerMET,metRef);}
      if (AcoString_=="Jet1Met" || AcoString_=="Jet1Jet2") {filterproduct.addObject(TriggerJet,ref1);}
      if (AcoString_=="Jet2Met" || AcoString_=="Jet1Jet2") {filterproduct.addObject(TriggerJet,ref2);}
      n++;
    }


  } // at least one jet


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

  return accept;
}