SelectHFMinBias.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:   BeamSplash
// Class:     BeamSPlash
//
//
// Original Author:  Luca Malgeri

#include <memory>
#include <vector>
#include <map>
#include <set>

// user include files

#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDFilter.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DPGAnalysis/Skims/interface/SelectHFMinBias.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"

using namespace edm;
using namespace std;

SelectHFMinBias::SelectHFMinBias(const edm::ParameterSet& iConfig)
{
}

SelectHFMinBias::~SelectHFMinBias()
{
}

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

edm::Handle<CaloTowerCollection> towers;
iEvent.getByLabel("towerMaker",towers);


int negTowers = 0;
int posTowers = 0;
for(CaloTowerCollection::const_iterator cal = towers->begin(); cal != towers->end(); ++cal) {
   for(unsigned int i = 0; i < cal->constituentsSize(); i++) {
      const DetId id = cal->constituent(i);
      if(id.det() == DetId::Hcal) {
        HcalSubdetector subdet=(HcalSubdetector(id.subdetId()));
        if(subdet == HcalForward) {
          if(cal->energy()>3. && cal->eta()<-3.)
            negTowers++;
          if(cal->energy()>3. && cal->eta()>3.)
            posTowers++;
        }
     }
   }
}
if(negTowers>0 && posTowers>0)
  return true;

return false;

}

//define this as a plug-in
DEFINE_FWK_MODULE(SelectHFMinBias);