EcalRecHitsFilter.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:   EcalRecHitsFilter
// Class:     EcalRecHitsFilter
//
//class EcalHighEnCosmicFilter EcalHighEnCosmicFilter.cc
//
// Original Author:  


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

// user include files
#include "CaloOnlineTools/EcalTools/plugins/EcalRecHitsFilter.h"

#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 "Geometry/CaloTopology/interface/CaloTopology.h"
#include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
#include "DataFormats/EgammaReco/interface/BasicCluster.h"
#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EcalElectronicsId.h"

using namespace edm;
using namespace std;
using namespace reco;

//
EcalRecHitsFilter::EcalRecHitsFilter(const edm::ParameterSet& iConfig) :
  NumBadXtalsThreshold_ (iConfig.getUntrackedParameter<int>("NumberXtalsThreshold")),
  EBRecHitCollection_ (iConfig.getParameter<edm::InputTag>("EcalRecHitCollectionEB"))

{
  EnergyCut = (iConfig.getUntrackedParameter<double>("energycut"));
}

EcalRecHitsFilter::~EcalRecHitsFilter()
{
}

bool EcalRecHitsFilter::filter( edm::Event& iEvent, const edm::EventSetup& iSetup)
{
  //int ievt = iEvent.id().event();
  Handle<EcalRecHitCollection> EBhits;
  iEvent.getByLabel(EBRecHitCollection_,EBhits);

  bool accepted = true;
  int nRecHitsGreater1GevPerEvent = 0;

  for (EcalRecHitCollection::const_iterator hitItr = EBhits->begin(); hitItr != EBhits->end(); ++hitItr)
    {

      EcalRecHit hit = (*hitItr);
      EBDetId det = hit.id();
  
      float ampli = hit.energy();
      if(ampli > EnergyCut /*1GeV*/)
    {
      nRecHitsGreater1GevPerEvent++;
      nRecHitsGreater1GevPerEvent_hist_MAP->Fill(det.iphi(),det.ieta());
    }
    }
  nRecHitsGreater1GevPerEvent_hist->Fill(nRecHitsGreater1GevPerEvent);
  if(nRecHitsGreater1GevPerEvent > NumBadXtalsThreshold_ ) accepted = false;
  return accepted;
}

void 
EcalRecHitsFilter::beginJob()
{
  nRecHitsGreater1GevPerEvent_hist = new TH1F("nRecHitsGreater1GevPerEvent_hist","nRecHitsGreater1GevPerEvent_hist",65000,0.,65000.);
  nRecHitsGreater1GevPerEvent_hist_MAP = new TH2F("nRecHitsGreater1GevPerEvent_hist_MAP","nRecHitsGreater1GevPerEvent_hist_MAP",360,1.,361.,171,-85.,86.);
}


void 
EcalRecHitsFilter::endJob()
{
  cout << "------EcalRecHitsFilter EndJob------>>>>>>>>>>" << endl;
  file = new TFile("RecHitFilter.root" , "RECREATE");
  file->cd();
  nRecHitsGreater1GevPerEvent_hist_MAP->Write();
  nRecHitsGreater1GevPerEvent_hist->Write();
  file->Close();
}