SiStripDetWithSomething.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripTools
// Class:      SiStripDetWithSomething
// 
//
// Original Author:  Andrea Venturi
//         Created:  Wed Oct 22 17:54:30 CEST 2008
//
//


// system include files
#include <memory>

// user include files
#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 <vector>

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

#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
#include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"

//
// class declaration
//

template <class T>
class SiStripDetWithSomething : public edm::EDFilter {
   public:
      explicit SiStripDetWithSomething(const edm::ParameterSet&);
      ~SiStripDetWithSomething();

   private:
      virtual void beginJob() ;
      virtual bool filter(edm::Event&, const edm::EventSetup&);
      virtual void endJob() ;
      
      // ----------member data ---------------------------

  edm::InputTag _digicollection;
  std::vector<unsigned int> _wantedmod;

};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
template <class T>
SiStripDetWithSomething<T>::SiStripDetWithSomething(const edm::ParameterSet& iConfig):
  _digicollection(iConfig.getParameter<edm::InputTag>("collectionName")),
  _wantedmod(iConfig.getUntrackedParameter<std::vector<unsigned int> >("selectedModules"))

{
   //now do what ever initialization is needed

  sort(_wantedmod.begin(),_wantedmod.end());

  edm::LogInfo("SelectedModules") << "Selected module list";
  for(std::vector<unsigned int>::const_iterator mod = _wantedmod.begin();mod!=_wantedmod.end();mod++) {
    edm::LogVerbatim("SelectedModules") << *mod ;
  }

}


template <class T>
SiStripDetWithSomething<T>::~SiStripDetWithSomething()
{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}


//
// member functions
//

// ------------ method called on each new Event  ------------
template <class T>
bool
SiStripDetWithSomething<T>::filter(edm::Event& iEvent, const edm::EventSetup& iSetup)
{
   using namespace edm;

   Handle<T> digis;
   iEvent.getByLabel(_digicollection,digis);

   for(typename T::const_iterator it = digis->begin();it!=digis->end();it++) {

     for(std::vector<unsigned int>::const_iterator mod=_wantedmod.begin();
     mod!=_wantedmod.end()&&it->detId()>=*mod;
     mod++) {
       if(*mod == it->detId()) {
     edm::LogInfo("ModuleFound") << " module " << *mod << " found with " 
                     << it->size() << " digis/clusters";
     return true;
       }
     }
   }

   return false;
}

// ------------ method called once each job just before starting event loop  ------------
template <class T>
void 
SiStripDetWithSomething<T>::beginJob()
{
}

// ------------ method called once each job just after ending the event loop  ------------
template <class T>
void 
SiStripDetWithSomething<T>::endJob() {
}

typedef SiStripDetWithSomething<edm::DetSetVector<SiStripDigi> > SiStripDetWithDigi;
typedef SiStripDetWithSomething<edmNew::DetSetVector<SiStripCluster> > SiStripDetWithCluster;

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