/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/HLTrigger/JetMET/src/HLT2jetGapFilter.cc

Go to the documentation of this file.

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

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

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

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

#include "DataFormats/JetReco/interface/CaloJetCollection.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"


//
// constructors and destructor
//
HLT2jetGapFilter::HLT2jetGapFilter(const edm::ParameterSet& iConfig)
{
   inputTag_ = iConfig.getParameter< edm::InputTag > ("inputTag");
   saveTag_    = iConfig.getUntrackedParameter<bool>("saveTag",false);
   minEt_   = iConfig.getParameter<double> ("minEt");
   minEta_= iConfig.getParameter<double> ("minEta"); 

   //register your products
   produces<trigger::TriggerFilterObjectWithRefs>();
}

HLT2jetGapFilter::~HLT2jetGapFilter(){}


// ------------ method called to produce the data  ------------
bool
HLT2jetGapFilter::filter(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  using namespace trigger;
  // The filter object
  std::auto_ptr<trigger::TriggerFilterObjectWithRefs> 
    filterobject (new trigger::TriggerFilterObjectWithRefs(path(),module()));
  if (saveTag_) filterobject->addCollectionTag(inputTag_);

  edm::Handle<reco::CaloJetCollection> recocalojets;
  iEvent.getByLabel(inputTag_,recocalojets);

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

//  std::cout << "HLT2jetGapFilter " << recocalojets->size() << " jets in this event" << std::endl;


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

    double etjet1=0.;
    double etjet2=0.;
    double etajet1=0.;
    double etajet2=0.;
    int countjets =0;

    for (reco::CaloJetCollection::const_iterator recocalojet = recocalojets->begin(); 
         recocalojet<=(recocalojets->begin()+1); recocalojet++) {
      
      if(countjets==0) {
        etjet1 = recocalojet->et();
        etajet1 = recocalojet->eta();
      }
      if(countjets==1) {
        etjet2 = recocalojet->et();
        etajet2 = recocalojet->eta();
      }
      countjets++;
    }

  if(etjet1>minEt_ && etjet2>minEt_ && (etajet1*etajet2)<0 && fabs(etajet1)>minEta_ && fabs(etajet2)>minEta_){
      for (reco::CaloJetCollection::const_iterator recocalojet = recocalojets->begin(); 
           recocalojet<=(recocalojets->begin()+1); recocalojet++) {
        reco::CaloJetRef ref(reco::CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet)));
        filterobject->addObject(TriggerJet,ref);
        n++;
      }
    }
    
  } // events with two or more jets
  
  
  
  // filter decision
  bool accept(n>=2);
  
  // put filter object into the Event
  iEvent.put(filterobject);
  
  return accept;
}