/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/HLTrigger/special/src/HLTEcalIsolationFilter.cc

Go to the documentation of this file.

#include "HLTrigger/special/interface/HLTEcalIsolationFilter.h"

#include "DataFormats/HcalIsolatedTrack/interface/IsolatedPixelTrackCandidate.h"

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

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

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

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


HLTEcalIsolationFilter::HLTEcalIsolationFilter(const edm::ParameterSet& iConfig)
{
  candTag_ = iConfig.getParameter<edm::InputTag> ("EcalIsolatedParticleSource");
  maxhitout = iConfig.getParameter<int> ("MaxNhitOuterCone");
  maxhitin  = iConfig.getParameter<int> ("MaxNhitInnerCone");
  maxenin = iConfig.getParameter<double> ("MaxEnergyInnerCone");
  maxenout = iConfig.getParameter<double> ("MaxEnergyOuterCone");
  maxetacand = iConfig.getParameter<double> ("MaxEtaCandidate");  

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

HLTEcalIsolationFilter::~HLTEcalIsolationFilter(){}

bool HLTEcalIsolationFilter::filter(edm::Event& iEvent, const edm::EventSetup& iSetup)
{

  // The Filter object
  std::auto_ptr<trigger::TriggerFilterObjectWithRefs> filterproduct (new trigger::TriggerFilterObjectWithRefs(path(),module()));

  // Ref to Candidate object to be recorded in filter object
  edm::Ref<reco::IsolatedPixelTrackCandidateCollection> candref;

  // get hold of filtered candidates
  edm::Handle<reco::IsolatedPixelTrackCandidateCollection> ecIsolCands;
  iEvent.getByLabel(candTag_,ecIsolCands);

  //Filtering

  unsigned int n=0;
  for (unsigned int i=0; i<ecIsolCands->size(); i++)
    {
      candref = edm::Ref<reco::IsolatedPixelTrackCandidateCollection>(ecIsolCands, i);
        
      if ((candref->nHitIn()<=maxhitin)&&(candref->nHitOut()<=maxhitout)&&(candref->energyOut()<maxenout)&&(candref->energyIn()<maxenin)&&fabs(candref->eta())<maxetacand)
        {
          filterproduct->addObject(trigger::TriggerTrack, candref);
          n++;
        }
    }
  
  
  bool accept(n>0);

  iEvent.put(filterproduct);

  return accept;

}