|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
|
// -*- C++ -*-
//
// Package: SiStripTools
// Class: SiStripDetWithSomething
//
/**\class SiStripDetWithSomething SiStripDetWithSomething.cc DPGAnalysis/SiStripTools/plugins/SiStripDetWithSomething.cc
Description: template EDFilter to select events with selected modules with SiStripDigis or SiStripClusters
Implementation:
*/
//
// 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/global/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::global::EDFilter<> {
public:
explicit SiStripDetWithSomething(const edm::ParameterSet&);
~SiStripDetWithSomething() override;
private:
void beginJob() override;
bool filter(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
void endJob() override;
// ----------member data ---------------------------
edm::EDGetTokenT<T> _digicollectionToken;
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)
: _digicollectionToken(consumes<T>(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::StreamID, edm::Event& iEvent, const edm::EventSetup& iSetup) const {
using namespace edm;
Handle<T> digis;
iEvent.getByToken(_digicollectionToken, 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);
|