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

Go to the documentation of this file.

#include "HLTrigger/HLTcore/interface/HLTFilter.h"

//
// class declaration
//
class HLTHcalTowerFilter : public HLTFilter 
{
public:
  explicit HLTHcalTowerFilter(const edm::ParameterSet &);
  ~HLTHcalTowerFilter();
  
private:
  virtual bool filter(edm::Event &, const edm::EventSetup &);
  
  edm::InputTag inputTag_;    // input tag identifying product
  bool saveTag_;              // whether to save this tag
  double min_E_HB_;           // energy threshold for HB in GeV
  double min_E_HE_;           // energy threshold for HE in GeV
  double min_E_HF_;           // energy threshold for HF in GeV
  int max_N_HB_;              // maximum number for HB
  int max_N_HE_;              // maximum number for HB
  int max_N_HF_;              // maximum number for HB
};

#include <memory>

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"

//
// constructors and destructor
//
HLTHcalTowerFilter::HLTHcalTowerFilter(const edm::ParameterSet& config) :
  inputTag_ (config.getParameter<edm::InputTag>("inputTag")),
  saveTag_  (config.getUntrackedParameter<bool>("saveTag",false)),
  min_E_HB_ (config.getParameter<double>       ("MinE_HB")),
  min_E_HE_ (config.getParameter<double>       ("MinE_HE")),
  min_E_HF_ (config.getParameter<double>       ("MinE_HF")),
  max_N_HB_ (config.getParameter<int>          ("MaxN_HB")),
  max_N_HE_ (config.getParameter<int>          ("MaxN_HE")),
  max_N_HF_ (config.getParameter<int>          ("MaxN_HF"))
{
  // register your products
  produces<trigger::TriggerFilterObjectWithRefs>();
}

HLTHcalTowerFilter::~HLTHcalTowerFilter()
{
}

//
// member functions
//

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

  // get hold of collection of objects
  Handle<CaloTowerCollection> towers;
  event.getByLabel(inputTag_, towers);

  // look at all objects, check cuts and add to filter object
  int n_HB = 0;
  int n_HE = 0;
  int n_HF = 0;
  double abseta = 0;
  for(CaloTowerCollection::const_iterator i = towers->begin(); i != towers->end(); ++i) 
    {
      abseta = fabs(i->eta());
      if(abseta<1.305)
        {
          if(i->hadEnergy() >= min_E_HB_) 
            {
              n_HB++;
              //edm::Ref<CaloTowerCollection> ref(towers, std::distance(towers->begin(), i));
              //filterobject->addObject(TriggerJet, ref);
            }
        }
      else if(abseta>=1.305 && abseta<3)
        {
          if(i->hadEnergy() >= min_E_HE_)
            {
              n_HE++;
              //edm::Ref<CaloTowerCollection> ref(towers, std::distance(towers->begin(), i));
              //filterobject->addObject(TriggerJet, ref);
            }
        }
      else
        {
          if(i->hadEnergy() >= min_E_HF_) 
            {
              n_HF++;
              //edm::Ref<CaloTowerCollection> ref(towers, std::distance(towers->begin(), i));
              //filterobject->addObject(TriggerJet, ref);
            }
        }
    }

  // filter decision
  bool accept(n_HB<max_N_HB_ && n_HE<max_N_HE_ && n_HF<max_N_HF_ );

  // put filter object into the Event
  event.put(filterobject);

  return accept;
}

// define as a framework module
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(HLTHcalTowerFilter);